fxlily2008

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3438.7
散金: 10
帖子: 141
在線: 210.9小時(shí)
蟲號(hào): 1680266
注冊(cè): 2012-03-10
專業(yè): 控制理論與方法

[求助] 關(guān)于matlab的LMI求助

各位學(xué)長學(xué)姐，我從來沒有用過matlab的LMI工具箱編寫過程序，現(xiàn)在我想用它的mincx實(shí)現(xiàn)一個(gè)極簡單的功能:
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

我試著編了一個(gè)小程序，但是結(jié)果好像不對(duì)，請(qǐng)高手幫助一下吧，不勝感激，小弟在線等。

最好回復(fù)一個(gè)完整的程序吧，謝謝��！

回復(fù)此樓

» 猜你喜歡

070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑已經(jīng)有9人回復(fù)
0703化學(xué)調(diào)劑，求各位老師收留已經(jīng)有9人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有11人回復(fù)
0703化學(xué) 305求調(diào)劑已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有5人回復(fù)
274求調(diào)劑已經(jīng)有3人回復(fù)
一志愿中海洋材料工程專碩330分求調(diào)劑已經(jīng)有3人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有9人回復(fù)
【同濟(jì)軟件】軟件（085405）考研求調(diào)劑已經(jīng)有3人回復(fù)
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑已經(jīng)有7人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助一個(gè)LMI問題已經(jīng)有10人回復(fù)
求助一個(gè)matlab關(guān)于分類討論的問題已經(jīng)有9人回復(fù)
求助關(guān)于Matlab中字符串的轉(zhuǎn)化問題已經(jīng)有4人回復(fù)
純新手，全方面求助關(guān)于MATLAB........ 已經(jīng)有9人回復(fù)
matlab畫圖求助已經(jīng)有11人回復(fù)
繼續(xù)求助Matlab繪圖！已經(jīng)有8人回復(fù)
求助matlab編程已經(jīng)有13人回復(fù)
求助svm在matlab中的應(yīng)用已經(jīng)有7人回復(fù)
【求助】matlab求解一個(gè)簡單的積分（急�。。。� 已經(jīng)有7人回復(fù)
求助�。￡P(guān)于MATLAB 已經(jīng)有6人回復(fù)
【求助】會(huì)用Matlab的LMI toolbox的請(qǐng)幫幫忙~~ 已經(jīng)有5人回復(fù)
【求助】時(shí)間同步的MATLAB仿真已經(jīng)有5人回復(fù)
【求助】matlab如何只讀取圖像的一部分？已經(jīng)有11人回復(fù)
【求助】matlab中如何將程序的結(jié)果輸出到一個(gè)文件中已經(jīng)有11人回復(fù)
求助一個(gè)差分方程的MATLAB實(shí)現(xiàn) 已經(jīng)有3人回復(fù)
【求助】關(guān)于matlab的高維數(shù)組已經(jīng)有6人回復(fù)
【求助】LMI求解-可行性解的問題已經(jīng)有18人回復(fù)
【求助】初學(xué)matlab 已經(jīng)有11人回復(fù)
【求助】matlab中保存計(jì)算數(shù)據(jù)的問題已經(jīng)有8人回復(fù)
【求助】用matlab編過有限差分法的請(qǐng)進(jìn) 已經(jīng)有9人回復(fù)

1樓 2013-05-15 19:35:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lijie169

銅蟲 (著名寫手)

應(yīng)助: 159 (高中生)
金幣: 161.5
散金: 2018
紅花: 7
帖子: 1338
在線: 491.5小時(shí)
蟲號(hào): 1180408
注冊(cè): 2010-12-30
專業(yè): 分離過程

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
xiegangmai: 金幣+1, 謝謝參與！ 2013-05-26 12:36:34

你這個(gè)條件就有問題，沒有解！

贊一下

回復(fù)此樓

2樓2013-05-15 20:07:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fxlily2008

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3438.7
散金: 10
帖子: 141
在線: 210.9小時(shí)
蟲號(hào): 1680266
注冊(cè): 2012-03-10
專業(yè): 控制理論與方法

