声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3807|回复: 2

[控制理论] 线性切换系统中LMI求解以及仿真

[复制链接]
发表于 2012-5-18 22:19 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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,得到

               

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

回复
分享到:

使用道具 举报

 楼主| 发表于 2012-5-18 22:20 | 显示全部楼层
我实在是找不出来错误啊~~~还有怎么设置~仿真模型啊~~~新手·新手··实在是不会啊·~~~谢谢大家啊··急救···
 楼主| 发表于 2012-5-21 22:47 | 显示全部楼层
哈哈··自己找到错误啦············
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 15:16 , Processed in 0.064974 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表