《实验三利用MATLAB求取状态空间模型的相似变换及其标准型、控制系统的不同状态模型实现》由会员分享，可在线阅读，更多相关《实验三利用MATLAB求取状态空间模型的相似变换及其标准型、控制系统的不同状态模型实现（13页珍藏版）》请在金锄头文库上搜索。
1、现代控制理论第一次上机实验报告实验三 利用 MATLAB 求取状态空间模型的相似变换及其标准型、控制系统的不同状态模型实现实验目的：1、通过实验掌握线性系统的对角线标准型、约旦标准型、模态标准型以及伴随矩阵标准型的表示及相应变换阵的求解；2、通过编程、上机调试，掌握系统可控性和可观测性的判别方法、系统的可控性和可观测性分解等；3、加深理解由控制系统传递函数建立能控、能观、约当标准型等不同状态模型的方法。实验要求： 1实现同一系统传递函数的状态模型是唯一的吗？ 2系统传递函数除上面三种不同状态模型实现外，常见的还有串连实现，对否？ 3对于上述系统传递函数，其输出稳态值与输入阶跃信号幅值有何关系？
2、实验步骤：1. 根据所给系统的已知条件（可自行参阅选择刘豹教材中的例题或习题） ，如传递函数、零极点模型或（A、B、C、D） ，实现状态空间模型之间的相似变换、写出其对角线标准型、约当标准型、模态标准型以及伴随矩阵标准型的表示及求解相应变换阵，采用MATLAB 的相关函数编写 m-文件。已知系统的传递函数如下： 3211()(2.5)8.50.5607.Gssss运行如下 m-文件，得到传递函数的状态空间模型：num=0 0 0 1;den=1 8.5 20 12.5;A,B,C,D=tf2ss(num,den)得到A =-8.5000 -20.0000 -12.50001.0000 0 00
3、 1.0000 0B =100C = 0 0 1D =0因此，传递函数的一个状态空间实现是123=8.520 12.51 0 00 1 0 123+100=0 0 1123G=ss(A,B,C,D);(1) 对角线标准型：计算矩阵 A 的特征值及与特征值对应的对角型变换矩阵 D 的 m-如下：V,D=eig(A)V,D=eig(A)V =-0.9798 0.9184 0.57740.1960 -0.3674 -0.5774-0.0392 0.1469 0.5774D =-5.0000 0 00 -2.5000 00 0 -1.0000由对角线标准型的变换阵 D，运行下列 m-文件的到对角线标准
4、型矩阵系数：G1=ss2ss(G,D)a = x1 x2 x3x1 -8.5 -40 -62.5x2 0.5 0 0x3 0 0.4 0b = u1x1 -5x2 0x3 0c = x1 x2 x3y1 0 0 -1d = u1y1 0Continuous-time model.由上可得，对角线标准型：123=8.540 62.50.5 0 00 0.4 0 123+500=0 0 1123对角型变换矩阵为：=5 0 00 2.5 00 0 1(2) 约旦标准型：计算矩阵 A 变换为约当标准型 J，并得到变换矩阵 V，运行下列 m-文件： V,J=jordan(A)V =2.5000 -1.6
5、667 0.1667-0.5000 0.6667 -0.16670.1000 -0.2667 0.1667J =-5.0000 0 00 -2.5000 00 0 -1.0000根据得到的约当标准型的变换矩阵 V，运行下列文件得到约当标准型的矩阵系数：G1=ss2ss(G,V)a = x1 x2 x3x1 -104 -613.6 -697.1x2 21 123.1 139.6x3 -4.2 -24.28 -27.58b = u1x1 2.5x2 -0.5x3 0.1c = x1 x2 x3y1 1 7.5 12.5d = u1y1 0Continuous-time model由上可得，约旦标准
6、型：123=104613.6697.121 123.1 139.64.224.2827.58123+2.50.50.1=1 7.5 12.5123约旦标准型的变换矩阵为：=2.5 1.66670.16670.5 0.66670,16670.1 0.26670.1667(3) 模态标准型运行以下 m-程序可得到模态标准型系数矩阵和其变换矩阵： G1,V=canon(G,modal)a = x1 x2 x3x1 -5 0 0x2 0 -2.5 0x3 0 0 -1b = u1x1 -0.825x2 -0.95x3 0.375c = x1 x2 x3y1 -0.1212 0.2807 0.4444d
7、 = u1y1 0Continuous-time model.V =-0.8250 -2.8875 -2.0625-0.9500 -5.7000 -4.75000.3750 2.8125 4.6875由上可得，模态标准型：123=5 0 00 2.5 00 0 1123+0.8250.950.375=0.12120.28070.4444123模态标准型的变换矩阵为：=0.8252.88752.06250.95 5.7 4.750.37 2.81254.6875(4) 伴随矩阵标准型运行以下 m-程序可得到伴随矩阵标准型系数矩阵和其变换矩阵： G1,V=canon(G,companion)a =
8、 x1 x2 x3x1 0 0 -12.5x2 1 0 -20x3 0 1 -8.5b = u1x1 1x2 0x3 0c = x1 x2 x3y1 0 0 1d = u1y1 0Continuous-time model.V =1.0000 8.5000 20.00000 1.0000 8.50000 0 1.0000由上可得，伴随矩阵标准型：123=0 0 12.51 0 200 1 8.5123+100=0 0 1123模态标准型的变换矩阵为：=1 8.5 200 1 8.50 0 12根据所给系统的已知条件（可自行参阅选择刘豹教材中的例题或习题） ，如（A、B、C、D）模型，判断其可控
9、性和可观测性并进行可控性和可观测性分解。判别可控、可观：(1) 构造系统的可控性判别矩阵 Tc 的 m-程序及结果如下： Tc=ctrb(A,B)Tc =1.0000 -8.5000 52.25000 1.0000 -8.50000 0 1.0000由 Tc 可得，系统可控。(2) 构造系统的可观测性判别矩阵 To 的 m-程序及结果如下： To=obsv(A,C)To =0 0 10 1 01 0 0由 To 可得，系统可观。运行以下 m-文件得到可控矩阵可观矩阵：可控矩阵： W=gram(G,c)W =0.0635 -0.0000 -0.0032-0.0000 0.0032 -0.0000
10、-0.0032 -0.0000 0.0022可观矩阵： W=gram(G,o)W =0.0022 0.0183 0.04000.0183 0.1591 0.36700.0400 0.3670 1.0294能控性分解 Ac,Bc,Cc,Tc,Kc=ctrbf(A,B,C)Ac =0 1.0000 00 0 -1.000012.5000 20.0000 -8.5000Bc =00 1Cc =-1 0 0Tc =0 0 -10 -1 01 0 0Kc =1 1 1 sum(Kc)ans =3由上可得，可控性分解子矩阵：123=0 1 00 0 112.520 8.5123+001=1 0 0123能
11、观测性分解 Ao,Bo,Co,To,Ko=obsvf(A,B,C)Ao =-8.5000 20.0000 -12.5000-1.0000 0 00 -1.0000 0Bo = -100Co =0 0 -1To =-1 0 00 1 00 0 -1Ko =1 1 1 sum(Ko)ans =3由上可得，可观性分解子矩阵：123=8.520 12.51 0 00 1 0 123+100=0 0 11233 按图 4.1 电路接线，输入阶跃信号，观察记录输出波形，观测稳态输出值(或稳态误差)和调整时间。(注意：电阻值可根据实际情况合理选取，但需尽量保证方框图中各环节的比例放大倍数。) 按图 4.2
12、图 4.3 分别接线，观察并记录两个电路相应的阶跃响应曲线，并与图4.1 所示系统阶跃响应曲线进行比较，它们是否一致？并简单解释其原因。 实验输出的参数要求及记录要求如下 4.1 仿真图4.1 仿真结果由 4.1 仿真结果图可知，稳态输出值为 0.08，调整时间为 64.2 仿真图4.2 仿真结果由 4.2 仿真结果图可知，稳态输出值为 0.08，调整时间为 6.34.3 仿真图4.3 仿真结果由 4.3 仿真结果图可知，稳态输出值为 0.078，调整时间为 7.7结论：由上可知，4.1 和 4.2、4.3 曲线变化趋势相同，但是稳态值和调节时间并不完全一致。实验要求： 1.实现同一系统传递函
13、数的状态模型是唯一的吗？ 答：不唯一。2. 系统传递函数除上面三种不同状态模型实现外，常见的还有串连实现，对否？ 答：对。3对于上述系统传递函数，其输出稳态值与输入阶跃信号幅值有何关系？答：成正比关系。输出稳态值变化比例和输入阶跃信号幅值变换比例相同。实验总结：通过此次实验，我更加深入地学习了状态空间模型的求解，及线性系统对角线标准型、约旦标准型、模态标准型、伴随矩阵标准型的表示方法，和相互之间进行变换的方法。还学习到了怎样根据系统矩阵系数构建能控、能观判别矩阵来判别系统的能控性、能观性；并进行能控性、能观性分解。同时也通过 simulink 仿真的形式对系统模型在单位阶跃输入情况下的输出有了更加形象的认识。在做实验的过程中，解决遇到的问题，不断地进行思考
