凝雪小小

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 421.8
散金: 10
帖子: 105
在線: 61.7小時(shí)
蟲號(hào): 1394131
注冊(cè): 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

運(yùn)行提示為：
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
麻煩各位指出其中有錯(cuò)誤的地方不勝感激！

回復(fù)此樓

» 猜你喜歡

被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有3人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工求調(diào)劑已經(jīng)有5人回復(fù)
277調(diào)劑已經(jīng)有4人回復(fù)
290求調(diào)劑已經(jīng)有4人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有3人回復(fù)
26考研求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以已經(jīng)有5人回復(fù)
311求調(diào)劑已經(jīng)有7人回復(fù)
303求調(diào)劑已經(jīng)有3人回復(fù)

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

matlab用多元線性最小二乘法求解多元線性方程的問題已經(jīng)有11人回復(fù)
求助有關(guān)Matlab有約束非線性最優(yōu)化問題已經(jīng)有9人回復(fù)
求matlab求解一個(gè)方程組的問題已經(jīng)有6人回復(fù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
怎樣求解約束條件為不等式的線性規(guī)劃問題？大蝦請(qǐng)進(jìn) 已經(jīng)有10人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
MATLAB非線性方程組求解~~~反應(yīng)釜模擬已經(jīng)有9人回復(fù)
請(qǐng)教！大家用哪種軟件求解約束非線性最優(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ùn)天使

銅蟲 (小有名氣)

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

樓主，你好！

我想請(qǐng)問你的問題解決了沒有？如果解決了，可否把代碼發(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小時(shí)
蟲號(hào): 1752935
注冊(cè): 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小時(shí)
蟲號(hào): 330854
注冊(cè): 2007-03-24
專業(yè): 控制理論與方法

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

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

1）最優(yōu)的γ²，需要加入一個(gè)循環(huán)語句。
2）無解原因：（1）首先仔細(xì)檢查你的理論推導(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小時(shí)
蟲號(hào): 330854
注冊(cè): 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小時(shí)
蟲號(hào): 1394131
注冊(cè): 2011-09-07
性別: GG
專業(yè): 控制理論與方法

引用回帖:

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

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

贊一下

回復(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小時(shí)
蟲號(hào): 1394131
注冊(cè): 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ù)是每個(gè)都一起調(diào)呢還是每次只調(diào)一個(gè)呢

贊一下

回復(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小時(shí)
蟲號(hào): 2022741
注冊(cè): 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小時(shí)
蟲號(hào): 3097413
注冊(cè): 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分鐘
蟲號(hào): 3974526
注冊(cè): 2015-07-16

我有一個(gè)問題，怎么才能傳上去

贊一下

回復(fù)此樓

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

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

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

[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 0703化學(xué)調(diào)劑 +6	妮妮ninicgb 2026-03-15	9/450	2026-03-16 16:40 by houyaoxu
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +7	曦熙兮 2026-03-15	7/350	2026-03-16 16:34 by houyaoxu
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 求調(diào)劑 +7	18880831720 2026-03-11	7/350	2026-03-13 16:10 by JourneyLucky
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 274求調(diào)劑 +3	S.H1 2026-03-12	3/150	2026-03-13 15:15 by JourneyLucky
[考研] 296求調(diào)劑 +3	大口吃飯身體健 2026-03-13	3/150	2026-03-13 10:31 by 學(xué)員8dgXkO
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655