系统的能控性能观测性稳定性分析

1、实验报告课程线性系统理论基础实验日期 年_月_日专业班级 _姓名 _ 学号 _同组人实验名称 系统的能控性、能观测性、稳定性分析及实现评分批阅教师签字一、实验目的加深理解能观测性、能控性、稳定性、最小实现等观念。掌握如何使用MATLABft行以下分析和实现。1、系统的能观测性、能控性分析；2、系统的稳定性分析；3、系统的最小实现。二、实验内容(1)能控性、能观测性及系统实现(a) 了解以下命令的功能；自选对象模型，进行运算，并写出结果gram, ctrb, obsv, lyap, ctrbf, obsvf, minreal(b)已知连续系统的传递函数模型,G(s)s as310s227s 18

2、当a分别取-1 , 0, 1时，判别系统的能控性与能观测性;6.666(c)已知系统矩阵为A 1010.66670.3333001，B1，121C 1 0 2，判别系统的能控性与能观测性;(d)求系统G(s)s10s2127s 18的最小实现(2)稳定性(a)代数法稳定性判据已知单位反馈系统的开环传递函数为:G(s)22。)，试对系统闭环判别其稳定性(b) 根轨迹法判断系统稳定性已知一个单位负反馈系统开环传递函数为G(s)，试在系统的闭环根轨迹图上选择一点，求k(s 3)s(s 5)( s 6)(s2 2s 2) 出该点的增益及其系统的闭环极点位置，并判断在该点系统闭环的稳定性。(c) Bod

3、e 图法判断系统稳定性 已知两个单位负反馈系统的开环传递函数分别为 用 Bode 图法判断系统闭环的稳定性。(d) 判断下列系统是否状态渐近稳定、是否 BIBO稳定。三、实验环境1、计算机 120 台；2、软件 1 套。四、实验原理(或程序框图)及步骤1、系统能控性、能观性分析 设系统的状态空间表达式如( 1-1 )所示。 系统的能控性、能观测性分析是多变量系统设计的基础，包括能控性、能观测性的定义和判别。 系统状态能控性定义的核心是：对于线性连续定常系统( 1-1) , 若存在一个分段连续的输入函数 u(t) ，在有限的时间( t1-t 0)内，能把任一给定的 初态 x(t 0)转移至预期的

4、终端 x(t 1) ，则称此状态是能控的。若系统所有的状 态都是能控的，则称该系统是状态完全能控的。能控性判别分为状态能控性判别和输出能控性判别。状态能控性分为一般判别和直接判别法，后者是针对系统的系数阵 A 是对角标准形或约当标准形的系统，状态能控性判别时不用计算，应用公式直接判断，是一种直接简易法；前者状态能控性分为一般判别是应用最广泛的 一种判别法。输出能控性判别式为：RankQcy Rank CB CABCAn 1Bp(2-1)状态能控性判别式为：RankQc Rank B ABAn 1 Bn(2-2)系统状态能观测性的定义：对于线性连续定常系统(2-1), 如果对 t 0时刻存在ta

