| 查看: 1496 | 回复: 1 | ||
huangmagudi铁虫 (小有名气)
|
[求助]
高手请进————matlab中的线性矩阵不等式问题
|
|
利用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')) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 此时能找到可行解，如果我把 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); 则不能找到可行解了。按理说，第二种情况具有更少保守性，因为此时N1和N2都是自由变量。第一种情况是第二种的特殊情形，为啥换了两行程序，就找不到可行解了？实在没搞明白 |
回复此楼» 猜你喜欢
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
基金申报
已经有3人回复
-
国自然申请面上模板最新2026版出了吗？
已经有9人回复
-
溴的反应液脱色
已经有6人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
推荐一本书
已经有10人回复
-
参与限项
已经有5人回复
-
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
huangmagudi
铁虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 3424.2
- 散金: 50
- 红花: 1
- 帖子: 210
- 在线: 944.8小时
- 虫号: 1693261
- 注册: 2012-03-15
- 专业: 控制理论与方法
2楼2015-12-03 15:48:02