LMI求助、论文的仿真程序

下面是我自己编的程序，之前一点matlab基础也没有，真的不知道该怎么编了 希望能够得到大家的帮助 这个程序运行之后W/M/F/E都是0，肯定是不对哒
A1=[1 0.2955;0 0.9704];
B0=[0.002 0.0198];
B1=[0.0025 0.0098];
pi11=[1 0;0 0];
pi13=[0 0;0 1];
pi21=[0 0;0 0];
pi23=[1 0;0 1];
pi31=[0 0;0 1];
pi33=[1 0;0 0];
pi41=[0 0;0 0];
pi43=[1 0;0 1];
b1=1.1361;b2=0.7247;
setlmis([]);
P=lmivar(1,[2 1]);
Q=lmivar(1,[2 1]);
R=lmivar(1,[1 1]);
J=lmivar(1,[2 1]);
W=lmivar(2,[2 2]);
M=lmivar(2,[2 2]);
F=lmivar(2,[1 2]);
E=lmivar(2,[1 2]);
lmiterm([1 1 1 P],-1,1);
lmiterm([1 1 5 P],A1,1);
lmiterm([1 1 6 E],B0',1);
lmiterm([1 1 7 R],B1',1);
lmiterm([1 1 8 F],B0',1);
lmiterm([1 2 2 Q],-1,1);
lmiterm([1 2 6 W],1,1);
lmiterm([1 2 8 M],1,1);
lmiterm([1 3 3 R],-1,1);
lmiterm([1 3 6 E],1,1);
lmiterm([1 3 8 F],1,1);
lmiterm([1 4 4 J],-1,1);
lmiterm([1 4 5 P],pi11*A1,1);
lmiterm([1 4 6 E],pi11*B0',1);
lmiterm([1 4 7 R],pi11*B1',1);
lmiterm([1 4 8 F],pi11*B0',1)
lmiterm([1 4 8 J],pi13,1)
lmiterm([1 5 5 P],-b1^(-2),1);
lmiterm([1 6 6 Q],-b1^(-2),1);
lmiterm([1 7 7 R],-b1^(-2),1);
lmiterm([1 8 8 J],-b1^(-2),1);
lmiterm([2 1 1 P],-1,1);
lmiterm([2 1 5 P],A1,1);
lmiterm([2 1 6 E],B0',1);
lmiterm([2 1 7 R],B1',1);
lmiterm([2 1 8 F],B0',1);
lmiterm([2 2 2 Q],-1,1);
lmiterm([2 2 6 W],1,1);
lmiterm([2 2 8 M],1,1);
lmiterm([2 3 3 R],-1,1);
lmiterm([2 3 6 E],1,1);
lmiterm([2 3 8 F],1,1);
lmiterm([2 4 4 J],-1,1);
lmiterm([2 4 5 P],pi21*A1,1);
lmiterm([2 4 6 E], pi21*B0',1);
lmiterm([2 4 7 R], pi21*B1',1);
lmiterm([2 4 8 F], pi21*B0',1)
lmiterm([2 4 8 J], pi23,1)
lmiterm([2 5 5 P],-b2^(-2),1);
lmiterm([2 6 6 Q],-b2^(-2),1);
lmiterm([2 7 7 R],-b2^(-2),1);
lmiterm([2 8 8 J],-b2^(-2),1);
lmiterm([3 1 1 P],-1,1);
lmiterm([3 1 5 P],A1,1);
lmiterm([3 1 6 E],B0',1);
lmiterm([3 1 7 R],B1',1);
lmiterm([3 1 8 F],B0',1);
lmiterm([3 2 2 Q],-1,1);
lmiterm([3 2 6 W],1,1);
lmiterm([3 2 8 M],1,1);
lmiterm([3 3 3 R],-1,1);
lmiterm([3 3 6 E],1,1);
lmiterm([3 3 8 F],1,1);
lmiterm([3 4 4 J],-1,1);
lmiterm([3 4 5 P], pi31*A1,1);
lmiterm([3 4 6 E], pi31*B0',1);
lmiterm([3 4 7 R], pi31*B1',1);
lmiterm([3 4 8 F], pi31*B0',1)
lmiterm([3 4 8 J], pi33,1)
lmiterm([3 5 5 P],-b1^(-2),1);
lmiterm([3 6 6 Q],-b1^(-2),1);
lmiterm([3 7 7 R],-b1^(-2),1);
lmiterm([3 8 8 J],-b1^(-2),1);
lmiterm([4 1 1 P],-1,1);
lmiterm([4 1 5 P],A1,1);
lmiterm([4 1 6 E],B0',1);
lmiterm([4 1 7 R],B1',1);
lmiterm([4 1 8 F],B0',1);
lmiterm([4 2 2 Q],-1,1);
lmiterm([4 2 6 W],1,1);
lmiterm([4 2 8 M],1,1);
lmiterm([4 3 3 R],-1,1);
lmiterm([4 3 6 E],1,1);
lmiterm([4 3 8 F],1,1);
lmiterm([4 4 4 J],-1,1);
lmiterm([4 4 5 P], pi41*A1,1);
lmiterm([4 4 6 E], pi41*B0',1);
lmiterm([4 4 7 R], pi41*B1',1);
lmiterm([4 4 8 F], pi41*B0',1)
lmiterm([4 4 8 J], pi43,1)
lmiterm([4 5 5 P],-b2^(-2),1);
lmiterm([4 6 6 Q],-b2^(-2),1);
lmiterm([4 7 7 R],-b2^(-2),1);
lmiterm([4 8 8 J],-b2^(-2),1);
lmiterm([-5 1 1 P],1,1);
lmiterm([-6 1 1 Q],1,1);
lmiterm([-7 1 1 R],1,1);
lmiterm([-8 1 1 J],1,1);
lmisys=getlmis;
[tmin,xfeas]=feasp(lmisys)
P=dec2mat(lmisys,xfeas,P);
Q=dec2mat(lmisys,xfeas,Q);
R=dec2mat(lmisys,xfeas,R);
J=dec2mat(lmisys,xfeas,J);
W=dec2mat(lmisys,xfeas,W);
M=dec2mat(lmisys,xfeas,M);
F=dec2mat(lmisys,xfeas,F);
E=dec2mat(lmisys,xfeas,E);