| 2 | 1/1 | 返回列表 |
| 查看: 1522 | 回復(fù): 1 | ||
huangmagudi鐵蟲(chóng) (小有名氣)
|
[求助]
高手請(qǐng)進(jìn)————matlab中的線性矩陣不等式問(wèn)題
|
|
利用matlab中Yalmip工具箱編的程序如下: %%%%%%%%%%%%%%%%%% d=57; A=[0.8 0;0 0.91]; Ad=[-0.1 0;-0.1 -0.1]; n=size(A,1); for i=1:3 e{i}=[]; for j=1:3 if j==i e{i}=[e{i} eye(n)]; else e{i}=[e{i} zeros(n)]; end end end f1=A*e{1}+Ad*e{2}; f2=f1-e{1}; f3=[e{1};(d+1)*e{3}-e{1}]; f4=[f1;(d+1)*e{3}-e{2}]; P1=sdpvar(2*n,2*n,'symmetric'); P2=sdpvar(n,n,'symmetric'); P3=sdpvar(n,n,'symmetric'); % N1=sdpvar(n,3*n); % N2=sdpvar(n,3*n); N1=1/d*[P2 -P2 zeros(n,n)]; N2=3/(d*(d+1))*[-P2 -P2 2*P2]; Delta=f4'*P1*f4-f3'*P1*f3+f2'*P2*f2-2/d*(e{1}'-e{2}')*N1+(2*d+2)/d*(e{1}'+e{2}'-2*e{3}')*N2+e{1}'*P3*e{1}-e{2}'*P3*e{2}; FAI=[Delta N1' (d+1)/sqrt(3)*N2';N1 -P2 zeros(n,n);(d+1)/sqrt(3)*N2 zeros(n,n) -P2]; con=[]; con1=[P1>=0]; con2=[P2>=0]; con3=[P3>=0]; con4=[FAI<=0]; con=[con,con1,con2, con3,con4]; sol=solvesdp(con,[],sdpsettings('solver','lmilab')) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 此時(shí)能找到可行解,如果我把 N1=1/d*[P2 -P2 zeros(n,n)]; N2=3/(d*(d+1))*[-P2 -P2 2*P2]; 這兩行程序換成如下程序 N1=sdpvar(n,3*n); N2=sdpvar(n,3*n); 則不能找到可行解了。按理說(shuō),第二種情況具有更少保守性,因?yàn)榇藭r(shí)N1和N2都是自由變量。第一種情況是第二種的特殊情形,為啥換了兩行程序,就找不到可行解了?實(shí)在沒(méi)搞明白 |
鐵蟲(chóng) (小有名氣)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085602 化工專(zhuān)碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 11/550 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書(shū)。 2026-03-28 | 3/150 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過(guò) 2026-03-27 | 6/300 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專(zhuān)業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲(chóng)不懶了 2026-03-27 | 5/250 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|