版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

凝雪小小

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 421.8
散金: 10
帖子: 105
在線: 61.7小時
蟲號: 1394131
注冊: 2011-09-07
性別: GG
專業(yè): 控制理論與方法

[求助] matlab線性矩陣不等式求解做優(yōu)化問題

在求解如下凸優(yōu)化問題中遇到了問題：
1.不知道如何編程得到最優(yōu)的γ²
2.求解提示為沒有可行解
是否是我退出的不等式有問題  或者其他方面有問題
A=[0,1,0,-1;-882,-28.4,0,28.4;0,0,0,1;1696.15,54.62,-657.12,-2680.39;];
B=[0,0,0,-0.039]';
Bw=[0,0,-1,2625.77]';
C1=[-882,-28.4,0,28.4];
C2=[0,0,10,0];
Ea=[0,0,-65.71,262.58]
%Eb=0
L=[0,0,0,1]'
I=eye(1)
gam=4.2; %γ=gam
a=0.01; %ρ=a
b=0.1; %ε=b
c=inv(b);
umax=3000

setlmis([]);
X=lmivar(1,[4 1]);    %定義決策變量
Z=lmivar(2,[1 4]);

lmiterm([1 1 1 X],A,1,'s');
lmiterm([1 1 1 Z],B,1,'s');
lmiterm([1 1 2 0],L);
lmiterm([1 1 3 X],1,Ea'); %Ea→Ea'
lmiterm([1 1 4 0],Bw);
lmiterm([1 1 5 X],1,C1');
lmiterm([1 2 2 0],-c);
lmiterm([1 3 3 0],-b);
lmiterm([1 4 4 0],-gam^2);
lmiterm([1 5 5 0],-1);
lmiterm([-2 1 1 X],1,1);
lmiterm([3 1 1 0],-1);
lmiterm([3 1 2 Z],0.95,1);
lmiterm([3 2 2 X],-3000,1);
lmiterm([4 1 1 0],-1);
lmiterm([4 1 2 X],C2,0.95);
lmiterm([4 2 2 X],-1,1);

lmisys=getlmis; %完成LTI框架的設(shè)設(shè)置
[tmin,xfeas]=feasp(lmisys);  %求解可行解問題
X=dec2mat(lmisys,xfeas,X);  %提取解矩陣  把決策變量轉(zhuǎn)化為矩陣形式
Z=dec2mat(lmisys,xfeas,Z);
P=inv(X);
K=Z*P

運行提示為：
Solver for LMI feasibility problems L(x) < R(x)
This solver minimizes  t  subject to  L(x) < R(x) + t*I
The best value of t should be negative for feasibility

Iteration : Best value of t so far

   1                   2487.161836
   2                   1661.789005
   3                   1200.565677
   4                   542.424422
   5                   311.999933
   6                   311.999933
   7                   279.917289
   8                   279.917289
***                new lower bound: 36.447556

Result:  best value of t: 279.917289
      f-radius saturation:  0.000% of R = 1.00e+009

These LMI constraints were found infeasible

K =

2.0861 0.1139  -16.8718 3.5340
麻煩各位指出其中有錯誤的地方不勝感激！

回復(fù)此樓

» 猜你喜歡

268求調(diào)劑已經(jīng)有8人回復(fù)
26申博已經(jīng)有4人回復(fù)
211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) 已經(jīng)有7人回復(fù)
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有7人回復(fù)
材料工程專碩274一志愿211求調(diào)劑已經(jīng)有6人回復(fù)
一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑已經(jīng)有4人回復(fù)
332求調(diào)劑已經(jīng)有5人回復(fù)
271求調(diào)劑已經(jīng)有14人回復(fù)
275求調(diào)劑已經(jīng)有4人回復(fù)

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

matlab用多元線性最小二乘法求解多元線性方程的問題已經(jīng)有11人回復(fù)
求助有關(guān)Matlab有約束非線性最優(yōu)化問題已經(jīng)有9人回復(fù)
求matlab求解一個方程組的問題已經(jīng)有6人回復(fù)
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
怎樣求解約束條件為不等式的線性規(guī)劃問題？大蝦請進已經(jīng)有10人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
MATLAB非線性方程組求解~~~反應(yīng)釜模擬已經(jīng)有9人回復(fù)
請教！大家用哪種軟件求解約束非線性最優(yōu)化問題比較好已經(jīng)有14人回復(fù)
【求助】matlab用lsqnonlin非線性擬合的問題【已解決】已經(jīng)有8人回復(fù)
【求助】matlab 二次規(guī)劃的優(yōu)化的問題已經(jīng)有4人回復(fù)
【分享】學(xué)習(xí)MATLAB的看過來，《MATLAB中文幫助》通過矩陣讓你更好的處理數(shù)據(jù)！已經(jīng)有16人回復(fù)

我是帝都小愣神

1樓 2012-02-26 22:12:31

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

回帖支持 ( 顯示支持度最高的前 50 名 )

幸運天使

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 378.2
散金: 38
紅花: 1
帖子: 112
在線: 70.3小時
蟲號: 1095076
注冊: 2010-09-10
專業(yè): 機器人學(xué)及機器人技術(shù)

樓主，你好！

我想請問你的問題解決了沒有？如果解決了，可否把代碼發(fā)給我看看，我也遇到了相同的問題，不甚感激！

贊一下(3人)

回復(fù)此樓

8樓2014-05-14 14:53:27

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

dc652037717

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 521.5
帖子: 42
在線: 22.7小時
蟲號: 1752935
注冊: 2012-04-13
性別: GG
專業(yè): 控制論中的數(shù)學(xué)方法

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
凝雪小小: 金幣+5 2012-04-16 10:05:47

求最優(yōu)的r值應(yīng)該是用mincx語句，如果只用feasp語句得加上循環(huán)語句。

贊一下(1人)

回復(fù)此樓

6樓2012-04-15 13:30:42

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

普通回帖

lhjyml

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1648.3
帖子: 73
在線: 193.5小時
蟲號: 330854
注冊: 2007-03-24
專業(yè): 控制理論與方法

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

感謝參與，應(yīng)助指數(shù) +1
凝雪小小(金幣+18): ★★★很有幫助 2012-02-27 14:54:55

1）最優(yōu)的γ²，需要加入一個循環(huán)語句。
2）無解原因：（1）首先仔細檢查你的理論推導(dǎo)，若沒有問題，檢查程序，若還沒有問題，在條件參數(shù)（比如b=0.1; %ε=b）

贊一下

回復(fù)此樓

2樓2012-02-27 12:15:32

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

lhjyml

木蟲 (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1648.3
帖子: 73
在線: 193.5小時
蟲號: 330854
注冊: 2007-03-24
專業(yè): 控制理論與方法

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

★ ★ ★ ★ ★ ★ ★
凝雪小小: 金幣+7 2012-05-09 17:29:00

無解原因：再調(diào)節(jié)相關(guān)參數(shù)（比如b=0.1; %ε=b）

贊一下

回復(fù)此樓

3樓2012-02-27 12:16:50

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

凝雪小小

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 421.8
散金: 10
帖子: 105
在線: 61.7小時
蟲號: 1394131
注冊: 2011-09-07
性別: GG
專業(yè): 控制理論與方法

引用回帖:

樓: Originally posted by lhjyml at 2012-02-27 12:15:32:
1）最優(yōu)的γ&sup2，需要加入一個循環(huán)語句。
2）無解原因：（1）首先仔細檢查你的理論推導(dǎo)，若沒有問題，檢查程序，若還沒有問題，在條件參數(shù)（比如b=0.1; %ε=b）

1）要加一個怎樣的循環(huán)語句呢之前沒有編過求高手指點呀
2）不確定的參數(shù)有什么試的經(jīng)驗么

贊一下

回復(fù)此樓

我是帝都小愣神

4樓2012-02-27 14:54:28

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

凝雪小小

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 421.8
散金: 10
帖子: 105
在線: 61.7小時
蟲號: 1394131
注冊: 2011-09-07
性別: GG
專業(yè): 控制理論與方法

引用回帖:

樓: Originally posted by lhjyml at 2012-02-27 12:16:50:
無解原因：再調(diào)節(jié)相關(guān)參數(shù)（比如b=0.1; %ε=b）

這些參數(shù)是每個都一起調(diào)呢還是每次只調(diào)一個呢

贊一下

回復(fù)此樓

我是帝都小愣神

5樓2012-02-29 20:36:17

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

ep22stone

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 148.3
紅花: 1
帖子: 64
在線: 55.8小時
蟲號: 2022741
注冊: 2012-09-23
性別: MM
專業(yè): 控制理論與方法

mincx語句如何用，樓主能不能介紹一下

贊一下

回復(fù)此樓

7樓2013-03-23 22:58:05

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

大母蟲子

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 422.8
帖子: 46
在線: 11小時
蟲號: 3097413
注冊: 2014-03-28
性別: GG
專業(yè): 控制理論與方法

LZ問題解決了嗎？過去好久了

贊一下

回復(fù)此樓

術(shù)業(yè)有專攻

9樓2015-01-17 14:58:18

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

guaixiaofu

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 2.5
帖子: 1
在線: 26分鐘
蟲號: 3974526
注冊: 2015-07-16

我有一個問題，怎么才能傳上去

贊一下

回復(fù)此樓

10樓2015-08-27 16:36:26

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主凝雪小小的主題更新

返回列表

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

[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 293求調(diào)劑 +5	世界首富 2026-03-11	5/250	2026-03-17 10:16 by Sammy2
[考研] 本人考085602 化學(xué)工程專碩 +15	不知道叫什么！ 2026-03-15	17/850	2026-03-17 10:13 by Sammy2
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 290求調(diào)劑 +5	孔志浩 2026-03-12	10/500	2026-03-16 09:01 by 余暉&
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 279求調(diào)劑 +3	抓著星星的女孩 2026-03-10	3/150	2026-03-13 23:47 by userper
[考研] 285 求調(diào)劑資源與環(huán)境一志愿北京化工大學(xué) +3	未名考生 2026-03-10	3/150	2026-03-13 23:04 by JourneyLucky
[考研] 279求調(diào)劑 +3	Dizzy123@ 2026-03-10	3/150	2026-03-13 23:02 by JourneyLucky
[考研] 290求調(diào)劑 +9	ADT 2026-03-11	9/450	2026-03-13 21:55 by JourneyLucky
[考研] 310求調(diào)劑 +3	【上上簽】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考研] 考研調(diào)劑 +4	芬達46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[論文投稿] 投稿問題 5+4	星光燦爛xt 2026-03-12	6/300	2026-03-13 14:17 by god_tian
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考研] 工科0856專碩化學(xué)工程269能調(diào)劑嗎 +10	我想讀研11 2026-03-10	10/500	2026-03-13 10:14 by Yuyi.

24小時熱門版塊排行榜

[求助] matlab線性矩陣不等式求解做優(yōu)化問題

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

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

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

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

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