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

Ch.5Ch.5Ch.5Ch.5李雅普诺夫稳定性分析目录(1/1)(1/1)(1/1)(1/1)目录�概述�5.15.15.15.1李雅普诺夫稳定性的定义�5.25.25.25.2李雅普诺夫稳定性的基本定理�5.35.35.35.3线性系统的稳定性分析�5.45.45.45.4非线性系统的稳定性分析�5.55.55.55.5MatlabMatlabMatlabMatlab问题�本章小结李雅普诺夫稳定性的基本定理(1/2)(1/2)(1/2)(1/2)5.2李雅普诺夫稳定性的基本定理�本节主要研究李雅普诺夫意义下各种稳定性的判定定理和判定方法。讨论的主要问题有:基本概念::::矩阵和函数的定号性(正定性、负定性等)基本方法::::非线性系统线性化方法李雅普诺夫第一法矩阵符号(正定性、负定性等)检验方法李雅普诺夫第二法难点喔!李雅普诺夫稳定性的基本定理(2/2)(2/2)(2/2)(2/2)�下面先讲述�李雅普诺夫第一法,然后讨论�李雅普诺夫第二法李雅普诺夫第一法(1/7)(1/7)(1/7)(1/7)5.2.1李雅普诺夫第一法�李雅普诺夫第一法又称间接法,它是研究动态系统的一次近似数学模型(线性化模型)稳定性的方法。它的基本思路是:�首先,对于非线性系统,可先将非线性状态方程在平衡态附近进行线性化,�即在平衡态求其一次Taylor展开式,�然后利用这一次展开式表示的线性化方程去分析系统稳定性。�其次,解出线性化状态方程组或线性状态方程组的特征值,然后根据全部特征值在复平面上的分布情况来判定系统在零输入情况下的稳定性。李雅普诺夫第一法(2/7)(2/7)(2/7)(2/7)�下面将讨论李雅普诺夫第一法的结论以及在判定系统的状态稳定性中的应用。�设所讨论的非线性动态系统的状态方程为xxxx’=ffff(xxxx)其中ffff(xxxx)为与状态向量xxxx同维的关于xxxx的非线性向量函数,其各元素对xxxx有连续的偏导数。李雅普诺夫第一法((((3/7)3/7)3/7)3/7)�欲讨论系统在平衡态xxxxe的稳定性,先必须将非线性向量函数ffff(xxxx)在平衡态附近展开成Taylor级数,即有其中A为n×n维的向量函数ffff(xxxx)与xxxx间的雅可比矩阵;R(xxxx-xxxxe)为Taylor展开式中包含xxxx-xxxxe的二次及二次以上的余项。雅可比矩阵A定义为eeeτee()()(-)(-)(-)(-)RAR==∂′=++∂=+eeeeeeeexxxxxxxxxxxxxxxxfxfxfxfxxfxxxxxxfxxxxxxfxxxxxxfxxxxxxxxxxxxxxxxxxxxxxxxx111τ1/.../()........./.../nnnnfxfxAfxfx==∂∂∂∂⎡⎤∂⎢⎥==⎢⎥∂⎢⎥∂∂∂∂⎣⎦eeeeeeeexxxxxxxxxxxxxxxxfxfxfxfxxxxx李雅普诺夫第一法((((4/7)4/7)4/7)4/7)�上述线性化方程的右边第一项A(xxxx-xxxxe)代表原非线性状态方程的一次近似式,如果用该一次近似式来表达原非线性方程的近似动态方程,即可得如下线性化的状态方程:xxxx’=A(xxxx-xxxxe)�由于对如上式所示的状态方程总可以通过n维状态空间中的坐标平移,将平衡态xxxxe移到原点。�因此,上式又可转换成如下原点平衡态的线性状态方程:xxxx’=Axxxx�判别非线性系统平衡态xxxxe稳定性的李雅普诺夫第一法的思想即为:�通过线性化,将讨论非线性系统平衡态稳定性问题转换到讨论线性系统xxxx’=Axxxx的稳定性问题。李雅普诺夫第一法((((5/7)5/7)5/7)5/7)�李雅普诺夫第一法的基本结论是:1.1.1.1.若线性化系统的状态方程的系统矩阵A的所有特征值都具有负实部,则原非线性系统的平衡态xxxxe渐近稳定,而且系统的稳定性与高阶项R(xxxx)无关。2.2.2.2.若线性化系统的系统矩阵A的特征值中至少有一个具有正实部,则原非线性系统的平衡态xxxxe不稳定,而且该平衡态的稳定性与高阶项R(xxxx)无关。3.3.3.3.若线性化系统的系统矩阵A除有实部为零的特征值外,其余特征值都具有负实部,则原非线性系统的平衡态xxxxe的稳定性由高阶项R(xxxx)决定。李雅普诺夫第一法((((6/7)6/7)6/7)6/7)�由上述李雅普诺夫第一法的结论可知,该方法与经典控制理论中稳定性判据的思路一致,需求解线性化状态方程或线性状态方程的特征值,根据特征值在复平面的分布来分析稳定性。�值得指出的区别是:�经典控制理论讨论的是输出稳定性问题,而李雅普诺夫方法讨论状态稳定性问题。�由于李雅普诺夫第一法需要求解线性化后系统的特征值,因此该方法也仅能适用于非线性定常系统或线性定常系统,而不能推广至时变系统。李雅普诺夫第一法((((7/7)7/7)7/7)7/7)————例5-15-15-15-1试确定系统在原点处的稳定性。�解1111：由状态方程知,原点为该系统的平衡态。�将系统在原点处线性化,则系统矩阵为⎥⎦⎤⎢⎣⎡−−=∂∂==1210)(KKAeeeexxxxxxxxxxxxxxxxffffτ0,)1(21122211221⎥⎦⎤⎢⎣⎡−−=⎥⎦⎤⎢⎣⎡′′KKxKxxKxxx�例5-15-15-15-1某装置的动力学特性用下列常微分方程组来描述:因此,系统的特征方程为|λI-A|=λ2+K1λ+K2=0李雅普诺夫第一法((((8/7)8/7)8/7)8/7)2.2.2.2.由李雅普诺夫第一法知,原非线性系统的原点为渐近稳定的充分条件为:K10和K20.李雅普诺夫第二法(1/3)(1/3)(1/3)(1/3)5.2.2李雅普诺夫第二法�由李雅普诺夫第一法的结论可知,该方法能解决部分弱非线性系统的稳定性判定问题,但对强非线性系统的稳定性判定则无能为力,而且该方法不易推广到时变系统。�下面我们讨论对所有动态系统的状态方程的稳定性分析都适用的李雅普诺夫第二法。李雅普诺夫第二法(2/3)(2/3)(2/3)(2/3)�李雅普诺夫第二法又称为直接法。�它是在用能量观点分析稳定性的基础上建立起来的。�若系统平衡态渐近稳定,则系统经激励后,其储存的能量将随着时间推移而衰减。当趋于平衡态时,其能量达到最小值。�反之,若平衡态不稳定,则系统将不断地从外界吸收能量,其储存的能量将越来越大。�基于这样的观点,只要能找出一个能合理描述动态系统的n维状态的某种形式的能量正性函数,通过考察该函数随时间推移是否衰减,就可判断系统平衡态的稳定性。李雅普诺夫第二法(3/3)(3/3)(3/3)(3/3)�在给出李雅普诺夫稳定性定理之前,下面先介绍一些�数学预备知识,然后介绍一些�李雅普诺夫稳定性定理的直观意义,最后介绍�李雅普诺夫稳定性定理数学预备知识(1/1)(1/1)(1/1)(1/1)1.1.1.1.数学预备知识�下面介绍在李雅普诺夫稳定性分析中需应用到的如下数学预备知识:�函数的正定性�二次型函数和对称矩阵的正定性�矩阵正定性的判别方法实函数的正定性(1/4)(1/4)(1/4)(1/4)————函数定号性定义((((1111))))实函数的正定性�实函数正定性问题亦称为函数定号性问题。�它主要讨论该函数的值在什么条件下恒为正,什么条件下恒为负的。�下面先给出n维向量xxxx的标量实函数V(xxxx)的正定性定义。�定义5-55-55-55-5设xxxx∈RRRRn,Ω是RRRRn中包含原点的一个区域,若实函数V(xxxx)对任意n维非零向量xxxx∈Ω都有V(xxxx)0;当且仅当xxxx=0时,才有V(xxxx)=0,�则称函数V(xxxx)为区域Ω上的正定函数。∆∆∆实函数的正定性(2/4)(2/4)(2/4)(2/4)————函数定号性定义�从定义可知,所谓正定函数,即指除零点外恒为正值的标量函数。由正定函数的定义,我们相应地可定义�负定函数、�非负定(又称半正定或正半定)函数、�非正定函数(又称半负定或负半定)和�不定函数。实函数的正定性(3/4)(3/4)(3/4)(3/4)————函数定号性定义�定义5-65-65-65-6设xxxx∈RRRRn,Ω是RRRRn中包含原点的一个区域,若实函数V(xxxx)对任意n维非零向量xxxx∈Ω,都有V(xxxx)0;当且仅当xxxx=0时,才有V(xxxx)=0,则称函数V(xxxx)为区域Ω上的负定函数。�若对任意n维非零向量xxxx∈Ω,都有V(xxxx)≥0,且V(0)=0,则称函数V(xxxx)为区域Ω上的非负定函数。�若对任意n维非零向量xxxx∈Ω,都有V(xxxx)≤0,且V(0)=0,则称函数V(xxxx)为区域Ω上的非正定函数。�若无论取多么小的原点的某个邻域,V(xxxx)可为正值也可为负值,则称函数V(xxxx)为不定函数。∆∆∆实函数的正定性(4/4)(4/4)(4/4)(4/4)�下面是几个在由变量x1和x2组成的2维线性空间中的正定函数、负定函数等的例子。1)1)1)1)正定函数222212221)2(2xxxxx+−+2)2)2)2)负定函数2122122215)2(2xxxxx−+−−−3)3)3)3)非负定函数22122)2(2xxx−4)4)4)4)非正定函数22121)2(3xxx+−−实函数的正定性(5/4)(5/4)(5/4)(5/4)5)5)5)5)不定函数�函数的定号性是一个相对概念,与其函数定义域有关。�如,函数对x1与x2组成的2维空间为非负定的,但对于1维空间x2则为正定的。�上面定义了时不变函数V(xxxx)的定号性,相应地可以定义标量时变函数V(xxxx,t)的定号性。2212212221)2()2(23xxxxxx+−−+−222x实函数的正定性((((6/4)6/4)6/4)6/4)�定义5-75-75-75-7设xxxx∈RRRRn,Ω是RRRRn中包含原点的一个封闭有限区域,实函数V(xxxx,t)是定义在[t0,∞)×Ω上的一个标量函数且V(0000,t)=0,标量连续函数α(||xxxx||)和β(||xxxx||)为非减(函数值单调增加)的且满足α(0)=β(0)=0,1)1)1)1)如果对任意t≥t0和xxxx≠0,V(xxxx,t)为有界正定的,即0α(||xxxx||)≤V(xxxx,t)≤β(||xxxx||)，称函数V(xxxx,t)为[t0,∞)×Ω上的(时变)正定函数。2)2)2)2)如果对任意t≥t0和xxxx≠0,分别为有界负定，即0-α(||xxxx||)≥V(xxxx,t)≥-β(||xxxx||)；有界非负定，即0≤V(xxxx,t)≤β(||xxxx||)；有界非正定，即0≥V(xxxx,t)≥-β(||xxxx||)，实函数的正定性((((6/4)6/4)6/4)6/4)分别称函数V(xxxx,t)为[t0,∞)×Ω上的(时变)负定函数、非负定函数和非正定函数。3)3)3)3)如果存在t≥t0,无论取多么小的原点的某个邻域,V(xxxx,t)可为正值也可为负值,则称函数V(xxxx,t)为不定函数。二次型函数和对称矩阵的正定性(1/4)(1/4)(1/4)(1/4)(2)(2)(2)(2)二次型函数和对称矩阵的正定性�二次型函数是一类特殊形式函数。�设V(xxxx)为关于n维变量向量xxxx的实二次型函数,则其可以表示为其中aij(i=1,2,…,n,j=i,…,n)为实常数。211112121122222221().........nnnnnnnnnijijijiVaxaxxaxxaxaxxaxaxx=≥=++++++++=∑∑xxxx二次型函数和对称矩阵的正定性(2/4)(2/4)(2/4)(2/4)�由线性代数知识知,实二次型函数V(x)又可表示为V(xxxx)=xxxxτPxxxx其中P称为二次型函数V(xxxx)的权矩阵,它为如下n×n维实对称矩阵:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=nnnnnnaaaaa