现代控制理论郑宏伟.docx

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,denTransfer function:s3 + 7 s2 + 24 s + 24-s4 + 10 s3 + 35 s2 + 50 s + 242、状态方程模型ss功能:生成状态方程,或者将零极点模型或传递函数模型转换成状态方程模型。格式:G=ss(A,B,C,D其中,

3、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 x4x1 -10 -2.188 -0.3906 -0.09375x2 16 0 0 0x3 0 8 0 0x4 0 0 2 0b =u1x1 1x2 0x3 0x4 0c =x1 x2 x3 x4y1 1 0.4375 0.1875 0.09375d =u1y1 0Continuous-time model

4、.3、零极点模型zpk功能:生成零极点模型,或将状态方程模型或传递函数模型转换成零极点模型。格式:G=zpk(z, p, K其中,z,p,K分别表示系统的零点、极点和增益。【例】:(自己举例并编程演示G=tf(-10 20 0,1 7 20 28 19 5sys=zpk(G;G=tf(-10 20 0,1 7 20 28 19 5Transfer function:-10 s2 + 20 s-s5 + 7 s4 + 20 s3 + 28 s2 + 19 s + 5 sys=zpk(GZero/pole/gain:-10 s (s-2-(s+13 (s2 + 4s + 5/(二连续时间系统离散化

5、函数名称:c2d格式:G=c2d(G1,Ts,其中Ts 为采样周期。 功能:连续时间系统离散化。要求:先进行理论求解,再与仿真结果相比较。【例】试写出连续时间系统u x x +-=102010 采样周期为T 的离散化状态方程。1、理论求解解:先求At e1111201(-+-=-=s s L A sI L e t At -=+=-t t e e s s s sL 22101(2112102(11 T t t T T G =(-=-T T e e 2201(211 d e e Bd T H T T -=-1001(211(0220d e e T -=-0221(21 021412122T e e

6、 -+=-+-=-T T e e T 222121414121 所以:(1(k u T H k x T G k x +=+(2121414121(01(2111(1(22212221k u e e T k x k x e e k x k x T T T T -+-+-=+- 2.MATLAB 仿真 A=0 1;0 -2; B=0;1; T=0.1; G1 H1=c2d(A,B,TG1 =1.0000 0.09060 0.8187H1 =0.00470.0906这里T=0.1;综上所述说明用MATLAB仿真与理论计算相同,MATLAB仿真是正确的。(三状态空间表达式的线性变换函数名称:ss2ss

7、功能:完成状态空间表达式的线性变换。格式: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(pa =x1 x2 x3x1 -1 1 0x2 1.665e-016 -1 0x3 0 0 2b =u1x1 -0.1111x2 -0.3333x3 0.1111c =x1 x2 x3y1 1 0 1d =u1y1 0Continuous-time model.五、思考题1 MATLAB中的函数其实都

8、是一些子程序,那么其ss2tf(函数是如何编写的? A=;B=;C=;D=;Sys=ss(A,B,C,D;G=tf(Sys2 在MATLAB中对连续系统进行离散化有何现实意义?用数字计算机求解连续系统方程或对连续的被控对象进行计算机控制时,由于数字计算机运算和处理均用数字量,这样就必须将连续系统方程离散化。在MA TLAB中对连续系统进行离散化,能够使得计算机能求解连续系统方程或对连续的被控对象进行控制。第二部分线性控制系统能控性、能观性和稳定性分析一、实验目的1 掌握线性控制系统能控性和能观测性的判别方法,了解不可控系统或不可观测系统的结构分解方法。2 掌握控制系统在李亚普诺夫意义下的稳定性

9、的分析方法。二、实验项目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(ob=n ,则系统状态完全可观测。(二 稳定性分析设系统的状态方程为:-=21211110x x x x 试确定系统在平衡状态处的稳定性。

10、五、实验报告要求将调试前的原程序及调试后的结果要一起写到实验报告上。举例如下: 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,BQ1 =1 0 1 -2 -2 -42 1 -5 -2 9 60 2 3 2 6 -4 Q2=obsv(A,CQ2 =1 0 00 -1 01 0 -11 2 0-2 0 -2-1 -4 -1 R1=rank(Q1R1 =3 R2=rank(Q2R2 =3从计算结果可以看出,系统能控性矩阵和能观测性矩阵的秩都是3,为满秩,因此该系统是能控的,也是能观测的。实验二 状态反馈控制系统的设计

11、第一部分 基于MATLAB 和极点配置法状态反馈控制系统的设计一、实验目的1 掌握极点配置法的基本思想。2利用MATLAB 中的函数设计状态反馈控制系统。二、实验项目运用MATLAB 和极点配置法设计状态反馈控制系统。三、实验设备与仪器1、计算机2、MATLAB 软件四、实验原理及内容1、SISO 系统极点配置 acker格式:k=acker(a,b,p说明:acker 函数可计算反馈增益矩阵K 。其中K 为行向量,p 为由期望极点构成的行向量。【例】:已知系统动态方程为ux x +-=100320100010 x y 0010=试用MATLAB 编程设计反馈增益矩阵K ,使闭环极点配置在-2

12、,-1+j ,-1-j 。 解:首先判断系统的能控性,输入以下语句A=0 1 0;0 0 1;0 -2 -3; B=0;0;1; R=rank(ctrb(A,BR =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,PK =4.0000 4.0000 1.0000(答案:k=4 4 12、MIMO系统极点配置place格式:k=place(A,B,p五、实验报告要求要从理论上分析极点配置的过程,并将调试前的原程序及调试后的结果要一起写到

13、实验报告上。六、思考题1 极点配置法的基本思想和设计思路是什么?状态反馈系统的稳定性和瞬态性能主要是由系统极点决定的。如果引入状态反馈将系统的极点配置在s左半平面的希望位置上,则可以得到满意的系统特性,一个系统引入状态反馈可以任意配置极点的条件是原系统能控。2 如何验证设计出的系统是否达到了设计要求?看配置的极点是否在s左半平面的希望位置上第二部分极点配置全状态反馈控制系统的设计一、实验目的1 学习并掌握用极点配置法来设计全状态反馈控制系统。2 用软件仿真方法研究参数对系统性能的影响。二、实验仪器与设备计算机一台、MATLAB 软件。 三、实验内容1 设计典型二阶系统的极点配置全状态反馈控制系

14、统,并进行软件仿真研究。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所示情况,软件仿真其阶跃响应。 五、实验原理1 典型二阶系统全状态反馈的极点配置设计方法

15、 (1被控对象状态方程与能控性若被控系统(A 、B 、C 完全能控,则通过状态反馈可以任意配置极点,取图5-1所示系统为实验系统。(s R (s C s110501+s .-1x 2x 图5-1 二阶实验系统结构图由图可见系统的开环传递函数为(10501+=s .s s G ,取图中21x ,x 为状态变量,将系统开环传递函数表示为被控对象状态方程S (A 、B 、C ,可以得:=+-=x y u x x 0110002020 故有:201200=Rank AB B Rank RankW c 可见状态完全能控。(2理想极点配置期望的性能指标为:超调量%p 25,峰值时间50.t p 秒。由经典

16、控制理论可知:%e/p 2521=-,选择阻尼比7070.=。s .t n p 5012-=,选择(s /n 110=于是可以得到系统的理想极点为:07707707707711.j .p ,.j .p -=+-=。系统的理想特征方程为:10014142222+=+s .s s s nn 。 (3状态反馈系数的确定加入全状态反馈后的系统特征方程为:(020*1=+=+-+=+-k k s k s k s k s BK A sI配置理想极点,则有:(100141420202021222+=+s .s k k s k s于是可以计算出:9591021.k k K -=按极点配置设计的具有全状态反馈的

17、系统结构如图5-2所示。(s R (s C s 110501+s .-1x 2x 图5-2 二阶全状态反馈实验系统结构图1k 2k系统的模拟电路图如图5-3 所示,图中的参数21x x R ,R 分别为k ,k 3318,接线时请注意反馈电路的连接。(t r k 2001x R 2x R k200k100u10k50k 50u1RR(t c +-+图5-3 二阶全状态反馈实验系统模拟电路图2 典型三阶系统全状态反馈的极点配置设计方法 (1典型三阶系统如图5-4所示。(s R (s C 55+s 22+s -s103x 2x 1x 图5-4 典型三阶实验系统结构图其开环传递函数为(25100+=

18、s s s s G闭环传递函数为(100107100123+=+=s s s s G s G s W 该闭环系统的模拟电路如图5-6所示。 图5-5 典型三阶闭环实验系统的阶跃响应曲线(t r k 200k200u1k500k 500u1R(t c +-+k200k200k200k100u1R-+图5-6 三阶输出反馈实验系统模拟电路图可以用劳斯判据判断该闭环系统是不稳定的。闭环系统的阶跃响应曲线如图5-5所示。选取图5-4中的321x ,x ,x 为状态变量,系统开环传递函数可以表示为被控对象状态方程S (A 、B 、C :=+=Cx y Bu Ax x其中0015005002200100=

19、-=C ,B ,A因为32=B A ABBRank RankW c ,所以系统状态完全能控。(2理想极点和理想闭环特征方程考虑到系统稳定性等要求,选择理想极点为:91-=S ,222j S +-=,223j S -=由此可得到理想的闭环特征方程为:072441323=+s s s(3全状态反馈系数设计取321x ,x ,x 为状态变量,带全状态的典型三阶系统结构如图5-7所示。求取加全状态反馈后的闭环特征方程,由图5-7可以得到:(010010101057132233=-+-+=+-k s k k s k s BK A sI令其与理想的闭环特征方程一致,可以求出全状态反馈系数为:2122720

20、321.k ,.k ,.k -=-=-=55+s 22+s s10(s R (s C 1x 2x 3x 1k 2k 3k 图5-7 带全状态反馈的典型三阶实验系统结构图(4全状态反馈的典型系统的模拟电路如图5-8所示,321x x x R ,R ,R 的阻值分别为k ,k ,k 15091270。(t r k 200k200u1k500k 500u1R(t c +-+k200k200k100u1R-+图5-8 三阶全状态反馈实验系统模拟电路图RR-+1x R 2x R 3x R六、思考题与实验报告要求1 思考题(1状态反馈控制器在模拟实验电路中是如何实现的? (2状态反馈控制为什么会优于输出反馈控制? 状态反馈系统方程为:X=Ax+B(V-Kx=(A-BKx+BVy=(C-DKx+DV 式(5-3输出反馈系统方程为:X=Ax+B(V-Hy=A-BH(I+DHCx+B-BH(I+DHDV y=(I+DHCx+(I+DHDV 式(5-5比较式（5-3）和式（5-5）可知，式（5-5）系数矩阵中的 HC 相当于状 态反馈系统中的 K 矩阵。由于 m=n 等原因，K 矩阵可以选择的自由度比较大， 而 H 矩阵可以选择的自由度相对 K 矩阵 来说要小些， 尤其是 HC 对改善系统性 能的效果同 K 矩阵相比要小得多，因此，输出反馈改善系统性能的能力要差些。 16