引用回帖:

2樓: Originally posted by lijie169 at 2013-05-15 20:07:30
你這個(gè)條件就有問題，沒有解！

不是啊，雖然x1,x2,x3,x4可以無限大，但min(……)本身就是約束啊，理論上正確的結(jié)果應(yīng)該是x1=10.0000001，x2=x3=x4=0對(duì)吧。但是用matlab的LMI怎么寫程序啊，求幫助

贊一下

回復(fù)此樓

3樓2013-05-15 21:09:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fxlily2008

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3438.7
散金: 10
帖子: 141
在線: 210.9小時(shí)
蟲號(hào): 1680266
注冊(cè): 2012-03-10
專業(yè): 控制理論與方法

另外，其中的等于號(hào)可以不考慮，沒關(guān)系的

贊一下

回復(fù)此樓

4樓2013-05-15 21:10:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

idnk

金蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 1158
散金: 11
帖子: 240
在線: 100.8小時(shí)
蟲號(hào): 150845
注冊(cè): 2005-12-29
性別: GG
專業(yè): 控制理論與方法

★ ★
感謝參與，應(yīng)助指數(shù) +1
1592203609: 金幣+2, 應(yīng)助指數(shù)-1, 非應(yīng)助帖不要選擇應(yīng)助 2013-05-15 23:00:23

把你程序貼出來，大家一起看看什么問題

贊一下

回復(fù)此樓

5樓2013-05-15 22:36:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

elvinsu

金蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 1000.1
帖子: 118
在線: 218小時(shí)
蟲號(hào): 2466860
注冊(cè): 2013-05-16
性別: GG
專業(yè): 控制理論與方法

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fxlily2008: 金幣+5, ★★★很有幫助 2013-05-18 14:45:27
xiegangmai: 金幣+1, 謝謝參與！ 2013-05-19 08:43:35

LMI是可以處理，但是用線性規(guī)劃不是更直接嗎？ x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options)

贊一下

回復(fù)此樓

科學(xué)的歷程就是在假設(shè)求證中不斷進(jìn)步--巴瑞•馬歇爾教授

6樓2013-05-18 11:56:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fxlily2008

木蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3438.7
散金: 10
帖子: 141
在線: 210.9小時(shí)
蟲號(hào): 1680266
注冊(cè): 2012-03-10
專業(yè): 控制理論與方法

引用回帖:

6樓: Originally posted by elvinsu at 2013-05-18 11:56:08
LMI是可以處理，但是用線性規(guī)劃不是更直接嗎？ x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options)

哦，這樣啊，這個(gè)函數(shù)我不熟，我先試試，謝謝前輩指導(dǎo)，我先給你5分吧

贊一下

回復(fù)此樓

7樓2013-05-18 14:44:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

elvinsu

金蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 1000.1
帖子: 118
在線: 218小時(shí)
蟲號(hào): 2466860
注冊(cè): 2013-05-16
性別: GG
專業(yè): 控制理論與方法

引用回帖:

7樓: Originally posted by fxlily2008 at 2013-05-18 14:44:53
哦，這樣啊，這個(gè)函數(shù)我不熟，我先試試，謝謝前輩指導(dǎo)，我先給你5分吧...

不敢當(dāng)，謝謝加分！

回復(fù)此樓

科學(xué)的歷程就是在假設(shè)求證中不斷進(jìn)步--巴瑞•馬歇爾教授

8樓2013-05-19 21:00:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

EVRYIBISC

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 6
帖子: 1
在線:
蟲號(hào): 2499131
注冊(cè): 2013-06-06

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
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)

運(yùn)行結(jié)果:
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

>>

贊一下(3人)

回復(fù)此樓

9樓2013-06-06 21:35:52

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考研] 268求調(diào)劑 +6	簡單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 326求調(diào)劑 +4	諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15	7/350	2026-03-16 17:11 by 諾貝爾化學(xué)獎(jiǎng)覬?/a>
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky