李雅普诺夫稳定性的基本定理

1、李雅普诺夫第一法(3/7)3.2.1李雅普诺夫第一法李雅普诺夫第一法又称间接法，它是研究动态系统的一次近似数学模型(线性化模型)稳定性的方法。它的基本思路是：首先，对于非线性系统，可先将非线性状态方程在平衡态 附近进行线性化，“即在平衡态求其一次Taylor展开式，/然后利用这一次展开式表示的线性化方程去分析系 统稳是程。其次，解出线性化状态方程组或线性状态方程组的特征值, 然后根据全部特征值在复平面上的分布情况来判定系统 在零输入情况下的稳定性。下面将讨论李雅普诺夫第一法的结论以及在判定系统的状态稳 定性中的应用。设所讨论的非线性动态系统的状态方程为 xr=f(x) 其中/(兀)为与状态向量

2、V同维的关于兀的非线性向量函数，其各元 素对兀有连续的偏导数。<>李雅普诺夫第一法(5/7)李雅普诺夫第一法的基本结论是:1. 若线性化系统的状态方程的系统矩阵A的所有特征值都 具有负实部，则原非线性系统的平衡态耳渐近稳定，而且系 统的稳定性与高阶项皿X)无关。2. 若线性化系统的系统矩阵A的特征值中至少有一个具有 正实部，则原非线性系统的平衡态冬不稳定，而且该平衡态 的稳定性与高阶项皿X)无关。3. 若线性化系统的系统矩阵A除有实部为零的特征值外，其 余特征值都具有负实部，则原非线性系统的平衡态眾的稳 定性由高阶项R(x)决定。由上述李雅普诺夫第_法的结论可知，该方法与经典控制理

3、论中稳定性判据的思路一致，需求解线性化状态方程或线性状态方程的特征值，根据特征值在复平面的分布来分析稳定性。值得指出的区别是:“经典控制理论讨论的是输出稳定性问题，而李雅普诺 夫方法讨论状态稳定性问题。由于李雅普诺夫第一法需要求解线性化后系统的特征值，因此该方法也仅能适用于非线性定常系统或线性定常系 统，而不能推广至时变系统。<>李雅普诺夫第一法(9/7)例3J某装置的动力学特性用下列常微分方程组来描述:试确定系统在原点处的稳定性。解1:由状态方程知，原点为该系统的平衡态。A将系统在原点处线性化,则系统矩阵为af(x) 01 _A =dxTX = XeL-2-kJ因此，系统的特征方

4、程为AI-A I=22+Ki2+K2=02.由李雅普诺夫第一法知，原非线性系统的原点为渐近稳定的充 分条件为：K>0 和 K2>Q.<>李雅普诺夫第二法(3/3)322李雅普诺夫第二法由李雅普诺夫第一法的结论可知，该方法能解决部分弱非线性 系统的稳定性判定问题，但对强非线性系统的稳定性判定则无 能为力，而且该方法不易推广到时变系统。下面我们讨论对所有动态系统的状态方程的稳定性分析 都适用的李雅普诺夫第二法。李雅普诺夫第二法又称为直接法。它是在用能量观点分析稳定性的基础上建立起来的。/若系统平衡态渐近稳定，则系统经激励后，其储存的能 量将随着时间推移而衰减。当趋于平衡态时

5、，其能量 达到最小值。/反之，若平衡态不稳定，则系统将不断地从外界吸收能量，其储存的能量将越来越大。基于这样的观点，只要能找出一个能合理描述动态系统的 维状态的某种形式的能量正性函数，通过考察该函数随 时间推移是否衰减，就可判断系统平衡态的稳定性。在给出李雅普诺夫稳定性定理之前，下面先介绍一些数学预备知识然后介绍一些李雅普诺夫稳定性定理的直观意义最后介绍 李雅普诺夫稳定性定理<>数学预备知识(1/1)1.数学预备知识下面介绍在李雅普诺夫稳定性分析中需应用到的如下数学预 备知识：函数的正定性二次型函数和对称矩阵的正定性矩阵正定性的判别方法<>实函数的正定性(3/4)函数定

