数学建模社区-数学中国

标题: 第一次用lingo，求大神帮忙改错，写论文急用，已经纠结好多天了，谢谢了 [打印本页]
作者: 童话王子    时间: 2015-5-21 22:25
标题: 第一次用lingo，求大神帮忙改错，写论文急用，已经纠结好多天了，谢谢了
各位大神，求帮忙查查错误，论文马上就要审查了，可是关键部分还是没有搞懂，第一次用lingo，看了好久都没看明白，现在很急用，我写的代码和原问题模型都在附件里了，求大神帮忙，谢谢了！

LINGO1.lg4

4.5 KB, 下载次数: 3, 下载积分: 体力 -2 点

这是我写的lingo代码，运行不出来

基于混合整数规划的铁路物流中心选址建模求解——与运筹学中的分枝定界类似，Lingo求解.pdf

255.78 KB, 下载次数: 4, 下载积分: 体力 -2 点

这是相关问题模型
作者: 童话王子    时间: 2015-5-21 22:36
跪求各位大神了，最好告诉我哪错了 我好能学习下
作者: 欧阳嫦娥    时间: 2015-7-22 15:26
为什么下载你给的程序竟然还要扣体力
作者: liwenhui    时间: 2015-7-23 10:54
欧阳嫦娥 发表于 2015-7-22 15:26
为什么下载你给的程序竟然还要扣体力

论坛把每个附件都作为一种资源，获取资源就会付出一定成本。
作者: liwenhui    时间: 2015-7-23 10:58
童话王子 发表于 2015-5-21 22:36
跪求各位大神了，最好告诉我哪错了 我好能学习下

发现第一个@For有问题，这个for针对的集合是factory，但后面的@sum被要求小于的值p不属于集合factory，这样当运行到这里时候Lingo不知道你这个约束到底是什么意思。
作者: liwenhui    时间: 2015-7-23 11:05
童话王子 发表于 2015-5-21 22:36
跪求各位大神了，最好告诉我哪错了 我好能学习下

做了一个简单的更改，能求解到局部最优解。全局最优解我的LINGO求解不出。
  1. model:
  2. sets:
  3. factory/F1..F9/:d;
  4. warhouse/C1..C5/:P,Q,E;
  5. tr/tr1..tr5/:Z;
  6. link1(factory,warhouse):S,X;
  7. endsets
  8. data:
  9. d=1500 1400 2000 1430 1220 2210 1810 1130 1350;
  10. P=12000 12000 21000 21000 12000;
  11. Q=2133 1890 2928 2815 1854;
  12. S=9.2 5.4 6.9 14.8 14.6 8.4 7.9 10.1 10.9
  13. 12.4 10.1 4.7 11.1 9.6 10.9 10.1 9.9 9.6
  14. 14.4 14.8 9.6 10.1 4.5 8.4 7.5 8.1 6.2
  15. 7.1 9.9 13.9 7.9 12.2 7.5 8.8 7.7 8.4
  16. 9.0 12.0 15.0 6.6 9.6 7.1 7.9 7.5 7.7;
  17. b=0.0728;
  18. a=100;
  19. c=3;
  20. E=11.6 7.5 8.8 15.6 15.2;
  21. h=0.0612;
  22. enddata
  23. min=@sum(link1(i,j):b*(S(i,j)*X(i,j))*Z(j))+@sum(link1(i,j):a*(X(i,j)-Q(j))*Z(j))+@sum(link1(i,j):C*X(i,j)*Z(j))+@sum(warhouse(j):h*E(j)*Z(j));
  24. @for(warhouse(j):@sum(factory(i):X(i,j))<P(j));
  25. @for(factory(i):@sum(warhouse(j):X(i,j))=d(i));
  26. @for(warhouse(j):9999*P(j)>@sum(link1(i,j):X(i,j)));
  27. @for(warhouse(j):
  28. @sum(tr(j):Z(j))<=3);
  29. @for(tr(j):@bin(Z));
  30. end
