数学建模社区-数学中国

标题: 求助，哪位高手帮我看看这道lingo程序哪错了 [打印本页]

作者: zhaoyunyaya 时间: 2010-5-22 17:12
标题: 求助，哪位高手帮我看看这道lingo程序哪错了
sets:
l/1..4/;
m(l,l):t,a,x;
mat/1..4/:b;
endsets
data:
a=825 4726 2436 528
  248 1045 763  183
  562 3263 1509 342
  437 2078 1125 325;
t=0.3 0.2 0.2 0.3
0.3 0.3 0.2 0.2
0.4 0.3 0.1 0.2
0.4 0.2 0.1 0.3;
b=1700 8000 3600 1200；
max=@sum(m(i,j):x(i,j)*t(i,j)/(m(i,j):a(i,j)));
@for(m(i,j)

gin(x(i,j);
@for(@sum(1（j):x(i,j))=b(i);
@for(m(i,j):x(i,j)/(m(i,j):a(i,j)))<1;
@for(m(1,j):x(1,j)/(m(1,j):a(1,j))>0.8;
@for(m(2,j):x(2,j)/(m(2,j):a(2,j))>0.7;
@for(m(3,j):x(i,j)/(m(i,j):a(3,j))>0.6;
@for(m(4,J):x(i,j)/(m(i,j):a(4,j))>0.8;
end
我想把b定义成一个一列四行的矩阵，可定义不了，谢谢大家啊

作者: olh2008 时间: 2010-5-22 21:11
Lingo中并不区分行向量与列向量，所以这个关系并不大。

@for(@sum(1（j):x(i,j))=b(i);5 c! K+ }9 }7

你这句@for怎么没有使用集，还有后面的@sum集怎么是数字1而不是字母L呢？？

作者: zhaoyunyaya 时间: 2010-5-23 08:46
我改成了@for(l(j)

sum(m(i,j):x(i,j))=b(i));可是还不对啊，怎么回事？

作者: zhaoyunyaya 时间: 2010-5-23 08:48
model:
sets:
l/1..4/;
m(l,l):t,a,x;
c/1..4/:b;
endsets
data:
a=825 4726 2436 528
  248 1045 763  183
  562 3263 1509 342
  437 2078 1125 325;
t=0.3 0.2 0.2 0.3
  0.3 0.3 0.2 0.2
  0.4 0.3 0.1 0.2
  0.4 0.2 0.1 0.3;
b=1700 8000 3600 1200;
enddata
max=@sum(m(i,j):x(i,j)*t(i,j)/a(i,j));
@for(m(i,j)

gin(x(i,j)));
@for(l(j)

sum(m(i,j):x(i,j))=b(i));
@for(m(i,j):x(i,j)/a(i,j)< 1);
@for(l(i):x(i,1)/a(i,1)>0.8);
@for(l(i):x(i,2)/a(i,2)>0.7);
@for(l(i):x(i,3)/a(i,3)>0.6);
@for(l(i):x(i,4)/a(i,4)>0.8);
end现在整体上是这样的，运行的时候老说有没定义的指数

作者: zhaoyunyaya 时间: 2010-5-23 08:49
回复 2# olh2008

能再帮我看一下吗，谢谢啊

作者: olh2008 时间: 2010-5-23 09:07
你代码里的类似于

@for(l(j)： @sum(m(i,j):x(i,j))=b(i));

这样的代码是有问题的。
你这里有@for和@sum，但@sum中的集却使用了二维的，还有最奇怪的是外层用@for，集为l(j),等号后面的表达式却是b(i)，这样是没有意义的啊，按理来说应该是b(j)才对。
建议LZ还是吧模型贴上来吧。

作者: zhaoyunyaya 时间: 2010-5-23 09:49
回复 6# olh2008

MAX=0.3/825*X11+0.3/4726*X12+0.2/2436*X13+0.3/528*X14+0.3/248*X21+0.3/1045*X22+0.2/763*X23+0.2/183*X24+0.4/562*X31+0.3/3263*X32+0.1/1509*X33+0.2/342*X34
+0.4/437*X41+0.2/2078*X42+0.1/1125*X43+0.3/325*X44;
X11+X21+X31+X41<1700;
X12+X22+X32+X42<8000;
X13+X23+X33+X43<3600;
X14+X24+X34+X44<1200;
X11/825<1;
X11/825>0.8;
X21/248<1;
X21/248>0.8;
X31/562<1;
X31/562>0.8;
X41/437<1;
X41/437>0.8;
X12/4726<1;
X12/4726>0.7;
X22/1045<1;
X22/1045>0.7;
X32/3263<1;
X32/3263>0.7;
X42/2078<1;
X42/2078>0.7;
X13/2436<1;
X13/2436>0.6;
X23/763<1;
X23/763>0.6;
X33/1509<1;
X33/1509>0.6;
X43/1125<1;
X43/1125>0.6;
X14/528<1;
X14/528>0.8;
X24/183<1;
X24/183>0.8;
X34/342<1;
X34/342>0.8;
X44/325<1;
X44/325>0.8;
@GIN(X11);@GIN(X12);@GIN(X13);@GIN(X14);
@GIN(X21);@GIN(X22);@GIN(X23);@GIN(X24);
@GIN(X31);@GIN(X32);@GIN(X33);@GIN(X34);
@GIN(X41);@GIN(X42);@GIN(X43);@GIN(X44);
END就是根据这个模型编的

作者: zhaoyunyaya 时间: 2010-5-23 10:09
回复 7# zhaoyunyaya
这句的意思就是X11+X21+X31+X41<b（1）;
X12+X22+X32+X42<b（2）;
X13+X23+X33+X43<b(3);
X14+X24+X34+X44<b(4);

作者: zhaoyunyaya 时间: 2010-5-23 10:09
回复 6# olh2008

这句的意思就是X11+X21+X31+X41<b（1）;
X12+X22+X32+X42<b（2）;
X13+X23+X33+X43<b(3);
X14+X24+X34+X44<b(4);谢谢你啊，呵呵

作者: alair006 时间: 2012-1-16 15:37
帮个忙，去下载下这个，我需要体力http://www.madio.net/thread-131123-1-1.html好人啊9504737谢谢啊97987

作者: zerodingying 时间: 2012-2-4 19:16
大神辛苦了。。。。

作者: zerodingying 时间: 2012-2-4 19:17
大神辛苦了。。。。

作者: alair002 时间: 2012-2-5 19:19
写的不错呀，收藏了先…谢谢啦1617437080777979

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