李雅普诺夫稳定性分析

1、第六章李雅普诺夫稳定性分析在反馈控制系统的分析设计中，系统的稳定性是首先需要考虑的问题之一。因为它关系到系统是 否能正常工作。经典控制理论中已经建立了劳斯判据、Huiwitz稳定判据、Nquist判据、对数判据、根轨迹判据等来判断线性定常系统的稳定性，但不适用于非线性和时变系统。分析非线性系统稳定性及自振的描述 函数法，则要求系统的线性部分具有良好的滤除谐波的性能；而相平面法则只适合于一阶、二阶非线 性系统。1892年俄国学者李雅普诺夫(Lyapunov)提出的稳定性理论是确定系统稳定性的更一般的理论， 它采用状态向量来描述，不仅适用于单变量、线性、定常系统，还适用于多变量、非线性、时变系统。

2、§ 6-1外部稳定性和内部稳定性系统的数学模型有输入输出描述(即外部描述)和状态空间描述(即内部描述)，相应的稳定性便分为外部稳定性和内部稳定性。一、外部稳定性1、定义(外部稳定性)：若系统对所有有界输入引起的零状态响应的输出是有界的，则称该系统是外部稳定的。(外部稳定性也称为 BIBO( Boun ded In put Bou nded Output)稳定性)说明：(1) 所谓有界是指如果一个函数h(t)，在时间区间0,：中，它的幅值不会增至无穷，即存在一个实常数k，使得对于所有的t e 0 co ,恒有h(t)兰k £ 成立。(2) 所谓零状态响应，是指零初始状态时非零

3、输入引起的响应。2、系统外部稳定性判据线性定常连续系统、(A,B,C)的传递函数矩阵为x 二 Ax Buy =CxsX =AX BUY =CX(sI -A)X =BUX =(sI -A)BUG(s) =C(sI _ A) / B当且仅当G(s)极点都在s的左半平面内时，系统才是外部稳定(或BIBO稳定)的。【例6.1.1】已知受控系统状态空间表达式为y 二'0 1x0 6 -2X = IX 十 I u'1 -t 11 一试分析系统的外部稳定性。一621 s+計 b J解：系统为SISO系统，传递函数为_ sG（S） =C（sl _ A）B = 0 1 1-1s 2= （s2）（

4、s 3）1一 s+3由于传递函数的极点位于s左平面，故系统是外部稳定的。二、内部稳定性对于线性定常系统x = Ax Bu， x（t0） = x0y = Cx如果外部输入U（t）二0，初始条件x0为任意，且由x0引起的零输入响应为x（t）二（t,to）Xo满足lim （t,to）xo 二 0t ）:则称系统是内部稳定的，或称为系统是渐近稳定的。说明：线性定常系统的渐近稳定与经典控制理论中的稳定性一致。【例6.1. 2】已知受控系统状态空间表达式为0 6 -2x=h _irii，试分析系统的内部稳定性。解：该系统为线性定常系统，其特征方程为：树 一 A =入（九 + 1） 6 =（九-2）（人 +

5、 3） = 0于是系统的特征值为'1 = 2 , ' 2 = 3，故系统不是内部稳定（渐近稳定）的。三、内部稳定性与外部稳定性的关系1、若系统是内部稳定（渐近稳定）的，则一定是外部稳定（BIB O稳定）的。2、 若系统是外部稳定（BIB O稳定）的，且又是可控可观测的，则系统是内部稳定（渐近稳定）的。此 时内部稳定和外部稳定是等价的。§ 6-2李雅普诺夫稳定性的基本概念、自治系统没有外界输入作用的系统叫自治系统。自治系统可用如下的显含时间t的状态方程来描述X = f (X,t) , x(to) = X。, t _to (6-1)其中X为n维状态向量。f (x,t)为线

