| 查看: 1516 | 回复: 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都是自由变量。第一种情况是第二种的特殊情形，为啥换了两行程序，就找不到可行解了？实在没搞明白 |
回复此楼» 猜你喜欢
材料专硕英一数二306
已经有6人回复
-
265求调剂
已经有4人回复
-
265求调剂
已经有7人回复
-
26调剂/材料/英一数二/总分289/已过A区线
已经有8人回复
-
281求调剂（0805）
已经有21人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有20人回复
-
279分求调剂 一志愿211
已经有6人回复
-
085600材料与化工 求调剂
已经有17人回复
-
能源材料化学课题组招收硕士研究生8-10名
已经有13人回复
-
317求调剂
已经有6人回复
huangmagudi
铁虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 3424.2
- 散金: 50
- 红花: 1
- 帖子: 210
- 在线: 944.8小时
- 虫号: 1693261
- 注册: 2012-03-15
- 专业: 控制理论与方法
2楼2015-12-03 15:48:02
