凝雪小小

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

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

在求解如下凸優(yōu)化問(wèn)題中遇到了問(wèn)題：
1.不知道如何編程得到最優(yōu)的γ²
2.求解提示為沒(méi)有可行解
是否是我退出的不等式有問(wèn)題  或者其他方面有問(wèn)題
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);  %求解可行解問(wèn)題
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)格式了，假專(zhuān)家開(kāi)始管格式了已經(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ù)
材料工程專(zhuān)碩調(diào)劑已經(jīng)有3人回復(fù)
26考研求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專(zhuān)業(yè)都可以已經(jīng)有5人回復(fù)
311求調(diào)劑已經(jīng)有7人回復(fù)
303求調(diào)劑已經(jīng)有3人回復(fù)

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

matlab用多元線性最小二乘法求解多元線性方程的問(wèn)題已經(jīng)有11人回復(fù)
求助有關(guān)Matlab有約束非線性最優(yōu)化問(wèn)題已經(jīng)有9人回復(fù)
求matlab求解一個(gè)方程組的問(wèn)題已經(jīng)有6人回復(fù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
怎樣求解約束條件為不等式的線性規(guī)劃問(wèn)題？大蝦請(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)化問(wèn)題比較好已經(jīng)有14人回復(fù)
【求助】matlab用lsqnonlin非線性擬合的問(wèn)題【已解決】已經(jīng)有8人回復(fù)
【求助】matlab 二次規(guī)劃的優(yōu)化的問(wèn)題已經(jīng)有4人回復(fù)
【分享】學(xué)習(xí)MATLAB的看過(guò)來(lái)，《MATLAB中文幫助》通過(guò)矩陣讓你更好的處理數(shù)據(jù)！已經(jīng)有16人回復(fù)

我是帝都小愣神

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

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

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

幸運(yùn)天使

銅蟲(chóng) (小有名氣)

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

樓主，你好！

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

贊一下(3人)

回復(fù)此樓

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

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

dc652037717

新蟲(chóng) (初入文壇)

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

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

普通回帖

lhjyml

木蟲(chóng) (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1648.3
帖子: 73
在線: 193.5小時(shí)
蟲(chóng)號(hào): 330854
注冊(cè): 2007-03-24
專(zhuān)業(yè): 控制理論與方法

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

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

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

贊一下

回復(fù)此樓

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

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

lhjyml

木蟲(chóng) (小有名氣)

應(yīng)助: 5 (幼兒園)
金幣: 1648.3
帖子: 73
在線: 193.5小時(shí)
蟲(chóng)號(hào): 330854
注冊(cè): 2007-03-24
專(zhuān)業(yè): 控制理論與方法

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

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

無(wú)解原因：再調(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的回帖

凝雪小小

新蟲(chóng) (小有名氣)

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

引用回帖:

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

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

贊一下

回復(fù)此樓

我是帝都小愣神

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

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

凝雪小小

新蟲(chóng) (小有名氣)

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

引用回帖:

樓: Originally posted by lhjyml at 2012-02-27 12:16:50:
無(wú)解原因：再調(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

銅蟲(chóng) (小有名氣)

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

mincx語(yǔ)句如何用，樓主能不能介紹一下

贊一下

回復(fù)此樓

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

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

大母蟲(chóng)子

銀蟲(chóng) (初入文壇)

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

LZ問(wèn)題解決了嗎？過(guò)去好久了

贊一下

回復(fù)此樓

術(shù)業(yè)有專(zhuān)攻

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

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

guaixiaofu

新蟲(chóng) (初入文壇)

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 311求調(diào)劑 +7	冬十三 2026-03-15	7/350	2026-03-17 15:26 by 霧散后相遇lc
[考研] 268求調(diào)劑 +6	好運(yùn)連綿不絕 2026-03-12	7/350	2026-03-17 14:56 by 呦呦憂郁
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9	柴郡貓_ 2026-03-12	9/450	2026-03-17 10:18 by Sammy2
[考研] 085601求調(diào)劑 +3	Du.11 2026-03-16	3/150	2026-03-16 20:42 by 無(wú)際的草原
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 344求調(diào)劑 +3	knight344 2026-03-16	3/150	2026-03-16 09:42 by 無(wú)際的草原
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 材料工程327求調(diào)劑 +3	xiaohe12w 2026-03-11	3/150	2026-03-14 20:20 by ms629
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 復(fù)試調(diào)劑 +9	Copy267 2026-03-10	9/450	2026-03-13 23:45 by userper
[考研] 材料專(zhuān)碩288分求調(diào)劑一志愿211 +4	在家想你 2026-03-11	4/200	2026-03-13 22:49 by JourneyLucky
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考博] 福州大學(xué)楊黃浩課題組招收2026年專(zhuān)業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655