6、性或非线性、定常或时变的n维向量函数。假定方程的解为X(t；Xo,t。),式中X。和to分别为初始状态向量和初始时刻，那么初始条件X。必满足X(t。； Xo,t。)=X。如果系统为线性系统，则(6-1)方程中的f(x,t)为X的线性向量函数，或按习惯表示为：X 二 A(t)x ,x(to)二 X。, t _to (6-2)二、平衡状态设控制系统的齐次状态方程为：X 二 f (x,t) , x(t。)= X。, t _t。对于所有t，如果存在某个状态xe ,满足：Xe 二 f(Xe，t)二。则称Xe为系统的一个平衡点或平衡状态。平衡状态的各分量相对时间不再发生变化。若已知系统状态方程，令X二。所

7、求得的解x，便是平衡状态。在大多数情况下，Xe =。(状态空间原点)为系统的一个平衡状态。当然，系统也可以有非零平衡状态。如果系统的平衡状态在状态空间中表现为彼此分隔的孤立点，则称其为孤立平衡状态。对于孤立平衡状态，总是可以通过移动坐标系而将其转换为状态空间的原点，所以在下面的讨论中，假定原点即Xe =。为平衡状态。所谓系统运动的稳定性，就是研究其平衡状态的稳定性，也即偏离平衡状态的受扰运动，能否只依靠系统内部的结构因素而返回到平衡状态，或者限制在平衡状态的附近。线性定常系统X = AX，其平衡状态满足 AXe二。，只要A非奇异，系统只有唯一的零解，即存 在一个位于状态空间原点的平衡状态；当

8、A为奇异矩阵时， Axe =。有无数解，也就是系统有无数个 平衡状态。对于非线性系统，f(xe,t)=。的解可能有多个，由系统状态方程决定。、李雅普诺夫意义下稳定设系统初始状态 X0位于以平衡状态 Xe为球心、半径为：的闭球域SG）内，即X0 -'xe 一 -；（ ；,to）t = to若能使系统方程的解X（t；Xo，to）在t-:的过程中，都位于以Xe为球心、任意规定的半径为:的闭球 域S（；）内，即X（t；X°,to） - Xe 乞；t -to则称该Xe是稳定的，通常称 Xe为李雅普诺夫意义下稳定的平衡状态。以二维系统为例，上述定义的平 面几何表示如图6-1所示。X2s(

9、；)XesC)XiXo-初始状态Xe-平衡状态图6-1二维空间李雅普诺夫意义下稳定性的几何解释示意图式中*称为向量的范数，其几何意义是空间距离的尺度。如Xo-Xe表示状态空间中Xo至Xe点之间的距离的尺度，其数学表达式为Xo _Xe、(Xio - Xie)2 (Xno - Xne)2在上述稳定性的定义中，如果、：只依赖于:而和初始时刻to的选取无关，则称平衡状态Xe是一致稳定的。对于定常系统，Xe的稳定等价于一致稳定。但对于时变系统，Xe的稳定并不意味着其为一致稳定。要注意到，按李雅普诺夫意义下的稳定性定义，当系统作不衰减的振荡运动时，将在平面描绘出一条封闭曲线，但只要不超过 S（），则认为稳

10、定，这同经典控制理论中线性定常系统稳定性的定义是有差异的。四、渐近稳定设Xe是系统X = f （X,t），x（to） =xo , t -to的一个孤立平衡状态，如果（1）Xe是李雅普诺夫意义下稳定的;则称此平衡状态是渐近稳定的。Xo-初始状态Xe-平衡状态图6-2二维空间渐近稳定性的几何解释示意图实际上，渐近稳定即为工程意义下的稳定，也就是经典控制理论中所讨论的稳定性。当关时，称平衡状态Xe是一致渐近稳定的。五、大范围（全局）渐近稳定X e是大范围渐近稳定的。这是因为线性系统稳定性，其:.总是有限的，故致渐近稳定。当初始条件扩展到整个状态空间，且具有渐近稳定性时，称此平衡状态对于严格线性系统，

11、如果它是渐近稳定的，必具有大范围渐近稳定性， 与初始条件的大小无关。一般非线性系统的稳定性与初始条件的大小密切相关通常只能在小范围内渐近稳定。当、：与t0无关时，称平衡状态 Xe是大范围六、不稳定不管把域s（j取得多么小，也不管把域s（；）取得如何的大，只要在s（j内存在一个非零初始状态Xo，使得有Xo出发的运动轨迹超出域 S（；）以外，则称平衡状态 Xe是不稳定的。线性系统的平衡状态不稳定，表征系统不稳定。非线性系统的平衡状态不稳定，只说明存在局部 发散的轨迹，至于是否趋于无穷远，要看S（；）域外是否存在其它平衡状态，若存在，如有极限环，则系统仍是李雅普诺夫意义下稳定的。下面介绍李雅普诺夫理

