数学建模社区-数学中国

标题: 向大神求教！用lingo12.0出错但用13.0可以，都是打开了global solve后出现的结果 [打印本页]
作者: zengyixuan401    时间: 2013-7-22 14:22
标题: 向大神求教！用lingo12.0出错但用13.0可以，都是打开了global solve后出现的结果
  1. MODEL:
  2. Title Location Problem;
  3. sets:
  4. demand/1..6/:a,b,d;
  5. supply/1,2/:x,y,e;
  6. link(demand,supply):c;
  7. endsets
  8. data:
  9. a=1.25,8.75,0.5,5.75,3,7.25;
  10. b=1.25 0.75 4.75 5 6.5 7.75;
  11. d=3,5,4,7,6,11;
  12. e=20,20;
  13. enddata
  14. init:
  15. x=5,2;y=1,7;
  16. endinit
  17. [obj]min=@sum(link(i,j):c(i,j)*((x(j)-a(i))^2+(y(j)-b(i))^2)^(1/2));
  18. @for(demand(i):@sum(supply(j):c(i,j))=d(i););
  19. @for(supply(j):@sum(demand(i):c(i,j))<=e(j););
  20. @for(supply:@free(x);@free(y));
  21. end
复制代码
12.0的没有结果，13.0的extended solver status栏目里没有数据
求解啊啊啊
感谢！！1
作者: wujianjack2    时间: 2013-7-22 20:43
本帖最后由 wujianjack2 于 2014-5-11 11:31 编辑

楼主的问题不错！是在看谢金星，薛毅老师的《优化建模与LINDO/LINGO软件》这本书吗？
首先，楼主的程序我个人觉得是没有什么问题，关键可能在于你使用的版本，我没有安装LINGO 12和LINGO 13，但如果你的LINGO软件的规模限制并不是"Unlimited"，则很可能在使用Global Solver时提示Liicense is too small for this problem这类的错误提示，你可以尝试用我的帖子中的LINGO 11运行下，不必选用Global Solver则可得到最优解。

程序代码如下：
MODEL:
TITLE LOCATION PROBLEM;
SETS:
DEMAND/1..6/:A,B,D;
SUPPLY/1..2/:X,Y,E;
LINK(DEMAND,SUPPLY):C;
ENDSETS
DATA:
A=1.25,8.75,0.5,5.75,3,7.25;
B=1.25,0.75,4.75,5,6.5,7.75;
D=3,5,4,7,6,11;
E=20,20;
ENDDATA
INIT:
X=5,2;Y=1,7;
ENDINIT
[OBJ] MIN=@SUM(LINK(I,J):C(I,J)*((X(J)-A(I))^2+(Y(J)-B(I))^2)^(1/2));
@FOR(DEMAND(I):@SUM(SUPPLY(J):C(I,J))=D(I););
@FOR(SUPPLY(J):@SUM(DEMAND(I):C(I,J))<=E(J););
@FOR(SUPPLY:@FREE(X);@FREE(Y););
END

运行的结果如下（仅给出重要的结果）：
  Local optimal solution found.
  Objective value:                              85.26604
  Infeasibilities:                              0.000000
  Total solver iterations:                            68


  Model Title: LOCATION PROBLEM
  X( 1)        3.254883            0.000000
  X( 2)        7.250000          -0.1853513E-05
  Y( 1)        5.652332            0.000000
  Y( 2)        7.750000          -0.1114154E-05
  C( 1, 1)        3.000000            0.000000
  C( 1, 2)        0.000000            4.008540
  C( 2, 1)        0.000000           0.2051358
  C( 2, 2)        5.000000            0.000000
  C( 3, 1)        4.000000            0.000000
  C( 3, 2)        0.000000            4.487750
  C( 4, 1)        7.000000            0.000000
  C( 4, 2)        0.000000           0.5535090
  C( 5, 1)        6.000000            0.000000
  C( 5, 2)        0.000000            3.544853
  C( 6, 1)        0.000000            4.512336
  C( 6, 2)        11.00000            0.000000

以上解答仅是我的个人意见，如有疑问，欢迎讨论交流！


Comment/May 11,2014

Global Optimal Solution Report by LINGO 14_Extended:
  Global optimal solution found.
  Objective value:                              85.26604
  Objective bound:                            85.26604
  Infeasibilities:                                    0.000000
  Extended solver steps:                        103
  Total solver iterations:                        22555
  Elapsed runtime seconds:                   2.05

  Model Class:                              NLP
  Total variables:                          17
  Nonlinear variables:                    16
  Integer variables:                       0

  Total constraints:                      10
  Nonlinear constraints:                1

  Total nonzeros:                         42
  Nonlinear nonzeros:                  16

  Model Title: LOCATION PROBLEM

                                Variable           Value        Reduced Cost
                                       Z        85.26604            0.000000
                                   X( 1)        3.254883            0.000000
                                   X( 2)        7.250000          -0.1024504E-05
                                   Y( 1)        5.652332            0.000000
                                   Y( 2)        7.750000           0.5438639E-06
                                C( 1, 1)        3.000000            0.000000
                                C( 1, 2)        0.000000            4.008540
                                C( 2, 1)        0.000000           0.2051358
                                C( 2, 2)        5.000000            0.000000
                                C( 3, 1)        4.000000            0.000000
                                C( 3, 2)        0.000000            4.487750
                                C( 4, 1)        7.000000            0.000000
                                C( 4, 2)        0.000000           0.5535090
                                C( 5, 1)        6.000000            0.000000
                                C( 5, 2)        0.000000            3.544853
                                C( 6, 1)        0.000000            4.512336
                                C( 6, 2)        11.00000            0.000000


作者: zengyixuan401    时间: 2013-7-23 10:19
谢谢，我再试试
作者: zengyixuan401    时间: 2013-7-23 10:31
我用你发的lingo11运行正确！！！谢谢啊！！
作者: wujianjack2    时间: 2013-7-23 11:36
zengyixuan401 发表于 2013-7-23 10:31
我用你发的lingo11运行正确！！！谢谢啊！！

哈哈，多谢多谢！共勉！



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