劳斯-霍尔维茨稳定性判据分解

第四节劳斯-霍尔维茨稳定性判据稳定性是控制系统最重要的问题，也是对系统最起码的要求。控制系统在实际运行中，总会受到外界和内部一些因素的扰动，例如负载或能源的波动、环境条件的改变、系统参数的变化等。如果系统不稳定，当它受到扰动时，系统中各物理量就会偏离其平衡工作点，并随时间推移而发散，即使扰动消失了，也不可能恢复原来的平衡状态。因此，如何分析系统的稳定性并提出保证系统稳定的措施，是控制理论的基本任务之一。常用的稳定性分析方法有：1.劳斯－赫尔维茨（Routh－Hurwitz）判据这是一种代数判据方法。它是根据系统特征方程式来判断特征根在S平面的位置，从而决定系统的稳定性.2.根轨迹法这是一种图解求特征根的方法。它是根据系统开环传递函数以某一（或某些）参数为变量作出闭环系统的特征根在S平面的轨迹，从而全面了解闭环系统特征根随该参数的变化情况。3.奈魁斯特（Nyquist）判据这是一种在复变函数理论基础上建立起来的方法。它根据系统的开环频率特性确定闭环系统的稳定性，同样避免了求解闭环系统特征根的困难。这一方法在工程上是得到了比较广泛的应用。4.李雅普诺夫方法上述几种方法主要适用于线性系统，而李雅普诺夫方法不仅适用于线性系统，更适用于非线性系统。该方法是根据李雅普诺夫函数的特征来决定系统的稳定性。一、稳定性的概念稳定性的概念可以通过图3-31所示的方法加以说明。考虑置于水平面上的圆锥体，其底部朝下时，若将它稍微倾斜，外作用力撤消后，经过若干次摆动，它仍会返回到原来状态。而当圆锥体尖部朝下放置时，由于只有一点能使圆锥体保持平衡，所以在受到任何极微小的扰动后，它就会倾倒，如果没有外力作用，就再也不能回到原来的状态了。(a)稳定的(b)不稳定的图3-31圆锥体的稳定性根据上述讨论，可以将系统的稳定性定义为，系统在受到外作用力后，偏离了正常工作点，而当外作用力消失后，系统能够返回到原来的工作点，则称系统是稳定的。瞬态响应项不外乎表现为衰减、临界和发散这三种情况之一，它是决定系统稳定性的关键。由于输入量只影响到稳态响应项，并且两者具有相同的特性，即如果输入量r(t)是有界的：|r(t)|∞，t≥0则稳态响应项也必定是有界的。这说明对于系统稳定性的讨论可以归结为，系统在任何一个有界输入的作用下，其输出是否有界的问题。一个稳定的系统定义为，在有界输入的作用下，其输出响应也是有界的。这叫做有界输入有界输出稳定，又简称为BIBO稳定。线性闭环系统的稳定性可以根据闭环极点在S平面内的位置予以确定。假如单输入单输出线性系统由下述的微分方程式来描述，即(3.58)则系统的稳定性由上式左端决定，或者说系统稳定性可按齐次微分方程式(3.59)来分析。这时，在任何初始条件下，若满足(3.60)()(1)(1)110()(1)(1)110nnnnmmmmacacacacbrbrbrbr()(1)(1)110nnnnacacacac0(1)(1)lim()lim()lim()ntttctctct0则称系统(3.58)是稳定的。为了决定系统的稳定性，可求出式(3.59)的解。由数学分析知道，式(3.59)的特征方程式为(3.61)设上式有k个实根pi（i=1，2，…，k），r对共轭复数根(si±jwi)(i=1，2，…，r)，k+2r=n，则齐次方程式(3.59)解的一般式为(3.62)式中系数Ai，Bi和Ci由初始条件决定。从式(3.62)可知：(1)若pi0，si0（即极点都具有负实部），则式(3.60)成立，系统最终能恢复至平衡状态，所以系统是稳定的。1110nnnnasasasa011()(cossin)iikrpttiiiiiiictCeeAtBtsww(3)若-pi或-si中有一个或一个以上是正数，则式(3.60)不满足。当t→∞时，c(t)将发散，系统是不稳定的。(4)只要-pi中有一个为零，或-si中有一个为零（即有一对虚根），则式(3.60)不满足。当t→∞时，系统输出或者为一常值，或者为等幅振荡，不能恢复原平衡状态，这时系统处于稳定的临界状态。总结上述，可以得出如下结论：线性系统稳定的充分必要条件是它的所有特征根均为负实数，或具有负的实数部分。根平面上是一个点，所以上述结论又可以这样说：线性系统稳定的充分必要条件是它的所有特征根，均在根平面的左半部分（见图3-32）。图3-32根平面表3.4列举了几个简单系统稳定性的例子。需要指出的是，对于线性定常系统，由于系统特征方程根是由特征方程的结构（即方程的阶数）和系数决定的，因此系统的稳定性与输入信号和初始条件无关，仅由系统的结构和参数决定。如果系统中每个部分都可用线性常系数微分方程描述，那么，当系统是稳定时，它在大偏差情况下也是稳定的。如果系统中有的元件或装置是非线性的，但经线性化处理后可用线性化方程来描述，则当系统是稳定时，我们只能说这个系统在小偏差情况下是稳定的，而在大偏差时不能保证系统仍是稳定的。程根，假如特征方程根能求得，系统稳定性自然就可断定。但是，要解四次或更高次的特征方程式，是相当麻烦的，往往需要求助于数字计算机。所以，就有人提出了在不解特征方程式的情况下，求解特征方程根在S平面上分布的方法。下面就介绍常用的劳斯判据和赫尔维茨判据。二、劳斯判据(一)系统稳定性的初步判别已知系统的闭环特征方程为(3.63)式中所有系数均为实数，且an0，则系统稳定的必要条件是上述系统特征方程的所有系数均为正数。证明如下：设式(3.63)有n个根，其中k个实根pj(j=1，2，…，k)，r对复根si±jwi(i=1，2，…，r)，n=k+2r。则特征方程式可写为1110()nnnnDsasasasa0111022221211()()()()[()][()]nnnnnkrrDsasasasaaspspspssswsw0假如所有的根均在左半平面，即pj0，si0，则pj0，si0。所以将各因子项相乘展开后，式（3.63）的所有系数都是正数。根据这一原则，在判别系统的稳定性时，可首先检查系统特征方程的系数是否都为正数，假如有任何系数为负数或等于零（缺项），则系统就是不稳定的。但是，假若特征方程的所有系数均为正数，并不能肯定系统是稳定的，还要做进一步的判别。因为上述所说的原则只是系统稳定性的必要条件，而不是充分必要条件。(二)劳斯判据这是1877年由劳斯(Routh)提出的代数判据。1.若系统特征方程式设an0，各项系数均为正数。2.按特征方程的系数列写劳斯阵列表：1110nnnnasasasa02411352123312341231101nnnnnnnnnnnsaaasaaasbbbscccsdddsfsg表中直至其余bi项均为零。2113142151311111nnnnnnnnnnnnnnnaabaaaaabaaaaabaaa67131121152131173141111nnnnnnaacbbbaacbbbaacbbb按此规律一直计算到n-1行为止。在上述计算过程中，为了简化数值运算，可将某一行中的各系数均乘一个正数，不会影响稳定性结论。3.考察阵列表第一列系数的符号。假若劳斯阵列表中第一列系数均为正数，则该系统是稳定的，即特征方程所有的根均位于根平面的左半平面。假若第一列系数有负数，则第一列系数符号的改变次数等于在右半平面上根的个数。例3.3系统特征方程为试用劳斯判据判别系统的稳定性。解从系统特征方程看出，它的所有系数均为正实数，满足系统稳定的必要条件。列写劳斯阵列表如下432ssss61211601126611061/66455/6106第一列系数均为正实数，故系统稳定。事实上，从因式分解可将特征方程写为其根为2，3，，均具有负实部，所以系统稳定。(s+2)(s+3)(s2+s+1)=013j220s1s2s3s4s例3.4已知系统特征方程式为解列写劳斯阵列表12531659（各系数均已乘3）-1115（各系数均已乘5/2）174（各系数均已乘11）15劳斯阵列表第一列有负数，所以系统是不稳定的。由于第一列系数的符号改变了两次(5→－11→174)，所以，系统特征方程有两个根的实部为正。54320sssss32565s4s3s2s1s0s4.两种特殊情况在劳斯阵列表的计算过程中，如果出现：(1)劳斯阵列表中某一行的第一个系数为零，其余各系数不为零（或没有其余项），这时可用一个很小的正数e来代替这个零，从而使劳斯阵列表可以继续运算下去（否则下一行将出现∞）。如果e的上下两个系数均为正数，则说明系统特征方程有一对虚根，系统处干临界状态；如果e的上下两个系数的符号不同，则说明这里有一个符号变化过程，则系统不稳定，不稳定根的个数由符号变化次数决定。例3.5设系统特征方程为s3+2s2+s+2=0解劳斯阵列表为由于e的上下两个系数(2和2)符号相同，则说明有一对虚根存在。上述特征方程可因式分解为(2)若劳斯阵列表中某一行（设为第k行）的所有系数均为零，则说明在根平面内存在一些大小相等，并且关于原点对称的根。在这种情况下可做如下处理：a.利用第k－1行的系数构成辅助多项式，它的次数总是偶数的；1122e22(1)()0ss20s1s2s3sb.求辅助多项式对s的导数，将其系数构成新行，代替第k行；c.继续计算劳斯阵列表；d.关于原点对称的根可通过令辅助多项式等于零求得。例3.6系统特征方程为解劳斯阵列表为11610160辅助多项式10+16000↓求导数200构成新行20s+016032sss101616000s1s2s3s2s从上表第一列可以看出，各系数均未变号，所以没有特征根位于右半平面。由辅助多项式知道10s2+160=0有一对共轭虚根为±j4。例3.7特征方程式为解劳斯阵列表如下：13-426-8辅助多项式2s4+6s2-8000↓求导数8120构成新行8s3+12s3-8100/3-84240sssss5323680s1s2s3s4s5s劳斯阵列表第一列变号一次，故有一个根在右半平面。由辅助多项式：可得s1,2=±１，s3,4=±j2，它们均关于原点对称，其中一个根在S平面的右半平面。(三)劳斯判据的应用应用劳斯判据不仅可以判别系统稳定不稳定，即系统的绝对稳定性，而且也可检验系统是否有一定的稳定裕量，即相对稳定性。另外劳斯判据还可用来分析系统参数对稳定性的影响和鉴别延滞系统的稳定性。2s4+6s2-8=01.稳定裕量的检验如图3-33所示，令(3.64)即把虚轴左移s1。将上式代入系统的特征方程式，得以z为变量的新特征方程式，然后再检验新特征方程式有几个根位于新虚轴(垂直线s=s1)的右边。如果所有根均在新虚轴的左边（新劳斯阵列式第一列均为正数），则说系统具有稳定裕量s1。s=z-s1图3-33稳定裕量s1例3.8检验特征方程式是否有根在右半平面，并检验有几个根在直线s=-１的右边。解劳斯阵列表为21310412.24第一列无符号改变，故没有根在S平面右半平面。再令s=z-1，代入特征方程式，得即320sss21013432(1)(1)(1)40zzz2101332410zzz23s2s1s0s则新的劳斯阵列表从表中可看出，第一列符号改变一次，故有一个根在直线s=-（即新座标虚轴）的右边，因此稳定裕量不到1。2.分析系统参数对稳定性的影响设一单位反馈控制系统如图3-34所示，其闭环传递函数为系统的特征方程式为z32-1z24-1z1-1/2z0-1()()()(1)()BCsKGsRssssK5320sssK65图3-34求K的范围列写劳斯阵列表：s315s26Ks1s0KK306若要使系统稳定，其充要条件是劳斯阵列表的第一列