救命！！！！！！！！2000年全国建模bB用LINGO求解程序

pphy1314

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

找不出来！

basketchch0503

不太懂

wujianjack2

程序本身其实问题不大，主要有如下两处：

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

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

弹你脑瓜崩

学习了