1、现代控制理论郑宏伟 评分: 现代控制理论实验院系:计算机与电子信息学院 班级: 电气07-2 姓名: 郑宏伟 学号: 07034020234 指导老师: 禹柳飞 第一个实验 日期:5月13日第二个实验 日期:6月3日茂名学院电信学院自动化实验室实验一 线性控制系统状态空间法分析第一部分 线性控制系统状态空间模型的建立及转换一、实验目的1 掌握线性控制系统状态空间模型的建立方法。2 掌握MATLAB中的各种模型转换函数。二、实验项目1 已知系统的传递函数求取其状态空间模型。2 MATLAB中各种模型转换函数的应用。3 连续时间系统的离散化。三、实验设备与仪器1、计算机2、MATLAB软件四、实验
2、原理及内容(一)系统数学模型的建立1、传递函数模型 tf功能:生成传递函数,或者将零极点模型或状态空间模型转换成传递函数模型。格式:G=tf(num,den)其中,(num,den)分别为系统的分子和分母多项式系数向量。返回的变量G为传递函数对象。【例】:(自己举例并编程演示) num=1 7 24 24; den=1 10 35 50 24 ;G=tf(num,den) Transfer function: s3 + 7 s2 + 24 s + 24-s4 + 10 s3 + 35 s2 + 50 s + 242、状态方程模型 ss功能:生成状态方程,或者将零极点模型或传递函数模型转换成状态
3、方程模型。格式:G=ss(A,B,C,D)其中,A,B,C,D分别为状态方程的系统矩阵、输入矩阵、输出矩阵和前馈矩阵。【例】:(自己举例并编程演示)num=1 7 24 24; den=1 10 35 50 24;%分子、分母多项式 G=tf(num,den); %获得系统的传递函数模型 sys=ss(G) 语句执行结果为a = x1 x2 x3 x4 x1 -10 -2.188 -0.3906 -0.09375 x2 16 0 0 0 x3 0 8 0 0 x4 0 0 2 0 b = u1 x1 1 x2 0 x3 0 x4 0 c = x1 x2 x3 x4 y1 1 0.4375 0.
4、1875 0.09375 d = u1 y1 0 Continuous-time model.3、零极点模型 zpk功能:生成零极点模型,或将状态方程模型或传递函数模型转换成零极点模型。格式:G=zpk(z, p, K)其中,z,p,K分别表示系统的零点、极点和增益。【例】:(自己举例并编程演示)G=tf(-10 20 0,1 7 20 28 19 5)sys=zpk(G);G=tf(-10 20 0,1 7 20 28 19 5)Transfer function: -10 s2 + 20 s-s5 + 7 s4 + 20 s3 + 28 s2 + 19 s + 5 sys=zpk(G)Ze
5、ro/pole/gain: -10 s (s-2)-(s+1)3 (s2 + 4s + 5)/(二)连续时间系统离散化 函数名称:c2d格式:G=c2d(G1,Ts),其中Ts为采样周期。功能:连续时间系统离散化。要求:先进行理论求解,再与仿真结果相比较。【例】试写出连续时间系统 采样周期为T的离散化状态方程。1、理论求解解:先求 所以: 2MATLAB仿真 A=0 1;0 -2; B=0;1; T=0.1; G1 H1=c2d(A,B,T)G1 = 1.0000 0.0906 0 0.8187H1 = 0.00470.0906这里T=0.1;综上所述说明用MATLAB仿真与理论计算相同,MA
6、TLAB仿真是正确的。(三)状态空间表达式的线性变换函数名称:ss2ss功能:完成状态空间表达式的线性变换。格式:G=ss2ss(G1,inv(P) 其中inv(p)为变换阵p的逆阵。例: a=0 1 0;0 0 1;2 3 0; b=0;0;1; c=1 0 0; p=1 0 1;-1 1 2;1 -2 4; G1=ss(a,b,c,0); G=ss2ss(G1,inv(p)a = x1 x2 x3 x1 -1 1 0 x2 1.665e-016 -1 0 x3 0 0 2 b = u1 x1 -0.1111 x2 -0.3333 x3 0.1111 c = x1 x2 x3 y1 1 0
7、1 d = u1 y1 0 Continuous-time model.五、思考题1 MATLAB中的函数其实都是一些子程序,那么其ss2tf()函数是如何编写的?A=;B=;C=;D=;Sys=ss(A,B,C,D);G=tf(Sys)2 在MATLAB中对连续系统进行离散化有何现实意义? 用数字计算机求解连续系统方程或对连续的被控对象进行计算机控制时,由于数字计算机运算和处理均用数字量,这样就必须将连续系统方程离散化。在MATLAB中对连续系统进行离散化,能够使得计算机能求解连续系统方程或对连续的被控对象进行控制。第二部分 线性控制系统能控性、能观性和稳定性分析一、实验目的1 掌握线性控制
8、系统能控性和能观测性的判别方法,了解不可控系统或不可观测系统的结构分解方法。2 掌握控制系统在李亚普诺夫意义下的稳定性的分析方法。二、实验项目1运用MATLAB分析给定系统的能控性和能观测性。2 系统的结构分解。3 运用MATLAB分析分析给定系统的稳定性。三、实验设备与仪器1、计算机2、MATLAB软件四、实验原理及内容(一) 系统可控性和可观测性判别1、可控性判别 (1)可控性判别矩阵co=ctrb(a,b)或co=ctrb(G) (2)如果rank(co)=n,则系统状态完全可控。2、可观测性判别 (1)可观测性判别矩阵ob=ctrb(a,c)或ob=ctrb(G) (2)如果rank(
9、ob)=n,则系统状态完全可观测。(二) 稳定性分析设系统的状态方程为: 试确定系统在平衡状态处的稳定性。五、实验报告要求将调试前的原程序及调试后的结果要一起写到实验报告上。举例如下: A=1 0 -1;-1 -2 0;3 0 1; B=1 0;2 1;0 2; C=1 0 0;0 -1 0; Q1=ctrb(A,B)Q1 = 1 0 1 -2 -2 -4 2 1 -5 -2 9 6 0 2 3 2 6 -4 Q2=obsv(A,C)Q2 = 1 0 0 0 -1 0 1 0 -1 1 2 0 -2 0 -2 -1 -4 -1 R1=rank(Q1)R1 = 3 R2=rank(Q2)R2 =
10、3从计算结果可以看出,系统能控性矩阵和能观测性矩阵的秩都是3,为满秩,因此该系统是能控的,也是能观测的。 实验二 状态反馈控制系统的设计第一部分 基于MATLAB和极点配置法状态反馈控制系统的设计一、实验目的1 掌握极点配置法的基本思想。2利用MATLAB中的函数设计状态反馈控制系统。二、实验项目运用MATLAB和极点配置法设计状态反馈控制系统。三、实验设备与仪器1、计算机2、MATLAB软件四、实验原理及内容1、SISO系统极点配置 acker格式:k=acker(a,b,p)说明:acker函数可计算反馈增益矩阵K。其中K为行向量,p为由期望极点构成的行向量。【例】:已知系统动态方程为 试
11、用MATLAB编程设计反馈增益矩阵K,使闭环极点配置在-2,-1+j,-1-j。 解:首先判断系统的能控性,输入以下语句A=0 1 0;0 0 1;0 -2 -3; B=0;0;1; R=rank(ctrb(A,B)R = 3这说明系统能控性矩阵满秩,系统能控,可以应用状态反馈,任意配置极点。A=0 1 0;0 0 1;0 -2 -3; B=0;0;1; C=10 0 0; P=-2 -1+j -1-j; K=place(A,B,P)K = 4.0000 4.0000 1.0000(答案:k=4 4 1) 2、MIMO系统极点配置 place格式:k=place(A,B,p)五、实验报告要求
12、要从理论上分析极点配置的过程,并将调试前的原程序及调试后的结果要一起写到实验报告上。六、思考题1 极点配置法的基本思想和设计思路是什么?状态反馈系统的稳定性和瞬态性能主要是由系统极点决定的。如果引入状态反馈将系统的极点配置在s左半平面的希望位置上,则可以得到满意的系统特性,一个系统引入状态反馈可以任意配置极点的条件是原系统能控。2 如何验证设计出的系统是否达到了设计要求?看配置的极点是否在s左半平面的希望位置上第二部分 极点配置全状态反馈控制系统的设计一、实验目的1 学习并掌握用极点配置法来设计全状态反馈控制系统。2 用软件仿真方法研究参数对系统性能的影响。二、实验仪器与设备计算机一台、MAT
13、LAB软件。三、实验内容1 设计典型二阶系统的极点配置全状态反馈控制系统,并进行软件仿真研究。2 设计典型三阶系统的极点配置全状态反馈控制系统,并进行软件仿真研究。四、实验步骤1 典型二阶系统(1)对一已知二阶系统(图5-1)用极点配置方法设计全反馈系数。(2)参照图5-2,图5-3,软件仿真其阶跃响应。(3)改变系统电路,使系统恢复到图5-1所示情况,软件仿真其阶跃响应。(4)对实验结果进行比较、分析,并完成实验报告。2 典型三阶系统(1)对一已知三阶系统(图5-4)用极点配置方法设计全反馈系数。(2)参照图5-7,图5-8,软件仿真其阶跃响应。(3)改变系统电路,使系统恢复到图5-6所示情
14、况,软件仿真其阶跃响应。五、实验原理1 典型二阶系统全状态反馈的极点配置设计方法(1)被控对象状态方程与能控性若被控系统(A、B、C)完全能控,则通过状态反馈可以任意配置极点,取图5-1所示系统为实验系统。由图可见系统的开环传递函数为,取图中为状态变量,将系统开环传递函数表示为被控对象状态方程S(A、B、C),可以得:故有: 可见状态完全能控。(2)理想极点配置期望的性能指标为:超调量,峰值时间秒。由经典控制理论可知:,选择阻尼比。,选择于是可以得到系统的理想极点为:。系统的理想特征方程为:。(3)状态反馈系数的确定加入全状态反馈后的系统特征方程为:配置理想极点,则有:于是可以计算出: 按极点
15、配置设计的具有全状态反馈的系统结构如图5-2所示。系统的模拟电路图如图5-3 所示,图中的参数分别为,接线时请注意反馈电路的连接。2 典型三阶系统全状态反馈的极点配置设计方法(1)典型三阶系统如图5-4所示。其开环传递函数为闭环传递函数为该闭环系统的模拟电路如图5-6所示。图5-5 典型三阶闭环实验系统的阶跃响应曲线可以用劳斯判据判断该闭环系统是不稳定的。闭环系统的阶跃响应曲线如图5-5所示。选取图5-4中的为状态变量,系统开环传递函数可以表示为被控对象状态方程S(A、B、C):其中因为,所以系统状态完全能控。(2)理想极点和理想闭环特征方程考虑到系统稳定性等要求,选择理想极点为:, 由此可得
16、到理想的闭环特征方程为:(3)全状态反馈系数设计取为状态变量,带全状态的典型三阶系统结构如图5-7所示。求取加全状态反馈后的闭环特征方程,由图5-7可以得到:令其与理想的闭环特征方程一致,可以求出全状态反馈系数为:(4)全状态反馈的典型系统的模拟电路如图5-8所示,的阻值分别为。六、思考题与实验报告要求1 思考题(1)状态反馈控制器在模拟实验电路中是如何实现的?(2)状态反馈控制为什么会优于输出反馈控制?状态反馈系统方程为:X=Ax+B(V-Kx)=(A-BK)x+BV y=(C-DK)x+DV 式(5-3)输出反馈系统方程为:X=Ax+B(V-Hy)=A-BH(I+DH)Cx+B-BH(I+DH)DVy=(I+DH)Cx+(I+DH)DV 式(5-5) 比较式(5-3)和式(5-5)可知,式(5-5)系数矩阵中的HC相当于状态反馈系统中的K矩阵。由于m=n等原因,K矩阵可以选择的自由度比较大,而H矩阵可以选择的自由度相对 K矩阵 来说要小些,尤其是HC对改善系统性能的效果同K矩阵相比要小得多,因此,输出反馈改善系统性能的能力要差些。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1