劳斯霍尔维茨稳定性判据.ppt

34劳斯-霍尔维茨稳定性判据,33反馈控制系统的稳态误差,31二阶系统的瞬态响应及性能指标,32增加零极点对二阶系统响应的影响,35控制系统灵敏度分析,第三章系统分析,综述,系统分析的概念。针对系统的稳定性(Stability)、稳态误差(Steady-State Error)和瞬态响应(Transient Response) 在时间域内，通过求解系统的微分方程得到 间接方法：稳定性判据(Stability Criterion)、根轨迹法(Root Locus Method)和频率法(Frequency Method) 发展,3.1. 二阶系统的瞬态响应及性能指标,基本概念 瞬态响应：系统的输出从输入信号r(t)作用时刻起，到稳定状态为止，随时间变化的过程。 分析系统的瞬态响应的方法： 1、直接求解法 2、间接评价法 3、计算机仿真法,3.1.1. 典型输入信号,控制系统的瞬态响应与输入信 号的形式有关,一个控制系统的实际输入信号往往具有多重形式，并且也常常难于事先确定,通常考虑某些典型输入信号对系统的影响,（一）阶跃信号 (Step Signal),图(3-1)阶跃信号,A=1时，称为单位阶跃信号，用1(t)表示,（二）斜坡信号(Ramp Signal),图(3-2)斜坡信号,当A=1时，则称为单位斜坡信号,（三）抛物线信号 (Parabolic Signal),图(3-3)抛物线信号,当A=1时，则称为单位抛物线信号,（四）脉冲信号 (Impulse Signal),图(3-4)脉冲信号,当e 趋于零时就得理想的单位脉冲信号(亦称d(t) 函数)。,（五）正弦信号(Sin Signal),图(3-4)脉冲信号,A为幅值,T为周期，w =2p/T为角频率,3.1.2. 系统的性能指标,性能指标用来衡量系统性能,常由系统在一定的典型输入信号作用下的具体性能指标来表示,性能指标有许多形式,性能指标,5、调整时间ts(settling time),主要包括：,1、最大超调量sp (percent overshot)：,2、延滞时间td (time delay)：响应曲线到达稳态值50%所需的时间,3、上升时间tr (rise time)：三种定义,4、峰值时间tp (peak time)：响应曲线到达第一个峰值所需的时间,性能指标图解,调整时间ts,峰值时间tp,上升时间tr,最大超调量sp,延滞时间td,其它性能指标,振荡次数，衰减比等,对于恒值控制系统，常以系统对单位扰动输入信号时的响应特性来衡量瞬态性能,图(3-7) 单位扰动输入,3.1.3. 瞬态响应 (Transient Response ),（一）一阶系统的瞬态响应,图(3-8) 一阶控制系统,代表一个电机的速度控制系统，其中t 是电机的时间常数,闭环传递函数,输出响应,一阶系统响应曲线,一阶系统的单位阶跃响应,图(3-9) 一阶系统的单位阶跃响应,系统时间常数 定义为系统响应达到稳态值的63.2%所需要的时间,响应的稳态值,(一)一阶系统的稳态误差,系统的稳态误差,（二）二阶系统的阶跃响应,闭环传递函数,图(3-10) 典型二阶系统结构,为阻尼比， wn为无阻尼自然振荡频率,以图(2-2)中R-L-C电路为例,传递函数,无阻尼自然振荡频率,当R=0时的谐振频率,阻尼比,电枢控制的直流电动机,输出w 与电枢电压ua之间传递函数为,典型二阶系统的阶跃响应特性,特征方程,解方程,0z1，欠阻尼情况 z =，临界阻尼情况 z ，过阻尼情况, 0z1 ，欠阻尼情况,系统传递函数,有阻尼振荡频率,输入r(t)=1(t),0z1 ，欠阻尼情况（续）,(a)根分布 (b)单位阶跃响应 图(3-11) 欠阻尼情况(0z 1),系统的误差为,当t时，稳态误差e ()。,z=，无阻尼情况,系统的特征根为一对共轭虚根s1,2= ±jwn,单位阶跃响应, z =，临界阻尼情况,两个相等的实数特征根：s1= s 2= -wn,(a)根分布 (b)单位阶跃响应 图(3-12) 临界阻尼情况(z 1),无超调，无振荡，单调响应过程, z ，过阻尼情况,两个不相等的实数根：,过阻尼情况（续）,(a)根分布 (b)单位阶跃响应 图(3-13) 过阻尼情况(z 1),无超调，过程拖得比 z =时长,不同z值下的二阶系统单位阶跃响应曲线族,二阶系统的脉冲响应特性,输入单位脉冲信号d(t)，即R(s)=,二阶系统单位脉冲响应的拉氏变换为,系统的单位脉冲响应,欠阻尼情况(0z 1),临界阻尼情况(z ),过阻尼情况(z 1),不同z 值时的单位脉冲响应曲线,（四）二阶系统的瞬态响应性能指标,主要的性能指标：上升时间 、峰值时间 、最大超调量 、调节时间,上升时间tr (rise time),c(tr)为1,峰值时间tp (peak time),一阶求导,到达第一个峰值时wd tp = p,最大超调量sp (percent overshot),代入t= tp,图(3-16) sp与z 的关系,调节时间ts (settling time),T=1/zwn,图(3-17) ts与z 的关系 图(3-18) z 稍微突变引起的ts突变,在z =0.69(或0.77)，ts有最小值，以后ts随z 的增大而近乎线性地上升,曲线的不连续性是由于z 在虚线附近稍微变化会引起ts突变造成,调节时间ts（近似方法）,或,小结,当wn一定，要减小tr和tp，必须减少z 值，要减少ts则应增大zwn值，而且z 值有一定范围，不能过大 增大wn，能使tr，tp和ts都减少 最大超调量sp只由z 决定，z 越小，sp越大。,3.1.4. 线性定常系统的重要特性,初始条件为零的线性定常系统,线性定常系统的重要特性(2),初始条件为零的线性定常系统,结论,单位脉冲响应是单位阶跃响应对时间的一阶导数。单位阶跃响应可以由单位斜坡响应和单位抛物线响应对时间的一阶导数和二阶导数求得 单位斜坡响应和单位抛物线响应是单位阶跃响应对时间的一重和二重积分 只要知道系统对某一种典型信号的响应，对其它典型信号的响应也可推知,3.2. 增加零极点对二阶系统响应的影响,高阶系统传递函数的一般形式,零极点的形式,高阶系统单位阶跃响应,假设没有重极点,高阶系统小结,主导极点,定义,假如高阶系统中距离虚轴最近的极点，其实数部分为其它极点的110或更小，并且附近又没有零点，则可认为系统的响应主要由该极点(或共轭复数极点)决定，这一分量衰减最慢。这种对系统瞬态响应起主要作用的极点，称为系统的主导极点,主导极点example,三阶系统传递函数,图(3-19) 三阶系统的零极点分布图,系统的性能指标可用二阶系统的曲线来表示,主导极点,主导极点example(续）,将三阶系统看成是由主导极点决定的二阶系统与一个惯性环节（一阶滤波器）串联而成,惯性环节的时间常数较大时，惯性环节的作用较强。二阶系统的输出c1(t)经过惯性环节的滤波后，震荡现象减弱很多,仿真结果,当z =0.45时，通过计算机仿真能够得到系统在单位阶跃输入下的响应,当t =2.25时，实数极点为-1/t = -0.444，而复数极点的实部为 -0.45，二者相差不大，系统是过阻尼的，响应没有超调,t 调整为0.9，即实数极点为-1.11，则计算得到的超调量为12%，调节时间为8.8秒,有零点情况,如果二阶系统包含有零点，且该零点位于主导极点附近，则会对系统的瞬态响应产生影响,标准二阶系统一个零点,有零点时的仿真结果,图(3-21) 含有一个零点二阶系统的阶跃响应,例3.1,解,1、从闭环传递函数可以看出，系统的传递系数(或静态增益)为1，所以系统对阶跃输入的稳态误差为零,零极点分布图,例3.1（续）,2、 应用MATLAB进行计算机仿真，得到单位阶跃响 应曲线,单位阶跃响应曲线,A: 原三阶系统，超调量sp%=37%，调节时间ts=1.6秒,B:忽略极点的系统 超调量sp%=54.5% ，调节时间ts=1.5秒,D:忽略零极点的系统 超调量sp%=9.5%，调节时间ts=1.2秒,C:忽略零点的系统 超调量sp%=5.5%，调节时间ts=1.4秒,不能忽略零极点的影响,一个不能忽略的零点对系统的影响是使超调量加大，响应速度加快,一个不能忽略的极点对系统的影响是使超调量减小，调节时间增加,3.3. 反馈控制系统的稳态误差,瞬态响应性能指标,稳态误差,稳态误差是对系统精度的一种衡量,系统在不同的输入信号作用下，会有不同的稳态误差。,系统静特性不稳定和参数变化等因素也会导致系统产生一定的稳态误差,3.3.1. 稳态误差的概念 (Steady-State Error),图(3-23) 单位反馈系统,图(3-24) 非单位反馈系统,误差信号e(t)与输入信号r(t)之间的传递函数,终值定理,稳态误差与输入信号和系统的结构、参数有关,(a) 阶跃输入 (b) 斜坡输入 (c) 抛物线输入 图(3-25) 不同典型信号作用下的稳态误差,3.3.2. 稳态误差的计算,控制系统的开环传递函数,系统有N个积分环节串联。因为系统的类型常按其开环传递函数中串联积分环节的数目分类，所以称此系统为N型系统增加型号数，可使系统精度提高，但对稳定性不利，实际系统中N,（一）单位阶跃输入时的稳态误差,位置误差系数,0型系统，N = 0,1型或1型以上的系统， N ,（二）单位斜坡输入时的稳态误差,稳态误差,速度误差系数,0型系统，N = 0,1型系统， N =,2型或高于2型系统， N 2,（三）单位抛物线信号(等加速度信号)输入时的稳态误差,加速度误差系数,，稳态误差,0型或1型系统，N = 0或1,2型系统， N =2,3型或高于3型系统， N 2,稳态误差小结,当输入信号是上述典型信号的组合时，为使系统满足稳态响应的要求，N值应按最复杂的输入信号来选定,例3.2,解：系统a为1型系统，其Ka = 0，不能紧跟r(t)的3t2分量，所以 ess= ； 系统b为2型系统，其Ka = K = 10/4，所以,3.3.3. 主扰动输入引起的稳态误差,定义,一般情况下，系统除受到输入信号的作用外，还可能承受各种扰动信号的作用，在扰动信号的作用下，系统也将产生稳态误差，称为扰动稳态误差,图(3-27) 主扰动的影响,主扰动,主扰动输入引起的稳态误差（续）,R(s),输出与扰动之间的传递函数,扰动为单位阶跃信号n(t)=1(t),3.3.4. 关于降低稳态误差问题,增大系统开环放大系数可以增强系统对参考输入的跟随能力；增大扰动作用点以前的前向通道放大系数可以降低扰动引起的稳态误差 增加前向通道中积分环节数，使系统型号提高，可以消除不同输入信号时的稳态误差 保证元件有一定的精度和稳定的性能，尤其是反馈通道元件 采用复合控制来降低系统误差，消除扰动影响,按输入反馈按扰动顺馈的复合控制系统,顺馈控制器的传递函数,按参考输入顺馈的复合控制系统,温度复合控制系统,3.4. 劳斯-霍尔维茨稳定性判据,劳斯赫尔维茨(RouthHurwitz)判据，代数判据方法 根轨迹法，图解求特征根的方法 奈魁斯特(Nyquist)判据，基于复变函数理论的方法 李雅普诺夫方法 ，适用于线性系统和非线性系统,常用的稳定性分析方法,3.4.1. 稳定性(Stability)的概念,(a) 稳定的 (b) 不稳定的,定义,系统在受到外作用力后，偏离了正常工作点，而当外作用力消失后，系统能够返回到原来的工作点，则称系统是稳定的,在有界输入的作用下，其输出响应也是有界的。这叫做有界输入有界输出稳定，又简称为BIBO稳定,单输入单输出线性系统,系统稳定,特征方程,设方程有k个实根-pi(i=1，2，k)，r对共轭复数根(-s i±jw i ) (i=1，2，r)，k+2r=n,讨论,-pi 0，-s i 0(即极点都具有负实部)，系统是稳定的。但由于存在复数根的wi，系统的运动是衰减振荡的；若wi =0，则系统的输出按指数曲线衰减 若-pi或-s i中有一个或一个以上是正数，系统是不稳定的 只要-pi中有一个为零，或-s i中有一个为零(即有一对虚根)，系统处于稳定的临界状态,线性系统稳定的充分必要条件,所有特征根均为负实数，或具有负的实数部分,所有特征根，均在根平面的左半部分,所有极点均位于S平面的左半部分,=,=,系统稳定性的简单例子,3.4.2. 劳斯判据(Routh Criterion),(一)系统稳定性的初步判别,证明,特征方程有n个根，其中k个实根- p j (j=1，2，k)，r对复根-s i±jw i (i=1，2，r)，n = k+2r。则特征方程式可写为,假如所有的根均在左半平面，即- p j 0，s i0。所以将各因子项相乘展开后，特征方程的所有系数都是正数,(二) 劳斯判据,系统特征方程,劳斯阵列表,第一列系数均为正数，系统稳定， 第一列系数有负数，则第一列系数符号的改变次数等于在右半平面上根的个数,例3.3,解,2、劳斯阵列表如下 s4 1 12 6 s3 6 11 0 s2 61/6 6 s1 455/61 0 s0 6,1、系统特征方程所有系数均为正实数，满足系统稳定的必要条件,3、第一列系数均为正实数，故系统稳定,解,1、劳斯阵列表如下 s5 1 2 5 s4 3 1 6 s3 5 9 (各系数均已乘3) s2 -11 15 (各系数均已乘5/2) s1 174 (各系数均已乘11) s0 15,2、第一列有负数，系统不稳定。由于第一列系数的符号改变了两次(511 174)，系统特征方程有两个根的实部为正,例3.4,特殊情况 (1),劳斯阵列表中某一行的第一个系数为零，其余各系数不为零(或没有其余项),解,1、劳斯阵列表 s 3 1 1 s 2 2 2 s 1 e s 0 2,2、e 的上下两个系数(2和2)符号相同，则说明有一对虚根存在 ,系统处于临界状态,用一个很小的正数e 来代替零,特殊情况 (2),劳斯阵列表中某一行(设为第k行)的所有系数均为零，则说明在根平面内存在一些大小相等，并且关于原点对称的根,处理步骤,利用第k1行的系数构成辅助多项式，它的次数总是偶数的,求辅助多项式对s的导数，将其系数构成新行，代替第k行,继续计算劳斯阵列表,关于原点对称的根可通过令辅助多项式等于零求得,例3.6,解,1、劳斯阵列表如下 s 3 1 16 s 2 10 160 辅助多项式10s 2 + 160 s 1 0 0 求导数 20 0 构成新行 20s + 0 s 0 160,2、第一列各系数均未变号，所以没有特征根位于右半平面,例3.7,解,1、劳斯阵列表如下 s 5 1 3 -4 s 4 2 6 -8 辅助多项式 2s 4 + 6s 2 - 8 s 3 0 0 0 求导数 8 12 0 构成新行 8s 3 + 12s s 2 3 -8 s 1 100/3 s 0 -8,2、第一列变号一次，故有一个根在右半平面 ,系统不稳定,(三) 劳斯判据的应用1：稳定裕量的检验,如果所有根均在新虚轴的左边，则说系统具有稳定裕量s 1,例3.8,解,1、劳斯阵列表 s 3 2 13 s 2 10 4 s 1 12.2 s 0 4,2、第一列无符号改变，故没有根在s平面右半平面。,例3.8（续）,3、令s= z-1，代入特征方程式，得 即,4、新的劳斯阵列表 z 3 2 -1 z 2 4 -1 z 1 -1/2 z 0 -1,5、第一列符号改变一次，故有一个根在直线s= -1的右边，因此稳定裕量不到1,2. 分析系统参数对稳定性的影响,劳斯阵列表：s 3 1 5 s 2 6 K s 1 s 0 K,例3.9,解,1、劳斯阵列表 s 4 1 T 100 s 3 2 10 s 2 T - 5 100 s 1 s 0 100,2、,3. 鉴别延滞系统的稳定性,注意,应用劳斯判据来判别延滞系统的稳定性，需要采用近似的方法处理,用有理函数近似,(1)用有限项简单有理函数的乘积近似,取n为有限值,为了保证一定的准确度，n 值往往较大，分析起来是相 当麻烦的,(2)用有理分式近似,泰勒级数,可用一个有理分式p(s)/q(s)来近似e-t s：派德(pade)近似,分子为m次，分母为n次 ，m和n的选取，应在满足近似准确度要求的前提下，尽可能少增加特征方程式的阶次,表3.5,n=1，m=3,K的临界值为1.13,3.4.3. 赫尔维茨判据,系统特征方程式为 ansnan-1sn-1a1sa0= 0,赫尔维茨判据：系统稳定的必要和充分条件是an0的情况下，对角线上所有子行列式i (i=1，2，,n)均大于零。,赫尔维茨判据(Hurwitz Criterion),n=1,a10，1= a0,n=2,n=3,n=4,例3.10,解,1、从特征方程式看出所有系数为正数，满足稳定的必要条件,2、赫尔维茨行列式,系统是稳定的,3.5. 控制系统灵敏度分析 (sensitivity),|G(s)H(s)|1,增加开环传递函数G(s)H(s)的大小会导致G(s)对输出影响减少,G(s)+G(s),C(s)=G(s)R(s),开环情况,闭环情况,G(s)H(s)G(s)H(s),控制系统灵敏度定义,定义,系统传递函数的变化率与被控过程传递函数变化率的比值,GB(s)=C(s)/R(s),开环系统的灵敏度等于1,闭环系统,反馈系统关于G (s)的灵敏度,闭环系统对反馈环节H(s) 灵敏度,控制系统引入反馈环节后能减少因参数变化而造成的影响，尤其是因被控过程参数变化所造成的影响,利用反馈减少灵敏度的简单例子,k = R1/Rf,A1时,当运算放大器处于开环状态(即无反馈电阻Rf )时，相对于增益A的开环灵敏度为1,闭环,A=104而且k = 0.1,闭环时相对于因子k的灵敏度,本章小结：,讨论了控制系统分析的基本内容；,稳定性、瞬态性能、稳态性能,开环系统：设备少、成本低，控制方法简单易行,反馈控制：成本大，系统变得复杂、原稳定的开环系统，由于反馈的引入，可能造成不稳定。,优点： 减少被控过程G(s)中参数变化时系统的灵敏度； 有力于控制系统的瞬态响应性能； 提高系统对于干扰的抑制力； 减少或消除系统的静态误差。,