合肥工业大学系统控制仿真综合实验报告

1、合肥工业大学电气与自动化工程学院综合实验报告实 验 名 称:系统仿真综合实验姓 名:学 号:专 业 班 级:实 验 地 点:逸夫楼807指 导 教 师:都海波 殷礼胜成 绩:日 期:实验一 MATLAB基本操作实验目的1熟悉MATLAB实验环境，练习MATLAB命令、m文件、Simulink的基本操作。2利用MATLAB编写程序进行矩阵运算、图形绘制、数据处理等。3利用Simulink建立系统的数学模型并仿真求解。实验原理MATLAB环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。MATLAB有3种窗口，即：命令窗口（The Command Window）、m-文件编辑窗口（The

2、 Edit Window）和图形窗口（The Figure Window），而Simulink另外又有Simulink模型编辑窗口。1命令窗口（The Command Window）当MATLAB启动后，出现的最大的窗口就是命令窗口。用户可以在提示符“”后面输入交互的命令，这些命令就立即被执行。在MATLAB中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。在命令窗口中输入该文件名，这一连串命令就被执行了。因为这样的文件都是以“.m”为后缀，所以称为m-文件。2m-文件编辑窗口（The Edit Window）我们可以用m-文件编辑窗口来产生新的m-文件，或者编辑已经存在的m

3、-文件。在MATLAB主界面上选择菜单“File/New/M-file”就打开了一个新的m-文件编辑窗口；选择菜单“File/Open”就可以打开一个已经存在的m-文件，并且可以在这个窗口中编辑这个m-文件。3图形窗口（The Figure Window）图形窗口用来显示MATLAB程序产生的图形。图形可以是2维的、3维的数据图形，也可以是照片等。 MATLAB中矩阵运算、绘图、数据处理等内容参见教材自动控制系统计算机仿真的相关章节。Simulink是MATLAB的一个部件，它为MATLAB用户提供了一种有效的对反馈控制系统进行建模、仿真和分析的方式。有两种方式启动Simulink:1在Com

4、mand window中，键入simulink，回车。2单击工具栏上Simulink图标。启动Simulink后，即打开了Simulink库浏览器（Simulink library browser）。在该浏览器的窗口中单击“Create a new model（创建新模型）”图标，这样就打开一个尚未命名的模型窗口。把Simulink库浏览器中的单元拖拽进入这个模型窗口，构造自己需要的模型。对各个单元部件的参数进行设定，可以双击该单元部件的图标，在弹出的对话框中设置参数。 实验内容1 用MATLAB可以识别的格式输入下面两个矩阵 再求出它们的乘积矩阵C，并将C矩阵的右下角23子矩阵赋给D矩阵。赋

5、值完成后，调用相应的命令查看MATLAB工作空间的占用情况。调用相应的命令查看MATLAB工作空间的占用情况。程序A=1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4B=1+4i,4,3,6,7,8;2,3,3,5,5,4+2i;2,6+7i,5,3,4,2;1,8,9,5,4,3C=A*BD=C(4:5,4:6)结果A = 1 2 3 3 2 3 5 7 1 3 5 7 3 2 3 9 1 8 9 4B= Columns 1 through 4 1.0000 + 4.0000i 4.0000 3.0000 6.0000 2.0000 3.0000 3.0000

6、5.0000 2.0000 6.0000 + 7.0000i 5.0000 3.0000 1.0000 8.0000 9.0000 5.0000 Columns 5 through 6 7.0000 8.0000 5.0000 4.0000 + 2.0000i 4.0000 2.0000 4.0000 3.0000 C = 1.0e+002 * Columns 1 through 4 0.1400 + 0.0400i 0.5200 + 0.2100i 0.5100 0.4000 0.2500 + 0.0800i 1.0300 + 0.3500i 1.0300 0.7700 0.2400 + 0

7、.0400i 0.9900 + 0.3500i 1.0000 0.7100 0.2200 + 0.1200i 1.0800 + 0.2100i 1.1100 0.8200 0.3900 + 0.0400i 1.1400 + 0.6300i 1.0800 0.9300 Columns 5 through 60.4100 0.3100 + 0.0400i 0.7700 0.5900 + 0.0600i 0.7000 0.5100 + 0.0600i 0.7900 0.6500 + 0.0400i 0.9900 0.7000 + 0.1600iD =82.0000 79.0000 65.0000 +

8、 4.0000i 93.0000 99.0000 70.0000 +16.0000i whoYour variables are:A B C D 2 分别用for和while循环结构编写程序，求出while循环K=0;i=0;while i B=0;0;0;1; C=1 0 0 0; D=0运行仿真，在Workspace中就得到tout和yout两个变量。其中输出变量yout的前两列为系统的输出信号，后4列为系统的状态变量。在Matlab的命令窗口输入plot(tout,yout(:,1:2)图片为8 建立下图所示非线性系统的Simulink模型，并观察在单位阶跃信号输入下系统的输出曲线和误差

9、曲线。输出曲线： 误差曲线：实验二 经典控制理论实验目的以MATLAB及Simulink为工具，对控制系统进行时域、频域及根轨迹分析。实验原理1、 时域分析法是根据系统的微分方程（或传递函数），利用拉普拉斯变换直接解出动态方程，并依据过程曲线及表达式分析系统的性能。时域响应指标如图所示。延迟时间td，指响应曲线第一次达到其终值一半所需要的时间。上升时间tr，指响应曲线从终值10%上升到终值90%所需要的时间；对于有振荡的系统，也可定义为响应从零第一次上升到终值所需要的时间。上升时间是系统响应速度的一种度量。峰值时间tp，指响应超过终值达到第一个峰值所需要的时间。调节时间ts，指响应达到并保持在

10、终值5%（或2%）内所需要的时间。超调量%，指响应的最大偏离量h(tp)与终值h()之差的百分比，即：稳态误差，描述系统稳态性能的一种性能指标。2、 频域分析法通常从频率特性出发对系统进行研究。在工程分析和设计中，通常把频率特性画成一些曲线，从频率特性曲线出发进行研究。这些曲线包括幅频特性和相频特性曲线，幅相频率特性曲线，对数频率特性曲线以及对数幅相曲线等，其中以幅相频率特性曲线，对数频率特性曲线应用最广。对于最小相位系统，幅频特性和相频特性之间存在着唯一的对于关系，故根据对数幅频特性，可以唯一地确定相应的相频特性和传递函数。根据系统的开环频率特性去判断闭环系统的性能，并能较方便地分析系统参量

11、对系统性能的影响，从而指出改善系统性能的途径。3、 根轨迹是求解闭环系统特征根的图解方法。由于控制系统的动态性能是由系统闭环零极点共同决定，控制系统的稳定性由闭环系统极点唯一确定，利用根轨迹确定闭环系统的零极点在s平面的位置，分析控制系统的动态性能。实验内容*1控制系统数学模型的转换自动控制系统计算机仿真教材第4章中的所有例题答：4-1：程序 num=5*2 0 3; den=conv(conv(conv(1 0 0,3,1),conv(1 2,1,2),5 0 3 8);prinsys(mun,den,s)printsys(num,den,s) 结果num/den = 10 s2 + 15

12、- 15 s8 + 65 s7 + 89 s6 + 83 s5 + 152 s4 + 140 s3 + 32 s2 程序 tf(num,den) 结果Transfer function: 10 s2 + 15-15 s8 + 65 s7 + 89 s6 + 83 s5 + 152 s4 + 140 s3 + 32 s24-2：程序 num=12 24 0 20; den=2 4 6 2 2; A B C D=tf2ss(num,den)结果A = -2 -3 -1 -1 1 0 0 0 0 1 0 0 0 0 1 0B = 1 0 0 0C = 6 12 0 10D = 04-3:程序 A=0

13、,1,1;0,0,1;-10,-17,-8;B=0;0;1;C=5,6,1;D=0; num,den=ss2tf(A,B,C,D); G=tf(num,den) 结果Transfer function: s2 + 11 s + 5-s3 + 8 s2 + 27 s + 104-4:程序 A=0,1,1;0,0,1;-10,-17,-8;B=0;0;1;C=5,6,1;D=0;z,p,k=ss2zp(A,B,C,D);Gn=zpk(z,p,k) 结果Zero/pole/gain: (s+10.52) (s+0.4751)-(s+0.4199) (s2 + 7.58s + 23.82)4-5:程序

14、num=20,10;den=1,15,74,120;R,P,H=residue(num,den)结果R = -55.0000 90.0000 -35.0000P = -6.0000 -5.0000 -4.0000H = 4-6:程序 A=0,1,0;0,0,1;-6,-11,-6;B=1,0;2,-1;0,2;C=1,-1,0;2,1,-1;D=zeros(2);T=0.1;G=ss(A,B,C,D);Gd=c2d(G,T) 结果a = x1 x2 x3 x1 0.9991 0.0984 0. x2 -0.02458 0.9541 0.07382 x3 -0.4429 -0.8366 0.51

15、12 b = u1 u2 x1 0.1099 -0. x2 0.1959 -0.0902 x3 -0.1164 0.1936 c = x1 x2 x3 y1 1 -1 0 y2 2 1 -1 d = u1 u2 y1 0 0 y2 0 0 Sampling time: 0.1Discrete-time model.4-7： 程序A1=-9,17;-1,3;B1=0,-1;-1,0;C1=-3,2;-13,18;D1=-1,0;-1,0;sys1=ss(A1,B1,C1,D1);sys2=tf(2,1,2);sys=feedback(sys1,sys2,2,2,-1) 结果a = x1 x2 x

16、3 x1 -9 17 1 x2 -1 3 0 x3 -26 36 -2 b = u1 u2 x1 0 -1 x2 -1 0 x3 -2 0 c = x1 x2 x3 y1 -3 2 0 y2 -13 18 0 d = u1 u2 y1 -1 0 y2 -1 0 Continuous-time model.4-8：程序 A=1,2,-1;0,2,1;1,-3,2;B=0;1;1;C=1,0,1; Qc=ctrb(A,B)结果Qc = 0 1 8 1 3 5 1 -1 -10程序 syms s;det(s*eye(3)-A) 结果ans =s3 - 5*s2 + 12*s - 11 程序 if r

17、ank(Qc)=3 disp(The system is controllable)else disp(The system is uncontrollable)endThe system is controllable Q=Qc*12,-5,1;-5,1,0;1,0,0,P=inv(Q)结果Q = 3 1 0 2 -2 1 7 -6 1P = 0.2353 -0.0588 0.0588 0.2941 0.1765 -0.1765 0.1176 1.4706 -0.4706 Ab=P*A*Q,Bb=P*B,Cb=C*QAb = 0 1.0000 0.0000 0 0 1.0000 11.000

18、0 -12.0000 5.0000Bb = 0 0.0000 1.0000Cb = 10 -5 14-9：程序 num=4,5,1;den=1,6,11,6;G=tf(num,den); A,B,C,D=tf2ss(num,den); G1=ss(A,B,C,D) 结果a = x1 x2 x3 x1 -6 -11 -6 x2 1 0 0 x3 0 1 0 b = u1 x1 1 x2 0 x3 0 c = x1 x2 x3 y1 4 5 1 d = u1 y1 0 Continuous-time model.程序 G2=ss(G) 结果a = x1 x2 x3 x1 -6 -2.75 -1.5

19、 x2 4 0 0 x3 0 1 0 b = u1 x1 2 x2 0 x3 0 c = x1 x2 x3 y1 2 0.625 0.125 d = u1 y1 0 Continuous-time model.程序 G3=canon(G,modal) 结果a = x1 x2 x3 x1 -3 0 0 x2 0 -2 0 x3 0 0 -1 b = u1 x1 -15.52 x2 -19.6 x3 5.745 c = x1 x2 x3 y1 -0.7086 0.3572 8.604e-016 d = u1 y1 0 Continuous-time model.程序G4=canon(G,compa

20、nion) 结果a = x1 x2 x3 x1 0 0 -6 x2 1 0 -11 x3 0 1 -6 b = u1 x1 1 x2 0 x3 0 c = x1 x2 x3 y1 4 -19 71 d = u1 y1 02已知二阶系统闭环传递函数为：（1） 编写程序求解系统的阶跃响应；计算系统的闭环根、阻尼比、无阻尼振荡频率；（2）编程计算并在同一幅图上绘出和的阶跃响应曲线；（3）编程计算并在同一幅图上绘出和的阶跃响应曲线。（4）绘出以下系统的阶跃响应曲线，并与原系统响应曲线比较（将5条曲线画在同一幅图上） ；记录各系统的峰值和峰值时间，计算超调量。（1）sys=tf(10,1 2 10);s

21、tep(sys)（2）num=10;den1=1 2*sqrt(10) 10;den2=1 4*sqrt(10) 10;sys1=tf(num,den1);sys2=tf(num,den2);subplot(1,2,1),step(sys1),subplot(1,2,2),step(sys2)num=10;den1=1 2*sqrt(10) 10;den2=1 4*sqrt(10) 10;sys1=tf(num,den1);sys2=tf(num,den2);step(sys1,sys2)（3）num1=2.5;den1=1 1 2.5;num2=40;den2=1 4 40;sys1=tf(

22、num1,den1);sys2=tf(num2,den2);subplot(1,2,1),step(sys1),subplot(1,2,2),step(sys2) num1=2.5;den1=1 1 2.5;num2=40;den2=1 4 40;sys1=tf(num1,den1);sys2=tf(num2,den2);step(sys1,sys2)（4）num1=10;num2=2 10;num3=1 0.5 10;num4=1 0.5 0;num5=1 0;den=1 2 10;sys1=tf(num1,den);sys2=tf(num2,den);sys3=tf(num3,den);s

23、ys4=tf(num4,den);sys5=tf(num5,den);step(sys1,sys2,sys3,sys4,sys5)num1=10;num2=2 10;num3=1 0.5 10;num4=1 0.5 0;num5=1 0;den=1 2 10;sys1=tf(num1,den);sys2=tf(num2,den);sys3=tf(num3,den);sys4=tf(num4,den);sys5=tf(num5,den);subplot(2,3,1),step(sys1),subplot(2,3,2),step(sys2),subplot(2,3,3),step(sys3),su

24、bplot(2,3,4),step(sys4),subplot(2,3,5),step(sys5)从图中可以读出其峰值时间，峰值，与超调量系统峰值峰值时间超调量G(s)1.35091.049235.0913G1(s)0.99652.4809-0.3465G2(s)0.99696.9156-0.3073G3(s)400.0079342.75843.9901e+004G4(s)24.898424.84342.3898e+0033已知某控制系统的开环传递函数试绘制系统的开环频率特性曲线（奈魁斯特曲线和伯德图），判断闭环系统是否稳定？并求出系统的幅值裕量与相位裕量。clcK=1.5;den=conv(

25、conv(1 0,1 1),1 2);roots(den) sys=tf(K,den)figure(1)nyquist(sys) figure(2)margin(sys) 答:从开环频率特性曲线可以看出，其相位裕量为正的，幅值裕量也是正的，则其为稳定的。幅值裕量为45.1deg, 相位裕量为12dB4. 已知开环传递函数令k1作伯特图，应用频域稳定判据确定闭环系统的稳定性，并确定使系统获得最大相位裕度的增益k值。num=1 1;den=conv(conv(1 0,1 0),0.1 1);sys=tf(num,den);figure(1),margin(sys)figure(2),rlocus(

26、sys)答:从频域稳定判据知当K=1时系统稳定，当获得最大相位的K=3.3时可获得最大裕度。5单位负反馈系统的开环传递函数为绘制根轨迹曲线，并求出使闭环系统稳定的K值范围。使用rltool（sys）函数观察不同K值的阶跃响应情况，分析K值变化对系统响应有何影响。clcnum=-0.5 1;den=conv(conv(0.5 1,0.2 1),0.1 1);sys=tf(num,den);rlocus(sys)K,poles=rlocfind(num,den)K = 1.1727poles = -16.5032 -0.2484 + 3.6199i -0.2484 - 3.6199i使系统稳定的K

27、范围是（0，1.1727）6. 非最小相位系统的开环传递函数如下 绘制频率特性曲线，并解释为什么这样的系统被称为“非最小相位”系统，试从其频率特性加以解释。在Simulink环境中建立对应闭环系统的仿真模型，观察它们的阶跃响应。num=-6 24;den=conv(0.5 1 0 0,0.1 1);sys=tf(num,den); nyquist(sys)margin(sys)num=10 -60 110 60; den=1 17 82 130 100; sys=tf(num,den);figure(1),nyquist(sys)figure(2),margin(sys)在Simulink环境

28、中建立对应闭环系统的仿真模型：从图中可以看出第二个不是稳定的.7单位负反馈系统的开环传递函数为：，其中 利用Simulink工具给出Simulink仿真框图及阶跃响应曲线。实验三 PID控制器的设计实验目的研究PID控制器对系统的影响。实验原理1模拟PID控制器控制系统结构如图所示模拟PID控制器的表达式为：式中，、和分别为比例系数、积分系数和微分系数。 对上式进行拉普拉斯变换，整理后得到连续PID控制器的传递函数 显然，、和三个参数一旦确定（注意：，），PID控制器的性能也就确定下来。为了避免纯微分运算，通常采用近似的PID控制器，其传递函数为 2. 离散PID控制器 如果采样周期为T，在第

29、k个采样周期的导数可近似表示为； 在k个采样周期内对的积分可近似表示为； 因此，离散PID控制器的表达式为； 离散PID控制器的表达式可简化为； 离散PID控制器的脉冲传递函数为；式中，、和分别为比例系数、积分系数和微分系数。实验内容1已知三阶对象模型，利用MATLAB编写程序，研究闭环系统在不同控制情况下的阶跃响应，并分析结果。(1) 时，在不同KP值下，闭环系统的阶跃响应；num=1;den=1 3 3 1;G=tf(num,den);Gc1=1;sys1=series(G,Gc1);G1=feedback(sys1,1);figure(1),step(G1)Gc2=3;sys2=seri

30、es(G,Gc2);G2=feedback(sys2,1); figure(2),step(G2) Gc3=5;sys3=series(G,Gc3);G3=feedback(sys3,1); figure(3),step(G3) Gc4=10;sys4=series(G,Gc4);G4=feedback(sys4,1); figure(4),step(G4) Kp=1 Kp=3 Kp=5 Kp=10(2) 时，在不同值下，闭环系统的阶跃响应；num=1;den=1 3 3 1;G=tf(num,den);Gc=tf(0.5 1,0.5 0);sys1=series(G,Gc);G1=feedb

31、ack(sys1,1);figure(1),step(G1)Gc2=tf(0.8 1,0.8 0);sys2=series(G,Gc2);G2=feedback(sys2,1); figure(2),step(G2) Gc3=tf(1.2 1,1.2 0);sys3=series(G,Gc3);G3=feedback(sys3,1); figure(3),step(G3) Gc4=tf(2 1,2 0);sys4=series(G,Gc4);G4=feedback(sys4,1); figure(4),step(G4) Ti=0.5 Ti=0.8 Ti=1.2 Ti=2(3) 时，在不同值下，

32、闭环系统的阶跃响应；num=1;den=1 3 3 1;G=tf(num,den);Gc1=tf(0.01 1 1,1 0);sys1=series(G,Gc1);G1=feedback(sys1,1);figure(1),step(G1)Gc2=tf(0.1 1 1,1 0);sys2=series(G,Gc2);G2=feedback(sys2,1);figure(2),step(G2)Gc3=tf(1 1 1,1 0);sys3=series(G,Gc3);G3=feedback(sys3,1);figure(3),step(G3)Gc4=tf(2 1 1,1 0);sys4=serie

33、s(G,Gc4);G4=feedback(sys4,1);figure(4),step(G4) Td=0.01 Td=0.1 Td=1 Td=2*2. 已知被控对象为一电机模型，传递函数为，输入信号为，采用PID控制方法设计控制器，利用MATLAB编程进行仿真，绘制具有PID控制器的控制系统正弦跟踪曲线。3. 完成自动控制系统计算机仿真教材第7章中的例7-4题。（1）使用期望特性法来设计PID控制器Kp=4.4156; Ki=119.34； Kd=0.067；（2）分析比例系数Kp对控制性能的影响，在Ki=119.34 Kd=0.067保持不变情况下，Kp分别取0.5、5、20，系统响应曲线如

34、下图所示Ki=119.34;Kd=0.067;Kp=0.5; Ki=119.34;Kd=0.067;Kp=5; Ki=119.34;Kd=0.067;Kp=20;（3）分析积分系数Ki对控制性能的影响，在Kp=4.4156 Kd=0.067保持不变情况下，Ki分别取20、120、300，系统响应曲线如下图所示 Ki=20;Kd=0.067;Kp=4.4156； Ki=120;Kd=0.067;Kp=4.4156; Ki=300;Kd=0.067;Kp=4.4156;（4）分析微分系数Kd对控制性能的影响，在Kp=4.4156 Ki=119.34保持不变情况下，Kd分别取0.01、0.07、0.

35、2，系统响应曲线如下图所示 Ki=119.34;Kd=0.01;Kp=4.4156Ki=119.34;Kd=0.07;Kp=4.4156Ki=119.34;Kd=0.2;Kp=4.41564. 完成自动控制系统计算机仿真教材第7章中的例7-5题。实验四 系统状态空间设计实验目的1 学习系统的能控性、能观测性判别计算方法；2 掌握状态反馈极点配置控制系统的设计方法。3 掌握状态观察器的设计方法实验原理如果给出了对象的状态方程模型，我们希望引入状态反馈控制，使得闭环系统的极点移动到指定位置，从而改善系统的性能，这就是极点配置。状态反馈与极点配置状态反馈是指从状态变量到控制端的反馈，如图所示。 设原

36、系统状态空间表达式为：引入状态反馈后，系统的状态空间表达式为：实验内容1已知对象模型要求将闭环系统的极点配置在-1，-2，-3 j 。判断能否通过状态反馈实现极点配置，若能请设计状态反馈阵。A=0 1 0 0;0 0 -1 0;0 0 0 1;0 0 11 0;B=0;1;0;-1;r=rank(ctrb(A,B)p=-1 -2 -3+j -3-j;K=place(A,B,p)运行结果：r = 4K = -2.0000 -4.2000 -43.0000 -13.20002已知对象模型判断能否通过状态反馈配置全部极点，并将能够配置的极点均配置到-2。A=0 1 0 0;0 5 0 0;0 0 -

37、7 0;0 0 0 -8;B=1;1;3;4;rc=rank(ctrb(A,B)p=-2 -2 -2 -2;K=acker(A,B,p)运行结果：rc = 4K = -0.0714 3.0925 2.4802 -3.11543已知对象模型(1) 要求将闭环系统的极点配置到-1,-2,-3，试利用MATLAB设计控制器，并绘出闭环系统的阶跃响应曲线。(说明：用两种函数place()和acker()配置极点)(2) 设计带有状态观测器的状态反馈系统，将闭环系统的所有极点均配置到-1，观测器极点配置到-4、-5、-6。(1)A=-0.3 0.1 -0.05;1 0.1 0;-1.5 -8.9 -0.05;B=2;0;4;C=1 2 3;rC=rank(ctrb(A,B)P=-1 -2 -3;K=place(A,B,P)运行结果：rC = 3K = 3.6496 5.7189 -0.3873方法一：A=-0.3 0.1 -0.05;1 0.1 0;-1.5 -8.9 -0.05;B=2;0;4; C=1