12、论中判断系统稳定性的方法。§ 6-3李雅普诺夫稳定性判别方法一、李雅普诺夫第一法（间接法）这是利用状态方程解的特性来判别系统稳定性的方法，它适用于线性定常、线性时变以及非线性 函数可线性化的情况。由于本章主要研究线性定常系统，所以在此仅介绍线性定常系统的特征值判据。线性定常系统的特征值判据：对于线性定常系统 x = Ax, x(0) = X。, t _0有(1 )系统的平衡状态是在李雅普诺夫意义下稳定的充分必要条件是，A的所有特征值均具有非正(负或零)实部，且具有零实部的特征值为A的最小多项式的单根。(2) 系统的惟一平衡状态 Xe=0是渐近稳定的充分必要条件是，A的所有特征值均具有

13、负实部。二、李雅普诺夫第二法(直接法)根据古典力学中的振动现象，若系统能量(含动能与位能)随时间推移而衰减，系统迟早回到达平衡状态，但要找到实际系统的能量函数表达式并非易事。李雅普诺夫提出，可虚构一个能量函数(后来被称为李雅普诺夫函数)，一般它与XX2，Xn及t有关，记为V(X,t)。若不显含t,则记为V(X)。它是一个标量函数，考虑到能量函数总是大于零，故为正定函数。能量衰减特性用 V(x,t)或V(x)表示。李雅普诺夫第二法利用 V及V的符号特征，直接对平衡状态稳定性作出判断，无需求出系统状态方程的解，故称直接法。用此方法解决了一些用其它稳定性判据难以解决的非线性系统的稳定性问题，遗 憾的

