eclatma0418 发表于 2007-8-18 16:38

[求助]LINGO!

<p>小弟编了一个lingo求解LP的程序,不知为何无法求解,请达人帮忙啊,急啊!!!!!!!!</p><p>model:<br/>sets:<br/>year/1..5/:y,s,f,smax,ymin,fmax;<br/>!y(i)-实施措施前预测产油量;<br/>!s(i)-实施措施前预测产水量;<br/>!f(i)-实施措施前预测总费用;<br/>!fmax(i)-第i年总费用上限;<br/>!smax(i)-第i年产水量上限;<br/>!ymin(i)-第i年产油量下限;<br/>cuoshi/1..7/:xx;<br/>!xx(j)-第j项措施总工作量上限;<br/>juzheng(year,cuoshi):x,yy,ss,ff;<br/>!x(i,j)-第i年采用第j项措施的井口数;<br/>!yy(i,j)-第i年第j项措施的单井产油量;<br/>!ss(i,j)-第i年第j项措施的单井产水量;<br/>!ff(i,j)-第i年第j项措施的单井费用;<br/>endsets<br/>data:<br/>y=1377.88 1237.76 1097.64 957.52 817.4;<br/>f=16861 20477 60851 80114 85014;<br/>smax=5706 5843 6152 6313 6720;<br/>s=5346.44 5373.32 5372.97 5346.68 5292.16;<br/>ymin=1500 1490 1485 1450 1423;<br/>fmax=86580 119100 171980 196132 199321;<br/>xx=426 390 430 374 445 534 476;<br/>yy= 0.0820&nbsp;&nbsp;&nbsp; 0.0720&nbsp;&nbsp;&nbsp; 0.0650&nbsp;&nbsp;&nbsp; 0.0540&nbsp;&nbsp;&nbsp; 0.0380&nbsp;&nbsp;&nbsp; 0.0570&nbsp;&nbsp;&nbsp; 0.0450<br/>&nbsp;&nbsp;&nbsp; 0.2160&nbsp;&nbsp;&nbsp; 0.0960&nbsp;&nbsp;&nbsp; 0.0860&nbsp;&nbsp;&nbsp; 0.0862&nbsp;&nbsp;&nbsp; 0.0760&nbsp;&nbsp;&nbsp; 0.0820&nbsp;&nbsp;&nbsp; 0.0750<br/>&nbsp;&nbsp;&nbsp; 0.2679&nbsp;&nbsp;&nbsp; 0.1479&nbsp;&nbsp;&nbsp; 0.1670&nbsp;&nbsp;&nbsp; 0.1509&nbsp;&nbsp;&nbsp; 0.0979&nbsp;&nbsp;&nbsp; 0.1009&nbsp;&nbsp;&nbsp; 0.1074<br/>&nbsp;&nbsp;&nbsp; 0.2304&nbsp;&nbsp;&nbsp; 0.1004&nbsp;&nbsp;&nbsp; 0.1344&nbsp;&nbsp;&nbsp; 0.1228&nbsp;&nbsp;&nbsp; 0.0704&nbsp;&nbsp;&nbsp; 0.0828&nbsp;&nbsp;&nbsp; 0.0893<br/>&nbsp;&nbsp;&nbsp; 0.1981&nbsp;&nbsp;&nbsp; 0.0881&nbsp;&nbsp;&nbsp; 0.0976&nbsp;&nbsp;&nbsp; 0.1006&nbsp;&nbsp;&nbsp; 0.0581&nbsp;&nbsp;&nbsp; 0.0686&nbsp;&nbsp;&nbsp; 0.0551;<br/>ss= 0.0720&nbsp;&nbsp;&nbsp; 0.0630&nbsp;&nbsp;&nbsp; 0.0570&nbsp;&nbsp;&nbsp; 0.0500&nbsp;&nbsp;&nbsp; 0.0320&nbsp;&nbsp;&nbsp; 0.0442&nbsp;&nbsp;&nbsp; 0.0675<br/>&nbsp;&nbsp;&nbsp; 0.2120&nbsp;&nbsp;&nbsp; 0.1880&nbsp;&nbsp;&nbsp; 0.1280&nbsp;&nbsp;&nbsp; 0.1760&nbsp;&nbsp;&nbsp; 0.1790&nbsp;&nbsp;&nbsp; 0.1262&nbsp;&nbsp;&nbsp; 0.1775<br/>&nbsp;&nbsp;&nbsp; 0.3750&nbsp;&nbsp;&nbsp; 0.2050&nbsp;&nbsp;&nbsp; 0.3120&nbsp;&nbsp;&nbsp; 0.3788&nbsp;&nbsp;&nbsp; 0.2131&nbsp;&nbsp;&nbsp; 0.2045&nbsp;&nbsp;&nbsp; 0.3196<br/>&nbsp;&nbsp;&nbsp; 0.5666&nbsp;&nbsp;&nbsp; 0.4086&nbsp;&nbsp;&nbsp; 0.5766&nbsp;&nbsp;&nbsp; 0.4986&nbsp;&nbsp;&nbsp; 0.4253&nbsp;&nbsp;&nbsp; 0.3584&nbsp;&nbsp;&nbsp; 0.5635<br/>&nbsp;&nbsp;&nbsp; 0.7585&nbsp;&nbsp;&nbsp; 0.6785&nbsp;&nbsp;&nbsp; 0.7846&nbsp;&nbsp;&nbsp; 0.6438&nbsp;&nbsp;&nbsp; 0.6653&nbsp;&nbsp;&nbsp; 0.4744&nbsp;&nbsp;&nbsp; 0.7874;<br/>ff=128.1000&nbsp;&nbsp; 78.1000&nbsp;&nbsp; 64.1000&nbsp;&nbsp; 43.0000&nbsp;&nbsp; 58.1000&nbsp;&nbsp; 36.9000&nbsp;&nbsp; 50.5000<br/>&nbsp;&nbsp; 107.0000&nbsp;&nbsp; 60.0000&nbsp;&nbsp; 51.0000&nbsp;&nbsp; 31.0500&nbsp;&nbsp; 49.0000&nbsp;&nbsp; 30.7500&nbsp;&nbsp; 43.8000<br/>&nbsp;&nbsp; 72.2500&nbsp;&nbsp; 52.2500&nbsp;&nbsp; 42.2500&nbsp;&nbsp; 22.2000&nbsp;&nbsp; 37.2500&nbsp;&nbsp; 23.7000&nbsp;&nbsp; 35.7000<br/>&nbsp;&nbsp; 56.5500&nbsp;&nbsp; 39.5500&nbsp;&nbsp; 33.5500&nbsp;&nbsp; 17.0500&nbsp;&nbsp; 30.5500&nbsp;&nbsp; 17.3000&nbsp;&nbsp; 27.3500<br/>&nbsp;&nbsp; 44.7000&nbsp;&nbsp; 30.7000&nbsp;&nbsp; 21.7000&nbsp;&nbsp; 13.1500&nbsp;&nbsp; 24.7000&nbsp;&nbsp; 12.0000&nbsp;&nbsp; 22.0500;</p><p>enddata<br/><a href="mailto:max=@sum(year(i):y(i)+@sum(cuoshi(j):x(i,j)*yy(i,j">max=@sum(year(i):y(i)+@sum(cuoshi(j):x(i,j)*yy(i,j</a>)));<br/>!目标一-产油总量最大;<br/>@for(year(i):@sum(cuoshi(j):x(i,j)*ss(i,j))+s(i)&lt;smax(i));<br/>!各年产水量不超过上限;<br/>@for(year(i):@sum(cuoshi(j):x(i,j)*yy(i,j))+y(i)&gt;ymin(i));<br/>!各年产油量不小于下限;<br/>@for(cuoshi(j):x(5,j)&lt;xx(j));<br/>!各项措施总工作量不超过上限;<br/>@for(year(i):@sum(cuoshi(j):x(i,j)*ff(i,j))+f(i)&lt;fmax(i));<br/>!各年总费用不超过上限;<br/>@for(juzheng(i,j)|i #LT# 4:x(i+1,j)&gt;x(i,j));<br/>!采取各措施的井口数逐年递增;<br/>@for(juzheng(i,j):@gin(x(i,j)));</p>

秦玉 发表于 2007-9-18 10:45

<p>我们前几天也在做这个题,没有解出来....帮不了你了</p>

bent2008 发表于 2007-10-20 19:11

暑假培训做过,不过都忘了,实在帮不了你<a href="http://www.abab123.com/bbs/down.asp?html=707637">http://www.abab123.com/bbs/down.asp?html=707637</a>

hubulwm 发表于 2009-5-21 21:09

你检查下这个语句代表的条件是否与其他约束冲突,因为把它去掉后有全局最优解.
@for(year(i):@sum(cuoshi(j):x(i,j)*yy(i,j))+y(i)>ymin(i));!各年产油量不小于下限;
页: [1]
查看完整版本: [求助]LINGO!