| 查看: 4172 | 回复: 1 | ||
[求助]
MATLAB中LMI程序求解H无穷可行解问题
|
|
MATLAB中LMI程序运行报错：Error using lmiterm (line 271)lhs of LMI #4, block (1,1): terms in diagonal blocks must be square； 程序如下： clc clear A=[0.0000000e+00 3.6301310e-21 -9.7026387e+00 -4.9849839e-05 0.0000000e+00 4.9698601e-07 0 0 0 0 0 0 0; 7.9257028e-24 0.0000000e+00 7.8831201e+03 0.0000000e+00 0.0000000e+00 0 0 0 0 0 0 0 0; 3.1226519e-07 1.9087724e-10 0.0000000e+00 1.3008627e-08 0.0000000e+00 1.5307304e-11 0.0000000e+00 0 0 0 0 0 0; -3.1226519e-07 -1.9087724e-10 0.0000000e+00 -1.3008627e-08 1.0000000e+00 -1.5307304e-11 0 0 0 0 0 0 0; -3.1895452e-23 5.2382432e-23 0.0000000e+00 -5.0419522e-07 -4.5526839e-08 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 1 0 0 0 0 0 0; 0 0 0 0 0 -0.49 -28 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 1 0; 0 0 0 0 0 0 0 0 0 0 0 0 1; 0 0 0 0 0 0 0 1 0 -1 0 -0.16 0; 0 0 0 0 0 0 0 0 1 0 -1 0 -0.16]; B=[0 0; 0 0; 0 0; 0 0; 0 1.5905035e-06; 0 0; 0.49 0; 0 0; 0 0; 0 0; 0 0; 0 0]; G=[1 0 0 0 0; 0 1 0 0 0; 0 0 1 0 0; 0 0 0 1 0; 0 0 0 0 1; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0; 0 0 0 0 0]; F=[1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; E1=[1 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; E2=[0 1; 0 0; 0 0; 0 0; 0 0; 0 0; 0 0; 0 0; 0 0; 0 0; 0 0; 0 0; 0 0]; C=[0 0 1 0 0 0 -1 0 0 0 0 0 0; 0 0 0 1 0 0 0 -1 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 0 0 ]; D=[0 0; 0 0; 1 0; 0 1]; setlmis([]) Q=lmivar(1,[13 1]); Y=lmivar(2,[2 13]); r=lmivar(1,[1 1]);% r1=lmivar(1,[1 1]);% r2=lmivar(1,[1 1]);% %strictly proper lmiterm([-1 1 1 Q],1,1); lmiterm([-2 1 1 r],1,1); lmiterm([-3 1 1 r2],1,1); lmiterm([4 1 1 Q],A,1,'s'); lmiterm([4 1 1 Y],B,1,'s'); lmiterm([4 1 1 r],1,F*F'); lmiterm([4 2 1 Q],E1,1); lmiterm([4 2 1 Y],E2,1); lmiterm([4 3 1 Q],C,1); lmiterm([4 3 1 Y],D,1); lmiterm([4 4 1 0],G'); lmiterm([4 2 2 r],-1,1); lmiterm([4 3 3 0],-1); lmiterm([4 4 4 r2],-1,1); lmisys=getlmis; [ropt xopt]=feasp(lmisys);% r=dec2mat(lmisys,xopt,r); r2=dec2mat(lmisys,xopt,r2); Q=dec2mat(lmisys,xopt,Q); Y=dec2mat(lmisys,xopt,Y); K=Y*inv(Q); 运行时报错，但是矩阵维数应该是没有问题，求大神帮助，LMI不等式如下  %@Y4C)Q1IUU4{KC@901WF(R.png |
回复此楼» 猜你喜欢
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
要不要辞职读博？
已经有6人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
理解，但是爱莫能助呀 
2楼2016-05-18 10:57:59