马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
这个是程序
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,得到
则由 得到
选择初始状态 ,闭环系统的状态轨线见怎么绘制仿真图的呐•。
|