6、号性定义(1) 实函数的正定性实函数正定性问题亦称为函数定号性问题。它主要讨论该函数的值在什么条件下恒为正，什么条件下 恒为负的。下面先给出维向量兀的标量实函数卩(兀)的正定性定义。定义35设灼卅q是卍中包含原点的一个区域，若实函数V(x) 对任意维非零向量r w Q都有Vd)o；当且仅当*0时，才有 V(x)=0,则称函数Vd)为区域。上的正定函数。从定义可知，所谓正定函数，即指除零点外恒为正值的标量函 数。由正定函数的定义，我们相应地可定义负定函数、非负定（又称半正定或正半定）函数、非正定函数（又称半负定或负半定）和-不定函数。定义36设兀是川中包含原点的一个区域，若实函数V(Q对任意维非

7、零向量都有V(QvO;当且仅当*0时才有V(x)=0,则称函数V(Q为区域Q上的负定函数。若对任意维非零向量vwQ，都有巾左0,且叫0)=0,则称函 数并兀)为区域G上的非负定函数。若对任意维非零向量兀wG，都有巾汇0,且V(0)=0,则称函 数只兀)为区域Q上的非正定函数。若无论取多么小的原点的某个邻域V(Q可为正值也可为 负值，则称函数V(Q为不定函数。A/A<>彳 + 2xf（%! - 2x9 ）2 + 兀；实函数的正定性（4/4） 下面是几个在由变量心和兀2组成的2维线性空间中的正定函数、 负定函数等的例子。1）正定函数2）负定函数22（Xj + 2 %2）_ 5兀3）非负

8、定函数（X - 2兀2）2<><>4）非正定函数亟5）不定函数3 Xj +2 %2（Xj 2兀2）（兀+ 2兀2）<>二次型函数和对称矩阵的正定性(3/4)(2) 二次型函数和对称矩阵的正定性二次型函数是一类特殊形式函数。示为设V(x)为关于维变量向量x的实二次型函数，则其可以表2V（X）=+ + a. x.x 丿1 111212In 1 n2+ + xx222In 2 n+2 + a x nn nn n=y1 V a.x 兀.厶厶 IJ I J心 1 J>i其中,nj=i.m)为实常数。由线性代数知识知，实二次型函数V(Q又可表示为V(x)=xTPx

9、其中P称为二次型函数V(x)的权矩阵，它为如下“X”维实对称矩阵:d 2 / 2 a. /2 In6Z12/2Q 22/2In an/2 r / 2 2n ann<>二次型函数和对称矩阵的正定性(5/4)二次型函数与一般函数一样，具有正定、负定、非负定.非正定和不定等定号性概念。二次型函数V(x)和它的对称权矩阵P是一一对应的。因此，由二次型函数的正定性同样可定义对称矩阵尸的正 定性。定义38设对称矩阵尸为二次型函数并兀)的权矩阵，当V(x)分 别为正定、负定、非负定、非正定与不定时，则称对称矩阵尸 相应为正定、负定、非负定、非正定与不定。<>二次型函数和对称矩阵的正定

10、性(#/4)矩阵定号性定义因此，由上述定义就可将判别二次型函数的正定性转换成为 判别对称矩阵的正定性。>对称矩阵尸为正定、负定、非负定与非正定时，并可分别 记为P>0, P<0, P>0, P<0o<A矩阵正定性的判别方法(2/5)“塞尔维斯特定理（3）矩阵正定性的判别方法判别矩阵的正定性（定号性）的方法主要有塞尔维斯特判别法、矩阵特征值判别法和合同变换法。 下面分别介绍。定理3-1(塞尔维斯特定理)实对称矩阵P为正定的充要条 件是尸的各阶顺序主子式均大于零，即L = pn>0 4 = Pi2 >0=1 P l> 0P21P22n其中p帀为

11、实对称矩阵p的第i行第/列元素。(2)实对称矩阵P为负定的充要条件是P的各阶顺序主子式满足:>0 i为偶数<0 i为奇数1,2,n矩阵正定性的判别方法(3/5)矩阵定号性判定定理定理32实对称矩阵尸为正定、负定、非负定与非正定的充 分必要条件是P的所有特征值分别大于零、小于零、大于等 于零与小于等于零；实对称矩阵P为不定的充分必要条件是P的特征值有正 有负。定理33实对称矩阵尸必定可经合同变换化成对角线矩阵，则P为正定、贪定、非负定与非正定葩充分必要棗件基的所肴 对角线元素分别大于零、小于零、大于等于零与小于等于尸为不定的充分必要条件是的对角线元素有正有负。定理33中的合同变换是指

12、对对称矩阵的同样序号的行和列 同时作同样的初等变换。上述三种判别实对称矩阵尸的定号性的方法，各有千秋。但总 的说来，基于塞尔维斯特定理的方法计算量较大，若将该方法推广 到判别非负定性和非正定性，则计算量成指数性地增加。特征值判别法需求解高阶特征方程以获得特征值，计算较 复杂，计算量也较大。合同变换法对矩阵只作初等变换，计算简单，便于应用。矩阵正定性的判别方法（#/5）例5-2口例3-2试用合同变换法判别下列实对称矩阵f的定号性：I 1-1-1P =3212 5解先对对称矩阵P作合同变换如下矩阵正定性的判别方法(5/5)例5-2VIA矩阵正定性的判别方法(#/5)例5-21-1-r行:(2)+(

13、1)t(2)_ 10-T-132=>列:(2)+(1)t(2)021-125-1150n行:(3)+ t(3)| 亠 I 列:(3)+(1)t(3)L°01n行:(3) -(2)/2t(3)列:(3) - (2)/2t(3)L°7/2因此，由定理33知，矩阵P为正定矩阵。VIA李雅普诺夫稳定性定理的直观意义(5/5)2.李雅普诺夫稳定性定理的直观意义从平衡态的定义可知，平衡态是使得系统静止不动(导数为零, 即运动变化的趋势为零)的状态。从能量的观点来说,静止不动即不存在运动变化所需要 的能量，即变化所需的能量为零。通过分析状态变化所反映的能量变化关系可以分析出 状态的

14、变迁或演变,可以分析出平衡态是否稳定或不稳 定。蔦嚴鑑雪能量变化来简介李雅普诺夫稳右图所示动力学系统的平衡态在一定范围内为渐近稳定的平衡态。对该平衡态的邻域，可定义其能量(动能+势能)函数如下:丄2V - mv + mgh21 卍=mxr + mg (xcos 0) > 02其中兀为位移#为速度，两者且选为状态变量。在图中所示状态，吟讥：由牛顿第二定律可知，其运动满足 如下方程：m(-x ,A) =mg cosQ -fmgsinQ其中伪摩擦阻尼系数。因此，有mx "=/wg(cos0 -/sinO)>因此，能量的变化趋势(导数)为V=mx'x"+mgx&

15、quot; cosd=-mgx7 (cosO-fsinQ) +mgxr cos0=/wgx7sin0当e取值为©90。，由于卩的方向与兀相反，疋为负，因此上式恒 小于零。/即渐近稳定的平衡态，其正定的能量函数的导数(变 化趋势)为负。对小球向上运动时亦可作同样分析。<>李雅普诺夫稳定性定理的直观意义(7/5)从直观物理意义的角度，也非常易于理解。由于物体运动所受到的摩擦力作负功，由能量守恒定律 可知，物体的能量将随物体运动减少，即其导数(变化趋势)为负。再如右图所示的动力学系统，其 平衡态在一定范围内为不稳定 的平衡态。对该平衡态的邻域，可定义其能量(动能+势能)函数如1

16、 2 1 2V = mvmv ° + mg (h - h()2 21 ,2 1 2=mxrmv o + mg (-x cos 0)2 2< 0<A李雅普诺夫稳定性定理的直观意义(9/5)由牛顿第二定律可知，其运动满足如 下方程：ma=/ngcos0 -fmgsinG因此，有nvc"=mg(cosQ -/sinO)因此，能量的变化趋势(导数)为Vf=mxxff-mgxcosQ=mgx(cosQ-fsirQymgxrcosQ=-mgxfsird当e取值为0,90。，由于为正，因此上式恒小于零。/即不稳定的平衡态，其负定的能量函数的导数(变化 趋势)为负。李雅普诺夫第

17、二法的基本思想就是通过定义和分析一个在平 衡态邻域的关于运动状态的广义能量函数来分析平衡态的稳 定性。通过考察该能量函数随时间变化是否衰减来判定平衡态 是渐近稳定,还是不稳定。基于上述关于函数的定号性的定义和上述物理意义解释, 下面阐述李雅普诺夫第二法关于/平衡态稳定、/渐近稳定、/大范围渐近稳定和/不稳定的几个定理。<1 >渐近稳定性定理(1/7)(1)渐近稳定性定理定理34设系统的状态方程为xz=f(x/) 其中冬=0为其平衡态。>若存在一个有连续一阶偏导数的正定函数V(x,0,满足下 述条件：1) 若F(x/)为负定的，则该系统在原点处的平衡态是一 致渐近稳定的；2)

18、更进一步，若随着11刘1->8,有V(x)too,那么该系统在原点处的平衡态是大范围一致渐近稳定的。对上述李雅普诺夫稳定性定理的使用有如下说明：1)此定理只为判别系统一致渐近稳定的充分条件,而非必 要条件。/也就是说，若找到满足上述条件的一个李雅普诺夫函数，则系统是一致渐近稳定或大范围一致渐近稳 定的。“瞬貓平廳薜聽貓诺夫函埶此时，我们或者0继续寻找满足条件的李雅普诺夫函数，或 者°蠶器養定理的结论来判别平衡态的<>渐近稳定性定理(5/7)2) 对于渐近稳定的平衡态，满足条件的李雅普诺夫函数总是 存在的，但并不唯一。3) 对于非线性系统，虽然具体的李雅普诺夫函数可证

