lingo一个奇怪的问题_跪求帮忙修改代码

Catnan

求高手帮忙修改代码，能让他正常运行
错误好像主要在@FOR(ZHILU(I):FR(I)=(FI(I)-FX(I))/FX(I));这句话附近
显示的错误如下：
Error 72
unable to solve for fixed variable
FR(ZHILU1)

in constraint
20

loosening  the variable's bounds may help
谢谢耐心解答

1. MODEL:
   
2. SETS:
   
3. JIZHU/JIZHU1..JIZHU8/:GI0,GI1,V;
   
4. DUAN/DUAN1..DUAN10/:D;
   
5. ZHILU/ZHILU1..ZHILU6/:CC,FI,FX,FD,FR;
   
6. LINKS(JIZHU,DUAN):DJ,DR;
   
7. LINK(ZHILU,JIZHU):C;
   
8. ENDSETS
   
9. DATA:
   
10. GI0=120,73,180,80,125,125,81.1,90;
    
11. V=2.2,1,3.2,1.3,1.8,2,1.4,1.8;
    
12. CC=1.111224786317819   1.318676323930283  -1.096597047551148   0.783854251638765   1.338126254232463   1.220559749861743;
    
13. FX=165     150      160     155     132     162;
    
14. FD= 186.4500  177.0000  174.4000  172.0500  151.8000  184.6800;
    
15. DR=
    
16. 70  0   50  0   0   30  0   0   0 40
    
17. 30  0   20  8   15  6   2   0   0   8
    
18. 110 0   40  0   30  0   20  40  0   40
    
19. 55  5   10  10  10  10  15  0   0   1
    
20. 75  5   15  0   15  15  0   10  10  10
    
21. 95  0   10  20  0   15  10  20  0   10
    
22. 50  15  5   15  10  10  5   10  3   2
    
23. 70  0   20  0   20  0   20  10  15  5;
    
24. DJ=
    
25. -505    0   124 168 210 252 312 330 363 489
    
26. -560    0   182 203 245 300 320 360 410 495
    
27. -610    0   152 189 233 258 308 356 415 500
    
28. -500    150 170 200 255 302 325 380 435 800
    
29. -590    0   116 146 188 215 250 310 396 510
    
30. -607    0   159 173 205 252 305 380 405 520
    
31. -500    120 180 251 260 306 315 335 348 548
    
32. -800    153 183 233 253 283 303 318 400 800;
    
33. 
    
34. C= 0.0814    0.0482    0.0523    0.1197   -0.0277    0.1187    0.1259   -0.0021
    
35.    -0.0550    0.1264   -0.0009    0.0331    0.0860   -0.1149   -0.0170    0.0988
    
36.    -0.0680    0.0627   -0.1565   -0.0099    0.1258    0.0040   -0.0039   -0.1994
    
37.    -0.0355   -0.1046    0.2048   -0.0207   -0.0140    0.0042    0.1472    0.0743
    
38.     0.0003    0.2408   -0.0662   -0.0414   -0.0662    0.0689   -0.0021   -0.0104
    
39.     0.2349   -0.0619   -0.0790    0.0930    0.0447   -0.0028    0.1677   -0.0011;
    
40. ENDDATA
    
41. MIN=@SUM(ZHILU(I):(FI(I)-FX(I))/FX(I));
    
42. 
    
43. @FOR(JIZHU(I):GI1(I)<=GI0(I)+GI0(I)*V(I));
    
44. @FOR(JIZHU(I):GI1(I)<=@SUM(DUAN(J):DR(I,J)));
    
45. @FOR(JIZHU(I):GI1(I)>=GI0(I)-GI0(I)*V(I));
    
46. @SUM(JIZHU(I):GI1(I))=1052.8;
    
47. 
    
48. @FOR(ZHILU(I):FI(I)=@SUM(JIZHU(J):C(I,J)*GI1(J)+CC(I)));
    
49. @FOR(ZHILU(I):FI(I)=@ABS(FI(I)));
    
50. @FOR(ZHILU(I):FI(I)<=FD(I));
    
51. !@FOR(ZHILU(I):FR(I)=(FI(I)-FX(I))/FX(I));
    
52. 
    
53. END

复制代码

madio

FR这个集的属性没有赋值，最后一个约束是判断FR(I)=(FI(I)-FX(I))/FX(I))一个等式是否成立，没有值也无法判断相等，这显然不符合要求。

Catnan

madio 发表于 2013-8-22 15:16
FR这个集的属性没有赋值，最后一个约束是判断FR(I)=(FI(I)-FX(I))/FX(I))一个等式是否成立，没有值也无法判 ...

FR(I)=(FI(I)-FX(I))/FX(I))这个地方我想表达的是：由自变量gi1得到变量fi   然后再通过FI得到FR。想单纯表达这么一个计算的过程  所以FR肯定没有值，FI也肯定没有值都是要通过lingo要求出的，所以我应该怎么写这种单纯的计算过程的式子呢？

lingogo

你这个问题根据错误提示,应该是约束条件里面设定了固定变量约束了,仔细分析一下,应该是第46行,@SUM(JIZHU(I):GI1(I))=1052.8;这里,删除这个约束即可得到最优值 -5.687719
Global optimal solution found.
  Objective value:                             -5.687719
  Objective bound:                             -5.687719
  Infeasibilities:                              0.000000
  Extended solver steps:                               0
  Total solver iterations:                             7

lingogo

进一步分析,如果,@SUM(JIZHU(I):GI1(I))=1052.8这个约束是必须的话,那么可以尝试,将这个条件换成,@SUM(JIZHU(I):GI1(I))>=1052.8和,@SUM(JIZHU(I):GI1(I))<=1052.8两个条件分别计算,,@SUM(JIZHU(I):GI1(I))>=1052.8这个条件依然报错,而@SUM(JIZHU(I):GI1(I))<=1052.8时,可以得到最优解,而且和删除@SUM(JIZHU(I):GI1(I))=1052.8约束得到的最优解一样,可以初步认定在@SUM(JIZHU(I):GI1(I))=1052.8约束下无法得到最优解,以上是我的一个简单分析,对此问题有更好处理的可以进一步交流

wujianjack2

lingogo 发表于 2013-8-22 19:29
进一步分析,如果,@SUM(JIZHU(I):GI1(I))=1052.8这个约束是必须的话,那么可以尝试,将这个条件换成,@SUM(JIZH ...

lingogo，你好！很感谢你的积极发言讨论，此题我已和楼主私聊，后来楼主自己发现了错误并改正，得到了自己想要的结果，欢迎以后继续积极发言讨论，视情况我将给你奖励！