Ch.5李雅普诺夫稳定性分析目录(1/1)目录概述5.1李雅普诺夫稳定性的定义5.2李雅普诺夫稳定性的基本定理5.3线性系统的稳定性分析5.4非线性系统的稳定性分析5.5Matlab问题本章小结李雅普诺夫稳定性的基本定理(1/2)5.2李雅普诺夫稳定性的基本定理本节主要研究李雅普诺夫意义下各种稳定性的判定定理和判定方法。讨论的主要问题有:基本概念:矩阵和函数的定号性(正定性、负定性等)基本方法:非线性系统线性化方法李雅普诺夫第一法矩阵符号(正定性、负定性等)检验方法李雅普诺夫第二法难点喔!李雅普诺夫稳定性的基本定理(2/2)下面先讲述李雅普诺夫第一法,然后讨论李雅普诺夫第二法李雅普诺夫第一法(1/7)5.2.1李雅普诺夫第一法李雅普诺夫第一法又称间接法,它是研究动态系统的一次近似数学模型(线性化模型)稳定性的方法。它的基本思路是:首先,对于非线性系统,可先将非线性状态方程在平衡态附近进行线性化,即在平衡态求其一次Taylor展开式,然后利用这一次展开式表示的线性化方程去分析系统稳定性。其次,解出线性化状态方程组或线性状态方程组的特征值,然后根据全部特征值在复平面上的分布情况来判定系统在零输入情况下的稳定性。李雅普诺夫第一法(2/7)下面将讨论李雅普诺夫第一法的结论以及在判定系统的状态稳定性中的应用。设所讨论的非线性动态系统的状态方程为x’=f(x)其中f(x)为与状态向量x同维的关于x的非线性向量函数,其各元素对x有连续的偏导数。李雅普诺夫第一法(3/7)欲讨论系统在平衡态xe的稳定性,先必须将非线性向量函数f(x)在平衡态附近展开成Taylor级数,即有其中A为nn维的向量函数f(x)与x间的雅可比矩阵;R(x-xe)为Taylor展开式中包含x-xe的二次及二次以上的余项。雅可比矩阵A定义为eeeτee()()(-)(-)(-)(-)RAReexxxxfxxfxxxxxxxxxx111τ1/.../()........./.../nnnnfxfxAfxfxeexxxxfxx李雅普诺夫第一法(4/7)上述线性化方程的右边第一项A(x-xe)代表原非线性状态方程的一次近似式,如果用该一次近似式来表达原非线性方程的近似动态方程,即可得如下线性化的状态方程:x’=A(x-xe)由于对如上式所示的状态方程总可以通过n维状态空间中的坐标平移,将平衡态xe移到原点。因此,上式又可转换成如下原点平衡态的线性状态方程:x’=Ax判别非线性系统平衡态xe稳定性的李雅普诺夫第一法的思想即为:通过线性化,将讨论非线性系统平衡态稳定性问题转换到讨论线性系统x’=Ax的稳定性问题。李雅普诺夫第一法(5/7)李雅普诺夫第一法的基本结论是:1.若线性化系统的状态方程的系统矩阵A的所有特征值都具有负实部,则原非线性系统的平衡态xe渐近稳定,而且系统的稳定性与高阶项R(x)无关。2.若线性化系统的系统矩阵A的特征值中至少有一个具有正实部,则原非线性系统的平衡态xe不稳定,而且该平衡态的稳定性与高阶项R(x)无关。3.若线性化系统的系统矩阵A除有实部为零的特征值外,其余特征值都具有负实部,则原非线性系统的平衡态xe的稳定性由高阶项R(x)决定。李雅普诺夫第一法(6/7)由上述李雅普诺夫第一法的结论可知,该方法与经典控制理论中稳定性判据的思路一致,需求解线性化状态方程或线性状态方程的特征值,根据特征值在复平面的分布来分析稳定性。值得指出的区别是:经典控制理论讨论的是输出稳定性问题,而李雅普诺夫方法讨论状态稳定性问题。由于李雅普诺夫第一法需要求解线性化后系统的特征值,因此该方法也仅能适用于非线性定常系统或线性定常系统,而不能推广至时变系统。李雅普诺夫第一法(7/7)—例5-1试确定系统在原点处的稳定性。解1:由状态方程知,原点为该系统的平衡态。将系统在原点处线性化,则系统矩阵为1210)(KKAexxxxf0,)1(21122211221KKxKxxKxxx例5-1某装置的动力学特性用下列常微分方程组来描述:因此,系统的特征方程为|I-A|=2+K1+K2=0李雅普诺夫第一法(8/7)2.由李雅普诺夫第一法知,原非线性系统的原点为渐近稳定的充分条件为:K10和K20.李雅普诺夫第二法(1/3)5.2.2李雅普诺夫第二法由李雅普诺夫第一法的结论可知,该方法能解决部分弱非线性系统的稳定性判定问题,但对强非线性系统的稳定性判定则无能为力,而且该方法不易推广到时变系统。下面我们讨论对所有动态系统的状态方程的稳定性分析都适用的李雅普诺夫第二法。李雅普诺夫第二法(2/3)李雅普诺夫第二法又称为直接法。它是在用能量观点分析稳定性的基础上建立起来的。若系统平衡态渐近稳定,则系统经激励后,其储存的能量将随着时间推移而衰减。当趋于平衡态时,其能量达到最小值。反之,若平衡态不稳定,则系统将不断地从外界吸收能量,其储存的能量将越来越大。基于这样的观点,只要能找出一个能合理描述动态系统的n维状态的某种形式的能量正性函数,通过考察该函数随时间推移是否衰减,就可判断系统平衡态的稳定性。李雅普诺夫第二法(3/3)在给出李雅普诺夫稳定性定理之前,下面先介绍一些数学预备知识,然后介绍一些李雅普诺夫稳定性定理的直观意义,最后介绍李雅普诺夫稳定性定理数学预备知识(1/1)1.数学预备知识下面介绍在李雅普诺夫稳定性分析中需应用到的如下数学预备知识:函数的正定性二次型函数和对称矩阵的正定性矩阵正定性的判别方法实函数的正定性(1/4)—函数定号性定义(1)实函数的正定性实函数正定性问题亦称为函数定号性问题。它主要讨论该函数的值在什么条件下恒为正,什么条件下恒为负的。下面先给出n维向量x的标量实函数V(x)的正定性定义。定义5-5设xRn,是Rn中包含原点的一个区域,若实函数V(x)对任意n维非零向量x都有V(x)0;当且仅当x=0时,才有V(x)=0,则称函数V(x)为区域上的正定函数。实函数的正定性(2/4)—函数定号性定义从定义可知,所谓正定函数,即指除零点外恒为正值的标量函数。由正定函数的定义,我们相应地可定义负定函数、非负定(又称半正定或正半定)函数、非正定函数(又称半负定或负半定)和不定函数。实函数的正定性(3/4)—函数定号性定义定义5-6设xRn,是Rn中包含原点的一个区域,若实函数V(x)对任意n维非零向量x,都有V(x)0;当且仅当x=0时,才有V(x)=0,则称函数V(x)为区域上的负定函数。若对任意n维非零向量x,都有V(x)≥0,且V(0)=0,则称函数V(x)为区域上的非负定函数。若对任意n维非零向量x,都有V(x)≤0,且V(0)=0,则称函数V(x)为区域上的非正定函数。若无论取多么小的原点的某个邻域,V(x)可为正值也可为负值,则称函数V(x)为不定函数。实函数的正定性(4/4)下面是几个在由变量x1和x2组成的2维线性空间中的正定函数、负定函数等的例子。1)正定函数222212221)2(2xxxxx2)负定函数2122122215)2(2xxxxx3)非负定函数22122)2(2xxx4)非正定函数22121)2(3xxx实函数的正定性(5/4)5)不定函数函数的定号性是一个相对概念,与其函数定义域有关。如,函数对x1与x2组成的2维空间为非负定的,但对于1维空间x2则为正定的。上面定义了时不变函数V(x)的定号性,相应地可以定义标量时变函数V(x,t)的定号性。2212212221)2()2(23xxxxxx222x实函数的正定性(6/4)定义5-7设xRn,是Rn中包含原点的一个封闭有限区域,实函数V(x,t)是定义在[t0,)上的一个标量函数且V(0,t)=0,标量连续函数(||x||)和(||x||)为非减(函数值单调增加)的且满足(0)=(0)=0,1)如果对任意tt0和x0,V(x,t)为有界正定的,即0(||x||)V(x,t)(||x||),称函数V(x,t)为[t0,)上的(时变)正定函数。2)如果对任意tt0和x0,分别为有界负定,即0-(||x||)V(x,t)-(||x||);有界非负定,即0V(x,t)(||x||);有界非正定,即0V(x,t)-(||x||),实函数的正定性(6/4)分别称函数V(x,t)为[t0,)上的(时变)负定函数、非负定函数和非正定函数。3)如果存在tt0,无论取多么小的原点的某个邻域,V(x,t)可为正值也可为负值,则称函数V(x,t)为不定函数。二次型函数和对称矩阵的正定性(1/4)(2)二次型函数和对称矩阵的正定性二次型函数是一类特殊形式函数。设V(x)为关于n维变量向量x的实二次型函数,则其可以表示为其中aij(i=1,2,…,n,j=i,…,n)为实常数。211112121122222221().........nnnnnnnnnijijijiVaxaxxaxxaxaxxaxaxxx二次型函数和对称矩阵的正定性(2/4)由线性代数知识知,实二次型函数V(x)又可表示为V(x)=xPx其中P称为二次型函数V(x)的权矩阵,它为如下nn维实对称矩阵:nnnnnnaaaaaaaaaP...2/2/............2/...2/2/...2/212221211211二次型函数和对称矩阵的正定性(3/4)二次型函数与一般函数一样,具有正定、负定、非负定、非正定和不定等定号性概念。二次型函数V(x)和它的对称权矩阵P是一一对应的。因此,由二次型函数的正定性同样可定义对称矩阵P的正定性。定义5-8设对称矩阵P为二次型函数V(x)的权矩阵,当V(x)分别为正定、负定、非负定、非正定与不定时,则称对称矩阵P相应为正定、负定、非负定、非正定与不定。□二次型函数和对称矩阵的正定性(4/4)--矩阵定号性定义因此,由上述定义就可将判别二次型函数的正定性转换成为判别对称矩阵的正定性。对称矩阵P为正定、负定、非负定与非正定时,并可分别记为P0,P0,P≥0,P≤0。矩阵正定性的判别方法(1/5)(3)矩阵正定性的判别方法判别矩阵的正定性(定号性)的方法主要有塞尔维斯特判别法、矩阵特征值判别法和合同变换法。下面分别介绍。矩阵正定性的判别方法(2/5)--塞尔维斯特定理定理5-1(塞尔维斯特定理)(1)实对称矩阵P为正定的充要条件是P的各阶顺序主子式均大于零,即0||...00222112112111PΔppppΔpΔn其中pij为实对称矩阵P的第i行第j列元素。(2)实对称矩阵P为负定的充要条件是P的各阶顺序主子式满足niiii,...,2,100为奇数为偶数矩阵正定性的判别方法(2/5)—矩阵定号性判定定理定理5-2实对称矩阵P为正定、负定、非负定与非正定的充分必要条件是P的所有特征值分别大于零、小于零、大于等于零与小于等于零;实对称矩阵P为不定的充分必要条件是P的特征值有正有负。□定理5-3实对称矩阵P必定可经合同变换化成对角线矩阵,则P为正定、负定、非负定与非正定的充分必要条件是的所有对角线元素分别大于零、小于零、大于等于零与小于等于零;P为不定的充分必要条件是的对角线元素有正有负。矩阵正定性的判别方法(3/5)—矩阵定号性判定定理定理5-3中的合同变换是指对对称矩阵的同样序号的行和列同时作同样的初等变换。上述三种判别实对称矩阵P的定号性的方法,各有千秋。但总的说来,基于塞尔维斯特定理的方法计算量较大,若将该方法推广到判别非负定性和非正定性,则计算量成指数性地