5、, t0taV ，根据t 0,t a上的y(t)的测量值，能够唯一地确定系统 在to时刻的任意初始状态Xo,则称系统在to时刻是状态完全能观测的，或简 称系统在t 0,t a区间上能观测。状态能观测性也分为一般判别和直接判别法，后者是针对系统的系数阵A 是对角标准形或约当标准形的系统，状态能观性判别时不用计算，应用公 式直接判断，是一种直接简易法；前者状态能观测性分为一般判别是应用最 广泛的一种判别法。状态能观测性判别式为：RankQo RankC CACAn 1 T n(2-3)系统的传递函数阵和状态空间表达式之间的有 (1-2) 式所示关系。已知系 统的传递函数阵表述， 求其满足 (1-2

6、) 式所示关系的状态空间表达式， 称为实 现。实现的方式不唯一，实现也不唯一。其中，当状态矩阵 A具有最小阶次的 实现称为最小实现，此时实现具有最简形式。五、程序源代码1.(a) 了解以下命令的功能；自选对象模型，进行运算，并写出结果。gram, ctrb, obsv, lyap, ctrbf, obsvf, minreal ；gram:求解用状态空间表示的系统的可控或客观 Gramia n矩阵 num=6 ;den=1 -1 ;H=tf(num,den, Ts ,Lc=gram(ss(H), c )H =6 zA2 - z -zA4 - zA3 + Z2 + z -Sample time:

7、secondsDiscrete-time transfer function.Lc =Ctrb ：计算矩阵可控性A= -1; 6 ; -2 ; -1 B=6 9;4 6;4 4;8 4;Tc=ctrb(A,B);rank(Tc)A = ans =3Obsv:计算可观察性矩阵A= -1; 6 ; -2 ; -1 B=6 9;4 6;4 4;8 4;C=1 2 3 4;Qo=obsv(A,C);Ro=rank(Qo)Ro =4Lyap: 解 lyapunov 方程A=0 0 -6;1 0 -11;0 1 -6;B=1 2 3;4 5 6;7 8 0;X=lyap(A,B)X =Ctrbf: 对线性

8、系统进行能控性分解A=0 0 -6;1 0 -11;0 1 -6;B=3;1;0;C=0 0 1;Abar,Bbar,Cbar,T,K=ctrbf(A,B,C)Abar =Bbar =Cbar = 0T =1 1 0Obsvf: 对线性系统进行能观性分解A=-2 1;1 -2;B=1;0;C=1 -1;AO,BO,CO,T,K=obsvf(A,B,C)AO = 0BO =CO =0T =Minreal 最小实现num=1 1;den=1 5 20;sys=tf(num,den)A B C D=tf2ss(num,den) sys=ss(A,B,C,D);sysr=minreal(sys)sys

9、 =s + 1sA2 + 5 s + 20Continuous-time transfer function.A = -5 -20B =10C =1 1D =0sysr =x2-200a = x1x1 -5x2 1b = u1x1 1x2 0 c = x1 x2y1 11 d = u1 y1 0Con ti nu ous-time state-space model.(b)已知连续系统的传递函数模型,G(s)s as3 10s227s 18当a分别取-1 , 0, 1时，判别系统的能控性与能观测性;a=-1 num=1,-1;den=1,10,27,18;a,b,c,d=tf2ss( nu m

10、,de n)n=len gth(a)Qc=ctrb(a,b)nc=rank(Qc) if n=nc,disp( 系统可控 ), else disp( 系统不可控 ),end Qo=obsv(a,c) no=rank(Qo)if n=no,disp( 系统可观 ), else disp( 系统不可观 ),end a=0num=1,0;den=1,10,27,18; a,b,c,d=tf2ss(num,den) n=length(a)Qc=ctrb(a,b) nc=rank(Qc)if n=nc,disp( 系统可控 ), else disp( 系统不可控 ),endQo=obsv(a,c)no=

11、rank(Qo)if n=no,disp( 系统可观 ), else disp( 系统不可观 ),end a=1num=1,1;den=1,10,27,18;a,b,c,d=tf2ss(num,den) n=length(a)Qc=ctrb(a,b)nc=rank(Qc) if n=nc,disp( 系统可控 ),else disp( 系统不可控 ),endQo=obsv(a,c)no=rank(Qo)if n=no,disp( 系统可观 ),else disp( 系统不可观 ),end6.66610.66670.33330矩阵为 A 101， B1 ， C 1 0 2 ，判别系统0121的能

12、控性与能观测性；a= ;1 0 1;0 1 2;b=0;1;1;c=1 0 2;d=0;n=length(a)Qc=ctrb(a,b)n c=ra nk(Qc)if n=nc,disp(系统可控)，else disp(系统不可控),end Qo=obsv(a,c)no=ran k(Qo)if n=no,disp(系统可观),else disp(系统不可观),end的最小实现(d)求系统G(s)num=1 1;den=1 10 27 18;G=tf( nu m,de n);Gs=ss(G);Gm=mi nreal(Gs);s 132s 10s27s 18Am=Bm=Cm=Dm=1 state r

13、emoved.Am =Bm =Cm =Dm =(2)稳定性(a) 代数法稳定性判据已知单位反馈系统的开环传递函数为：G(s)10(S 2)，试对系统s(s 1)(s 20)闭环判别其稳定性num=0 0 100 200; den=1 21 20 0;z,p,k=tf2zp( nu m,de n) z =-2P =0-20100(b) 根轨迹法判断系统稳定性已知一个单位负反馈系统开环传递函数为G(s)k(s 3)s(s 5)( s 6)(s2 2s 2)，试在系统的闭环根轨迹图上选择一点，求出该点的增益及其系统的闭环极点位置，并判断在该点系统闭环的稳定性。n1= 1,3;d仁 con v(1,0

14、,co nv(1,5,co nv(1,6,1,2,2);s1=tf(n 1,d1);rlocus(s1);k,poles=rlocfi nd(s1)(c) Bode图法判断系统稳定性已知两个单位负反馈系统的开环传递函数分别为用Bode图法判断系统闭环的稳定性。G1(s) num=;den=1,5,4,0;w=logspace(-1,2,47); mag,pha=bode(num,den,w); magdB=20*log10(mag);subplot(211);semilogx(w,magdB);grid on;title(Bode Diagram); xlabel(Frequency(rad/

15、sec); ylabel(Gain dB);subplot(212);semilogx(w,pha);grid on;xlabel(Frequency(rad/sec);ylabel(phase deg)G2(s)num=;den=1,5,-4,0;w=logspace(-1,2,47);mag,pha=bode(num,den,w);magdB=20*log10(mag);subplot(211);semilogx(w,magdB);grid on;title(Bode Diagram);xlabel(Frequency(rad/sec);ylabel(Gain dB);subplot(21

16、2);semilogx(w,pha);grid on;xlabel(Frequency(rad/sec);ylabel(phase deg)BIBO稳定。(d) 判断下列系统是否状态渐近稳定、是否A=0 1 0;0 0 1;250 0 -5;B=0;0;10;C=-25 5 0;D=0;z,p,k=ss2zp(A,B,C,D)六、实验数据、结果分析( b) a=-1-10 -27 -181 0 00 1 0 b =100c =0 1 -1d =0n =Qc =1 -10 730 1 -100 0 1 nc =3系统可控Qo =0 1 -11 -1 0-11 -27 -18 no =3系统可观

17、a=0-27 -180 0-1010b =100c =0d =0n =31 0Qc =1 -10 730 1 -100 0 1 nc =3系统可控Qo =0 1 01 0 0-10 -27 -18 no =系统可观a=1-27 -180 01 0a =-1010b =100c =0d =0n =Qc =1 -10 730 1 -100 0 1nc =3系统可控Qo =0 1 11 1 0-9 -27 -18no =6.666 10.6667 0.3333 0（C）已知系统矩阵为A 101, B 10 1 2 1C 1 0 2 ，判别系统的能控性与能观测性； n =3QC =0nC =3系统可控Qo =no =系统可观的最小实现(d)求系统G