1、x第四章系统稳定性及其李雅普诺夫稳定第四章系统稳定性及其李雅普诺夫稳定4- 1稳定性一般概念4-2李雅普诺夫稳定性定义4-3李雅普诺夫第一法(间接法)4-4李雅普诺夫第二法(直接法)4-5线性定常系统的李雅普诺夫稳定分析及系统参数优化4- 1稳定性一般概念对于一个实际的控制系统，其工作的稳定性无疑是一个极其重要的问题， 因为一个不稳定的系统在实际应用中是很难有效地发挥作用的。 从直观上看，系统的稳定性就是一个处于稳态的系统，在某一干扰信号的作用下，其状态偏离了原有平衡位置，如果该系统是稳定的， 那么当干扰取消后有限的时间内，系统会在自身作用下回到平衡状态；反之若系统不稳定， 则系统永远不会回到
2、原来的平衡位置。系统的稳定一般有外部稳定和内部稳定两种。 外部稳定又称作输出稳定， 也就是当系统在干扰取消后，在一定时间内，其输出会恢复到原来的稳态输出。 输出稳定有时描述为系统的BIBO稳定，即有限的系统输入只能产生有限的系统输出。 系统内部稳定主要针对系统内部状态，反映的是系统内部状态受干扰信号的影响。 当扰动信号取消后， 系统的内部状态会在一定时间内恢复到原来的平衡状态，则称系统状态稳定。在经典控制论中，研究对象都是用高阶微分方程或传递函数描述的单输入单输出(SISO)系统，反映的仅是输入输出的关系，不会涉及系统内部的状态。因此经典控制论 中只讨论系统的输出稳定问题。系统的稳定性是系统本
3、身的特性，与系统的外部输入(控制)无关。在经典控制论中， 我们通过研究线性定常系统的特征根的情况来判断系统的输出稳定性： 如果系统的特征根都有负的实部(即都在复平面的左部)，则系统输出稳定。对于n阶线性连续系统，其特征方程为：+爲-丄+ . +內尹十。(4 1)当n4时，要求出其所有特征根是非常困难的，从而要想通过解出高阶系统的特征根 来判别系统稳定性也是不现实的。 所以1877年劳斯(Routh )和1895年霍尔维茨(Hurwitz )分别提出了有名的劳斯一霍尔维茨稳定判据， 它可以通过线性定常系统特征方程的系数的简单代数运算来判别系统输出稳定性，而不必求出各个特征根。有关 Routh-H
4、urwitz判据的详细内容请参阅有关经典控制论教材。当系统不是线性定常系统时，或者对于系统内部状态稳定问题， 经典控制论中的方法就不好解决了，这就需要下面介绍的李雅普诺夫( Lyapu nov )稳定性的理论。4-2李雅普诺夫稳定性定义4-2 1系统的平衡状态设控制系统的齐次状态方程为:其中，X(t)为系统的n维状态向量, 不一定是线性定常的。如果对于所有f(Xe,t)三 0 (4 3)则我们称Xe为系统的平衡状态。对于一般控制系统，它可能没有，也可能有一个或多个平 衡状态。如果系统是一个线性定常系统，即：X = AX那么当A为非奇异时，Xe=0是系统的唯一平衡状态； 当A为奇异矩阵时，AX=
5、0有无数解,也就是系统有无数个平衡状态。系统的状态稳定性是针对系统的平衡状态的， 当系统有多个平衡状态时， 需要对每个平衡状态分别进行讨论。对系统矩阵A非奇异的线性定常系统，Xe =0是系统的唯一平衡状态， 所以对线性定常(LTI)系统，我们一般笼统用 Xe的稳定性代表系统稳定性。42 2李雅普诺夫稳定假设(4 2)所示一般控制系统的解为:(4 4)它是与初始时间t0及其初始状态 X0有关的，体现系统状态从(t0, X0)出发的一条状态轨 迹。设Xe为系统的一个平衡点， 如果给定一个以Xe为球心，以为半径的n维球域S( e ), 总能找到一个同样以 Xe为球心，S (e, t0)为半径的n维球
6、域S(3 )，使得从S(3 )球域出 发的任意一条系统状态轨迹 0 (t; X。，t。)在tt0的所有时间内，都不会跑出 S(e )球域，则 称系统的平衡状态 Xe是李雅普诺夫稳定 的(Lyapunov Stability) o般来说，S的大小不但与e有关，而且与系统的初始时间 t有关。当S仅与e有关时，称X e是一致稳定 的平衡状态。进一步的，如果 Xe不仅是李雅普诺夫稳定的平衡状态，而且当时间 t无限增加时，从S(S )球域出发的任一条状态轨迹 0 (t； X0, t0)都最终收敛于球心平衡点 Xe，那么称Xe是渐进稳定 的(Asymptotic Stability)。更进一步，如果从s(
7、g),即整个系统状态空间的任一点出发的任一条状态轨迹 o (t； X0,t0),当tis时，都收敛到平衡点 Xe，那么称Xe是大范围渐进稳定 的。很明显，这时的 Xe 是系统的唯一的平衡点。反之，对于给定 S(e )，不论S 0取得多么小，从 S(S )球域出发的状态轨迹 0 (t； X0, t0)，至少有一条跑出 S( e )球域，那么称平衡点 Xe是不稳定的。4-3李雅普诺夫第一法（间接法）李雅普诺夫第一法通过分析系统微分方程的显式解来分析系统的稳定性， 对线性定常系统，它可以直接通过系统的特征根情况来分析。 李雅普诺夫第一法的基本思路与经典控制论中的稳定性判别思路基本一致。设线性定常系统
8、的动态方程为：(4 - 5)在讨论系统状态稳定性（内部稳定）时，可以不考虑系统的输入结构和输入信号 u,只 从系统的齐次状态方程或矩阵 A出发。很明显，当-时，Xe= 0是系统的唯一平衡点。对于Xe= 0的稳定性，我们有如下判据（Xe大范围渐进稳定的充要条件）：当线性定常系统的系统矩阵 A的所有特征根都有负的实部时，其唯一的状态平衡点 Xe=0是渐进稳定的，而且是大范围渐进稳定。对于（4 5）所示系统，其输入输出的传递函数为：当W（S）的极点全部都有负实部时，该系统有界的输入将引起有界的输出（ BIBO ）,也就是说系统是输出稳定的。可以证明，当（4 5）式所示系统的传递函数 W（s）没有零极
9、点对消时，系统的状态稳定性和系统的输出稳定性是一致的，因为这时系统矩阵的特征根就是系统传递函数的极点。【例4 1】 某系统的动态方程为:请分析系统的状态稳定性和输出稳定性。解:Xe=0是系统的唯一平衡点。 系统特征方程为：円_卜几；1 二=（2 + 1）（1-1）-0因为有大于零的特征根，所以系统状态不稳定。系统传递函数为:极点S= -10 ,则称V(x)是正定的。2.V(x) 0,则称V(x)是半正定的(非负定)。3.V(x)0 ,称V(x)是负定的。4.V(x) 04 42李雅普诺夫第二法设系统状态方程为：(4 8)左二f仅上) 那么根据其中Xe=0为系统的一个平衡状态。如果存在一个正定的
10、标量函数 V(x),并且具有连续的一阶偏导数, 符号性质，我们有：1若，：，则j -不稳定；2若则李雅普诺夫稳定；3若一，或者但时.不恒为零，则渐进稳定；4若厂 渐进稳定，并且当 时，- ，则大范围渐进稳定。应当指出，上述稳定性判据只是一个充分条件，并不是必要条件。如果给定的 V(x)满足上述四个条件之 上述任何一个条件，那么只能说明所选的 V(x)对(4 8)所示系统失效，必须重新构造 V(x)。【例4 - 2】控制系统齐次状态方程为：1 -X(4 9)- 请研究其状态稳定性。解: -雇系统的唯一平衡状态当取 V (x)二 2x1 x2dv/ dv/很明显，V(x)0，并且 连续。沪艾)三4
11、巧右+ 2x2 =4工工2 - 2工尹2工f = 2x2 - 2xa3可以看出，.是不定的，所以上述 V(x)不是一个李雅普诺夫函数，它不能用来判断(4 9)式所表示的系统的稳定性。我们重新选择很明显，它是正定的，并且有连续的一阶偏导数。矶瓦)=细缶+ 希=细西-2?存i TjQf 级 兰0 (410)果T -,则由(4 10)可知，也即，由状态方程； - 可知，,), 所以-是渐进稳定的。更进一步，因为当 时，小；二；：丄：二二，所以- - 1是大范围渐进稳定的。为了避免上述复杂的补充说明，我们再取下列二次型函数：咚”*坷+伺+习尸+巧勺巩r)= 2昭1十(心+七X右+坊)诙=(xi +工J
12、) 0。= ifpx + 丈 P文=(AXfPX + jfP(AX) = + fPAX=疋(且9 + PA)X式(4 13)称作 李雅普诺夫方程。因此得:V(x)=-丈 QN其中；二为对称矩阵。如果 Q0，那么V(x) ： 0 ,所以Xe=0为渐进稳定,而且是大范围渐进稳定。在实际应用中，我们是先给定一正定矩阵 Q，然后通过李雅普诺夫方程 (4 - 13)求出实对称矩阵P，最后通过赛尔维斯特准则判别 P的正定性。如果 P0,则系统稳定。在应用李雅普诺夫方程时，应注意下列几点：(1)由李雅普诺夫方程求得的 P为正定是Xe=0渐进稳定的充要条件。(2)Q的选取是任意的，只要满足对称且正定(在一定条
13、件下可以是半正定) ，Q的选取不会影响系统稳定性判别的结果。(3)当XM 0时.不恒为零，那么 Q可以取半正定阵，即 Q 0。(4)当Q取为单位阵I时，李雅普诺夫方程变为：- / ( 414)这是一个比较简单的李雅普诺夫方程。【例4 3】 设系统状态方程为： L1 一1，试分析其稳定性。解：很明显Xe=0为系统唯一平衡点，用李雅普诺夫方程判断其稳定性。取Q=I,贝李氏方程为： 0 1 T11Pn+11Pl2_ 0 1 _-1 0_-1 -1P2P22_卩22 一-1 -10 -1考虑到P12二P21，所以化简得-212 Pu-Pu 011 P12 P22 （.Pl2P22）J 0 一1f- 2
14、去埠=_1lAi-Pu - = 02（阳-戸33-I解得，Pll = % 戸22=13 r2 21 1所以,2根据赛尔维斯特准则，纠 了 A D Aa = - 0可知P0，所以Xe=0渐进稳定。0 0如果取e4o J-0,那么(x) = = -X 益 0且当一： H时，厂二，由状态方程込=一孔一xt可知，： 。所以只有原点才使得 $弐J，所以Q可以取半正定。因此代入李雅普一 2 Pit = 0-Fu - -Pu - -Pii = 目皿-屉) _1求得：1加=2Pu 01所以Xe=0渐进稳定。图4 1控制系统方块图解:图。取图42各积分器输出为状态变量，由图可写出齐次状态方程为:首先求出图41方
15、块图所示系统的状态方程。图 4 1可等效变换为图42所示方块= -2卷 4- 绻=一 g -利用李雅普诺夫方程，并取: 求得： 0 1 0_I耳2刃刃 1 Pl2 Pl3 0 1 0_0 0 0_0 -2 1+及戈忽0 -2 10 0 0一疋 0 -1_P31 P32 P53_戶1 P32刃屯一上 0 -10 0-1由式（4 15）可知,xi=0, X2=0 ,所以Q 0取法合理。将Q,A代入李雅普诺夫方程:因为J -时X3=0 ,求得:要使上-渐进稳定，P必须正定，即:【例4- 5】利用李雅普诺夫方程求系统参数优化问题。求使J最小的系统参数K的值。解:-为系统唯一的平衡点。根据题意，一-是大范围渐进稳定的。根据李雅普诺 夫方程，并取题中的 Q0代入李雅普诺夫方程，得:0-1解得:12 加211 +戸22k考虑到大于零，所以：鸟吶冷弓L】m所以，P0,即稳定。所以 .一是系统的一个李雅普诺夫函数，而且 。再看系统性能指标J:J =QXdt = V（x）dt =-叭x）|：=-疋刃韦 =x1 p）px（q）- 疋（b）px（b）由于系统在-处是大范围渐进稳定的，所以 ：二； ，所以：(OjPXfO) = I丄+戸+ 一2A112丄+戶2A -2所以系统的最优参数不但与性能泛函中的参数 有关，而且与系统的初始状态Xo有关。如果取：-，则:
