新人求解惑,LINGO求解时有个约束条件用不上,在线等!
这是代码:model:
sets:
t/1..50/:v;
endsets
data:
l=5;
w=2;
c=50;
a=20;
b=40;
va=0.64;
vb=0.25;
!x=3.013;
enddata
max=3600*x*(1-(@pow(x,c)*@pow(l,c)*@prod(t(i):1/(i*v(i))))/(1+@sum(t(i):@pow(x,i)*@pow(l,i)*@prod(t(j)|j#le#i:1/(j*v(j))))));
@bnd(0,x,c);!问题出在这,x的范围用不上,还提示解不可行,最优解应该是3左右;
v(1)=1.5;
r=@log(@log(va/v(1))/@log(vb/v(1)))/@log((a-1)/(b-1));
beta=(a-1)/@pow(@log(v(1)/v(a)),1/r);
@for(t(i):v(i)=v(1)*@exp(-1*@pow((i-1)/beta,r)));
end
占个沙发,这是非线性规划,请懂的人进来
页:
[1]