数学建模社区-数学中国

标题: 救命！！！！！！！！2000年全国建模bB用LINGO求解程序 [打印本页]

作者: pphy1314 时间: 2006-12-23 14:47
标题: 救命！！！！！！！！2000年全国建模bB用LINGO求解程序

sets:
GCH/S1 S2 S3 S4 S5 S6 S7 /;SI:
ZHD/A2,A3,A5,A6,A7,A8,A9,A10,A11,A12,A13,A14,A15/:HM,YJ,ZJ,AJ;
YL(GCH,ZHD):C,X
ENDSETS
DATA:
SI=800,800,1000,2000,2000,2000,3000;
HM=2,3,5,6,7,8,9,10,11,12,13,14,15;
AJ=301,1356,194,205,201,680,480,300,220,210,420,500,0;
C=
320.3,300.2,198,180.5,163.1,181.2,224.2,252,256,266,281.2,288,302
360.3,345.2,266,250.5,243.1,226.2,269.2,297,301,311,326.2,333,347
375.3,355.2,276,260.5,253.1,241.2,203.2,237,241,251,266.2,273,287
410.3,395.2,316,300.5,293.1,276.2,244.2,222,211,221,236.2,243,257
400.3,380.2,301,285.5,278.1,266.2,234.2,212,188,206,226.2,228,242
405.3,385.2,306,290.5,283.1,271.2,234.2,212,201,195,176.2,161,178
425.3,405.2,326,310.5,303.1,291.2,259.2,237,226,216,198.2,186,162;
ENDDATA
!The objective:
MIN=@SUM[YL(I,J):C(I,J)*X(I,J)]+0.05*@SUM(ZHD:YJ+YJ*YJ+ZJ+ZJ*ZJ);
!The constraints;
[Y1]@FOR(GCH(I)SUM(ZHD(J):X(I,J))<=SI(I));
[Y2]@SUM(ZHD(J):X(7,J))=0;
[Y3]@FOR(ZHD(J)SUM(GCH(I):X(I,J))=YJ(J)+ZJ(J));
[Y4]@FOR(ZHD(J)|HM(J)#LT#15:ZJ(J)+YJ(J+1)=A(J));
[Y5]YJ(1)=104;
[Y6]ZJ(13)=0;
END

运行不出来！找错误！！！！

作者: qazxcv11 时间: 2007-1-10 22:41
找不出来！

作者: basketchch0503 时间: 2007-4-5 13:14

不太懂

作者: wujianjack2 时间: 2013-7-22 13:47
程序本身其实问题不大，主要有如下两处：

1.集合段"ZHD/A2,A3,A5..A15/:HM,YJ,ZJ,AJ;"看上去没什么问题，但我通过实际操作，发现LINGO并不支持这样的定义，因此我建议改定义为"ZHD/A3,A4,A5..A15/:HM,YJ,ZJ,AJ;"然后在注释里说明”这里的A3表示实际的A2，A4表示实际的A3“，可能这是LINGO功能的一个缺陷吧！

2.目标函数中@SUM()函数的圆括号误写为"[]".

按照我的方案，LINGO程序代码如下：
MODEL:
SETS:
GCH/S1..S7/:SI;
ZHD/A3,A4,A5..A15/:HM,YJ,ZJ,AJ;
YL(GCH,ZHD):C,X;
ENDSETS
DATA:
SI=800,800,1000,2000,2000,2000,3000;
HM=2,3,5,6,7,8,9,10,11,12,13,14,15;
AJ=301,1356,194,205,201,680,480,300,220,210,420,500,0;
C=
320.3 300.2 198 180.5 163.1 181.2 224.2 252 256 266 281.2 288 302
360.3 345.2 266 250.5 241 226.2 269.2 297 301 311 326.2 333 347
375.3 355.2 276 260.5 251 241.2 203.2 237 241 251 266.2 273 287
410.3 395.2 316 300.5 291 276.2 244.2 222 211 221 236.2 243 257
400.3 380.2 301 285.5 276 266.2 234.2 212 188 206 226.2 228 242
405.3 385.2 306 290.5 281 271.2 234.2 212 201 195 176.2 161 178
425.3 405.2 326 310.5 301 291.2 259.2 236 226 216 198.2 186 162;
ENDDATA
MIN=@SUM(YL(I,J):C(I,J)*X(I,J))+0.05*@SUM(ZHD:YJ+YJ*YJ+ZJ+ZJ*ZJ);
@FOR(GCH(I):@SUM(ZHD(J):X(I,J))<=SI(I));
@SUM(ZHD(J):X(7,J))=0;
@FOR(ZHD(J):@SUM(GCH(I):X(I,J))=YJ(J)+ZJ(J));
@FOR(ZHD(J)|HM(J)#LT#15:ZJ(J)+YJ(J+1)=AJ(J));
YJ(1)=104;ZJ(13)=0;
END

运行结果为(仅列出X(I,J)的非零值)：
  Local optimal solution found.
  Objective value:                            1278632.
  Infeasibilities:                            0.1136868E-12
  Total solver iterations:                         38

                     Variable          Value       Reduced Cost
                  X( S1, A5)       334.5000          0.000000
                  X( S1, A6)       200.0000          0.000000
                  X( S1, A7)       265.5000          0.000000
                  X( S2, A5)       179.0000          0.000000
                  X( S2, A4)       96.96236          0.000000
                  X( S2, A5)       224.0376          0.000000
                  X( S2, A8)       300.0000          0.000000
                  X( S3, A4)       63.08959          0.000000
                  X( S3, A5)       272.9104          0.000000
                  X( S3, A9)       664.0000          0.000000
                  X( S5, A4)       347.9480          0.000000
                  X( S5, A5)       252.0520          0.000000
                  X( S5, A10)       312.9032          0.000000
                  X( S5, A11)       415.0000          0.000000
                  X( S6, A10)       38.09683          0.000000
                  X( S6, A12)       86.00000          0.000000
                  X( S6, A13)       333.0000          0.000000
                  X( S6, A14)       621.0000          0.000000
                  X( S6, A15)       165.0000          0.000000

以上结果仅供参考，欢迎大家指出错误并交流，谢谢！

作者: 弹你脑瓜崩 时间: 2019-7-4 15:54
学习了

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)