数学建模社区-数学中国
标题:
求助,哪位高手帮我看看这道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/)
Powered by Discuz! X2.5