数学建模社区-数学中国

标题: [求助]高手帮忙看下~感谢! [打印本页]

作者: vera20    时间: 2008-5-26 23:36
标题: [求助]高手帮忙看下~感谢!

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

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=∑∑∑cijxijk +∑∑∑(si-ai)xijk+∑∑∑(si+cij-di)xijk,

约束:jk(xijk)=1

ixijk<=1,

(si+cij)xij,k-1<=skxijk,

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





欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5