目标规划lingo

!求第一级目标;
model:

sets:
value/1..2/:x;  !x为变量;
num/1..4/:g,dz,df; !g为柔性约束等号右侧目标值，dz为正偏差变量，df为负偏差变量;
cost(num,value):c;   !c为柔性约束系数值;
endsets

data:
g = 1500 0 16 15;
c = 200 300
    2 -1 
    4 0
    0 5;
enddata

min = df(1);
2*x(1) + 2*x(2) <=12;!刚性约束;
@for(num(i):@sum(value(j):c(i,j)*x(j))+df(i)-dz(i)=g(i));!柔性约束;
@for(value(i):@gin(x(i)));!限制x为整数;
end

!求第二级目标;
model:

sets:
value/1..2/:x;  !x为变量;
num/1..4/:g,dz,df; !g为柔性约束等号右侧目标值，dz为正偏差变量，df为负偏差变量;
cost(num,value):c;   !c为柔性约束系数值;
endsets

data:
g = 1500 0 16 15;
c = 200 300
    2 -1 
    4 0
    0 5;
enddata

min = df(2)+dz(2);
2*x(1) + 2*x(2) <=12;!刚性约束;
@for(num(i):@sum(value(j):c(i,j)*x(j))+df(i)-dz(i)=g(i));!柔性约束;
df(1)=0;!一级目标约束;
@for(value(i):@gin(x(i))); !限制x为整数;
end

!求第三级目标;
model:

sets:
value/1..2/:x;  !x为变量;
num/1..4/:g,dz,df; !g为柔性约束等号右侧目标值，dz为正偏差变量，df为负偏差变量;
cost(num,value):c;   !c为柔性约束系数值;
endsets

data:
g = 1500 0 16 15;
c = 200 300
    2 -1 
    4 0
    0 5;
enddata

min = 3*dz(3)+3*df(3)+dz(4);
2*x(1) + 2*x(2) <=12;!刚性约束;
@for(num(i):@sum(value(j):c(i,j)*x(j))+df(i)-dz(i)=g(i));!柔性约束;
df(1)=0;!一级目标约束;
df(2)+dz(2)=0;!二级目标约束;
@for(value(i):@gin(x(i))); !限制x为整数;
end

结果：

生产甲2件，乙4件；
最大利润为：2200+4300=1600元