阳光一点 发表于 2012-5-18 22:19

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

这个是程序

clc;clear;
setlmis([])
A=;
B=;
C=;
D=rand(4,2)
E1=;
E2=rand(2,2)
X=lmivar(1,);
W=lmivar(2,);
lmiterm(,A,1,'s');
lmiterm(,B,1,'s');
lmiterm(,D');
lmiterm(,E1,1);
lmiterm(,E2,1);
lmiterm(,-1);
lmiterm(,-1);
lmiterm(,-1,1);
LMIs=getlmis;
options=;
=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=
   tt=0;
   for i=1;1550
       t=0.01*I;tt=;
   F=;
   xl=x0+0.01*(A+B*K+D*F*(E1+E2*K))*x0;
   x=;
   x0=xl;
   end
plot(tt,x)
end



这个是在MATLAB中出现错误
>> setlmis([])
A=;
B=;
C=;
D=rand(4,2)
E1=;
E2=rand(2,2)
X=lmivar(1,);
W=lmivar(2,);

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(,A,1,'s');
lmiterm(,B,1,'s');
lmiterm(,D');
lmiterm(,E1,1);
lmiterm(,E2,1);
lmiterm(,-1);
lmiterm(,-1);
lmiterm(,-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,得到

               

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

阳光一点 发表于 2012-5-18 22:20

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

阳光一点 发表于 2012-5-21 22:47

哈哈··自己找到错误啦············
页: [1]
查看完整版本: 线性切换系统中LMI求解以及仿真