线性切换系统中LMI求解以及仿真

阳光一点

这个是程序

clc;clear;
setlmis([])
A=[0.21-10;1100;-1000;000.50];
B=[00;10;00;0-1];
C=[1000;0010;0001];
D=rand(4,2)
E1=[102-1;00-1];
E2=rand(2,2)
X=lmivar(1,[4 1]);
W=lmivar(2,[2 4]);
lmiterm([1 1 1 X],A,1,'s');
lmiterm([1 1 1 W],B,1,'s');
lmiterm([1 2 1 0],D');
lmiterm([1 3 1 X],E1,1);
lmiterm([1 3 1 W],E2,1);
lmiterm([1 2 2 0],-1);
lmiterm([1 3 3 0],-1);
lmiterm([2 1 1 X],-1,1);
LMIs=getlmis;
options=[le-5,0,0,0,0];
[tmin,xfeas]=feasp(LMIs,options);
if tmin>0 stop
else
XX=dec2mat(LMIs,xfeas,X)
   WW=dec2mat(LMIs,xfeas,W)
   K=WW*inv(XX)
   x0=rand(4,1);
   x=[x0]
   tt=0;
   for i=1;1550
       t=0.01*I;tt=[tt t];
   F=[sin(t)0;0 cos(t)];
   xl=x0+0.01*(A+B*K+D*F*(E1+E2*K))*x0;
   x=[x xl];
   x0=xl;
   end
plot(tt,x)
end



这个是在MATLAB中出现错误
>> setlmis([])
A=[0.21-10;1100;-1000;000.50];
B=[00;10;00;0-1];
C=[1000;0010;0001];
D=rand(4,2)
E1=[102-1;00-1];
E2=rand(2,2)
X=lmivar(1,[4 1]);
W=lmivar(2,[2 4]);

D =

    0.8147    0.6324
    0.9058    0.0975
    0.1270    0.2785
    0.9134    0.5469


E2 =

    0.9575    0.1576
    0.9649    0.9706

>> lmiterm([1 1 1 X],A,1,'s');
lmiterm([1 1 1 W],B,1,'s');
lmiterm([1 2 1 0],D');
lmiterm([1 3 1 X],E1,1);
lmiterm([1 3 1 W],E2,1);
lmiterm([1 2 2 0],-1);
lmiterm([1 3 3 0],-1);
lmiterm([2 1 1 X],-1,1);
LMIs=getlmis;
Error using lmiterm (line 233)
lhs of LMI #1, block (1,1): incompatible dimensions in A*X*B（就是这个••••）
【例3.1】        考虑不确定线性系统（3-1-1），其中
,
,
, ，
利用MATLAB中LMI工具箱，求解LMI，得到

               

则由 得到
               
选择初始状态 ，闭环系统的状态轨线见怎么绘制仿真图的呐•。

阳光一点

我实在是找不出来错误啊~~~还有怎么设置~仿真模型啊~~~新手·新手··实在是不会啊·~~~谢谢大家啊··急救···

阳光一点

哈哈··自己找到错误啦············