24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1022  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

XZQFDU

金虫 (小有名气)

[求助] 求以下非线性Lingo代码,真的不会编已有1人参与

求以下非线性Lingo代码,真的不会编

求以下非线性Lingo代码,真的不会编
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TRAMS

新虫 (初入文坛)

【答案】应助回帖

本人尝试了LINGO编码,未必正确,仅供参考,求解模型代码如下,亦可见附件task2.lg4。希望将公式背景介绍得尽量详细,截图更清晰点,方便编码。

model:
!帖子
求以下非线性Lingo代码,真的不会编
http://muchong.com/bbs/viewthread.php?tid=10575262&fpage=1&target=blank
;
sets:
    D/1..15/value;!通过修改15为 16/17/18/19/20 依次求解模型;!本人推测n应该为正整数;
endsets
data:
enddata

n=@size(D);
min=@sum(D(i):280-Dvalue(i)) + 245 + @sum(D(i)|i#ge#2Dvalue(i-1)-Dvalue(i))^2*0.005) + (280-DValue(1))^2*0.005 + (DValue(n)-35)^2*0.005;

@for(D(i):
    Dvalue(i)<=280;
);
@for(D(i):
    Dvalue(i)>=35;
);

@for(D(i)|i#ge#2:
    (Dvalue(i-1)-Dvalue(i)) / Dvalue(i-1) <=0.13;
);
@for(D(i)|i#ge#2:
    (Dvalue(i-1)-Dvalue(i)) / Dvalue(i-1) >=0.1;
);

(Dvalue(n)-35) / Dvalue(n) <= 0.13;

@for(D(i)|i#ge#2:
    Dvalue(i-1)>=Dvalue(i);
);
!@bnd(15,n+1,20);!本人推测n应该为正整数;
end

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : task2.lg4
  • 2017-02-01 16:24:32, 32 K
6楼2017-02-01 16:25:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

rsd014

木虫 (正式写手)

lingo不会 不过这玩意用MATLAB倒是很简单

发自小木虫IOS客户端
2楼2016-08-19 09:17:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

XZQFDU

金虫 (小有名气)

引用回帖:
2楼: Originally posted by rsd014 at 2016-08-19 09:17:46
lingo不会 不过这玩意用MATLAB倒是很简单

Matlab也可以,只是我不知道N也是参数,怎么处理

发自小木虫IOS客户端
3楼2016-08-19 10:24:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rsd014

木虫 (正式写手)

引用回帖:
3楼: Originally posted by XZQFDU at 2016-08-19 10:24:58
Matlab也可以,只是我不知道N也是参数,怎么处理
...

先弄成正常的数学规划吧

发自小木虫IOS客户端
4楼2016-08-19 13:14:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见