24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2646  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fxlily2008

木虫 (小有名气)

[求助] 关于matlab的LMI求助

各位学长学姐,我从来没有用过matlab的LMI工具箱编写过程序,现在我想用它的mincx实现一个极简单的功能:
target:  min(0.1*x1+0.2*x2+0.3*x3+0.4*x4)
   s.t. :  x1,x2,x3,x4>=0
            x1+x2+x3+x4>10

我试着编了一个小程序,但是结果好像不对,请高手帮助一下吧,不胜感激,小弟在线等。

最好回复一个完整的程序吧,谢谢!!
回复此楼

» 猜你喜欢

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

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

EVRYIBISC

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
fxlily2008: 金币+5, ★★★很有帮助 2013-06-07 09:04:54
xiegangmai: 金币+2, 谢谢参与! 2013-06-07 20:36:25
clear all;
close all;
clc

setlmis([]);

x1 = lmivar(1,[1,1]);
x2 = lmivar(1,[1,1]);
x3 = lmivar(1,[1,1]);
x4 = lmivar(1,[1,1]);

br1 = newlmi;

lmiterm([-br1 1 1 x1],1,1);

br2 = newlmi;

lmiterm([-br2 1 1 x2],1,1);

br3 = newlmi;

lmiterm([-br3 1 1 x3],1,1);

br4 = newlmi;

lmiterm([-br4 1 1 x4],1,1);

br5 = newlmi;

lmiterm([-br5 1 1 x1],1,1);
lmiterm([-br5 1 1 x2],1,1);
lmiterm([-br5 1 1 x3],1,1);
lmiterm([-br5 1 1 x4],1,1);
lmiterm([br5 1 1 0],10);

lmisys = getlmis;

n = decnbr(lmisys);
c = zeros(n,1);
c(n,1)=1;

for jj = 1:n
    [x1j,x2j,x3j,x4j] = defcx(lmisys,jj,x1,x2,x3,x4);
    c(jj) = 0.1*x1j+0.2*x2j+0.3*x3j+0.4*x4j
end
%c = mat2dec(lmisys,zeros(2,2),zeros(1,2),1);

options = [1e-5,0,0,0,0];

[copt,xopt] = mincx(lmisys,c,options);

x1opt = dec2mat(lmisys,xopt,x1)
x2opt = dec2mat(lmisys,xopt,x2)
x3opt = dec2mat(lmisys,xopt,x3)
x4opt = dec2mat(lmisys,xopt,x4)

运行结果:
c =

    0.1000
         0
         0
    1.0000


c =

    0.1000
    0.2000
         0
    1.0000


c =

    0.1000
    0.2000
    0.3000
    1.0000


c =

    0.1000
    0.2000
    0.3000
    0.4000


Solver for linear objective minimization under LMI constraints

Iterations   :    Best objective value so far

     1                   1.734323
     2                   1.133169
***                 new lower bound:     0.563924
     3                   1.022152
***                 new lower bound:     0.834201
     4                   1.013815
***                 new lower bound:     0.977789
     5                   1.000993
***                 new lower bound:     0.997951
     6                   1.000010
***                 new lower bound:     0.998732
     7                   1.000010
***                 new lower bound:     0.999878
     8                   1.000001
***                 new lower bound:     0.999973
     9                   1.000001
***                 new lower bound:     0.999994

Result:  feasible solution of required accuracy
          best objective value:     1.000001
          guaranteed relative accuracy:  6.36e-06
          f-radius saturation:  0.000% of R =  1.00e+09


x1opt =

   10.0000


x2opt =

   1.2620e-06


x3opt =

   6.3082e-07


x4opt =

   4.2051e-07

>>
9楼2013-06-06 21:35:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

lijie169

铜虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与! 2013-05-26 12:36:34
你这个条件就有问题,没有解!
2楼2013-05-15 20:07:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxlily2008

木虫 (小有名气)

引用回帖:
2楼: Originally posted by lijie169 at 2013-05-15 20:07:30
你这个条件就有问题,没有解!

不是啊,虽然x1,x2,x3,x4可以无限大,但min(……)本身就是约束啊,理论上正确的结果应该是x1=10.0000001,x2=x3=x4=0对吧。但是用matlab的LMI怎么写程序啊,求帮助
3楼2013-05-15 21:09:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fxlily2008

木虫 (小有名气)

另外,其中的等于号可以不考虑,没关系的
4楼2013-05-15 21:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 材料085601调剂 +25 何润采123 2026-04-10 27/1350 2026-04-10 23:17 by Ftglcn90
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 一志愿中南大学物理学,英一66,求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 调剂 +24 不逢春 2026-04-07 25/1250 2026-04-10 09:39 by 690616278
[考研] 277求调剂 +19 倪建设 2026-04-06 19/950 2026-04-10 09:24 by guosr9609
[考研] 工科273调剂 +5 X1999 2026-04-09 6/300 2026-04-10 07:52 by 1753564080
[考研] 293分求调剂,外语为俄语 +7 加一一九 2026-04-07 10/500 2026-04-08 20:14 by yutian743
[考研] 一志愿郑州大学085600求调剂 +21 吃的不少 2026-04-05 24/1200 2026-04-08 16:47 by sunhuadong
[考研] 331求调剂 +5 张元一 2026-04-07 6/300 2026-04-07 22:13 by hemengdong
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 301求调剂 +18 121. 2026-04-04 18/900 2026-04-07 17:49 by 蓝云思雨
[考研] 085602调剂 初试总分335 +10 19123253302 2026-04-05 10/500 2026-04-07 15:23 by 小乔同学ya
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 一志愿河北工业大学材料工程,初试344求专硕调剂 +6 15933906766 2026-04-05 6/300 2026-04-06 13:21 by 无际的草原
[考研] 332求调剂 +17 小小孟... 2026-04-05 18/900 2026-04-06 09:51 by 蓝云思雨
[考研] 302分 085601求调剂推荐 +11 zyx上岸! 2026-04-05 11/550 2026-04-05 22:13 by dongzh2009
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
信息提示
请填处理意见