查看: 2119|回复: 4

求助：食品油生产问题，程序报错。

字体大小: 正常放大

23 主题	8 听众	372 积分

升级 24%

TA的每日心情

	开心 2014-10-4 18:19

签到天数: 102 天

[LV.6]常住居民II

自我介绍: 本人是在校学生，即将参加有史以来的第一场国赛，希望能在这里获得帮助。

电梯直达

1^#

发表于 2013-8-7 18:36 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

本帖最后由 wujianjack2 于 2013-8-11 13:57 编辑

程序运行，然后报错，我找了半天还是不知道哪里错了，麻烦各位帮帮眼

model:
title:食品油生产;
sets:
class/v1,v2,o1,o2,o3/:r;
!为第i个月第j种油的采购价格;
month/1,2..6/:monthz,monthr;
link(class,month):l,lw,p,b;
!l1为第i个月第j种油的耗用量，l2为第i个月第j种油的储存量
p为第i个月第j种油的采购量,b为第i个月第j种油的采购价格;
endsets
data:
b = 1100 1200 1300 1100 1150
1300 1300 1100 900 1150
1100 1400 1300 1000 950
1200 1100 1200 1200 1250
1000 1200 1500 1100 1050
900 1000 1400 800 1350;

r = 8.8 6.1 2.0 4.2 5.0;
enddata
@for(month(i):l(i,1)+l(i,2)<=200);
@for(month(i):l(i,3)+l(i,4)+l(i,5)<=250);
@for(month(i):monthr(i)=@sum(class(j):l(i,j)*r(j))/@sum(class(j):l(i,j)));
@for(class(j):500+p(1,j)-l(1,j) = lw(1,j));
@for(link(i,j)|i#ge#2:p(i,j)+lw(i-1,j)-l(i,j) = lw(i,j));
@for(month(i):l(i,6)=500);
@for(link:lw<=1000);
@for(month(i):@free(monthz(i)));
@for(month:monthr>=3);
@for(month:monthr<=6);@for(month(i):
monthz(i)=1800*@sum(class(j):l(i,j))-270*(l(i,1)+l(i,2))-310*(l(i,3)+l(i,4)+l(i,5))-50*@sum(class(j):lw(i,j))-@sum(class(j):b(i,j)*p(i,j)));

[obj] max = @sum(month(i):monthz(i));
end

zan

lingo, 运行

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

相关帖子

使用道具举报

madio

3万主题	1311 听众	5万积分

TA的每日心情

	奋斗 2024-7-1 22:21

签到天数: 2014 天

[LV.Master]伴坛终老

自我介绍: 数学中国站长

群组: 数学建模培训课堂1

群组: 数学中国美赛辅助报名

群组: Matlab讨论组

群组: 2013认证赛A题讨论群组

群组: 2013认证赛C题讨论群组

2^#

发表于 2013-8-8 10:09 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

给你改了一下，主要是集定义和使用的问题，出现了下标越界，你看看吧！

model:
title:食品油生产;
sets:
class/v1,v2,o1,o2,o3/:r;
!为第i个月第j种油的采购价格;
month/1,2..6/:monthz,monthr;
link(month,class):l,lw,p,b;
!l1为第i个月第j种油的耗用量，l2为第i个月第j种油的储存量
p为第i个月第j种油的采购量,b为第i个月第j种油的采购价格;
endsets
data:
b = 1100 1200 1300 1100 1150
1300 1300 1100 900 1150
1100 1400 1300 1000 950
1200 1100 1200 1200 1250
1000 1200 1500 1100 1050
900 1000 1400 800 1350;
r = 8.8 6.1 2.0 4.2 5.0;
enddata
[obj] max = @sum(month(i):monthz(i));
@for(month(i):l(i,1)+l(i,2)<=200);
@for(month(i):l(i,3)+l(i,4)+l(i,5)<=250);
@for(month(i):monthr(i)=@sum(class(j):l(i,j)*r(j))/@sum(class(j):l(i,j)));
@for(class(j):500+p(1,j)-l(1,j) = lw(1,j));
@for(link(i,j)|i#ge#2:p(i,j)+lw(i-1,j)-l(i,j) = lw(i,j));
@for(month(i):l(i,5)=500);
@for(link:lw<=1000);
@for(month(i):@free(monthz(i)));
@for(month:monthr>=3);
@for(month:monthr<=6);@for(month(i):
monthz(i)=1800*@sum(class(j):l(i,j))-270*(l(i,1)+l(i,2))-310*(l(i,3)+l(i,4)+l(i,5))-50*@sum(class(j):lw(i,j))-@sum(class(j):b(i,j)*p(i,j)));
end