| 查看: 1025 | 回复: 6 | ||
XZQFDU金虫 (小有名气)
|
[求助]
求以下非线性Lingo代码,真的不会编已有1人参与
|
2楼2016-08-19 09:17:46
XZQFDU
金虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 1700.1
- 散金: 99
- 帖子: 161
- 在线: 92.5小时
- 虫号: 3148354
- 注册: 2014-04-19
- 专业: 对流层大气物理学
3楼2016-08-19 10:24:58
4楼2016-08-19 13:14:35
5楼2016-08-19 14:02:29
【答案】应助回帖
|
本人尝试了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#2 Dvalue(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
7楼2017-02-01 16:27:55













回复此楼

value;!通过修改15为 16/17/18/19/20 依次求解模型;!本人推测n应该为正整数;
Dvalue(i-1)-Dvalue(i))^2*0.005) + (280-DValue(1))^2*0.005 + (DValue(n)-35)^2*0.005;