sets:
MI/1..17/:qi,Hi,pi,Si,ni;
endsets
data:
cf=20000;
cu=1000;
cp=1;
f=4;
L1=22;
L2=27;
V=60;
h0=0.06;
e=0.01;
m=400;
ef=0.1;
ew=2.0;
qi=1610,2111,2532,2623,2712,2703,3018,3282,3417,3540,4101,4365,3663,3675,2869,2679,1743;
enddata
!目标函数1;
Z1=TR-TC;
TR=f*Q;
TC=Cf+Cu1+Cp1;
Cu1=@sum(MI(i):(L1*cu)/(V*Hi(i)));
Cp1=cp*Q;
!目标2;
@for(MI(i):ni=qi/si);
@for(MI(i):Si=60*x*m/Hi);
z2=@sum(MI(i):ni(i));
!目标;
max=0.5*z1+0.5*z2
!约束条件;
!3;
P=@sum(MI(i):pi(i));
Q=@sum(MI(i):qi(i));
Q=P*(1-ef*f)-0.5*ew*@sum(MI(i):pi(i)*Hi(i));
@for(MI(i):qi(i)=pi(i)*(1-ew*0.5*Hi(i)-ef*f);
!2;
@for(MI(i):(L1*m*x)/(V*Hi(i))>=qi(i));
@for(MI(i):Hi(i)<=(1-ef*f-@sqrt(@sqr(1-ef*f)-(2*ew*x*m*L1)/pi(i)*V))/ew);
!1;
@for(MI(i):Hi(i)>=(h0+e));
!0;
@gin(x);