机控-系统的时间响应

1、中北大学机械工程与自动化学院中北大学机械工程与自动化学院主讲：马维金主讲：马维金 教授教授2021-12-21 1机械控制工程基础机械控制工程基础杨叔子杨叔子 主编主编 华中科技大学出版社华中科技大学出版社 20112011年年 3.1 时间响应及其组成 3.2 典型输入信号 3.3 一阶系统 3.4 二阶系统 3.5 高阶系统 3.6 系统误差分析与计算 3.7 函数在时间响应中的应用 3.8利用MATLAB分析时间响应2第3章 系统的时间响应33.0 系统的时间响应分析q建立系统数学模型之后，可以用不同的方法，通过系统的数学模型来分析系统的特性，时间响应分析是重要的系统特性分析方法之一。q

2、时间响应：在输入作用下，系统的输出在时域的表现形式，在数学上就是系统的动力学方程在一定初始条件下的解。q对系统的传统分析方法主要三种：时间响应分析 (第三章)、 频域响应分析 (第四章)和稳定性分析 (第五章)43.1 时间响应及其组成 举例说明53.1 时间响应及其组成 一般情况63.1 时间响应及其组成 讨论73.1 时间响应及其组成83.1 系统的特征根影响系统自由响应的收敛性iiijs数的特征根系统的特征根即传递函93.1 系统的特征根影响系统自由响应的收敛性iiijs数的特征根系统的特征根即传递函103.1 系统的特征根影响系统自由响应的收敛性iiijs数的特征根系统的特征根即传递函

3、113.1 系统的特征根影响系统自由响应的收敛性iiijs数的特征根系统的特征根即传递函123.1 系统的特征根影响系统自由响应的收敛性iiijs数的特征根系统的特征根即传递函133.1 系统的特征根影响系统自由响应的收敛性iiijs数的特征根系统的特征根即传递函143.1 系统的特征根影响系统自由响应的收敛性iiijs数的特征根系统的特征根即传递函153.1 系统的特征根影响系统自由响应的收敛性163.2 控制系统的典型输入信号173.3 一阶系统q可用一阶微分方程描述的系统称为一阶系统，其微分方程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） qT称为一阶系统的时间常数，也是一

4、阶系统的特征参数，表征系统的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo183.3.1 一阶系统的单位脉冲响应193.3.2 一阶系统的单位阶跃响应203.3.3 一阶系统的单位斜坡响应213.3.4 一阶系统的性能指标223.3.5 一阶系统不同T的响应情况11)()()()()()(TssXsXsGtxtxdttdxTioioo233.4 二阶系统q一般控制系统均系高阶系统，但在一定条件下，可忽略某些次要因素近似的用一个二阶系统来表示。因此研究二阶系统有较大的实际意义。 二阶系统的动力学方程及传递函数分别为： 22222222)()()()()()(

5、2)(nnnioinononosssXsXsGtxtxdttdxdttxd243.4 二阶系统 式中： 称为无阻尼固有频率； 称为阻尼比。它们是二阶系统的特征参数。 由上式的分母可以得到二阶系统的特征方程： 此方程的两个特征根是 上式可知，随着阻尼比 取值不同，二阶系统的特征根也不同。 0222nnss122, 1nns253.4 二阶系统 （1）当0 1时，两特征根为共轭复数，即：此时，系统为欠阻尼系统，二阶系统的传递函数的极点是一对位于复数平面左半面内的共轭复数极点。 22, 11nnjsdjs2, 1263.4 二阶系统 （2）当 0时，两特征根为共轭纯虚根，即：此时系统为无阻尼系统。n

6、js2, 1273.4 二阶系统 （3）当 1时，特征方程有两个相等的实根，即： 此时系统称为临界阻尼系统。 ns2, 1283.4 二阶系统 （4）当 1时，特征方程有两个不等实根，即： 此时系统称为过阻尼系统，它可看作是两个一阶惯性环节的组合，又可以视为两个一阶环节的并联，还可以看成两个一阶环节的串联。 122, 1nns293.4.1 二阶系统的单位脉冲响应11)()()()()()(TssXsXsGtxtxdttdxTioioo303.4.1 二阶系统的单位脉冲响应11)()()()()()(TssXsXsGtxtxdttdxTioioo313.4.1 二阶系统的单位脉冲响应11)()

