[求助]LINGO!

eclatma0418

小弟编了一个lingo求解LP的程序,不知为何无法求解,请达人帮忙啊,急啊!!!!!!!!

model:
sets:
year/1..5/:y,s,f,smax,ymin,fmax;
!y(i)-实施措施前预测产油量;
!s(i)-实施措施前预测产水量;
!f(i)-实施措施前预测总费用;
!fmax(i)-第i年总费用上限;
!smax(i)-第i年产水量上限;
!ymin(i)-第i年产油量下限;
cuoshi/1..7/:xx;
!xx(j)-第j项措施总工作量上限;
juzheng(year,cuoshi):x,yy,ss,ff;
!x(i,j)-第i年采用第j项措施的井口数;
!yy(i,j)-第i年第j项措施的单井产油量;
!ss(i,j)-第i年第j项措施的单井产水量;
!ff(i,j)-第i年第j项措施的单井费用;
endsets
data:
y=1377.88 1237.76 1097.64 957.52 817.4;
f=16861 20477 60851 80114 85014;
smax=5706 5843 6152 6313 6720;
s=5346.44 5373.32 5372.97 5346.68 5292.16;
ymin=1500 1490 1485 1450 1423;
fmax=86580 119100 171980 196132 199321;
xx=426 390 430 374 445 534 476;
yy= 0.0820    0.0720    0.0650    0.0540    0.0380    0.0570    0.0450
    0.2160    0.0960    0.0860    0.0862    0.0760    0.0820    0.0750
    0.2679    0.1479    0.1670    0.1509    0.0979    0.1009    0.1074
    0.2304    0.1004    0.1344    0.1228    0.0704    0.0828    0.0893
    0.1981    0.0881    0.0976    0.1006    0.0581    0.0686    0.0551;
ss= 0.0720    0.0630    0.0570    0.0500    0.0320    0.0442    0.0675
    0.2120    0.1880    0.1280    0.1760    0.1790    0.1262    0.1775
    0.3750    0.2050    0.3120    0.3788    0.2131    0.2045    0.3196
    0.5666    0.4086    0.5766    0.4986    0.4253    0.3584    0.5635
    0.7585    0.6785    0.7846    0.6438    0.6653    0.4744    0.7874;
ff=128.1000   78.1000   64.1000   43.0000   58.1000   36.9000   50.5000
   107.0000   60.0000   51.0000   31.0500   49.0000   30.7500   43.8000
   72.2500   52.2500   42.2500   22.2000   37.2500   23.7000   35.7000
   56.5500   39.5500   33.5500   17.0500   30.5500   17.3000   27.3500
   44.7000   30.7000   21.7000   13.1500   24.7000   12.0000   22.0500;

enddata
max=@sum(year(i):y(i)+@sum(cuoshi(j):x(i,j)*yy(i,j)));
!目标一-产油总量最大;
@for(year(i)sum(cuoshi(j):x(i,j)*ss(i,j))+s(i)<smax(i));
!各年产水量不超过上限;
@for(year(i)sum(cuoshi(j):x(i,j)*yy(i,j))+y(i)>ymin(i));
!各年产油量不小于下限;
@for(cuoshi(j):x(5,j)<xx(j));
!各项措施总工作量不超过上限;
@for(year(i)sum(cuoshi(j):x(i,j)*ff(i,j))+f(i)<fmax(i));
!各年总费用不超过上限;
@for(juzheng(i,j)|i #LT# 4:x(i+1,j)>x(i,j));
!采取各措施的井口数逐年递增;
@for(juzheng(i,j)gin(x(i,j)));

秦玉

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

[em06]

bent2008

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

hubulwm

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