一个小小的LINGO建模技巧

wujianjack2

  本帖将以一个看似简单，实则求解很难的优化测试题来讨论LINGO的求解能力（LINDO API），个人见识有限，敬请谅解！
  参考帖子：
  http://forum.simwe.com/thread-806242-1-1.html     （原问题提出者为shamohu版主，16楼为我的回复）

1.问题陈述：
   参考文献：
   1stOpt and Global Optimization Platform-Comparison and Case Study.caj
   
   
   文献给出参考结果：（取较优解）
   

2.LINGO程序代码：（以LINGO 14及【LINGO学习杂谈 二】中Patch的LINDO API版本为求解引擎）
   方案一：直接编写
  （1）源代码一：
MIN=(X1^X2+X2^X1-5*X1*X2*X3-85)^2+(X1^3*X2^X3*X3^X2-60)^2+(X1^X3+X3^X1-X2-0.55)^2;            @FREE(X1);@FREE(X2);@FREE(X3);
   
   运行结果截图：
   运行状态：
   

   运行结果：
   

（2）源代码二：
X1^X2+X2^X1-5*X1*X2*X3-85=0;X1^3*X2^X3*X3^X2-60=0;
X1^X3+X3^X1-X2-0.55=0;
@FREE(X1);@FREE(X2);@FREE(X3);

  运行结果截图：
  Global Solver Status:
  
   
  Local Nonlinear Solver Status:
  
  
   Local Nonlinear Solver运行结果：
   
  综上，可以看出，如果直接编写，LINGO对于此题的求解结果并不理想，但这并不意味着此题LINGO无法求解，只是，我们需要改变下策略。

   方案二：把复杂的约束简化（前提：把X1,X2,X3的定义域缩小为非负）
   对于方案一的代码，LINDO Systems公司的Mark Wiley先生给出了评价，并给出了合理的建议：
   

   如下代码基本沿用Mark Wiley先生给出的思路，稍作调整：
[C1]  X1^X2+X2^X1-5*X1*X2*X3-85=0;
[LC2] 3*@LOG(X1)+X3*@LOG(X2)+X2*@LOG(X3)-@LOG(60)=0;
[C3]  X1^X3+X3^X1-X2-0.55=0;

   使用Global Solver与否均可较快得到较好的结果：（以下结果为Global Solver求解结果）
   
   有兴趣的可尝试其它版本，包括各种Crack，由于LINDO API版本号的不同，求解效果也可能出现差异。

xiangbeihai

不错不错！！！！！！！

myron_lee

这论坛怎么没人讨论lingo api在第三方高级编译语言的一些知识点啊，难道都是直接在lingo上写然后求解看结果吗？我觉得真正发挥价值，提高系统友好性还得跟业务系统集成整合。

拿着菜刀砍蛆

求体力啊！！！！！！

colinwang

赞，，，，，，，，，，

我叫纪坑坑

我是在限定条件中，老师说lingo没有定义不等于，l例如，要使a*b不等于0在lingo不能直接实现，而如果a*b等于一个极小的数，lingo就认定a*b=0了，为防止这种情况，可以令@abs(a*b)>=10^（-100），转个弯解决不等于0.    不知道这种情况是不是类似。。

我叫纪坑坑

因为解非常非常小吗，今天上课也遇到类似问题，不过不是在求解中，是限定条件。

buaaflylc

线性优化用它的确不错

wangrui2014

学习了，lingo确实强大

屋外的雪

感谢楼主精彩演示，楼主大材啊

z919953051

我想知道咋获得体力值。虽然我看不太懂那个外国先生说的话。。。