线性系统理论Matlab实践仿真报告.doc

《线性系统理论Matlab实践仿真报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、线性系统理论Matlab实验报告1、本题目是在已知状态空间描述的情况下要求设计一个状态反馈控制器，从而使得系统具有实数特征根，并要求要有一个根的模值要大于5，而特征根是正数是系统不稳定，这样的设计是无意义的，故而不妨设采用状态反馈后的两个期望特征根为-7，-9，这样满足题目中所需的要求。（1）要对系统进行状态反馈的设计首先要判断其是否能控，即求出该系统的能控性判别矩阵，然后判断其秩，从而得出其是否可控；判断能控程序设计如下：>>A=[-0.80.02;-0.020];B=[0.051;0.0010]

2、;Qc=ctrb(A,B)Qc=0.05001.0000-0.0400-0.80000.00100-0.0010-0.0200Rc=rank(Qc)Rc=2Qc=0.05001.0000-0.0400-0.80000.00100-0.0010-0.0200得出结果能控型判别矩阵的秩为2，故而该系统是完全可控的，故可以对其进行状态反馈设计。（2）求取状态反馈器中的K，设的期望特征根为-7，-9；其设计程序如下：>>A=[-0.80.02;-0.020];B=[0.051;0.0010];P=[-7-9

3、];k=place(A,B,P)k=1.0e+003*-0.02009.00000.0072-0.4500程序中所求出的k即为所求状态反馈控制器的状态反馈矩阵，即由该状态反馈矩阵所构成的状态反馈控制器能够满足题目要求。2、（a）要求求该系统的能控型矩阵，并验证该系统是不能控的。设计程序：>>A=[01000;-0.1-0.5000;0.50000;001000;0.51000];>>B=[0;1;0;0;0];>>C=[00010];>>Qc=ctrb(A,B)Qc=01.0000-0.50000.

4、1500-0.02501.0000-0.50000.1500-0.0250-0.0025000.5000-0.25000.07500005.0000-2.500001.00000-0.10000.0500>>Rc=rank(Qc)Rc=4从程序运行的结果可得，系统能控型判别矩阵的秩为4，而系统为5阶系统，故而就验证了该系统为不可控的。而该系统的能控型矩阵就为程序中的Qc矩阵。（b）要求求出其传递函数，并建立新的状态变量模型。设计程序：>>A=[01000;-0.1-0.5000;0.50000;00

5、1000;0.51000];B=[0;1;0;0;0];C=[00010];I=eye(5);symss;E=s*I-A;F=collect(inv(E));G=C*F*BG=50/(10*s^4+5*s^3+s^2)>>num=[50];>>den=[105100];>>Gf=tf(num,den)Transferfunction:50--------------------10s^4+5s^3+s^2>>sys=ss(Gf)a=x1x2x3x4x1-0.5-0.400x20.25000x300.

6、12500x4000.250b=u1x132x20x30x40c=x1x2x3x4y100020d=u1y10Continuous-timemodel.结果分析：可以看到所求得的新系统为4阶系统，其原因就在于原系统为不可控的，并且其能控型判别矩阵的秩为4，从而其传递函数表现的也只是可控部分的4个变量，从而导致了新系统为4阶系统。（c）要求证明其上所求的新系统为可控的，只需求其新系统能控型判别矩阵的秩，看是否为4即可。设计程序：>>A1=[-0.5-0.400;0.25000;00.12500;000

7、.250];>>B1=[32;0;0;0];>>C1=[00020];>>Qc1=ctrb(A1,B1)Qc1=32.0000-16.00004.8000-0.800008.0000-4.00001.2000001.0000-0.50000000.2500>>Rc1=rank(Qc1)Rc1=4其程序运行结果为能控型判别矩阵的秩为4，从而说明该新系统是完全能控的。其实在求得其能控矩阵时就可以判断出其为能控的，因为该矩阵为上三角矩阵，即为满秩矩阵，即秩为4，也即说明该新系统为完全可控。（d）要求判断系

8、统的稳定性，可以采用李雅普洛夫特征值法进行判定；设计程序：A=[01000;-0.1-0.5000;0.50000;001000;0.51000];>>d=eig(A)d=000-0.2500+0.1936i-0.2500-0.1936i从求得的结果中可以看出其特征值的根的实部都不是正数，从而就说明了该系统是李雅普洛夫意义下稳定的。（c）讨论状态变量模型的能控性元复杂性的关系。很直观地讲，一个系统要能控，必须要其能控型判别矩阵的秩等于系统的阶数也即就是状态变量的数目，