[讨论]在Lingo里怎样表示偏微分？

Hisence

&#160;&#160;&#160; 今天我们老师布置了一道题，我想用lingo来编，可是在求某个变量的偏微分的时候卡住了，因为是第一次接触到lingo，比较不清楚，请大虾指点指点。谢了<br/>

包不同

<p>MODEL:<br/>SETS:<br/>venders1/v1,v2/:c1,x1;<br/>goods1/g1,g2/:h1,i1,m1,o1;<br/>vg1(venders1,goods1):q1,d1,z1,g1,b1,p1;<br/>venders2/s1,s2/:c2,x2;<br/>goods2/w1,w2/:h2,i2,m2,o2;<br/>vg2(Venders2,goods2):q2,d2,z2,g2,b2,p2;<br/>ENDSETS<br/>DATA:<br/>c1= 1200,2000;<br/>h1= 1,1.5;<br/>m1= 3500,4000;<br/>q1= 0.04 0.02<br/>&nbsp;&nbsp;&nbsp; 0.03 0.04;<br/>d1= 0.18 0.09<br/>&nbsp;&nbsp;&nbsp; 0.2 0.1;<br/>g1= 2 1.8<br/>&nbsp;&nbsp;&nbsp; 1 0.8;<br/>b1= 3000 5000<br/>&nbsp;&nbsp;&nbsp; 3000 4000;<br/>p1= 5 5.1<br/>&nbsp;&nbsp;&nbsp; 2.2 2.1;<br/>c2= 1200,2000;<br/>h2= 1,1.5;<br/>m2= 3200,3500;<br/>q2= 0.04 0.02<br/>&nbsp;&nbsp;&nbsp; 0.03 0.04;<br/>d2= 0.18 0.09<br/>&nbsp;&nbsp;&nbsp; 0.2 0.1;<br/>G2= 2 1.8<br/>&nbsp;&nbsp;&nbsp; 1 0.8;<br/>b2= 3500 4500<br/>&nbsp;&nbsp;&nbsp; 3000 4000;<br/>p2= 5 5.1<br/>&nbsp;&nbsp;&nbsp; 2.2 2.1;</p><p>o1=100,80;</p><p>o2=100,80;<br/>i1= 0,0;<br/>ENDDATA<br/><a href="mailto:min=@sum(vg1(I,J)1(I,J)*z1(I,J))+@sum(venders1:x1*c1)+@sum(vg1(i,j):g1(i,j)*z1(i,j))+@sum(goods1:h1*i1)+@sum(vg2(i,j)2(i,j)*z2(i,j))+@sum(venders2:x2*c2)+@sum(vg2(i,j):g2(i,j)*z2(i,j)+@sum(goods2:h2*i2" target="_blank"><font color="#000000">min=@sum(vg1(I,J)1(I,J)*z1(I,J))+@sum(venders1:x1*c1)+@sum(vg1(i,j):g1(i,j)*z1(i,j))+@sum(goods1:h1*i1)+@sum(vg2(i,j):P2(i,j)*z2(i,j))+@sum(venders2:x2*c2)+@sum(vg2(i,j):g2(i,j)*z2(i,j)+@sum(goods2:h2*i2</font></a>));<br/>@FOR(vg1(i,j):z1(i,j)&lt;=b1(i,j));<br/>@for(vg2(i,j):z2(i,j)&lt;=b2(i,j));</p><p>@for(goods1(j):i1(j)+@sum(goods1(j):z1(i,j)&gt;=m1(j)+1.28o1(j)));<br/>venders1(i):x1(i)=@if(vg1(i,j):z1(i,j)#gt#0,1,0);</p><p>venders2(i):x2(i)=@if(vg2(i,j):z2(i,j)#gt#0,1,0);</p><p>goods2:i2(j)=goods1:i1(j)+@sum(goods2(j):z2(i,j)-m2(j));<br/>end</p><p>这是一个求两阶段规划问题,其中goods2:i2(j)=goods1:i1(j)+@sum(goods2(j):z2(i,j)-m2(j));是状态转移方程,由于本人刚接触lingo&nbsp;所以用的方法比较笨拙,出现错误也不知如何修改,哪位高手能指点一下吗?谢谢</p><p>也不知道为什么我不能发新帖子</p>