12 / 2 页下一页

查看: 3047|回复: 12

求助，哪位高手帮我看看这道lingo程序哪错了

[复制链接]

字体大小: 正常放大

zhaoyunyaya

5 主题	2 听众	33 积分

升级 29.47%

TA的每日心情

	开心 2011-9-1 09:55

签到天数: 1 天

[LV.1]初来乍到

自我介绍: 200 字节以内

不支持自定义 Discuz! 代码

群组: 东北三省联盟

电梯直达

1^#

发表于 2010-5-22 17:12 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

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定义成一个一列四行的矩阵，可定义不了，谢谢大家啊

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

olh2008

88 主题	42 听众	1万积分

船长

TA的每日心情

	开心 2018-9-1 14:36

签到天数: 86 天

[LV.6]常住居民II

群组: Latex研学群

群组: 数学建模

群组: Mathematica研究小组

群组: LINGO

群组: Matlab讨论组

2^#

发表于 2010-5-22 21:11 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

Lingo中并不区分行向量与列向量，所以这个关系并不大。

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

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

我也说一句

发表

生命，到最后总能成诗……

使用道具举报

zhaoyunyaya

5 主题	2 听众	33 积分

升级 29.47%

TA的每日心情

	开心 2011-9-1 09:55

签到天数: 1 天

[LV.1]初来乍到

自我介绍: 200 字节以内

不支持自定义 Discuz! 代码

群组: 东北三省联盟

3^#

发表于 2010-5-23 08:46 |只看该作者 |招呼Ta 关注Ta

我改成了@for(l(j)

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

我也说一句

发表

使用道具举报

zhaoyunyaya

5 主题	2 听众	33 积分

升级 29.47%

TA的每日心情

	开心 2011-9-1 09:55

签到天数: 1 天

[LV.1]初来乍到

自我介绍: 200 字节以内

不支持自定义 Discuz! 代码

群组: 东北三省联盟

4^#

发表于 2010-5-23 08:48 |只看该作者 |招呼Ta 关注Ta

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

5 主题	2 听众	33 积分

升级 29.47%

TA的每日心情

	开心 2011-9-1 09:55

签到天数: 1 天

[LV.1]初来乍到

自我介绍: 200 字节以内

不支持自定义 Discuz! 代码

群组: 东北三省联盟

5^#

发表于 2010-5-23 08:49 |只看该作者 |招呼Ta 关注Ta

回复 2# olh2008

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

我也说一句

发表

使用道具举报

olh2008

88 主题	42 听众	1万积分

船长

TA的每日心情

	开心 2018-9-1 14:36

签到天数: 86 天

[LV.6]常住居民II

群组: Latex研学群

群组: 数学建模

群组: Mathematica研究小组

群组: LINGO

群组: Matlab讨论组

6^#

发表于 2010-5-23 09:07 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

你代码里的类似于

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

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

我也说一句

发表

生命，到最后总能成诗……

使用道具举报

zhaoyunyaya

5 主题	2 听众	33 积分

升级 29.47%

TA的每日心情

	开心 2011-9-1 09:55

签到天数: 1 天

[LV.1]初来乍到

自我介绍: 200 字节以内

不支持自定义 Discuz! 代码

群组: 东北三省联盟

7^#

发表于 2010-5-23 09:49 |只看该作者 |招呼Ta 关注Ta

回复 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就是根据这个模型编的