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

Word格式完美整理实验报告课程线性系统理论基础实验日期年月日专业班级姓名学号同组人实验名称系统的能控性、能观测性、稳定性分析及实现评分批阅教师签字一、实验目的加深理解能观测性、能控性、稳定性、最小实现等观念。掌握如何使用MATLAB进行以下分析和实现。1、系统的能观测性、能控性分析；2、系统的稳定性分析；3、系统的最小实现。二、实验内容（1）能控性、能观测性及系统实现（a）了解以下命令的功能；自选对象模型，进行运算，并写出结果。gram,ctrb,obsv,lyap,ctrbf,obsvf,minreal；（b）已知连续系统的传递函数模型，182710)(23sssassG，Word格式完美整理当a分别取-1，0，1时，判别系统的能控性与能观测性；（c）已知系统矩阵为2101013333.06667.10666.6A，110B，201C，判别系统的能控性与能观测性；（d）求系统1827101)(23sssssG的最小实现。（2）稳定性（a）代数法稳定性判据已知单位反馈系统的开环传递函数为：)20)(1()2(100)(sssssG，试对系统闭环判别其稳定性（b）根轨迹法判断系统稳定性已知一个单位负反馈系统开环传递函数为)22)(6)(5()3()(2ssssssksG，试在系统的闭环根轨迹图上选择一点，求出该点的增益及其系统的闭环极点位置，并判断在该点系统闭环的稳定性。（c）Bode图法判断系统稳定性已知两个单位负反馈系统的开环传递函数分别为ssssGssssG457.2)(,457.2)(232231用Bode图法判断系统闭环的稳定性。（d）判断下列系统是否状态渐近稳定、是否BIBO稳定。Word格式完美整理xyuxx0525,100050250100010三、实验环境1、计算机120台；2、MATLAB6.X软件1套。四、实验原理（或程序框图）及步骤1、系统能控性、能观性分析设系统的状态空间表达式如（1-1）所示。系统的能控性、能观测性分析是多变量系统设计的基础，包括能控性、能观测性的定义和判别。系统状态能控性定义的核心是：对于线性连续定常系统（1-1）,若存在一个分段连续的输入函数u(t)，在有限的时间（t1-t0）内，能把任一给定的初态x(t0)转移至预期的终端x(t1)，则称此状态是能控的。若系统所有的状态都是能控的，则称该系统是状态完全能控的。能控性判别分为状态能控性判别和输出能控性判别。状态能控性分为一般判别和直接判别法，后者是针对系统的系数阵A是对角标准形或约当标准形的系统，状态能控性判别时不用计算，应用公式直接判断，是一种直接简易法；前者状态能控性分为一般判别是应用最广泛的一种判别法。输出能控性判别式为：pBCACABCBRankRankQncy1（2-1）Word格式完美整理状态能控性判别式为：nBAABBRankRankQnc1（2-2）系统状态能观测性的定义：对于线性连续定常系统（2-1）,如果对t0时刻存在ta，t0ta，根据[t0,ta]上的y(t)的测量值，能够唯一地确定系统在t0时刻的任意初始状态x0，则称系统在t0时刻是状态完全能观测的，或简称系统在[t0,ta]区间上能观测。状态能观测性也分为一般判别和直接判别法，后者是针对系统的系数阵A是对角标准形或约当标准形的系统，状态能观性判别时不用计算，应用公式直接判断，是一种直接简易法；前者状态能观测性分为一般判别是应用最广泛的一种判别法。状态能观测性判别式为：nCACACRankRankQTno1（2-3）系统的传递函数阵和状态空间表达式之间的有(1-2)式所示关系。已知系统的传递函数阵表述，求其满足(1-2)式所示关系的状态空间表达式，称为实现。实现的方式不唯一，实现也不唯一。其中，当状态矩阵A具有最小阶次的实现称为最小实现，此时实现具有最简形式。五、程序源代码1.(a)了解以下命令的功能；自选对象模型，进行运算，并写出结果。gram,ctrb,obsv,lyap,ctrbf,obsvf,minreal；gram:求解用状态空间表示的系统的可控或客观Gramian矩阵num=[6-0.6-0.12];den=[1-10.250.25-0.125];Word格式完美整理H=tf(num,den,'Ts',0.1)Lc=gram(ss(H),'c')H=6z^2-0.6z-0.12-------------------------------------z^4-z^3+0.25z^2+0.25z-0.125Sampletime:0.1secondsDiscrete-timetransferfunction.Lc=10.76517.87693.6759-0.00007.876910.76517.87691.83793.67597.876910.76513.9385-0.00001.83793.93852.6913Ctrb：计算矩阵可控性A=[-2.2-0.71.5-1;0.2-6.36-1.5;0.6-0.9-2-0.5;1.4-0.1-1-3.5]B=[69;46;44;84];Tc=ctrb(A,B);rank(Tc)A=-2.2000-0.70001.5000-1.00000.2000-6.30006.0000-1.50000.6000-0.9000-2.0000-0.50001.4000-0.1000-1.0000-3.5000ans=3Word格式完美整理Obsv:计算可观察性矩阵A=[-2.2-0.71.5-1;0.2-6.36-1.5;0.6-0.9-2-0.5;1.4-0.1-1-3.5]B=[69;46;44;84];C=[1234];Qo=obsv(A,C);Ro=rank(Qo)A=-2.2000-0.70001.5000-1.00000.2000-6.30006.0000-1.50000.6000-0.9000-2.0000-0.50001.4000-0.1000-1.0000-3.5000Ro=4Lyap:解lyapunov方程A=[00-6;10-11;01-6];B=[123;456;780];X=lyap(A,B)X=-3.2833-3.9000-0.1167-5.5000-8.6500-0.40000.2833-0.0000-0.0333Ctrbf:对线性系统进行能控性分解A=[00-6;10-11;01-6];B=[3;1;0];Word格式完美整理C=[001];[Abar,Bbar,Cbar,T,K]=ctrbf(A,B,C)Abar=-3.00000.0000-0.00009.4868-3.30000.95398.6189-3.13440.3000Bbar=-0.0000-0.00003.1623Cbar=-0.94350.33150T=-0.10480.3145-0.9435-0.29830.89500.33150.94870.31620K=110Obsvf:对线性系统进行能观性分解A=[-21;1-2];B=[1;0];C=[1-1];[AO,BO,CO,T,K]=obsvf(A,B,C)AO=-1.00000Word格式完美整理0.0000-3.0000BO=0.70710.7071CO=01.4142T=0.70710.70710.7071-0.7071K=10Minreal最小实现num=[11];den=[1520];sys=tf(num,den)[ABCD]=tf2ss(num,den)sys=ss(A,B,C,D);sysr=minreal(sys)sys=s+1--------------s^2+5s+20Continuous-timetransferfunction.A=-5-2010B=Word格式完美整理10C=11D=0sysr=a=x1x2x1-5-20x210b=u1x11x20c=x1x2y111d=u1y10Continuous-timestate-spacemodel.（b）已知连续系统的传递函数模型，182710)(23sssassG，当a分别取-1，0，1时，判别系统的能控性与能观测性；a=-1num=[1,-1];den=[1,10,27,18];[a,b,c,d]=tf2ss(num,den)Word格式完美整理n=length(a)Qc=ctrb(a,b)nc=rank(Qc)ifn==nc,disp('系统可控'),elsedisp('系统不可控'),endQo=obsv(a,c)no=rank(Qo)ifn==no,disp('系统可观'),elsedisp('系统不可观'),enda=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)ifn==nc,disp('系统可控'),elsedisp('系统不可控'),endQo=obsv(a,c)no=rank(Qo)ifn==no,disp('系统可观'),elsedisp('系统不可观'),enda=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)ifn==nc,disp('系统可控'),elsedisp('系统不可控'),endQo=obsv(a,c)no=rank(Qo)ifn==no,disp('系统可观'),elsedisp('系统不可观'),endWord格式完美整理矩阵为2101013333.06667.10666.6A，110B，201C，判别系统的能控性与能观测性；a=[6.666-10.6667-0.3333;101;012];b=[0;1;1];c=[102];d=0;n=length(a)Qc=ctrb(a,b)nc=rank(Qc)ifn==nc,disp('系统可控'),elsedisp('系统不可控'),endQo=obsv(a,c)no=rank(Qo)ifn==no,disp('系统可观'),elsedisp('系统不可观'),end（d）求系统1827101)(23sssssG的最小实现。num=[11];den=[1102718];G=tf(num,den);Gs=ss(G);Gm=minreal(Gs);Am=Gm.aBm=Gm.bCm=Gm.cDm=Gm.d1stateremoved.Am=3.5391-12.15405.1323-12.5391Bm=Word格式完美整理0.0606-0.2425Cm=0.25000.0625Dm=0（2）稳定性（a）代数法稳定性判据已知单位反馈系统的开环传递函数为：)20)(1()2(100)(sssssG，试对系统闭环判别其稳定性num=[00100200];den=[121200];[z,p,k]=tf2zp(num,den)z=-2p=0-20-1k=100（b）根轨迹法判断系统稳定性已知一个单位负反馈系统开环传递函数为)22)(6)(5()3()(2ssssssksG，试在系统的闭环根轨迹图上选择一点，求出该点的增益及其系统的闭环极点位置，并判断在该点系统闭环的稳定性。n1=[1,3];Word格式完美整理d1=conv([1,0],conv([1,5],conv([1,6],[1,2,2])));s1=tf(n1,d1);rlocus(s1);[k,poles]=rlocfind(s1)（c）Bode图法判断系统稳定性已知两个单位负反馈系统的开环传递函数分别为ssssGss