二阶系统的阶跃响应.pdf

《二阶系统的阶跃响应.pdf》由会员分享，可在线阅读，更多相关《二阶系统的阶跃响应.pdf（14页珍藏版）》请在三一文库上搜索。

1、1 第 3 章 辅导 控制系统典型的输入信号 1. 阶跃函数 阶跃函数的定义是 0,0 0, )( t tA r tx 式中 A 为常数。 A 等于 1的阶跃函数称为单位阶跃函数，如图所示。它表示为 xr(t) l(t)，或 xr(t)=u(t) 单位阶跃函数的拉氏变换为 Xr(s)=L1(t)=1/s 在 t 0 处的阶跃信号，相当于一个不变的信号突然加到系统上；对于恒值系统，相当 于给定值突然变化或者突然变化的扰动量；对于随动系统， 相当于加一突变的给定位置信号。 2. 斜坡函数 这种函数的定义是 0,0 0, )( t ttA txr 式中 A 为常数。该函数的拉氏变换是 Xr(s)=L

2、At=A/s 2 这种函数相当于随动系统中加入一按恒速变化的位置信号，该恒速度为A。当 Al 时， 称为单位斜坡函数，如图所示。 2 3. 抛物线函数 如图 所示，这种函数的定义是 0,0 0, t )( 2 t tA txr 式中 A 为常数。这种函数相当于随动系统中加入一按照恒加速变化的位置信号，该恒 加速度为A。抛物线函数的拉氏变换是 Xr(s)=LAt 2=2A/s3 当 A1/2 时，称为单位抛物线函数，即Xr(s)=1/s3。 4. 脉冲函数 这种函数的定义是 0)(0, )0(,0,0 )( t A tt txr 式中 A 为常数， 为趋于零的正数。脉冲函数的拉氏变换是 A A

3、LsXr lim 0 )( 当 A1，0 时，称为单位脉冲函数(t)，如图所示。单位脉冲函数的面积等于l， 即 3 1)( dtt 在 tt0处的单位脉冲函数用(t-t0)来表示，它满足如下条件 幅值为无穷大、持续时间为零的脉冲纯属数学上的假设，但在系统分析中却很有用处。 单位脉冲函数(t)可认为是在间断点上单位阶跃函数对时间的导数，即 反之，单位脉冲函数(t)的积分就是单位阶跃函数。 控制系统的时域性能指标 对控制系统的一般要求归纳为稳、准、快。工程上为了定量评价系统性能好坏，必须给 出控制系统的性能指标的准确定义和定量计算方法。 1 动态性能指标 动态性能指标通常有如下几项： 延迟时间 d

4、 t阶跃响应第一次达到终值)(h的 50所需的时间。 上升时间 r t阶跃响应从终值的10上升到终值的90所需的时间； 对有振荡的系统， 也可定义为从0 到第一次达到终值所需的时间。 峰值时间 p t阶跃响应越过稳态值)(h达到第一个峰值所需的时间。 调节时间 s t阶跃响到达并保持在终值)(h5误差带内所需的最短时间；有时也用 终值的2误差带来定义调节时间。 超调量峰值)( p th超出终值)(h的百分比，即 100 )( )()( h hthp 在上述动态性能指标中，工程上最常用的是调节时间 s t（描述 “快”） ，超调量（描 述“匀”）以及峰值时间 p t。 2 稳态性能指标 稳态误差

5、是时间趋于无穷时系统实际输出与理想输出之间的误差，是系统控制精度或抗 干扰能力的一种度量。稳态误差有不同定义，通常在典型输入下进行测定或计算。 4 一阶系统的阶跃响应 一. 一阶系统的数学模型 由一阶微分方程描述的系统，称为一阶系统。一些控制元部件及简单系统如RC 网络、 发电机、空气加热器、液面控制系统等都是一阶系统。 因为单位阶跃函数的拉氏变换为R(s)=1/s，故输出的拉氏变换式为 1 11 1 1 )()()( Ts T ssTs sRssC 取 C(s)的拉氏反变换得 t T ec(t) 1 1 或写成 ttss ccc(t) 式中， css=1，代表稳态分量； t T tt ec

6、1 代表暂态分量。当时间t 趋于无穷，暂态分 量衰减为零。 显然， 一阶系统的单位阶跃响应曲线是一条由零开始，按指数规律上升并最终 趋于 1 的曲线，如图所示。响应曲线具有非振荡特征，故又称为非周期响应。 一阶系统的单位阶跃响应 二阶系统的阶跃响应 典型二阶系统方框图，其闭环传递函数为： vm v mv mv KssT K sTsK sTsK sR sC s 2 ) 1(/1 )1(/ 22 2 2 nn n ss 式中 Kv-开环增益； n-无阻尼自然频率或固有频率， m v n T K ； -阻尼比， mnT 2 1 。 二阶系统的闭环特征方程为 s 2+2 ns+ 2 n=0 其特征根为

7、 5 n s1 2 2,1 1. 临界阻尼 (=1) 其时域响应为 )1(1tetc n t n 上式包含一个衰减指数项。c(t)为一无超调的单调上升曲线，如图3-8b 所示。 (a) (b) (c) 1时二阶系统的特征根的分布与单位阶跃响应 2. 过阻尼 (1) 具有两个不同负实根)1(, 2 21n ss的惯性环节单位阶跃响应拉氏变换 式。其时域响应必然包含二个衰减的指数项，其动态过程呈现非周期性，没有超调和振荡。 图为其特征根分布图。 3. 欠阻尼（ 01） 图 3-9 01时二阶系统特征根的分布图 3-10 欠阻尼时二阶系统的单位阶跃响应 4. 无阻尼 (0) )( 22 2 n n

8、ss sC 其时域响应为 ttc n cos1 在这种情况下，系统的响应为等幅(不衰减 )振荡， 6 图0 时特征根的分布图0 时二阶系统的阶跃响应 5. 负阻尼（ 0） 当0 时，特征根将位于复平面的虚轴之右，其时域响应中的e的指数将是正的时间 函数，因而 t n e 为发散的，系统是不稳定的。 显然， 0 时的二阶系统都是不稳定的，而在 1 时，系统动态响应的速度又太慢， 所以对二阶系统而言，欠阻尼情况是最有实际意义的。下面讨论这种情况下的二阶系统的动 态性能指标。 欠阻尼二阶系统的动态性能指标 1. 上升时间 tr 上升时间 tr是指瞬态响应第一次到达稳态值所需的时间。 2 1 n d

9、r t 由此式可见，阻尼比越小，上升时间tr则越小； 越大则 tr越大。固有频率n越大， tr越小，反之则 tr越大。 2. 峰值时间 tp及最大超调量 Mp 2 1 n d p t 最大超调量 )1/( max 2 )(eccM p 最大超调百分数%100. )( )( % )1/(max 2 e c cc c 3. 调整时间 ts 707.00 4 )1ln( 2 1 4 1 %)2( 707. 00 3 )1ln( 2 1 3 1 %)5( 2 2 ， ， nn s nn s t t 图 3-13 二阶系统单位阶跃响应的一对包络线图 3-14 调节时间和阻尼比的近似关系 7 根据以上分析

10、， 二阶振荡系统特征参数和n与瞬态性能指标( 4. 振荡次数 在调整时问ts之内，输出c(t)波动的次数称为振荡次数 ，显然 f s t t 式中 2 1 22 n d f t，称为阻尼振荡的周期时间。 122 1 22 TSST s 这一系统的单位阶跃响应瞬态特性指标为： 最大超调百分数 %3.4%100 )1/( % 2 e 上升时间 Tt n r 7.4 1 2 调整时间 Tts43.8%2(用近似式求得为8T) Tts14. 4%5(用近似式求得为6T) 有一位置随动系统其中Kk4。求该系统的 ()固有频率； (2)阻尼比； (3)超调量和调 整时间； (4)如果要求实现工程最佳参数l

11、2，开环放大系数 k k值应是多少？ 【解】系统的闭环传递函数为 k k Kss K s 2 4 k K 与二阶系统标准形式的传递函数 22 2 2 nn n ss s 对比得： (1) 固有频率24 kn K (2) 阻尼比由12 n 得25.0 2 1 n (3) 超调%47%100% )1/( 2n e (4) 调整时间st n s 6 3 %5 当要求 2 1 时，由12 n 得5.0, 2 12 nkn K 可见该系统要满足工程最佳参数的要求，须降低开环放大系数 k K的值。但是，降低 k K 值将增大系统的误差。 8 劳斯稳定判据 将系统的特征方程式写成如下标准式 0 1 2 2

12、1 10nn nnn asasasasa 将各系数组成如下排列的劳斯表 1 1 1 21 2 4321 3 4321 2 7531 1 6420 gs fs ees ccccs bbbbs aaaas aaaas o n n n n 表中的有关系数为 1 3021 1 a aaaa b 1 5041 2 a aaaa b 1 7061 3 a aaaa b 系数 i b的计算，一直进行到其余的b 值全部等于零为止。 1 2131 1 b baab c 1 3151 2 b baab c 1 4171 3 b baab c 这一计算过程，一直进行到n 行为止。为了简化数值运算，可以用一个正整数去

13、除或 乘某一行的各项，这时并不改变稳定性的结论。 (l) 第一列所有系数均不为零的情况第一列所有系数均不为零时，劳斯判据指出， 特征方程式的实部为正实数根的数目等于劳斯表中第一列的系数符号改变的次数。方程式的 根全部在复平面的左半平面的充分必要条件是，方程式的各项系数全部为正值，并且劳斯表 9 的第一列都具有正号。 例如，三阶系统的特征方程式为 0 32 2 1 3 0 asasasa 列出劳斯表为 3 0 1 30211 31 2 20 3 as a aaaa s aas aas 则系统稳定的充分必要条件是 0 0 a，0 1 a，0 2 a，0 3 a，0)( 3021 aaaa 系统的特

14、征方程为 05432 2345 sssss 试用劳斯判据判断系统的稳定性。 解 计算劳斯表中各元素的数值，并排列成下表 5 32 059 031 532 411 0 1 2 3 4 5 s s s s s s 由上表可以看出，第一列各数值的符号改变了两次，由+2 变成 -1，又由 -1 改变成 +9。 因此该系统有两个正实部的根，系统是不稳定的。 (2) 某行第一列的系数等于零而其余项中某些项不等于零的情况在计算劳斯表中的 各元素的数值时，如果某行的第一列的数值等于零，而其余的项中某些项不等于零，那么可 以用一有限小的数值来代替为零的那一项，然后按照通常方法计算阵列中其余各项。如果 零()上面

15、的系数符号与零()下面的系数符号相反，则表明这里有一个符号变化。 例如，对于下列特征方程式 0122 234 ssss 劳斯表为 1 2 2 1)0( 022 111 0 1 2 3 4 s s s s s 现在观察第一列中的各项数值。当趋近于零时， 2 2的值是一很大的负值，因此可 以认为第一列中的各项数值的符号改变了两次。由此得出结论， 该系统特征方程式有两个根 具有正实部，系统是不稳定的。 如果零 ( )上面的系数符号与零()下面的系数符号不变，则表示系统有纯虚根。例如， 10 对下列特征方程式 022 23 sss 劳斯表为 2 22 11 0 1 2 3 s s s s 可以看出，

16、第一列各项中 的上面和下面的系数符号不变，故有一对虚根。 将特征方程 式分解，有 0)2)(1( 2 ss 解得根为 1 2, 1 jp，2 3 p (3) 某行所有各项系数均为零的情况如果劳斯表中某一行的各项均为零，或只有等于 零的一项，这表示在s 平面内存在一些大小相等但符号相反的特征根。在这种情况下， 可利用全零行的上一行各系数构造一个辅助方程，式中s 均为偶次。将辅助方程对s求导， 用所得的导数方程系数代替全零行，然后继续计算下去。至于这些大小相等，符号相反的根， 可以通过解辅助方程得到。 系统特征方程式为 01616201282 23456 ssssss 试用劳斯判据判断系统的稳定性

17、。 解劳斯表中的 6 s 3 s各项为 162081 6 s 000 861 016122 3 4 5 s s s 由上表可以看出， 3 s行的各项全部为零。为了求出 3 s- 0 s各项，将 4 s行的各项组成辅 助方程为 86)( 24 sssA 将辅助方程)(sA对 s 求导数得 ss ds sdA 124 )(3 用上式中的各项系数作为 3 s行的各项系数，并计算以下各行的各项系数，得劳斯表为 124 861 016122 162081 3 4 6 s s s s s 8 3 4 83 0 1 2 s s s 11 从上表的第一列可以看出，各项符号没有改变，因此可以确定在右半平面没有特

18、征方程 式的根。另外，由于 3 s行的各项皆为零，这表示有共轭虚根。这些根可由辅助方程求出。 本例中的辅助方程式是 086 24 ss 由之求得特征方程式的大小相等符号相反的虚根为 2 2, 1 jp，2 4,3 jp，21 6, 5 jp 稳态误差及其计算 误差本身是时间t 的函数，在时域中以te表示。稳定系统误差的终值称为稳态误差 ss e， 即为误差信号的稳态分量，则稳态误差为 )()( limlim 0 ssEtee st ss 系统的误差传递函数 )()(1 1 sHsGsR sE 故 )()(1sHsG sR sE 将系统误差的拉氏变换E(s)代入 (3-38)，得稳态误差的计算公

19、式为 )()(1 )( lim 0sHsG ssR e s ss 控制系统的型别 控制系统的一般开环传递函数可以写成 Nn j j N m i ik sTs sTK sHsG 1 1 )1( ) 1( )()( 式中 k K为开环放大系数或称为开环传递系数； i T、 j T为时间常数； N 表示开环传递函 数中串联的积分环节个数。这是一个很重要的结构参数。根据N 的数值，可将系统分为几 种不同类型。N0 的系统称为0 型系统； N1 的系统称为I 型系统； N=2 的系统称为II 型系统。当N2 时，要使系统稳定是很困难的。因此，一般采用的是0 型、 I 型和 II 型系 统。 典型输入下系

20、统的稳态误差 对于不同输入函数，下面分析系统的稳态误差。 1. 单位阶跃输入下的稳态误差 单位阶跃输入( s sR 1 )下的系统稳态误差，由式(3-40)得 )()(1 lim 0sHsG s e s ss )()(1 11 sHsGs 12 定义 )()( lim 0 sHsGk s p p k称为位置误差系数，则 p ss K e 1 1 0 型系统的稳态误差为 k ss K e 1 1 I 型或高于I 型的系统的位置稳态误差为 0 ss e 2. 单位斜坡输入下的稳态误差 单位斜坡输入( 2 1 s sR)的系统稳态误差 )()( 11 )()(1 lim lim 0 2 0sHssG

21、ssHsG s e s s ss 定义 )()( lim 0 sHssGK s K称为速度误差系数。 则 K ess 1 对于 0 型系统 0 )1( ) 1( 1 0 1 0 lim n j j m i ik s sTs sTK k 所以 ss e 对于 I 型系统 k n j j m i ik s K sTs sTK sk 1 1 0 )1( )1( lim 所以 k ss Kk e 11 对于 II 型或更高型系统 13 n j j m i ik s sTs sTK sk 1 2 1 0 ) 1( )1( lim 所以 0 ss e ，0 型系统不能跟踪斜坡输入；单位反馈的型系统能跟踪斜

22、坡输入，但总有一 定误差， 3. 单位抛物线输入下的稳态误差 定义 )()( 2 0 lim sHsGsK s a a K称为加速度误差系数 则 a ss K e 1 对于 0 型系统 0 ) 1( ) 1( 1 0 12 0 lim n j j m i ik s a sTs sTK sk 所以 ss e 对于型系统 0 )1( ) 1( 1 1 12 0 lim n j j m i ik s a sTs sTK sk 所以 ss e 对于型系统 k n j j m i ik s a K sTs sTK sk 1 1 12 0 )1( )1( lim 所以 k ss k e 1 由此可知， 0 型和型系统都不能跟踪抛物线输入，型系统能跟踪抛物线输入，但存 在稳态误差。 典型输入信号作用下的稳态误差和误差系数 系 统 型 别 静态误差系数 阶跃输入 )( 1)(tAtr 斜坡输入 tAtr)( 加速度输入 2 )( 2 tA tr p K v K a K 位置误差 p ss K A e 1 速度误差 v ss K A e 加速度误差 a ss K A e 14 0 K 0 0 K A 1 I K 0 0 K A II K 0 0 K A 。 减小稳态误差的方法 1. 引入给定量顺馈 2. 引入扰动量顺馈