7、()()()()(TssXsXsGtxtxdttdxTioioo323.4.1 二阶系统的单位脉冲响应333.4.2 二阶系统的单位阶跃响应q可用一阶微分方程描述的系统称为一阶系统，其微分方程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） T称为一阶系统的时常数，也是一阶系统的特征参数，表征系统的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo343.4.2 二阶系统的单位阶跃响应q可用一阶微分方程描述的系统称为一阶系统，其微分方程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） T称为一阶系统的时常数，也是一阶系统的特征参数，表征系统

8、的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo353.4.2 二阶系统的单位阶跃响应11)()()()()()(TssXsXsGtxtxdttdxTioioo363.4.3 二阶系统的单位阶跃响应373.4.3 二阶系统响应的性能指标q可用一阶微分方程描述的系统称为一阶系统，其微分方程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） T称为一阶系统的时常数，也是一阶系统的特征参数，表征系统的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo383.4.3 二阶系统响应的性能指标q可用一阶微分方程描述的系统称为一

9、阶系统，其微分方程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） T称为一阶系统的时常数，也是一阶系统的特征参数，表征系统的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo393.4.3 二阶系统响应的性能指标q可用一阶微分方程描述的系统称为一阶系统，其微分方程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） T称为一阶系统的时常数，也是一阶系统的特征参数，表征系统的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo403.4.3 二阶系统响应的性能指标q可用一阶微分方程描述的系统称为一阶系统，其微分方

10、程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） T称为一阶系统的时常数，也是一阶系统的特征参数，表征系统的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo413.4.3 二阶系统响应的性能指标423.4.3 二阶系统响应的性能指标q可用一阶微分方程描述的系统称为一阶系统，其微分方程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） T称为一阶系统的时常数，也是一阶系统的特征参数，表征系统的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo433.4.3 二阶系统响应的性能指标443.4.3 二阶系统响应

11、的性能指标 举例q可用一阶微分方程描述的系统称为一阶系统，其微分方程和传递函数的一般形式分别为：（所对应的典型环节是惯性环节） T称为一阶系统的时常数，也是一阶系统的特征参数，表征系统的固有特性。11)()()()()()(TssXsXsGtxtxdttdxTioioo453.4.3 二阶系统响应的性能指标 举例463.4.3 二阶系统响应的性能指标 举例473.5 高阶系统483.5 高阶系统493.5 高阶系统503.5 高阶系统513.5 高阶系统523.6 系统误差分析与计算q系统误差：期望输出与实际输出的差。q误差产生的原因：存在随机干扰、元件的性能不完善等。引起误差的内因是系统本身

12、的结构，引起误差的外因是输入量及其导数的变化。q系统的输出由瞬态分量和稳态分量所组成,因而系统的误差也由瞬态误差和稳态误差所组成。在过渡过程开始时，瞬态误差是主要部分，但它随着时间而逐渐衰减，稳态误差将成为主要部分。533.6.1 系统误差e(t)与偏差(t)系统误差：理想输出与实际输出之差543.6.1 系统误差e(t)与偏差(t)系统误差：理想输出与实际输出之差553.6.2 系统误差e(t)的一般计算563.6.3 系统稳态误差ess(t)与稳态偏差ss (t)系统的稳态性能描述：稳态误差与稳态偏差573.6.4 与输入有关的稳态偏差ss (t)583.6.4 与输入有关的稳态偏差ss

13、(t)593.6.4 与输入有关的稳态偏差ss (t)603.6.5 与干扰有关的稳态偏差ss (t)系统的稳态偏差：与结构特征和输入信号有关613.6.5 与干扰有关的稳态偏差ss (t)系统的稳态偏差：与结构特征和输入信号有关623.7 函数在时间响应中的作用q 函数定义：-1d)(:0)(ttttt且633.7 函数在时间响应中的作用q单位脉冲函数与传递函数互为拉普拉斯变换对。 q如果对系统输入一单位脉冲函数(t)，则系统的单位脉冲响应函数为)()()()(),()(1)(, )()()()(tgttgtwsWtwLtLsGtLsGsW643.7 函数在时间响应中的作用653.7 函数在

14、时间响应中的作用663.7 函数在时间响应中的作用673.7 函数在时间响应中的作用683.7 函数在时间响应中的作用693.8 利用MATLAB分析时间响应3.8.1 用MATLAB求系统时间响应q impulse 函数：用于生成单位脉冲响应703.8 利用MATLAB分析时间响应3.8.1 用MATLAB求系统时间响应q step 函数：用于生成单位阶跃响应713.8 利用MATLAB分析时间响应3.8.1 用MATLAB求系统时间响应q lsim 函数：用于生成对任意输入的时间响应723.8 利用MATLAB分析时间响应3.8.1 用MATLAB求系统时间响应q 举例：图示系统的传递函数

15、如下，求系统在时间常数取不同值时的单位脉冲响应、单位阶跃响应和正弦函数响应。73q 解：令时间常数分别取不同值时，分析单位脉冲响应和单位阶跃响应的MATLAB文本如下。t=0:0.01:0.8; nG=50;tao=0;dG=0.05 1+50*tao 50;G1=tf(nG,dG);tao=0.0125;dG=0.05 1+50*tao 50;G2=tf(nG,dG);tao=0.025;dG=0.05 1+50*tao 50;G3=tf(nG,dG)y1,T=impulse(G1,t);y1a,T=step(G1,t);y2,T=impulse(G2,t);y2a,T=step(G2,t)

16、;y3,T=impulse(G3,t);y3a,T=step(G3,t);subplot(121),plot(T,y1,-,T,y2,-.,T,y3,-)legend(tao=0,tao=0.0125,tao=0.025)xlabel(t(sec),ylabel(x(t);grid on;subplot(122),plot(T,y1a,-,T,y2a,-.,T,y3a,-)legend(tao=0,tao=0.0125,tao=0.025)grid on;xlabel(t(sec),ylabel(x(t);74q 解：令时间常数分别取不同值时，用MATLAB分析单位脉冲响应和单位阶跃响应的结果

17、如图。3.8 利用MATLAB分析时间响应3.8.1 用MATLAB求系统时间响应00.10.20.30.40.50.60.70.8-10-50510152025t(sec)x(t) 00.10.20.30.40.50.60.70.800.20.40.60.811.21.4t(sec)x(t) tao=0tao=0.0125tao=0.025tao=0tao=0.0125tao=0.02575q 对于任意输入，例如正弦输入，应用lsim函数求时间常数=0.025时的时间相应及误差曲线，MATLAB文本如下。t=0:0.01:1;u=sin(2*pi.*t);%tao=0.025;nG=50;d

18、G=0.05 1+50*tao 50;G=tf(nG,dG);%y=lsim(G,u,t);%plot(t,u,-.,t,y,-,t,u-y,-.,linewidth,1)legend(u(t),xo(t),e(t)grid;xlabel(t(sec),ylabel(x(t);3.8 利用MATLAB分析时间响应3.8.1 用MATLAB求系统时间响应76q 解：令时间常数=0.025时，用MATLAB分析正弦响应的结果如图。3.8 利用MATLAB分析时间响应3.8.1 用MATLAB求系统时间响应00.10.20.30.40.50.60.70.80.91-1-0.8-0.6-0.4-0.2

19、00.20.40.60.81t(sec)x(t) u(t)xo(t)e(t)77q举例：对于上述例题， 取不同值时，求出单位阶跃响应以后，用MATLAB计算的瞬态性能指标如下表。q从上表可见，系统引入速度负反馈以后，系统的调整时间和最大超调量都得到减小。3.8 利用MATLAB分析时间响应3.8.2 用MATLAB求系统的瞬态性能指标trtPMPts00.0640 0.1050 0.3509 0.35300.0125 0.0780 0.1160 0.1523 0.25000.025 0.1070 0.1410 0.0415 0.188078q举例：用MATLAB计算瞬态性能指标的文本如下。t=

20、0:0.001:1; yss=1; dta=0.02; nG=50; tao=0;dG=0.05 1+50*tao 50;G1=tf(nG,dG); tao=0.0125;dG=0.05 1+50*tao 50;G2=tf(nG,dG); tao=0.025;dG=0.05 1+50*tao 50;G3=tf(nG,dG); y1=step(G1,t);y2=step(G2,t);y3=step(G3,t); r=1;while y1(r)1-dta & y1(s)1+dta;s=s-1;end ts1=(s-1)*0.001; r=1;while y2(r)1-dta &y3(s)1+dta;s=s-1;end ts2=(s-1)*0.001; r=1;while y3(r)1-dta & y3(s)1+dta;s=s-1;end ts3=(s-1)*0.001 tr1 tp1 mp1 t