《现代控制理论》实验指导书

《《现代控制理论》实验指导书》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、《现代控制理论》实验指导书授课教师：龙威昆明理工大学机电学院2011年4月实验1系统的传递函数阵和状态空间表达式的转换[实验设备]-17-PC计算机1台（要求P4-1.8G以上），MATLAB6.X软件1套。[实验目的]1学习多变量系统状态空间表达式的建立方法、了解系统状态空间表达式与传递函数相互转换的方法；2通过编程、上机调试，掌握多变量系统状态空间表达式与传递函数的转换方法。[实验内容]1设系统的模型如式(1.1)示。(1.1)其中A为n×n维系数矩阵、B为n×m维输入矩阵C为p×n维输出矩阵，D为传递

2、阵，一般情况下为0，只有n和m维数相同时，D=1。系统的传递函数阵和状态空间表达式之间的关系如式(1.2)示。(1.2)式(1.2)中，表示传递函数阵的分子阵，其维数是p×m；表示传递函数阵的按s降幂排列的分母。2实验步骤①根据所给系统的传递函数或（A、B、C阵），依据系统的传递函数阵和状态空间表达式之间的关系如式(1.2)，采用MATLA的file.m编程。注意：ss2tf和tf2ss是互为逆转换的指令；②在MATLA界面下调试程序，并检查是否运行正确。③[例1.1]已知SISO系统的状态空间表达式为(1

3、.3)，求系统的传递函数。(1.3)★程序:%首先给A、B、C阵赋值；A=[010;001;-4-3-2];B=[1;3;-6];-17-C=[100];D=0;%状态空间表达式转换成传递函数阵的格式为[num,den]=ss2tf(a,b,c,d,u)[num,den]=ss2tf(A,B,C,D,1)☆程序运行结果:num=01.00005.00003.0000den=1.00002.00003.00004.0000从程序运行结果得到:系统的传递函数为:……………………..(1.4)①[例1.2]从系统

4、的传递函数(1.4)式求状态空间表达式。★程序:num=[0153];%在给num赋值时，在系数前补0，使num和den赋值的个数相同；den=[1234];[A,B,C,D]=tf2ss(num,den)☆程序运行结果:A=-2-3-4100010B=100C=153-17-D=0由于一个系统的状态空间表达式并不唯一,[例1.2]程序运行结果虽然不等于式(1.3)中的A、B、C阵，但该结果与式(1.3)是等效的。不防对上述结果进行验证。⑤[例1.3]对上述结果进行验证编程★程序:%将[例1.2]上述结果赋

5、值给A、B、C、D阵；A=[-2-3-4;100;010];B=[1;0;0];C=[153];D=0；[num,den]=ss2tf(A，B，C，D,1)程序运行结果与[例1.1]完全相同。[实验要求]在运行以上[例]程序的基础上，应用MATLAB对(1.5)系统仿照[例1.2]编程，求系统的A、B、C、阵；然后再仿照[例1.3]进行验证。并写出实验报告。(1.5)提示：num=[0012；0153]；实验2状态空间控制模型系统仿真及状态方程求解-17-[实验设备]PC计算机1台（要求P4-1.8G以上）

6、，MATLAB6.X软件1套。[实验目的]1、熟悉线性定常离散与连续系统的状态空间控制模型的各种表示方法。2、熟悉系统模型之间的转换功能。3、利用MATLAB对线性定常系统进行动态分析[实验内容]1、考虑由以下状态方程描述的系统求该系统状态对初始状态的时间响应。★程序:得到如下图2.1所示系统状态对初始条件的响应曲线：-17-图2.1系统状态对初始条件的相应2、考虑以下系统试给出该系统的单位阶跃响应曲线。★程序:可以得到如图2.2所示的四条单位阶跃响应：-17-[实验要求]1、进行模型间的相互转换。2、试求

7、以下系统在余弦输入信号和初始状态下的状态响应。编程实现，并给出响应曲线。-17-实验3能控能观判据及稳定性判据[实验设备]PC计算机1台（要求P4-1.8G以上），MATLAB6.X软件1套。[实验目的]1、利用MATLAB分析线性定常及离散系统的可控性与可观性。2、利用MATLAB进行线性定常及离散系统的李雅普诺夫稳定性判据。[实验内容]1、已知系统状态空间方程如下所示，对系统进行可控性、可观性分析。（1），（2）（3）★程序:(以第一题为例)（1）a=[-1-22;0-11;10-1]a=-1-220-

8、1110-1>>b=[201]'b=201-17->>c=[120]c=120>>Qc=ctrb(a,b)Qc=20001011-1rank(Qc)ans=3，系统满秩，故系统能控。rank(obsv(a,c))ans=3，系统满秩，故系统能观。（2）、（3）两题计算方法相同。2、已知系统状态空间方程描述如下：，，试判定其稳定性，并绘制出时间响应曲线来验证上述判断。★程序:A=[-10-35-50-24;1000