关于线性矩阵不等式的出错问题

fang_317

哪位高手能帮我看看这个问题!
clear
close
clc
A=[-2.73 4.26;2.08 -3.12];
B=[1.27 0.72;0.45 1.18];
C=[1.20 2.50;0.80 1.80];
G=[1 2;0 5];
E1=[0.16 0 0.1 0 0 0;0 0.16 0 0.1 0 0];
E2=[0 0 0 0 0.12 0;0 0 0 0 0 0.15];
F1=[1 0;0 1;0 0;0 0;1 0;0 1];
F2=[0 0;0 0;1 0;0 1;0 0;0 0];
L=[0 0;0 1];
I=[1 0;0 1];
I1=[1 0 0 0 0 0;
    0 1 0 0 0 0;
    0 0 1 0 0 0;
    0 0 0 1 0 0;
    0 0 0 0 1 0;
    0 0 0 0 0 1];

setlmis([])
X=lmivar(1,[2 1]);
Q=lmivar(1,[2 1]);
r=lmivar(1,[1 0]);
e1=lmivar(1,[1 0]);
e2=lmivar(1,[1 0]);

lmiterm([1 1 1 X],1,A','s');
lmiterm([1 1 1 Q],B*L,1,'s');
lmiterm([1 1 1 e1],1,E1*E1');
lmiterm([1 1 2 0],G);
lmiterm([1 2 2 r],-1,I);
lmiterm([1 3 1 X],C,1);
lmiterm([1 3 2 0],0);
lmiterm([1 3 3 r],-1,I);
lmiterm([1 3 3 e2],1,E2*E2');
lmiterm([1 4 1 X],F1,1);
lmiterm([1 4 1 Q],F2*L,1);
lmiterm([1 4 2 0],0);
lmiterm([1 4 3 0],0);
lmiterm([1 4 4 e2],-1,I1);
lmiterm([1 5 1 X],F1,1);
lmiterm([1 5 1 Q],F2*L,1);
lmiterm([1 5 2 0],0);
lmiterm([1 5 3 0],0);
lmiterm([1 5 4 0],0);
lmiterm([1 5 5 e1],-1,I1);

Xpos=newlmi;
lmiterm([-Xpos 1 1 X],1,1);
rpos=newlmi;
lmiterm([-rpos 1 1 r],1,1);
e1pos=newlmi;
lmiterm([-rpos 1 1 e1],1,1);
e2pos=newlmi;
lmiterm([-rpos 1 1 e2],1,1);
uncertainHinf2=getlmis;

c=mat2dec(uncertainHinf2,eye(2),eye(2),1,1,1);
options=[1e-5,0,0,0,0];
[copt,xopt]=mincx(uncertainHinf2,c,options);
X=dec2mat(uncertainHinf2,xopt,X)
Q=dec2mat(uncertainHinf2,xopt,Q)
r=dec2mat(uncertainHinf2,xopt,r)
e1=dec2mat(uncertainHinf2,xopt,e1)
e2=dec2mat(uncertainHinf2,xopt,e2)

K=Q*inv(X)
最后就是显示不了结果!显示
Result:  could not establish feasibility nor infeasibility

The LMI constraints were found INFEASIBLE

??? Error using ==> dec2mat
DECVARS must be a vector of length 9

Error in ==> D:\Matlab\work\uncertainHinf2.m
On line 94  ==> X=dec2mat(uncertainHinf2,xopt,X)
不知道什么原因,用feasp的时候就好用为什么换成mincx的时候就出现错误!请高手指点指点,多谢

[ 本帖最后由 eight 于 2007-9-2 21:14 编辑 ]

eight

原帖由 fang_317 于 2007-9-2 11:02 发表

                               
登录/注册后可看大图

哪位高手能帮我看看这个问题!
clear
close
clc
A=[-2.73 4.26;2.08 -3.12];
B=[1.27 0.72;0.45 1.18];
C=[1.20 2.50;0.80 1.80];
G=[1 2;0 5];
E1=[0.16 0 0.1 0 0 0;0 0.16 0 0.1 0 0];
E2=[0 0 0  ...

根据出错地方自己调试，详细描述请看看置顶帖：聚宝盆

Angel-wing

我也遇到同样的问题
不知道你解决了没有.
请与我联系
QQ:15709027

草三

楼上你好 我也遇到了相同的问题 请问怎么解决的？ 特别着急  加我qq好吗？649013663

happy

矩阵不等式没有可行解

相关可能原因及解决办法如下：
1）、先检查以下LMI有没有问题，比如推导过程、输入到MATLAB，也许是打字打错了。如果还不行，可能是理论有问题。不是所有正确的推导都有解出现；
2）、程序可以运行，说明语法没有问题。有可能需要检查一下LMI，检查矩阵的主对角线元素，若对角线元素有严格大于0存在，肯定是无可行解的；
3）、一般可行性问题的程序都有多个LMI，可逐个检查，一个LMI一个的检查是否可行，找到到底是哪个LMI不可行，再想办法解决。