在lmi中计算增广矩阵时，计算结果为0，lmi的语法应该时没有问题的，而且当我改变常数矩阵A,B,C时候，LMI中会出现复数矩阵，导致无法运行

#相关代码出错

% %%%modle
% fai=[-0.9671,0.1484;-0.2968,0.5219];
% B1=[1.4142;2.2361];
% B2=[-1.4142;-2.2361];
% D1=[0.5;0];
% D2=[-0.4836;0.1484];
% E=[-1.9221,0.2236;1.9221,-0.4472];
% epsilon=0.014;
% al=1.1361;
% a2=0.6524;
% a3=0.724;
% a4=0.7288;
% I=[1,0;0,1];
%%%%model
A=[0,1;-2,-3] ;
B=[0;1] ;
C=[0,1] ;
D=0 ;
I=[1,0;0,1]; 
B1=[1.4142;2.2361];
B2=[-1.4142;-2.2361];
al=1.1361;
a2=0.6524;
a3=0.724;
a4=0.7288;
T=0.2;
tao=0.007;
epsilon=0.014;
fai=expm(A*T);
fai0=expm(A*(T/2));
syms t
N=int((expm(A*t))*B,0,T);
N1=double(N);
syms t
D1=int((expm(A*t))*B,0,T/2);
D11=double(D1);
syms t
D2=int((expm(A*t))*B,T/2,T);
D22=double(D2);
delta=norm(D11);
E=delta*fai0;
%%%%lmi model
setlmis([]) ;
Y=lmivar(1,[2 1]) ;
S=lmivar(1,[2 1]) ;
X=lmivar(1,[1 0]) ;
J=lmivar(2,[2 1]) ;
%%%%%%lmi1
lmiterm([1 1 1 Y],-al^(-2),1);
lmiterm([1 1 4 Y],1,fai');
lmiterm([1 1 5 Y],1,fai');
lmiterm([1 2 2 S],-al^(-2),1);
lmiterm([1 2 4 J],1,D1');
lmiterm([1 2 5 J],1,D1');
lmiterm([1 2 6 J],1,1);
lmiterm([1 2 7 J],1,B1');
lmiterm([1 3 3 X],-al^(-2),1);
lmiterm([1 3 4 X],1,D2');
lmiterm([1 3 5 X],1,D2');
lmiterm([1 3 7 X],1,B2');
lmiterm([1 4 4 S],-1,1);
lmiterm([1 4 4 0],epsilon*E*(E'));
lmiterm([1 4 5 0],epsilon*E*(E'));
lmiterm([1 5 5 Y],-1,1);
lmiterm([1 5 5 0],epsilon*E*(E'));
lmiterm([1 6 6 X],-1,1);
lmiterm([1 7 7 0],-epsilon*I);
%%%%%%lmi2
lmiterm([2 1 1 Y],-a2^(-2),1);
lmiterm([2 1 4 Y],1,fai');
lmiterm([2 2 2 S],1-a2^(-2),1);
lmiterm([2 3 3 X],1-a2^(-2),1);
lmiterm([2 3 4 X],1,(D1+D2)');
lmiterm([2 3 5 X],1,(B1+B2)');
lmiterm([2 4 4 Y],-1,1);
lmiterm([2 4 4 0],epsilon*E*(E'));
lmiterm([2 5 5 0],-epsilon*I);
%%%%%%lmi3
lmiterm([3 1 1 Y],-a3^(-2),1);
lmiterm([3 1 4 Y],1,fai');
lmiterm([3 2 2 S],1-a3^(-2),1);
lmiterm([3 2 4 J],1,D1');
lmiterm([3 2 5 J],1,1);
lmiterm([3 2 6 J],1,B1');
lmiterm([3 3 3 X],-al^(-2),1);
lmiterm([3 3 4 X],1,D2');
lmiterm([3 3 6 X],1,B2');
lmiterm([3 4 4 Y],-1,1);
lmiterm([3 4 4 0],epsilon*E*(E'));
lmiterm([3 5 5 X],-1,1);
lmiterm([3 6 6 0],-epsilon*I);
%%%%%%lmi4
lmiterm([4 1 1 Y],-a4^(-2),1);
lmiterm([4 1 4 Y],1,fai');
lmiterm([4 1 5 Y],1,fai');
lmiterm([4 2 2 S],-a4^(-2),1);
lmiterm([4 3 3 X],1-a4^(-2),1);
lmiterm([4 3 4 X],1,(D1+D2)');
lmiterm([4 3 6 X],1,(B1+B2)');
lmiterm([4 4 4 S],-1,1);
lmiterm([4 4 4 0],epsilon*E*(E'));
lmiterm([4 4 5 0],epsilon*E*(E'));
lmiterm([4 5 5 Y],-1,1);
lmiterm([4 5 5 0],epsilon*E*(E'));
lmiterm([4 6 6 0],-epsilon*I);
lmisys=getlmis; 
[tmin,xfeasp]=feasp(lmisys); 
YY=dec2mat(lmisys,xfeasp,Y);

SS=dec2mat(lmisys,xfeasp,S);

XX=dec2mat(lmisys,xfeasp,X);

JJ=dec2mat(lmisys,xfeasp,J);

目前有两个问题
1：算出来的增广矩阵J为0
2：当改变系数矩阵时，LMI将会出现复数矩阵，导致运行失败