复制代码
结果显示如下：

  Local optimal solution found.
  Objective value:                             -9004142.
  Objective bound:                             -9004142.
  Infeasibilities:                              0.000000
  Extended solver steps:                               0
  Total solver iterations:                            29

  Model Class:                                     MINLP

  Total variables:                     50
  Nonlinear variables:                 50
  Integer variables:                    5

  Total constraints:                   25
  Nonlinear constraints:                1

  Total nonzeros:                     190
  Nonlinear nonzeros:                  50

                                                    Variable           Value
                                                           B       0.7280000E-01
                                                           A        100.0000
                                                           C        3.000000
                                                           H       0.6120000E-01
                                                      D( F1)        1500.000
                                                      D( F2)        1400.000
                                                      D( F3)        2000.000
                                                      D( F4)        1430.000
                                                      D( F5)        1220.000
                                                      D( F6)        2210.000
                                                      D( F7)        1810.000
                                                      D( F8)        1130.000
                                                      D( F9)        1350.000
                                                      P( C1)        12000.00
                                                      P( C2)        12000.00
                                                      P( C3)        21000.00
                                                      P( C4)        21000.00
                                                      P( C5)        12000.00
                                                      Q( C1)        2133.000
                                                      Q( C2)        1890.000
                                                      Q( C3)        2928.000
                                                      Q( C4)        2815.000
                                                      Q( C5)        1854.000
                                                      E( C1)        11.60000
                                                      E( C2)        7.500000
                                                      E( C3)        8.800000
                                                      E( C4)        15.60000
                                                      E( C5)        15.20000
                                                     Z( TR1)        1.000000
                                                     Z( TR2)        1.000000
                                                     Z( TR3)        1.000000
                                                     Z( TR4)        1.000000
                                                     Z( TR5)        1.000000
                                                  S( F1, C1)        9.200000
                                                  S( F1, C2)        5.400000
                                                  S( F1, C3)        6.900000
                                                  S( F1, C4)        14.80000
                                                  S( F1, C5)        14.60000
                                                  S( F2, C1)        8.400000
                                                  S( F2, C2)        7.900000
                                                  S( F2, C3)        10.10000
                                                  S( F2, C4)        10.90000
                                                  S( F2, C5)        12.40000
                                                  S( F3, C1)        10.10000
                                                  S( F3, C2)        4.700000
                                                  S( F3, C3)        11.10000
                                                  S( F3, C4)        9.600000
                                                  S( F3, C5)        10.90000
                                                  S( F4, C1)        10.10000
                                                  S( F4, C2)        9.900000
                                                  S( F4, C3)        9.600000
                                                  S( F4, C4)        14.40000
                                                  S( F4, C5)        14.80000
                                                  S( F5, C1)        9.600000
                                                  S( F5, C2)        10.10000
                                                  S( F5, C3)        4.500000
                                                  S( F5, C4)        8.400000
                                                  S( F5, C5)        7.500000
                                                  S( F6, C1)        8.100000
                                                  S( F6, C2)        6.200000
                                                  S( F6, C3)        7.100000
                                                  S( F6, C4)        9.900000
                                                  S( F6, C5)        13.90000
                                                  S( F7, C1)        7.900000
                                                  S( F7, C2)        12.20000
                                                  S( F7, C3)        7.500000
                                                  S( F7, C4)        8.800000
                                                  S( F7, C5)        7.700000
                                                  S( F8, C1)        8.400000
                                                  S( F8, C2)        9.000000
                                                  S( F8, C3)        12.00000
                                                  S( F8, C4)        15.00000
                                                  S( F8, C5)        6.600000
                                                  S( F9, C1)        9.600000
                                                  S( F9, C2)        7.100000
                                                  S( F9, C3)        7.900000
                                                  S( F9, C4)        7.500000
                                                  S( F9, C5)        7.700000
                                                  X( F1, C1)        0.000000
                                                  X( F1, C2)        1500.000
                                                  X( F1, C3)        0.000000
                                                  X( F1, C4)        0.000000
                                                  X( F1, C5)        0.000000
                                                  X( F2, C1)        0.000000
                                                  X( F2, C2)        1400.000
                                                  X( F2, C3)        0.000000
                                                  X( F2, C4)        0.000000
                                                  X( F2, C5)        0.000000
                                                  X( F3, C1)        0.000000
                                                  X( F3, C2)        2000.000
                                                  X( F3, C3)        0.000000
                                                  X( F3, C4)        0.000000
                                                  X( F3, C5)        0.000000
                                                  X( F4, C1)        0.000000
                                                  X( F4, C2)        0.000000
                                                  X( F4, C3)        1430.000
                                                  X( F4, C4)        0.000000
                                                  X( F4, C5)        0.000000
                                                  X( F5, C1)        0.000000
                                                  X( F5, C2)        0.000000
                                                  X( F5, C3)        1220.000
                                                  X( F5, C4)        0.000000
                                                  X( F5, C5)        0.000000
                                                  X( F6, C1)        0.000000
                                                  X( F6, C2)        2210.000
                                                  X( F6, C3)        0.000000
                                                  X( F6, C4)        0.000000
                                                  X( F6, C5)        0.000000
                                                  X( F7, C1)        0.000000
                                                  X( F7, C2)        0.000000
                                                  X( F7, C3)        1810.000
                                                  X( F7, C4)        0.000000
                                                  X( F7, C5)        0.000000
                                                  X( F8, C1)        0.000000
                                                  X( F8, C2)        0.000000
                                                  X( F8, C3)        0.000000
                                                  X( F8, C4)        0.000000
                                                  X( F8, C5)        1130.000
                                                  X( F9, C1)        0.000000
                                                  X( F9, C2)        1350.000
                                                  X( F9, C3)        0.000000
                                                  X( F9, C4)        0.000000
                                                  X( F9, C5)        0.000000

                                     Row    Slack or Surplus      Dual Price
                                       1       -9004142.           -1.000000
                                       2        12000.00            0.000000
                                       3        3540.000            0.000000
                                       4        16540.00            0.000000
                                       5        21000.00            0.000000
                                       6        10870.00            0.000000
                                       7        0.000000            0.000000
                                       8        0.000000            0.000000
                                       9        0.000000            0.000000
                                      10        0.000000            0.000000
                                      11        0.000000            0.000000
                                      12        0.000000            0.000000
                                      13        0.000000            0.000000
                                      14        0.000000            0.000000
                                      15        0.000000            0.000000
                                      16       0.1199880E+09        0.000000
                                      17       0.1199795E+09        0.000000
                                      18       0.2099745E+09        0.000000
                                      19       0.2099790E+09        0.000000
                                      20       0.1199869E+09        0.000000
                                      21        2.000000            0.000000
                                      22        2.000000            0.000000
                                      23        2.000000            0.000000
                                      24        2.000000            0.000000
                                      25        2.000000            0.000000




作者: 欧阳嫦娥    时间: 2015-7-23 18:41
liwenhui 发表于 2015-7-23 10:54
论坛把每个附件都作为一种资源，获取资源就会付出一定成本。

这样   




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5