19、明所讨 论的系统在平衡态的邻域内是渐近稳定的，但并不意味 着在其他的区域系统是或不是渐近稳定的；/对于线性系统，如果存在着渐近稳定的平衡态，则它 必是大范围渐近稳定的。4) 此定理不仅适用于线性系统，同样适用于非线性系统;既 适用于定常系统，同样也适用于时变系统。/因此李雅普诺夫第二法是判别平衡态稳定性的具有普遍性的方法。5) 李雅普诺夫第二法的结论并没有指明寻找李雅普诺夫函 数的方法。/寻找李雅普诺夫函数的方法将依具体的系统和状态 方程而具体分析。对于二阶系统，容易给出上述定理的直观几何解释(右图为李雅普诺夫函数V(x/)为欧氏距离的一个二维系统的帀立2相平面图)。李雅普诺夫函数V(x,f)

20、相当于定义 为表征系统的某种广义能量的一 神正定鹵薮。令为不同的常数，则相当于 在维状态空间上定义了一簇以 原点为中心,形状相似的同心超球 面。V随f增大，V随f减少, 则平衡态 则平衡态不稳泄 渐近稳定导函数wg)表征系统的广义能量函数的变化速率。<>渐近稳定性定理（#/7）V随f增大，V随f减少, 则平衡态 则平衡态 不稳定 渐近稳定卩（对）为负定同时也表示系统状态 将从现在所处于的在该封闭超球面 簇中超球面向原点方向（向内）运动， 最后逐渐趋向原点。渐近稳定性定理(7/7)“例3<>渐近稳定性定理(#/7)“例3例3-3试确定用如下状态方程描述的系统的平衡态稳定性

