查看: 1682|回复: 7

求助！关于Matlab线性规划的问题。。。请大家帮帮忙，谢谢！

字体大小: 正常放大

JimOu

1 主题	7 听众	118 积分

升级 9%

TA的每日心情

	开心 2015-3-17 00:22

签到天数: 30 天

[LV.5]常住居民I

自我介绍: 数模人

电梯直达

1^#

发表于 2013-1-22 23:09 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

5体力

                              题目 max f=0.05x0+0.27x1+0.19x2+0.185x3+0.185x4
                                       st.  x0+1.01x1+1.02x2+1.045x3+1.065x4=1
                                             0.025x1<=a
                                             0.015x2<=a
                                             0.055x3<=a
                                             0.026x4<=a
                                             x0,x1,x2,x3,x4>=0
                  a=0，从步长为0.001进行循环搜索，到a=0.05为止。

编程如下：
                  clc,clear a=0;
                  hold on
                  while a<0.05
                     c=[-0.05,-0.27,-0.19,-0.185,-0.185];
                     A=[zeros(4,1),diag([0.025,0.015,0.055,0.026])];
                     b=a*ones(4,1);
                     Aeq=[1,1.01,1.02,1.045,1.065];
                     beq=1;
                     LB=zeros(5,1);
                     [x,Q]=linprog(c,A,b,Aeq,beq,LB);
                     Q=-Q;
                     plot(a,Q,'.r');
                     a=a+0.001;
                  end
                     xlabel('a'),ylabel('Q')
结果输出的图像为空白（就是图像是什么也没有），想问下大家这个编程代码哪里错了？
图片：

最佳答案

新手ttt 查看完整内容

你把clear后面加一个分号就可以了，你把a给clear了！！！ clc,clear ;a=0; hold on while a

zan

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

使用道具举报

新手ttt

3 主题	7 听众	63 积分

升级 61.05%

TA的每日心情

	开心 2013-2-1 20:05

签到天数: 13 天

[LV.3]偶尔看看II

自我介绍: 擅长数学

2^#

发表于 2013-1-22 23:09 |只看该作者 |招呼Ta 关注Ta

你把clear后面加一个分号就可以了，你把a给clear了！！！
clc,clear ;a=0;
                  hold on
                  while a<0.05
                     c=[-0.05,-0.27,-0.19,-0.185,-0.185];
                     A=[zeros(4,1),diag([0.025,0.015,0.055,0.026])];
                     b=a*ones(4,1);
                     Aeq=[1,1.01,1.02,1.045,1.065];
                     beq=1;
                     LB=zeros(5,1);
                     [x,Q]=linprog(c,A,b,Aeq,beq,LB);
                     Q=-Q;
                     plot(a,Q,'.r');
                     a=a+0.001;
                  end
                     xlabel('a'),ylabel('Q')

我也说一句

发表

使用道具举报

JimOu

1 主题	7 听众	118 积分

升级 9%

TA的每日心情

	开心 2015-3-17 00:22

签到天数: 30 天

[LV.5]常住居民I

自我介绍: 数模人

3^#

发表于 2013-1-23 09:29 |只看该作者 |招呼Ta 关注Ta

自顶。。。

我也说一句

发表

使用道具举报

JimOu

1 主题	7 听众	118 积分

升级 9%

TA的每日心情

	开心 2015-3-17 00:22

签到天数: 30 天

[LV.5]常住居民I

自我介绍: 数模人

4^#

发表于 2013-1-24 11:15 |只看该作者 |招呼Ta 关注Ta

自顶，希望有能力的人可以答复。

我也说一句

发表

使用道具举报

新手ttt

3 主题	7 听众	63 积分

升级 61.05%

TA的每日心情

	开心 2013-2-1 20:05

签到天数: 13 天

[LV.3]偶尔看看II

自我介绍: 擅长数学

5^#

发表于 2013-1-24 11:21 |只看该作者 |招呼Ta 关注Ta

你把clear后面加一个分号就可以了，你把a给clear了！！！
clc,clear ;a=0;
                  hold on
                  while a<0.05
                     c=[-0.05,-0.27,-0.19,-0.185,-0.185];
                     A=[zeros(4,1),diag([0.025,0.015,0.055,0.026])];
                     b=a*ones(4,1);
                     Aeq=[1,1.01,1.02,1.045,1.065];
                     beq=1;
                     LB=zeros(5,1);
                     [x,Q]=linprog(c,A,b,Aeq,beq,LB);
                     Q=-Q;
                     plot(a,Q,'.r');
                     a=a+0.001;
                  end
                     xlabel('a'),ylabel('Q')

点评

JimOu 哎啊~太粗心了。。。详情回复发表于 2013-1-25 09:52

我也说一句

发表

使用道具举报

新手ttt

3 主题	7 听众	63 积分

升级 61.05%

TA的每日心情

	开心 2013-2-1 20:05

签到天数: 13 天

[LV.3]偶尔看看II

自我介绍: 擅长数学

6^#

发表于 2013-1-24 11:26 |只看该作者 |招呼Ta 关注Ta

附件是你的图

我也说一句

发表

使用道具举报

新手ttt

3 主题	7 听众	63 积分

升级 61.05%

TA的每日心情

	开心 2013-2-1 20:05

签到天数: 13 天

[LV.3]偶尔看看II

自我介绍: 擅长数学

7^#

发表于 2013-1-24 11:34 |只看该作者 |招呼Ta 关注Ta

你把clear后面加一个分号就可以了，你把a给clear了！！！ clc,clear ;a=0; hold on while a<0.05 c=[-0.05,-0.27,-0.19,-0.185,-0.185]; A=[zeros(4,1),diag([0.025,0.015,0.055,0.026])]; b=a*ones(4,1); Aeq=[1,1.01,1.02,1.045,1.065]; beq=1; LB=zeros(5,1); [x,Q]=linprog(c,A,b,Aeq,beq,LB); Q=-Q; plot(a,Q,'.r'); a=a+0.001; end xlabel('a'),ylabel('Q')