向大神求教！用lingo12.0出错但用13.0可以，都是打开了global solve后出现的结果

zengyixuan401

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

楼主的问题不错！是在看谢金星，薛毅老师的《优化建模与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

谢谢，我再试试

zengyixuan401

我用你发的lingo11运行正确！！！谢谢啊！！

wujianjack2

zengyixuan401 发表于 2013-7-23 10:31
我用你发的lingo11运行正确！！！谢谢啊！！

哈哈，多谢多谢！共勉！