model:
sets:
hang/1..3/:n,f;
zhou/1..15/;
zong(hang,zhou):x,fmin,fmax,k,piao;
endsets
data:
n=250 450 500;
fmin=1500 1500 1600 1650 1700 1700 1800 1800 1800 1800 1800 1750 1700 1650 1600
1000 1000 1100 1100 1150 1150 1300 1300 1300 1300 1250 1250 1150 1150 1150
680 680 740 740 780 780 950 950 950 950 950 950 860 860 820;
fmax=1650 1650 1750 1800 1900 1900 2000 2000 2000 2000 2000 1950 1850 1850 1800
1100 1100 1250 1250 1450 1450 1450 1450 1450 1450 1450 1450 1350 1350 1230
750 750 800 800 850 850 990 990 990 990 990 990 900 950 900;
k=0 0 10 3 7 48 113 147 60 30 42 37 19 8 0
7 37 15 23 25 189 34 72 120 150 63 54 32 10 3
11 36 64 25 63 48 34 111 84 18 107 51 28 9 3;
enddata
max=@sum(zong:(fmax-x)/(fmax-fmin)*k*x);
@for(hang(i):@for(zhou(j):x(i,j)<fmax(i,j)));
@for(hang(i):@for(zhou(j):x(i,j)>fmax(i,j)));
@for(hang(j):@for(zhou(i)|i#ge#2:@abs(x(j,i)-x(j,i-1))/x(j,i)<0.2));
@for(hang(j):@sum(zhou(i):(fmax(j,i)-x(j,i))/(fmax(j,i)-fmin(j,i))*k(j,i))<n(j));
@for(hang(j):@for(zhou(i):piao(j,i)=(fmax(j,i)-x(j,i))/(fmax(j,i)-fmin(j,i))*k(j,i)));
@for(zong(i,j):@gin(piao(i,j)));
@for(hang(j):f(j)=@sum(zhou(i):(fmax(j,i)-x(j,i))/(fmax(j,i)-fmin(j,i))*k(j,i)*x(j,i)));
end
sets:
hang/1..3/:n,f;
zhou/1..15/;
zong(hang,zhou):x,fmin,fmax,k,piao;
endsets
data:
n=250 450 500;
fmin=1500 1500 1600 1650 1700 1700 1800 1800 1800 1800 1800 1750 1700 1650 1600
1000 1000 1100 1100 1150 1150 1300 1300 1300 1300 1250 1250 1150 1150 1150
680 680 740 740 780 780 950 950 950 950 950 950 860 860 820;
fmax=1650 1650 1750 1800 1900 1900 2000 2000 2000 2000 2000 1950 1850 1850 1800
1100 1100 1250 1250 1450 1450 1450 1450 1450 1450 1450 1450 1350 1350 1230
750 750 800 800 850 850 990 990 990 990 990 990 900 950 900;
k=0 0 10 3 7 48 113 147 60 30 42 37 19 8 0
7 37 15 23 25 189 34 72 120 150 63 54 32 10 3
11 36 64 25 63 48 34 111 84 18 107 51 28 9 3;
enddata
max=@sum(zong:(fmax-x)/(fmax-fmin)*k*x);
@for(hang(i):@for(zhou(j):x(i,j)<fmax(i,j)));
@for(hang(i):@for(zhou(j):x(i,j)>fmax(i,j)));
@for(hang(j):@for(zhou(i)|i#ge#2:@abs(x(j,i)-x(j,i-1))/x(j,i)<0.2));
@for(hang(j):@sum(zhou(i):(fmax(j,i)-x(j,i))/(fmax(j,i)-fmin(j,i))*k(j,i))<n(j));
@for(hang(j):@for(zhou(i):piao(j,i)=(fmax(j,i)-x(j,i))/(fmax(j,i)-fmin(j,i))*k(j,i)));
@for(zong(i,j):@gin(piao(i,j)));
@for(hang(j):f(j)=@sum(zhou(i):(fmax(j,i)-x(j,i))/(fmax(j,i)-fmin(j,i))*k(j,i)*x(j,i)));
end