[求助]高手帮忙看下~感谢!

vera20

在同学的帮助下,程序改进了不少,但约束部分还是不懂,程序还是运行不出来~~~只有继续请教高手啦!

MODEL:
sets:
ship/1..11/:s,a,d,l,m;
berth/1..4/:r;
time:;
link(ship,berth):c;
contact(ship,berth,time):x;
endsets
data:
m=428 45 259 172 684 356 435 350 150 250 333;
a=9 9 0.5 21 0.5 8.5 7 11.5 21.5 22 9;
l=20 21 13 23.83 23.83 21 20.5 23.83 23.83 23.83 23.83;
r=90 60 30 30;
enddata
!Objective function;
[OBJ]min=@sum(contact(i,j,k):c(i,j)*x(i,j,k))+@sum(contact(i,j,k)s(i)-a(i))*x(i,j,k)+@sum(contact(i,j,k)s(i)+c(i,j)-d(i))*x(i,j,k));

@for(link(i,j):c(i,j)*r(j)=m(i));
@sum(link(i,j):x(i,j,k))=1;
@for(ship(i)sum(ship(i):x(i,j,k));
@for(link(i,j,k):s(i)*x(i,k-1)+c(i,j)*x(i,j,k-1)<s(k)*x(i,j,k);
@for(k<11,k>0);
end

数学模型:Min z=∑∑∑c_ijx_ijk +∑∑∑(s_i-a_i)x_ijk+∑∑∑(s_i+c_ij-d_i)x_ijk,

_约束:∑_j∑_k(x_ijk)=1

∑_ix_ijk<=1,

(s_i+c_ij)x_ij,k-1<=s_kx_ijk,

[此贴子已经被作者于2008-5-26 23:39:14编辑过]