数学建模社区-数学中国

标题: 跪求大神帮忙看看这个lingo程序问题在哪儿 [打印本页]

作者: xiangbeihai 时间: 2017-3-10 13:42
标题: 跪求大神帮忙看看这个lingo程序问题在哪儿
model:
sets:
!1发点5收点4种道路类型;
notes/1..6/:cap,dem;
road/1,2,3,4/:p,a,v;
links(notes,notes,road):q,d,vel,x;
endsets
!目标函数;
submodel obj1:
min=z1;
z1=@sum(links:q*d);
endsubmodel
submodel obj2:
max=z2;
z2=@sum(links:vel*x);
endsubmodel
submodel obj3:
min=z3;
z3=@sum(links:p*x);
endsubmodel
!需求约束;
@SUM(notes(i):@SUM(notes(j):@sum(road(r):x(i,j,r))))=1;!4;
@for(notes(l):DEM(l)=@sum(links|i#ne#l:Q(i,l,r))-@sum(links|j#ne#l:Q(j,l,r))); !6;
@for(links:@sum(links:Q(i,j,r))>@sum(links:DEM(j)*x(i,j,r)));!7;
@for(links:@sum(links:Q(i,j,r))<@sum(links:CAP*x(i,j,r)-DEM(i)*x(i,j,r)));!8;
@for(links:VEL(i,j,r)<V(r)*x(i,j,r));!9;
@for(links:Pr=VEL(i,j,r)*a(r));!10;
@FOR(links(m,l,k):@sum(links(m,l,k)|m#ne#l:x(m,l,k))=@sum(links(l,n,k)|n#ne#l:x(1,n,k)));!5;
@bin(x(i,j,r));
@for(links:d(i,j,r)=@rand(seed)*50+10);
@for(links:VEL(i,j,r)=@rand(seed)*120+5);
data:
DEM=0  1000  1002  1300  1200  1007;
CAP=1000;
a=0.25  0.1  0.05  0.15;
V=30  55  80  120;
enddata
end

作者: xiangbeihai 时间: 2017-3-10 13:45
谁能排查成功，找我微信18066817737，鄙人悬赏100+ RMB！！！

作者: xiangbeihai 时间: 2017-3-10 13:46
谁能排查成功，找我微信18066817737，鄙人悬赏100+ RMB！！！

作者: TGFzhou 时间: 2017-3-11 13:25

xiangbeihai 发表于 2017-3-10 13:45
谁能排查成功，找我微信18066817737，鄙人悬赏100+ RMB！！！

我修改后运行了一下，但好像没有可行解

作者: fox 时间: 2017-3-14 09:26
重赏之下必有勇夫啊！顶起来

作者: 数字中国meme 时间: 2017-3-16 19:01
表示本白白还欠火候哈

作者: 我是大白菜 时间: 2020-2-16 09:46
谢谢分享。，，，，

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)