matlab的lmi仿真例子,LMI工具箱编程求助

matlab的lmi仿真例子,LMI工具箱编程求助

本人初学MATLAB的LMI工具箱，请问哪位大神指导下面的程序为什么没有解呢？是编程出问题了吗？这个是一篇文章里的例子，应该是有解的啊。谢谢大家了！

A1=[0 1.0000 0 0; 1.9600 0 -0.4000 -0.6000; 0 0 0 1.0000; 1.9600 0 -2.4000 -3.6000];

A2=[0 1.0000 0 0; 0 0 -0.4000 -0.6000; 0 0 0 1.0000; 0 0 -2.4000 -3.6000];

B1=[0; 0; 0; -1];

B2=[0; 0; 0; -1];

x1=[1.0000 0 0 0]';

setlmis([]);%初始化

%定义矩阵变量

X=lmivar(2,[4,4]);

M1=lmivar(2,[1,4]);

M2=lmivar(2,[1,4]);

K1=lmivar(2,[1,4]);

K2=lmivar(2,[1,4]);

P=lmivar(2,[4,4]);

%描述第一个矩阵

lmiterm([1 1 1 X],A1,1,'s');

lmiterm([1 1 1 M1],B1,1,'s');

%描述第二个矩阵

lmiterm([2 1 1 X],A2,1,'s');

lmiterm([2 1 1 M2],B2,1,'s');

%描述第三个矩阵

lmiterm([3 1 1 X],A1,1,'s');

lmiterm([3 1 1 M2],B1,1,'s');

lmiterm([3 1 1 X],A2,1,'s');

lmiterm([3 1 1 M1],B2,1,'s');

lmisys=getlmis;

[r,xfeas]=gevp(lmisys,3,[0, 0, 0, 0, 1]);

X=dec2mat(lmisys, xfeas, X);

M1=dec2mat(lmisys, xfeas, M1);

M2=dec2mat(lmisys, xfeas, M2);

P=inv(X);

K1=M1*P;

K2=M2*P;