14、是对一般非线性系统仍未形成构造李雅普诺夫函数的通用方法。对于线性系统，通常用二次型函 数XT Px作为李雅普诺夫函数。1、标量函数V(x)符号性质的几个定义(1) 正定性标量函数V(x)在域s中对所有非零状态(x = 0)有V(x) 0且V(0) = 0，则称V(x)在域s内Xo-初始状态2 2 正定。如 V (x)=捲x2是正定的。(2) 负定性标量函数V(x)在域s中对所有非零状态(x = 0)有V(x) ： 0且V(0) = 0，则称V(x)在域s内负定。如V(x)二-(X； ' x；)是正定的。(3) 正半定性V(0) = 0，且标量函数 V(x)在域s内某些非零状态处有V(x

15、) = 0，而在其它非零状态处有V(x) 0且，则称V(x)在域s内正半定。女口 V(x) = (Xi 2x2)2，当Xi = -2X2时有V(x)=0 ；当X1 - 2x2时有V ( X) 0，故V(x)为正半定。(4) 负半定性V(0) =0，且标量函数V(x)在域s内某些状态处有 V(x) =0，而在其它状态处有 V(x) :： 0且,则称V(x)在域s内负半定。女口 V(x -(x1 2x2)2是负半定的。(5) 不定性标量函数V(x)在域s内可正可负，则称 V(x)不定。如VajnXjX?是不定的。2、标量函数V(x)取二次型时的符号V(x) =xTPx - lx1 x2-PiiPi

16、2Pin 1_Xi 1Xn】p21P22P2naX2-_PniPn2Pnn 一X 一实二次型V(x)是正定的充要条件是矩阵式中P为对称矩阵，有Pij二Pji 。显然满足V(0) = 0。当P阵的每一个元都为实数时， 称作实二次型。P的各顺序主子行列式均大于零(赛尔维斯特准则)，即PiiPi2>0，PiiPi nP2iP22PniPnnPii0，则V(x)正定，且称P为正定矩阵。当矩阵P的各顺序主子行列式负、正相间时，即PiiPi2>0，, (-i)nPiiaPinaP2iP22PniPnnPii ：°，则V(x)负定，且称P为负定矩阵。若矩阵P的各顺序主子行列式含有等于零

17、的情况,则V( x)为正半定或负半定。 不属于以上所有情况的，V(x)为不定。【例6.3.1】证明下列的二次型是正定的2 2 2V(x) = 10捲 4x2 x3 2xx2 - 2x2x3 - 4XX3证明：上式用矩阵形式表示为101-2切V(x)=xTPx =X2X3】14-1 X22- 11LX3_j由于P11=10 0P11p12101p21p2214P11p12p13p21p22p23=p31p32p3310=39 0-2-1所以，V(x)是正定的。3、李雅普诺夫第二法定理1设系统的状态方程为X = f(x)f (0) =0(2) V(x)是所以该系统如果存在一个标量函数 V(x)，它

18、有连续的一阶偏导数，而且满足：(1)V(x)是正定的;负定的，则系统在原点处的平衡状态是渐近稳定的。【例6.3.2】设系统的状态方程是X1 = X2 - X1 (x1xf)x2 二-xx2(x.|2 X；)试分析系统在原点处的平衡状态是否为渐近稳定的。解：令x1 = 0及x2 = 0 ,解得x1 =0 , x2 = 0，故原点为平衡状态，且只有一个平衡状态。2 2 2设V(x)x；，则 V(x) =2X1/ 2X2X2 二-2(X1X2)显然，对于x=0存在V(x) ：0以及V(0)=0,故V(x)是负定的。又由于 V(x)是正定的, 在原点处的平衡状态是渐进稳定的。因为只有一个平衡状态，故该

19、系统是渐进稳定的。【例6.3.4】设系统的状态方程是xj01丁捲1丨=丨I-X2 -T 1X2-试分析系统在原点处的平衡状态是否为渐近稳定的。 解：这个系统的平衡点是 X| = X2 = 0。(1) 先猜想一个能代表系统能量的正定函数V(x) =2x2 x2则:V(x)工4洛為 2x2x2由原状态方程知： x, = x2， x2- x - x2，代入上式得V (x)二 2x,x 2x； 它是不定的，因而无法判断系统的稳定性。(2) 重新构造正定函数 V( X):2 2V(x)=xi X2 2因此：V (x 2x1x1 2x2x2 - -2x2它是负半定的。但是因为 V(x)不恒为零，所以可知该

20、系统是渐近稳定的。V(x)是李亚普诺夫函数。(3) 再取函数V(x)二1 (x, x2)2 2x； x； 12因此：V(x) = -(x,2 - x；)是负定的。且：|x|t 旳，V(X)t 血，故是大范围渐近稳定。从上述例子可以看出，找到一个系统的李亚普诺夫函数证明系统是稳定时，则系统必定是稳定的；反 之，如果用李亚普诺夫第二法不能判断系统是稳定时，并不能得出系统是不稳定的。4.3 线性定常系统的稳定性分析线性定常系统可以应用各种方法，诸如劳斯一霍尔维茨，奈奎斯特法来判断系统的稳定性。李亚普诺夫直接法也提供了一种稳定性判据的方法。在上一节的例4-3中，我们用猜想和试探、验证的方法来找李亚普诺

21、夫函数，这是很费时的，对于高阶的系统尤其费时，有时甚至是困难的。在本节中， 采用二次型函数来判定系统稳定性，这种方法较为简捷，并且还可以利用这种方法作为基础来解参数 最优问题以及某些控制系统的设计问题。设线性定常系统为式中x是n维状态向量。如果采用二次型函数作为李亚普诺夫函数，即（4-19）（4-20）（4-21）（4-22）（4-23）V x 二xTPx式中P是正定的实对称矩阵。那么对上式求导得V x =xTPx xTPx将式（4-18 ）代入（4-20 ）式得V x = xTPAx xTATPx 二 XT PA AtP x令_Q = PA AtP代入式（4-21 ）得V x = -xTQx从上式可以看出，为了判定V（X）是否负定，只需要判定 Q是否正定。因此，判定线性定常系统的稳定性可如下进行：选定一个正定实对称矩阵P,按式（4-