程序如下:
MODEL:
sets:
num_i/1..12/:y;
num_j/1..2/:a,c;
link_ij(num_i,num_j):x;
endsets
data:
y=7,8,8,12,9,10,7,10,11,12,12,13;
x=1,1,1,2,1,3,1,4,1,5,1,6,1,7,1,8,1,9,1,10,1,11,1,12;
enddata
[OBJ]min=@sum(num_i(i):(@sum(num_j(j):a(j)*x(i,j)+c(j)*@abs(x(i,j))-y(i))^2+(@sum(num_j(j):a(j)*x(i,j)-c(j)*@abs(x(i,j))-y(i))^2));
@for(num_i(i):@sum(num_j(j):a(j)*x(i,j)+c(j)*@abs(x(i,j)))>=y(i));
@for(num_i(i):@sum(num_j(j):a(j)*x(i,j)-c(j)*@abs(x(i,j)))<=y(i));
@for(num_i(i):@sum(num_j(j):c(j)*@abs(x(i,j)))>=0);
END