出版一本关于LINGO8.0的书，不知大家会买否？

wanbaocheng

<>我已贴到网上一部分了</P>
<><a href="http://www.shumo.com/bbs/dispbbs.asp?boardid=10&amp;ID=7115&amp;replyID=48081" target="_blank" >http://www.shumo.com/bbs/dispbbs.asp?boardid=10&amp;ID=7115&amp;replyID=48081</A></P>

legence

随着建模软件的应用日益广泛，会对不同的建模软件有需求，建议还是早点处为好。

legence

<>MODEL:
   TITLE LOCATION PROBLEMS;</P>
<>   SETS:
   DEMAND/1..6/:A,B,D;
   SUPPLY/1..2/:X,Y,E;
   LINK/DEMAND,SUPPLY/:C;
   ENDSETS</P>
<>   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  </P>
<P>  INIT:
   X,Y=5,1,2,7;
   ENDINIT
   !objective function;

   [OBJ] MIN = @SUM(LINK(I,J):C(I,J) * ((X(J)-A(I))^2 + (Y(J)-B(I))^2)^(1/2) );
   !demand constraints;
   @FOR(DEMAND(I):[DEMAND_CON] @SUM(SUPPLY(J):C(I,J)) = D(I););
   !supply contraints;
   @FOR(SUPPLY(I):[SUPPLY_CON] @SUM(DEMAND{J):C(J,I) &lt;= E(I););
   @FOR(SUPPLY: @FREE(X); @FREE(Y););
END</P>
<P>编译时出现以下错误：</P>
<IMG src="http://qflt.512j.com/upload/forum7_f_2_0_1098778763.jpg" border=0>

wanbaocheng

<><FONT face=宋体>MODEL:
   TITLE LOCATION PROBLEMS;</FONT></P>
<><FONT face=宋体>   SETS:
   DEMAND/1..6/:A,B,D;
   SUPPLY/1..2/:X,Y,E;
   LINK<FONT color=#ff0000>(</FONT>DEMAND,SUPPLY<FONT color=#f70909>)</FONT>:C;
   ENDSETS</FONT></P>
<><FONT face=宋体>   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  </FONT></P>
<P><FONT face=宋体>  INIT:
   X,Y=5,1,2,7;
   ENDINIT
   !objective function;</FONT></P>
<P><FONT face=宋体>   [OBJ] MIN = @SUM(LINK(I,J):C(I,J) * ((X(J)-A(I))^2 + (Y(J)-B(I))^2)^(1/2) );
   !demand constraints;
   @FOR(DEMAND(I):[DEMAND_CON] @SUM(SUPPLY(J):C(I,J)) = D(I););
   !supply contraints;
   @FOR(SUPPLY(I):[SUPPLY_CON] @SUM(DEMAND<FONT color=#ff0000>(</FONT>J):C(J,I)<FONT color=#ff3300>)</FONT> &lt;= E(I););
   @FOR(SUPPLY: @FREE(X); @FREE(Y););
END</FONT></P>
<P><FONT face=宋体>红色字符部分是修改或加的，我已用LINGO8.0运行通过。有什么问题在联系。</FONT></P>

legence

<>为什么</P><>LINK/DEMAND,SUPPLY/:C;不对呢，我见很多资料都是这么写的，请指教，谢谢</P><>LINK<FONT color=#ff0000>(</FONT>DEMAND,SUPPLY<FONT color=#f70909>)</FONT>:C;就对呀？</P>

legence

谢谢，知道我错在那里了！

敏夫

我肯定会买

legence

<>总是出问题，请问斑竹这是什么问题</P>
<>MODEL:</P>
<>   sets:
   roter/1..10/:r,q1r,q2r;
   endsets</P>
<P>   data:
   r = 1,2,3,4,5,6,7,8,9,10;
   p = 2000;u = 0.27;w = 6283.2;
   k = 2.198;m1 = 1.198;m2 = -3.198;
   enddata</P>
<P>   init:
   a,b = 0.075,0.112;
   endinit</P>
<P>   q1a = (a^2) * p * (w^2) * (3+u)/(2 * k * (-3 + k));
   q1b = (b^2) * p * (w^2) * (3+u)/(2 * k * (-3 + k));
   q2a = (a^2) * p * (w^2) * (3+u)/(2 * k * (-3 - k));
   q2b = (b^2) * p * (w^2) * (3+u)/(2 * k * (-3 - k));
   c1 = ((a^(m1+1)) *(q2a - q1a) - (b^(1+k)) * (q2b - q1b))/(a^(2*k) - b^(2*k));
   c2 = (a^k) * (b^k) * (-a * (b^k) * (q2a - q1a) + b * (a^k) * (q2b - q1b));
  
   @for(roter(i):q1r(i) = ((a+(b-a)*r(i)/10)^2)*p*(w^2)*(3+u)/(2 * k * (-3+k)););   
   @for(roter(i):q2r(i) = ((a+(b-a)*r(i)/10)^2)*p*(w^2)*(3+u)/(2 * k * (-3-k)););
  
   @max(roter(i):c1 *(m1+1)* ((a+(b-a)*r(i)/10)^m1) + c2 *(m2+1)* ((a+(b-a)*r(i)/10)^m2) + (m1+1) * q1r(i)
        - (m2+1) * q2r(i) + p*(w^2)*(a+(b-a)* r(i)/10)) &lt;= 1.28e+9;
   @max(roter(i):c1 * ((a+(b-a)*r(i)/10)^m1) + c2 * ((a+(b-a)*r(i)/10)^m2) + q1r(i) + q2r(i)) &lt;= 1.58e+8;</P>
<P>END</P>

[此贴子已经被作者于2004-10-27 21:43:30编辑过]

legence

怎么没有人解答呢

wanbaocheng

你的问题不是语法错误，而是算法在计算中遇到像0作分母等无意义的情况，对你的优化问题不清楚，故不能乱改。