1、1第 3 章 辅导控制系统典型的输入信号1. 阶跃函数阶跃函数的定义是 0 , ,)(ttArtx式中 A 为常数。A 等于 1 的阶跃函数称为单位阶跃函数，如图所示。它表示为xr(t)l(t) ，或 xr(t)=u(t)单位阶跃函数的拉氏变换为Xr(s)=L1(t)=1/s在 t0 处的阶跃信号，相当于一个不变的信号突然加到系统上；对于恒值系统，相当于给定值突然变化或者突然变化的扰动量；对于随动系统，相当于加一突变的给定位置信号。2. 斜坡函数这种函数的定义是 0 ,)(tAtxr式中 A 为常数。该函数的拉氏变换是Xr(s)=LAt=A/s2这种函数相当于随动系统中加入一按恒速变化的位置信
2、号，该恒速度为 A。当 Al 时，称为单位斜坡函数，如图所示。23. 抛物线函数如图 所示，这种函数的定义是 0 ,t)(2Atxr式中 A 为常数。这种函数相当于随动系统中加入一按照恒加速变化的位置信号，该恒加速度为 A。抛物线函数的拉氏变换是Xr(s)=LAt2=2A/s3当 A1/2 时，称为单位抛物线函数，即 Xr(s)=1/s3。4. 脉冲函数这种函数的定义是 0)( 0 , , ,)(tAtxr式中 A 为常数， 为趋于零的正数。脉冲函数的拉氏变换是 ALsXrlim0)(当 A1，0 时，称为单位脉冲函数 (t)，如图 所示。单位脉冲函数的面积等于l，即31)(dt在 tt 0
3、处的单位脉冲函数用 (t-t 0)来表示，它满足如下条件幅值为无穷大、持续时间为零的脉冲纯属数学上的假设，但在系统分析中却很有用处。单位脉冲函数 (t)可认为是在间断点上单位阶跃函数对时间的导数，即反之，单位脉冲函数 (t)的积分就是单位阶跃函数。控制系统的时域性能指标对控制系统的一般要求归纳为稳、准、快。工程上为了定量评价系统性能好坏，必须给出控制系统的性能指标的准确定义和定量计算方法。1 动态性能指标动态性能指标通常有如下几项：延迟时间 阶跃响应第一次达到终值 的 50所需的时间。dt )(h上升时间 阶跃响应从终值的 10上升到终值的 90所需的时间；对有振荡的系统，r也可定义为从 0
4、到第一次达到终值所需的时间。峰值时间 阶跃响应越过稳态值 达到第一个峰值所需的时间。pt )(调节时间 阶跃响到达并保持在终值 误差带内所需的最短时间；有时也s h5用终值的 误差带来定义调节时间。2超调量 峰值 超出终值 的百分比，即)(pth)( 10htp在上述动态性能指标中，工程上最常用的是调节时间 （描述“快” ） ，超调量 st （描述“匀” ）以及峰值时间 。 pt2 稳态性能指标稳态误差是时间趋于无穷时系统实际输出与理想输出之间的误差，是系统控制精度或抗干扰能力的一种度量。稳态误差有不同定义，通常在典型输入下进行测定或计算。4一阶系统的阶跃响应一. 一阶系统的数学模型由一阶微分
5、方程描述的系统，称为一阶系统。一些控制元部件及简单系统如 RC 网络、发电机、空气加热器、液面控制系统等都是一阶系统。因为单位阶跃函数的拉氏变换为 R(s)=1/s，故输出的拉氏变换式为 11)()( TsTsRsC取 C(s)的拉氏反变换得tec(t)或写成 tst式中，c ss=1，代表稳态分量； 代表暂态分量。当时间 t 趋于无穷，暂态分Ttec1量衰减为零。显然，一阶系统的单位阶跃响应曲线是一条由零开始，按指数规律上升并最终趋于 1 的曲线，如图所示。响应曲线具有非振荡特征，故又称为非周期响应。一阶系统的单位阶跃响应二阶系统的阶跃响应典型二阶系统方框图，其闭环传递函数为： vmvmKs
6、TsKsRC2)1(/1 22ns式中Kv-开环增益； n-无阻尼自然频率或固有频率， ；mvn-阻尼比， 。 mnT21二阶系统的闭环特征方程为s2+2 ns+ 2n=0 其特征根为 12,151. 临界阻尼(=1)其时域响应为)1(tetcntn上式包含一个衰减指数项。c(t)为一无超调的单调上升曲线，如图 3-8b 所示。(a) (b) (c)1 时二阶系统的特征根的分布与单位阶跃响应2. 过阻尼(1)具有两个不同负实根 的惯性环节单位阶跃响应拉氏变换)1(,221ns式。其时域响应必然包含二个衰减的指数项，其动态过程呈现非周期性，没有超调和振荡。图为其特征根分布图。3. 欠阻尼（01）
7、图 3-9 01 时二阶系统特征根的分布 图 3-10 欠阻尼时二阶系统的单位阶跃响应4. 无阻尼(0) )(2nsC其时域响应为 ttco1在这种情况下，系统的响应为等幅(不衰减) 振荡， 6图 0 时特征根的分布 图 0 时二阶系统的阶跃响应5. 负阻尼（0）当 0 时，特征根将位于复平面的虚轴之右，其时域响应中的 e 的指数将是正的时间函数，因而 为发散的，系统是不稳定的。tne显然，0 时的二阶系统都是不稳定的，而在 1 时，系统动态响应的速度又太慢，所以对二阶系统而言，欠阻尼情况是最有实际意义的。下面讨论这种情况下的二阶系统的动态性能指标。欠阻尼二阶系统的动态性能指标1. 上升时间
8、tr上升时间 tr 是指瞬态响应第一次到达稳态值所需的时间。21ndrt由此式可见，阻尼比 越小，上升时间 tr 则越小； 越大则 tr 越大。固有频率 n 越大，t r 越小，反之则 tr 越大。2. 峰值时间 tp 及最大超调量 Mp21ndpt最大超调量 )/(max2)ecp最大超调百分数 %10.(%)/(2c3. 调整时间 ts 70. 4)1ln(24)2( . 3352，nnst图 3-13 二阶系统单位阶跃响应的一对包络线 图 3-14 调节时间和阻尼比的近似关系7根据以上分析，二阶振荡系统特征参数 和 n 与瞬态性能指标(4. 振荡次数 在调整时问 ts 之内，输出 c(t
9、)波动的次数称为振荡次数 ，显然fst式中 ，称为阻尼振荡的周期时间。212ndft 12TSs这一系统的单位阶跃响应瞬态特性指标为：最大超调百分数 %3.40)1/(%2e上升时间 Ttnr 7.2调整时间(用近似式求得为 8T)ts43.8(用近似式求得为 6T)1%5有一位置随动系统其中 Kk4。求该系统的( )固有频率；(2)阻尼比；(3) 超调量和调整时间；(4) 如果要求实现工程最佳参数 l ，开环放大系数 值应是多少？2k【解】系统的闭环传递函数为kKs24k与二阶系统标准形式的传递函数 22ns对比得：(1) 固有频率 4kn(2) 阻尼比 由 得 125.02n(3) 超调
10、%7%)/(2ne(4) 调整时间 stns635当要求 时，由 得 2115.0,212nkK可见该系统要满足工程最佳参数的要求，须降低开环放大系数 的值。但是，降低k8值将增大系统的误差。kK劳斯稳定判据将系统的特征方程式写成如下标准式01210 nnn assasa将各系数组成如下排列的劳斯表 122432132756420gsfeccsbbaon 表中的有关系数为 13021ab54217063ab系数 的计算，一直进行到其余的 b 值全部等于零为止。ib1231ac52b1473ac这一计算过程，一直进行到 n 行为止。为了简化数值运算，可以用一个正整数去除或乘某一行的各项，这时并不
11、改变稳定性的结论。(l) 第一列所有系数均不为零的情况 第一列所有系数均不为零时，劳斯判据指出，特征方程式的实部为正实数根的数目等于劳斯表中第一列的系数符号改变的次数。方程式9的根全部在复平面的左半平面的充分必要条件是，方程式的各项系数全部为正值，并且劳斯表的第一列都具有正号。例如， 三阶系统的特征方程式为 032130assa列出劳斯表为 3010213as则系统稳定的充分必要条件是， ， ， ，0a1200)(321a系统的特征方程为54345 sss试用劳斯判据判断系统的稳定性。解 计算劳斯表中各元素的数值，并排列成下表 53209150345ss由上表可以看出，第一列各数值的符号改变了
12、两次，由+2 变成-1，又由-1 改变成+9 。因此该系统有两个正实部的根，系统是不稳定的。(2) 某行第一列的系数等于零而其余项中某些项不等于零的情况 在计算劳斯表中的各元素的数值时，如果某行的第一列的数值等于零，而其余的项中某些项不等于零，那么可以用一有限小的数值 来代替为零的那一项，然后按照通常方法计算阵列中其余各项。如果零( ) 上面的系数符号与零( ) 下面的系数符号相反，则表明这里有一个符号变化。例如，对于下列特征方程式 01234ss劳斯表为 12)0(03ss现在观察第一列中的各项数值。当 趋近于零时， 的值是一很大的负值，因此可以认为第一列中的各项数值的符号改变了两次。由此得
13、出结论，该系统特征方程式有两个根具有正实部，系统是不稳定的。10如果零() 上面的系数符号与零() 下面的系数符号不变，则表示系统有纯虚根。例如，对下列特征方程式 023ss劳斯表为 2102s可以看出，第一列各项中 的上面和下面的系数符号不变，故有一对虚根。将特征方程式分解，有)(12s解得根为， 2,1jp23p(3) 某行所有各项系数均为零的情况 如果劳斯表中某一行的各项均为零，或只有等于零的一项， 这表示在 s 平面内存在一些大小相等但符号相反的特征根。 在这种情况下，可利用全零行的上一行各系数构造一个辅助方程，式中 s 均为偶次。将辅助方程对 s 求导，用所得的导数方程系数代替全零行
14、，然后继续计算下去。至于这些大小相等，符号相反的根，可以通过解辅助方程得到。系统特征方程式为 016201823456 sss试用劳斯判据判断系统的稳定性。解 劳斯表中的 各项为6s36s0812345s由上表可以看出， 行的各项全部为零。为了求出 - 各项，将 行的各项组成辅3s 3s04s助方程为 86)(24sA将辅助方程 对 s 求导数得)(Asds13用上式中的各项系数作为 行的各项系数，并计算以下各行的各项系数，得劳斯表为3s12486036ss11834012s从上表的第一列可以看出，各项符号没有改变，因此可以确定在右半平面没有特征方程式的根。另外，由于 行的各项皆为零，这表示有
15、共轭虚根。这些根可由辅助方程求出。3s本例中的辅助方程式是 08624s由之求得特征方程式的大小相等符号相反的虚根为， ， 2,1jp4,3jp216,5jp稳态误差及其计算误差本身是时间 t 的函数，在时域中以 表示。稳定系统误差的终值称为稳态误差te，即为误差信号的稳态分量，则稳态误差为se)()(limli0sEtsts系统的误差传递函数)(1HGsRE故 s将系统误差的拉氏变换 E(s)代入(3-38)，得稳态误差的计算公式为)(1li0ses控制系统的型别控制系统的一般开环传递函数可以写成NnjjmiiksTKsHG1)()(式中 为开环放大系数或称为开环传递系数； 、 为时间常数；
16、N 表示开环传递kKij函数中串联的积分环节个数。这是一个很重要的结构参数。根据 N 的数值，可将系统分为几种不同类型。N0 的系统称为 0 型系统；N 1 的系统称为 I 型系统；N=2 的系统称为II 型系统。当 N2 时，要使系统稳定是很困难的。因此，一般采用的是 0 型、I 型和 II 型系统。典型输入下系统的稳态误差对于不同输入函数，下面分析系统的稳态误差。121. 单位阶跃输入下的稳态误差单位阶跃输入( )下的系统稳态误差，由式(3-40) 得sR1)(lim0sHGes)(1s定义0ksp称为位置误差系数，则pkpsKe10 型系统的稳态误差为ksI 型或高于 I 型的系统的位置
17、稳态误差为0se2. 单位斜坡输入下的稳态误差单位斜坡输入( )的系统稳态误差21sR )(1)(limli020 sHGsHGess 定义)(0Ks称为速度误差系数。K则 es1对于 0 型系统 0)1(0limnjjiikssTKk所以 se对于 I 型系统 knjjmiiks KsTk10)(l13所以 ksKe1对于 II 型或更高型系统njjmiikssTk120)(l所以 se，0 型系统不能跟踪斜坡输入；单位反馈的型系统能跟踪斜坡输入，但总有一定误差， 3. 单位抛物线输入下的稳态误差定义)(20limsHGKsa称为加速度误差系数aK则 ase1对于 0 型系统 0)1(02linjjmiiksasTKk所以 se对于型系统 0)1(120linjjmiiksasTKk所以 se对于型系统 knjjmiiksa KsTk120)(l所以 kse由此可知，0 型和型系统都不能跟踪抛物线输入，型系统能跟踪抛物线输入，但存在稳态误差。典型输入信号作用下的稳态误差和误差系数14静态误差系数 阶跃输入 )(1)(tAtr斜坡输入 tAr)(加速度输入 2)(tAr系统型别 pKva位置误差 psKe速度误差 vsKe加速度误差 asKe0 0 0 1I 0 0AII K0 0 KA。减小稳态误差的方法1. 引入给定量顺馈2. 引入扰动量顺馈