21、。解显然，原点(0,0)是给定系统的唯一平衡态，如果我们选择正 定函数为李雅普诺夫函数,那么沿任意轨迹x(/),V(x)对时间的全导数V r(x) = 2兀X： + 2xx - 一2(兀:+ 兀；)2 < 0是负定函数。此外，当11x11>8时，必有V(X)>00o因此，由定理34知，在原点处的平衡态是大范围一致渐近稳定的。<>例3-4试确定用如下状态方程描述的系统的平衡态稳定性。解显然，原点(0,0)是给定系统的唯一平衡态，如果我们选择正 定函数V (x) = %.2 + 兀；为李雅普诺夫函数,那么沿任意轨迹x(/),V(x)对时间的全导数V x) = 2xjXj += -2%2 § 0是负定函数，故由定理34知，根据所选的李雅普诺夫函数分析不 出该平衡态是否渐近稳定或稳定。但这也并不意味着该平衡态就并不渐近稳定。<渐近稳定性定理(9/7)定理34中严格要求选择的李雅普诺夫函数为正定函数，其导数为负定函数。这给该定理的应用，特别是寻找适宜的李雅普诺夫函数 带来一定困难。下面给出一个定理对上述定理34作一补充，以减弱判别 条件。<>（2）稳定性定理定理35设系统的状态方程为其中冬=0为其平衡态。 若存在一个有连续一阶偏导数的正定函数V（xQ满足下述条 件：1）卩（X,°为非正定（半负定）的，则