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è)極簡(jiǎn)單的功能:
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ù)此樓

» 猜你喜歡

一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有4人回復(fù)
一志愿武理材料305分求調(diào)劑已經(jīng)有6人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有3人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有7人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線已經(jīng)有7人回復(fù)
295求調(diào)劑已經(jīng)有5人回復(fù)
354求調(diào)劑已經(jīng)有5人回復(fù)
材料專業(yè)求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(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è)簡(jiǎn)單的積分（急�。。。� 已經(jīng)有7人回復(fù)
求助��！關(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ā)表

[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 求b區(qū)學(xué)校調(diào)劑 +3	周56 2026-03-11	3/150	2026-03-13 16:20 by JourneyLucky
[考研] 求調(diào)劑 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-12	10/500	2026-03-13 15:42 by ms629
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu