控制系统的时域分析.ppt

《控制系统的时域分析.ppt》由会员分享，可在线阅读，更多相关《控制系统的时域分析.ppt（49页珍藏版）》请在三一办公上搜索。

1、控制系统的时域分析,本章主要内容,原理要点12.1控制系统的动态性能指标分析12.1.1 控制系统的动态性能指标控制系统动态性能指标MATLAB求取实例12.2控制系统的稳态性能指标分析系统的稳态性能指标12.2.2 控制系统稳态性能指标MATLAB求取示例12.3MATLAB时域响应仿真的典型函数应用MATLAB时域响应仿真的典型函数MATLAB时域响应仿真的典型函数应用实例12.4MATLAB图形化时域分析MATLAB LTI Viewer应用实例,通常在系统阶跃响应曲线上来定义系统动态性能指标。因此，在用MATLAB求取系统动态性能指标之前，首先给出单位阶跃响应函数step的用法。对给定

2、系统，可使用表12.2所列函数调用方式得到系统阶跃响应。,控制系统动态性能指标MATLAB求取实例,表12.2系统阶跃响应函数用法表,表12.3 状态空间方程表示的系统阶跃响应函数用法表,对于状态空间方程表示的系统，也可直接使用step函数，其用法见表12.3。,例1：设单位负反馈系统的开环传递函数为：试求系统单位阶跃响应。,num=0.3 1;den=1 0.5 0;G=tf(num,den);G0=feedback(G,1)%得到反馈系统 Transfer function:0.3 s+1-s2+0.8 s+1 step(G0)%直接得到系统单位阶跃响应曲线 y,t=step(G0);%返

3、回系统单位阶跃响应曲线参数 plot(t,y)%由plot函数绘制单位阶跃响应曲线,(a)直接绘制结果(b)返回参数调用其它函数绘制结果,图1 例1运行结果,例2：求如下系统的单位阶跃响应。,A=0 1;-25-4;B=1 1;0 1;C=eye(2);D=zeros(2);%给定系统 step(A,B,C,D)%求系统阶跃响应曲线,图2 系统阶跃响应曲线,step(A,B,C,D,1)%得出第一路输入的响应曲线 title(输入=u_1的阶跃响应)step(A,B,C,D,2)%得到第二路输入的响应曲线 title(输入=u_2的阶跃响应),图3 第1路输入的响应曲线,图4 第2路响应曲线,

4、例3：单位负反馈系统的开环传递函数为：，试求系统动态性能指标。Gk=tf(10,2 1 0);G0=feedback(Gk,1)Transfer function:10-2 s2+s+10 step(G0)%得到系统单位阶跃响应曲线 title(系统10/(2s2+s+10)的单位阶跃响应,Fontsize,16)%设置属性,图5 例3运行结果,得到系统的单位阶跃响应曲线后，在图形窗口上点击右键，在 Characeristics下的子菜单中可以选择Peak Response（峰值）、Settling Time(调整时间)、Rise Time(上升时间)和Steady State（稳态值）等参数

5、进行显示，操作如图5，其显示参数的系统响应曲线如图6。其它属性如title、x-label、y-label，调整时间的值等也可进入Properties子菜单设置。用户还可以在曲线上任选一点并用鼠标拖动之，系统将同时显示这点的时间及幅值。,图6 例3显示参数的系统响应,例4：单位负反馈系统的开环传递函数为：，编写程序求系统动态性能指标。程序如下：s=tf(s);Gk=7/s/(s+1);G0=feedback(Gk,1,-1)y,t=step(G0);%返回系统时域响应曲线值C=dcgain(G0);%得到系统终值%峰值时间计算max_y,k=max(y);peak_time=t(k)%超调量计

6、算max_overshoot=100*(max_y-C)/C%上升时间计算，以从稳态值的10%上升到90%定义r1=1;while(y(r1)0.1*C)r1=r1+1;endr2=r1;,while(y(r2)0.98*Cendsettling_time=t(s)step(G0)%求取系统响应曲线,图7 例4运行结果，思考为什么Characeristics产生的结果和编程获得结果不同？如何解决此问题？写到实验报告中！,12.2控制系统的稳态性能指标分析,系统的稳态性能指标 稳态误差：系统误差为，而稳态误差即当时间t趋于无穷时，系统输出响应的期望值与实际值之差。这种定义被称为是在输出端定义的稳

7、态误差。表12.4 给出不同输入信号下系统的稳态误差计算方式。,表12.4 给出不同输入信号下系统的稳态误差计算方式。,表12.4 不同输入信号下系统的稳态误差计算,在MATLAB中，各稳态误差系数可由以下命令求取：Kp=dcgain(numk,denk)Kv=dcgain(numk 0,denk)Ka=dcgain(numk 0 0,denk),12.2.2 控制系统稳态性能指标MATLAB求取示例,例6：单位负反馈系统的开环传递函数为：，试求单位阶跃输入下 的稳态误差。手工计算：由题知，系统为0型系统。查表可知，系统在单位阶跃输入下存在稳态误差，且稳态误差为1/11。MATLAB程序：,s

8、=tf(s);Gk=10/(0.1*s+1)/(0.5*s+1);G=feedback(Gk,1)%得到闭环系统 Transfer function:10-0.05 s2+0.6 s+11 step(G)%得到系统阶跃响应曲线 ess=1-dcgain(G)%得到稳态误差 ess=0.0909,图8 直接绘图求稳态值请用误差系数法求解，并把程序和结论写在实验报告中！,分析：手工计算和MATLAB程序得出的结果比较是一致的。可见由MATLAB程序很容易得到稳态误差。使用 Simulink求取稳态误差更方便，因为在 Simulink下可以直接将误差信号引出到示波器观察。,12.3MATLAB时域响

9、应仿真的典型函数应用,MATLAB时域响应仿真的典型函数MATLAB时域响应仿真的典型输入函数除step（单位阶跃函数）外，还有impulse（单位脉冲函数），lsim（求任意函数作用下系统响应的函数）等。各函数的用法如表12.5。,表12.5 求取时域响应函数及用法,MATLAB时域响应仿真的典型函数应用实例例7：求一阶惯性环节的脉冲响应曲线，观察T变化对系统性能的影响。程序如下：t=0:0.1:50;for T=1 5 10 G=tf(1,T 1);impulse(G,t);hold on end title(系统1/(Ts+1)脉冲响应曲线.T取1,5,10,Fontsize,16);,

10、图9例7一阶惯性环节脉冲响应曲线（为便于观察，图形进行了放大）,对于如上曲线，也可以采用多图绘制的方法完成，程序如下：t=0:0.1:100;T=1 5 10for n=1:3 G=tf(1,T(n)1);y(:,n)=impulse(G,t);%得到系统响应返回参数endplot(t,y)title(系统1/(Ts+1)脉冲响应曲线.T取1,5,10,Fontsize,16);figure(2);%另外生成一个图subplot(2,2,1)plot(t,y(:,1);title(T=1);subplot(2,2,2)plot(t,y(:,2);title(T=5);subplot(2,2,3

11、)plot(t,y(:,3);title(T=10);,图10 例7运行结果,例8：已知某控制系统的闭环传递函数为 120/(s2+12s+120)（1）求在单位斜坡输入作用下，系统的输出响应曲线。（2）求在带偏置正弦信号(2+sint)作用下，系统的输出响应曲线。,系统的单位斜坡输入响应曲线求取方式1：t=0:0.1:10;num=120;den=1 12 120;G=tf(num,den);u=t;y=lsim(G,u,t);plot(t,y,g,t,u,b-);axis(0 2.5 0 2.5);title(Unit-Ramp Response,-by use of lsim funct

12、ion);,图12 方式1运行结果,系统的单位斜坡输入响应曲线求取方式2：t=0:0.1:10;num=120;den=1 12 120 0;y=step(num,den,t);plot(t,y,g,t,t,b-);axis(0 2.5 0 2.5);title(Unit-Ramp Response-by use of step function);,图11 方式2运行结果，为什么可以用阶跃函数来产生斜坡信号的响应？请回答到实验报告中！,（2）系统在输入信号作用下的响应程序 t=0:0.01:10;num=120;den=1 12 120;G=tf(num,den);u=2+sin(t);y=

13、lsim(G,u,t);plot(t,y,t,u,b-);title(Response to 2+sin(t)-by use of lsim function);,图13,输入信号,作用下的响应曲线,12.4MATLAB图形化时域分析,除应用函数直接进行时域分析之外，也可以利用MATLAB的图形工具，得到系统的响应曲线及性能指标，以供进一步分析。,MATLAB LTI Viewer应用实例,有关MATLAB LTI Viewer的用法，在系统稳定性分析一章中已有初步介绍。,例9：当取0.2,0.4,0.6时，通过LTI Viewer工具观察二阶系统 的阶跃响应曲线和脉冲响应曲线。,1.编写MA

14、TLAB程序，求取不同值时各系统传递函数。for i=1:3 zeta(i)=0.2*i;ss(i)=tf(1,1 2*zeta(i)1);end,2.打开MATLAB LTI Viewer（用命令ltiview）。3.导入已经建立的系统 ss（Import菜单），如图14。,图14 导入已建立系统窗口,值时的阶跃响应曲线,图15 不同,4.用户可以使用快捷菜单在默认窗口上观察系统的阶跃响应性能指标。也可以选择Plot TypesImpulse项以显示单位脉冲响应曲线。如图16。,图16 系统单位脉冲响应曲线,5.用户还可改变显示方式。如通过EditPlot Configurations选项打开如图17窗口：,图17 显示曲线配置窗口 图18 上下分区方式显示响应曲线,在如图17的选择中，可以上下分区的方式同时显示系统阶跃、脉冲响应曲线。如图18所示。通过右键菜单IO Selector选项打开如图19窗口，在图19的选择下，只显示第一种情况下的响应曲线，如图20。,图19 IO Selector窗口 图20 只显示输入1的响应曲线,