| 查看: 1526 | 回复: 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都是自由变量。第一种情况是第二种的特殊情形，为啥换了两行程序，就找不到可行解了？实在没搞明白 |
回复此楼» 猜你喜欢
材料工程281还有调剂机会吗
已经有15人回复
-
一志愿厦大0856，306求调剂
已经有14人回复
-
271求调剂
已经有18人回复
-
山东高校教师考核超级无底线，员工过不下去啦
已经有8人回复
-
有没有学校材料专业收跨调(一志愿085410)
已经有9人回复
-
280求调剂
已经有12人回复
-
344 材料专业 求调剂211 无地域要求
已经有3人回复
-
调剂
已经有6人回复
-
290求调剂
已经有11人回复
-
0854调剂
已经有9人回复
huangmagudi
铁虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 3424.2
- 散金: 50
- 红花: 1
- 帖子: 210
- 在线: 944.8小时
- 虫号: 1693261
- 注册: 2012-03-15
- 专业: 控制理论与方法
2楼2015-12-03 15:48:02
