12 / 2 页下一页

查看: 7643|回复: 11

[问题求助] help！！！！！

[复制链接]

字体大小: 正常放大

nuaahappy

473 主题	6 听众	61 积分

升级 58.95%

该用户从未签到

自我介绍: 学生

电梯直达

1^#

发表于 2010-6-9 15:29 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

各位xdjm们，我用lingo编的程序老是报错，但是如果把红色标记的数字改小（如100）就能得出结果了，有谁能帮我解决下啊，很紧急啊，不胜感谢！！！
下面是程序：
model:
max=a1+a2+b1+b2+c1+c2+d1+d2;
38*a1+950*b1+7.8*c1+195*d1<=3537830-yu1;
35*a2+900*b2+7.2*c2+190*d2<=3537830-yu2;
yu1=1238*x1+1393*y1+1548*z1;
yu2=1518*x2+1651*y2+1783*z2;
a1+a2>=7500;
b1+b2>=3300;
c1+c2>=8000;
d1+d2>=1860;
31.2*x1+43.2*y1+55.2*z1>=0.098*a1+0.06*b1+0.0245*c1+0.015*d1;
8000*x1+10000*y1+12000*z1>=1.5*a1+15*b1+2.5*c1+50*d1;
31.2*x2+43.2*y2+55.2*z2>=0.098*a2+0.06*b2+0.0245*c2+0.015*d2;
8000*x2+10000*y2+12000*z2>=1.5*a2+15*b2+2.5*c2+50*d2;
@gin(x1);@gin(y1);@gin(z1);@gin(x2);@gin(y2);@gin(z2);
@gin(a1);@gin(b1);@gin(c1);@gin(d1);
@gin(a2);@gin(b2);@gin(c2);@gin(d2);

zan

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

使用道具举报

nuaahappy

473 主题	6 听众	61 积分

升级 58.95%

该用户从未签到

自我介绍: 学生

2^#

发表于 2010-6-9 15:51 |只看该作者 |招呼Ta 关注Ta

对不起，一激动发错了，应该是
model:
max=a1+a2+b1+b2+c1+c2+d1+d2;
38*a1+950*b1+7.8*c1+195*d1+35*a2+900*b2+7.2*c2+190*d2<=3537830-yu1-yu2;
yu1=1238*x1+1393*y1+1548*z1;
yu2=1518*x2+1651*y2+1783*z2;
a1+a2>=7500;
b1+b2>=3300;
c1+c2>=8000;
d1+d2>=1860;
31.2*x1+43.2*y1+55.2*z1>=0.098*a1+0.06*b1+0.0245*c1+0.015*d1;
8000*x1+10000*y1+12000*z1>=1.5*a1+15*b1+2.5*c1+50*d1;
31.2*x2+43.2*y2+55.2*z2>=0.098*a2+0.06*b2+0.0245*c2+0.015*d2;
8000*x2+10000*y2+12000*z2>=1.5*a2+15*b2+2.5*c2+50*d2;
@gin(x1);@gin(y1);@gin(z1);@gin(x2);@gin(y2);@gin(z2);
@gin(a1);@gin(b1);@gin(c1);@gin(d1);
@gin(a2);@gin(b2);@gin(c2);@gin(d2);

我也说一句

发表

使用道具举报

nuaahappy

473 主题	6 听众	61 积分

升级 58.95%

该用户从未签到

自我介绍: 学生

3^#

发表于 2010-6-9 15:53 |只看该作者 |招呼Ta 关注Ta

第一次发帖，出状况了
model:
max=a1+a2+b1+b2+c1+c2+d1+d2;
38*a1+950*b1+7.8*c1+195*d1+35*a2+900*b2+7.2*c2+190*d2<=3537830-yu1-yu2;
yu1=1238*x1+1393*y1+1548*z1;
yu2=1518*x2+1651*y2+1783*z2;
a1+a2>=7500;
b1+b2>=3300;
c1+c2>=8000;
d1+d2>=1860;
31.2*x1+43.2*y1+55.2*z1>=0.098*a1+0.06*b1+0.0245*c1+0.015*d1;
8000*x1+10000*y1+12000*z1>=1.5*a1+15*b1+2.5*c1+50*d1;
31.2*x2+43.2*y2+55.2*z2>=0.098*a2+0.06*b2+0.0245*c2+0.015*d2;
8000*x2+10000*y2+12000*z2>=1.5*a2+15*b2+2.5*c2+50*d2;
@gin(x1);@gin(y1);@gin(z1);@gin(x2);@gin(y2);@gin(z2);
@gin(a1);@gin(b1);@gin(c1);@gin(d1);
@gin(a2);@gin(b2);@gin(c2);@gin(d2);

我也说一句

发表

使用道具举报

linmatsas

53 主题	13 听众	3592 积分

逍遥游

TA的每日心情

	奋斗 2014-12-2 09:53

签到天数: 54 天

[LV.5]常住居民I

自我介绍: 额。。。。世界上最讨厌的事情就是自我介绍。。。

群组: Matlab讨论组

群组: 数学建模

群组: 小草的客厅

群组: 2012数学一考研交流

群组: C 语言讨论组

4^#

发表于 2010-6-9 16:53 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

错误代码是什么啊？我运行怎么没有出错呢？运行结果如下：
Global optimal solution found.
  Objective value:                            34941.00
  Objective bound:                            34941.00
  Infeasibilities:                            0.000000
  Extended solver steps:                            0
  Total solver iterations:                         14

                     Variable          Value       Reduced Cost
                           A1       0.000000          -1.000000
                           A2       0.000000          -1.000000
                           B1       0.000000          -1.000000
                           B2       3300.000          -1.000000
                           C1       0.000000          -1.000000
                           C2       29781.00          -1.000000
                           D1       0.000000          -1.000000
                           D2       1860.000          -1.000000
                        YU1       0.000000          0.000000
                        YU2       0.000000          0.000000
                           X1       0.000000          0.000000
                           Y1       0.000000          0.000000
                           Z1       0.000000          0.000000
                        YYU2       586745.0          0.000000
                           X2       158.0000          0.000000
                           Y2       114.0000          0.000000
                           Z2       89.00000          0.000000
                        DA1       7500.000          0.000000

                        Row Slack or Surplus    Dual Price
                           1       34941.00          1.000000
                           2       6.800000          0.000000
                           3       0.000000          0.000000
                           4       0.000000          0.000000
                           5       0.000000          0.000000
                           6       0.000000          0.000000
                           7       21781.00          0.000000
                           8       0.000000          0.000000
                           9       0.000000          0.000000
                           10       0.000000          0.000000
                           11       13811.67          0.000000
                           12       3255048.          0.000000

还有提点建议，程序后面要写 end
而且lingo内部运行线性条件比较快例如：
x+y<=0；这样子，还有就是变量设的有点乱…………