24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2001  |  回复: 4

小龙_1217

木虫 (初入文坛)

[交流] 【求助】matlab 二次规划的优化的问题 已有1人参与

我们MATLAB里的二次规划算法去优化一个方程 我编的程序如下:
H=eye(12,12);
c=zeros(12,1);
A=[1 0 -0.5 0 0 0 0 0 0 0 0 0;
      0 0 0 1 0 -0.5 0 0 0 0 0 0;
      0 0 0 0 0 0 1 0 -0.5 0 0 0;
      0 0 0 0 0 0 0 0 0 1 0 -0.5];
b=[0;0;0;0];
Aeq=[eye(3) eye(3) eye(3) eye(3);L1 L2 L3 L4];beq=[80;0;200;0;12;0];
        x_1=-0.15;y_1=0.4;z_1=0.35;
        x_2=0.15;y_2=0.4;z_2=0.35;
        x_3=-0.15;y_3=0;z_3=0.35;
        x_4=0.15;y_4=0;z_4=0.35;
        L1=[0 z_1 y_1;-z_1 0 -x_1;-y_1 x_1 0];
        L2=[0 z_2 y_2;-z_2 0 -x_2;-y_2 x_2 0];
        L3=[0 z_3 y_3;-z_3 0 -x_3;-y_3 x_3 0];
        L4=[0 z_4 y_4;-z_4 0 -x_4;-y_4 x_4 0];
VLB=[];VUB=[];
[x,z]=quadprog(H,c,A,b,Aeq,beq,VLB,VUB)
在程序上我个人感觉好像没什么大错呀 得出的 结果是:
x =

   32.3333     (x1)
    0.0000      (x2)
   64.6667     (x3)
  -32.3333     (x4)
   -0.0000      (x5)
  -64.6667     (x6)
   64.3333     (x7)
    0.0000      (x8)
  168.6667    (x9)
   15.6667     (x10)
   -0.0000      (x11)
   31.3333     (x12)
这个结果肯定是不对的  做过很多修改 其结果与上面的都差不多  这里 首先就是要保证x3,x6,x9,x12是大于零的 (貌似这个限制条件我也不会加呀)x6就是小于零的  再者 x7 x9这个数都过于的大 我怎么弄都不知道如何让这数减小  这些理想的值应该在【20 0 50 20 0 50 20 0 50 20 0 50】这个范围左右   我试过家范围 VLB=[0;-100;0;0;-100;0;0;-100;0;0;-100;0];VUB=[110;100;100;110;100;100;110;100;100;110;100;100]; 但是出来的结果还是差不多的 x6还是负的  x7 x9结果过于大 尤其是x9  给了范围 结果还是在范围之外  规定在50左右 每次都150 160的 太离谱了  貌似范围没有作用  我尝试修改了很多 但是每次都是失败 求求哪位大侠能帮帮忙我将万分感谢 小弟不胜感激呀
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

包公

新虫 (初入文坛)

会不会
2楼2010-09-17 11:21:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dreamred8087

金虫 (正式写手)

小龙_1217(金币+10): 谢谢 2011-03-16 16:57:04
建议楼主加上变量的上下界,不然Matlab可能找到的是其他的局部解,不知道对不对,呵呵...
3楼2010-09-17 17:46:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

winjade

银虫 (小有名气)

小龙_1217(金币+10): 谢谢 2011-03-16 16:57:17
从函数的使用来看,quadprog()是求最小值的,
x=quadprog (H,f,A,b,Aeq,beq,lb,ub)
其中f为规划函数,而你的程序中是一个零矩阵,再确认一下。
4楼2010-10-09 12:04:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ab77977

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by winjade at 2010-10-09 12:04:17
从函数的使用来看,quadprog()是求最小值的,
x=quadprog (H,f,A,b,Aeq,beq,lb,ub)
其中f为规划函数,而你的程序中是一个零矩阵,再确认一下。

f是规划函数吗?
5楼2013-04-19 10:09:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小龙_1217 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 复试调剂 +3 asdasdassda 2026-04-05 3/150 2026-04-05 17:26 by zhousanduo
[考研] 318求调剂 +11 ykyhsa 2026-04-05 13/650 2026-04-05 12:44 by aidndnjck
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 一志愿西北农林畜牧专硕336分求调剂 +3 5ourr 2026-04-03 3/150 2026-04-05 10:40 by JOKER0401
[考研] 295求调剂 +4 A你好研究生 2026-04-04 5/250 2026-04-04 22:46 by yu221
[考研] 321求调剂 +6 认真求上学 2026-04-03 6/300 2026-04-04 19:51 by dongzh2009
[考研] 321求调剂 +13 认真求上学 2026-04-02 13/650 2026-04-04 18:23 by macy2011
[考研] 0856调剂 +8 曲听筠 2026-03-30 8/400 2026-04-04 08:46 by tianyyysss
[考研] 考研调剂 +3 15615482637 2026-04-03 3/150 2026-04-03 22:50 by ms629
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
[考研] 320求调剂 +5 振—TZ 2026-04-02 5/250 2026-04-03 14:42 by fxue1114
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦! 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 07生物学求调剂 一志愿同济大学359分 +3 LAMC. 2026-03-30 3/150 2026-04-02 10:26 by 18828373951
[考研] 0817化工学硕调剂 +11 努力上岸中! 2026-03-31 11/550 2026-04-01 20:30 by 赖春艳
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
[考研] 262求调剂 +7 ZZ..000 2026-03-30 8/400 2026-03-31 10:05 by cal0306
[考研] 083000环境科学与工程调剂,总分281 +4 橙子(胜意) 2026-03-30 4/200 2026-03-31 00:44 by Linzejun
信息提示
请填处理意见