数学建模社区-数学中国

标题: lingo 错误代码11，求改正，谢谢！ [打印本页]

作者: 0061gc 时间: 2013-8-28 11:53
标题: lingo 错误代码11，求改正，谢谢！
本帖最后由 wujianjack2 于 2013-8-28 12:21 编辑

model:
sets:
var/1..10/;
col/1..21/:x,p;
link(var,var):Ml,Md,Ld;
endsets
data:
Ld=0 26 31 17 41 15 20 45 35 19
10 0 29 19 59 54 39 47 17 23
20 28 0 27 36 31 39 12 37 20
10 27 45 0 15 24 30 38 25 25
42 36 21 39 0 33 56 24 44 55
50 53 14 37 31 0 33 45 35 30
31 37 38 49 12 30 0 25 40 36
48 58 27 44 50 11 50 0 20 40
39 51 12 12 38 26 31 36 0 33
46 45 47 19 49 13 22 58 51 0;
enddata
min=@sum(col(i):p*x);
@for(var(i):@sum(var(j):ML(i,j))=1);
@for(var(j):@sum(var(i):ML(i,j))=1);
@for(link(i,j):@bin(ML));
@for(link(i,k)|i#ne#k:Md(i,k)>0);
@for(link(i,k):@for(link(j,l):Md(i,k)=@if((ML(i,j)*ML(k,l))#eq#1,Ld(j,l),Md(i,k)));
p(1)=282*(Md(1,2)+Md(2,3)+Md(3,5)+Md(5,7)+Md(7,8)+Md(8,10));
p(2)=282*(Md(3,4)+Md(4,5)+Md(5,6)+Md(6,7)+Md(7,8)+Md(8,9)+Md(9,10));
p(3)=261*(Md(1,4)+Md(4,10));
p(4)=261*(Md(1,4)+Md(4,8));
p(5)=261*(Md(1,4)+Md(4,9));
p(6)=277*(Md(4,5)+Md(5,6)+Md(6,8));
p(7)=277*(Md(4,5)+Md(5,6)+Md(6,7));
p(8)=277*(Md(4,5)+Md(5,7)+Md(7,9)+Md(9,10));
p(9)=218*(Md(2,3)+Md(3,4)+Md(4,5)+Md(5,6));
p(10)=218*(Md(2,3)+Md(3,4)+Md(4,5)+Md(5,7)+Md(7,8));
p(11)=218*(Md(2,3)+Md(3,4)+Md(4,5)+Md(5,7)+Md(7,10));
p(12)=397*(Md(4,5)+Md(5,7)+Md(7,8)+Md(8,10));
p(13)=397*(Md(4,5)+Md(5,8)+Md(8,9)+Md(9,10));
p(14)=397*(Md(4,5)+Md(5,8)+Md(8,10));
p(15)=451*(Md(1,2)+Md(2,3)+Md(3,4)+Md(4,5)+Md(5,6)+Md(6,7)+Md(7,8)+Md(8,9)+Md(9,10));
p(16)=388*(Md(4,5)+Md(5,7)+Md(7,9));
p(17)=388*(Md(5,6)+Md(6,8)+Md(8,9)+Md(9,10));
p(18)=388*(Md(5,6)+Md(6,7)+Md(7,8)+Md(8,9));
p(19)=389*(Md(8,9)+Md(9,10));
p(20)=389*(Md(6,7)+Md(7,8));
p(21)=389*(Md(3,4)+Md(4,5));
@for(col(i):@bin(x));
x(1)+x(2)=1;
x(3)+x(4)+x(5)=1;
x(6)+x(7)+x(8)=1;
x(9)+x(10)+x(11)=1;
x(12)+x(13)+x(14)=1;
x(15)=1;
x(16)+x(17)+x(18)=1;
x(19)+x(20)+x(21)=1;
end

作者: wujianjack2 时间: 2013-8-28 12:28
楼主你好，很复杂的问题呢！
有一个错误是："@for(link(i,k):@for(link(j,l):Md(i,k)=@if((ML(i,j)*ML(k,l))#eq#1,Ld(j,l),Md(i,k)));"似乎缺少了一个括号，应作：
"@for(link(i,k):@for(link(j,l):Md(i,k)=@if((ML(i,j)*ML(k,l))#eq#1,Ld(j,l),Md(i,k))));"
改正后程序是可以运行的，但是好像比较复杂，一时难以找到可行解，楼主也可以再检查下模型。
个人一点拙见，仅供楼主参考。

作者: 0061gc 时间: 2013-8-28 13:16

wujianjack2 发表于 2013-8-28 12:28
楼主你好，很复杂的问题呢！
有一个错误是："@for(link(i,k)for(link(j,l):Md(i,k)=@if((ML(i,j)* ...

谢谢你的解答，确实少了一个括号，但是改正后，仍然运行不出来啊，，按理应该是可以的有可行解的，或许是哪约束的不全面吧，有错误的那句话，是想表达：当ML(i，j)=1且ML（k，l）=1时，Md（i，k）=Ld（j，l）这个意思，你有什么好的建议吗

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