1、- 79 -5 控制系统的李雅普诺夫稳定性分析内容提要稳定性是系统的又一重要特性。所谓系统的稳定性,就是系统在受到小的外界扰动后,被调量与规定量之间的偏差值的过渡过程的收敛性。显然,稳定性是系统的一个动态属性。在控制理论和控制工程中,无论是调节器理论、观测器理论还是滤波预测、自适应理论,都不可避免的要遇到系统稳定性问题。稳定性问题一直是一个最基本的和最重要的问题。随着控制理论与控制工程所涉及的领域由线性时不变系统扩展为时变系统和非线性系统,稳定性分析的复杂程度也在急剧的增长。直到目前,虽然有许多判据可应用于线性时不变系统或其它各自相应类型的问题中,以判断系统稳定情况,但能同时有效地适用于线性、
2、非线性、定常、时变等各类系统的方法,则是俄国数学家李雅普诺夫(Lyaponov) 在 19世纪所提出的方法。这就是控制系统稳定性分析的李雅普诺夫方法。李雅普诺夫稳定性理论是稳定性分析、应用与研究的最重要基础。习题与解答5.1 判断下列函数的正定性1) 221313()Vxxx2) 22288x3) 1313()xx4) 221231044Vx5) 123()xx解 1) , 因为顺序主子式 10()3TVxAxx210,50,32130所以 , 为正定函数。A()Vx- 80 -2) , 因为主子式 841()2TVxAxx8218,210,0,70,0,416261所以 不定, 为不定函数。
3、A()Vx3) , 因为顺序主子式 120()TAx120,0,1204所以 为不定矩阵, 为不定函数。A()Vx4) , 因为顺序主子式 102()4Tx10,39,41041029所以 , 为正定函数。0A()Vx5) , 因为顺序主子式 1()32TAx10,0,31321402所以 , 为正定函数。 A()Vx- 81 -5.2 用李雅普诺夫第一方法判定下列系统在平衡状态的稳定性。 21212()xx解 解方程组 得三个孤立平衡点(0,0) , (1,1)和2121()xx(1,1) 。在(0,0)处将系统近似线性化,得 ,由于原系统为定常系统,且矩阵12xx的特征根 均具有负实部,于
4、是根据李雅普诺夫定理可知系统在原1212i点(0,0)附近一致渐近稳定。在(1,1)和(1,1)处将系统近似线性化,得 ,由于矩阵31xx的特征根 ,根据李雅普诺夫定理可知系统在点(1,1)附近不330稳定。在(1,1)处将系统近似线性化,得 ,由于原系统为定常系统,且31xx矩阵 的特征根 ,根据李雅普诺夫定理可知系统在点(1,1)和330点(1,1)附近不稳定。该题求解时往往容易忽略平衡点(1,1)和(1,1) 。 5.3 试用李雅普诺夫稳定性定理判断下列系统在平衡状态的稳定性。 123xx解 由于题中未限定利用哪一种方法,且系统为线性定常系统,所以利用第一方法比- 82 -较合适。经计算
5、知矩阵 的特征根为 ,所以系统在原点是大范围渐近123230稳定的。对于线性系统关于稳定性的结果是大范围的全局性结果。5.4 设线性离散时间系统为 01()() m02xkxk试求在平衡状态系统渐近稳定的 值范围。m解 令 由方程 得,QITGPQ121312133 30001 12Pmm解此方程得2210841Pm若要 应有 。 025.5 试用李雅普诺夫方法求系统 12axx在平衡状态 为大范围渐近稳定的条件。0x解 用李雅普诺夫第一方法。首先求系统矩阵的特征方程 121212212()0aIAaa- 83 -由韦达定理,两个特征值同时具有负实部的充要条件为 , 。120a1212a5.6
6、 系统的状态方程为 10xx试计算相轨迹从 点出发,到达 区域内所需要的时间。1(0)x221(.)解 由于 ,该系统发散, 单调增加。注意到12,()A21x,所以此题无解。 21(0)0.x5.7 给定线性时变系统, 01xxt0t判定其原点 是否是大范围渐近稳定。e解 取 ,则221(,)()Vxtt22x1 12()0)xtxt29(0t因为 ,所以系统在原点处大范围渐近稳定。 lim)xV5.8 考虑四阶线性自治系统, ,1423321400xbxAb 0i 4,321i应用李雅普诺夫的稳定判据,试以 , 表示这个系统的平衡点 渐近稳定i,0x- 84 -的充要条件。解 在李雅普诺夫
7、矩阵方程式WVA中,令 为W210b显然, 是半正定矩阵。求矩阵方程式的解 , 是对称矩阵。V434213214312vv将方程左边的 行 列元素记成 元素,可求得下面的一系列等式:ij),(ji(1,1)元素 0214vb(1,2)元素 243(1,3)元素 12vv(1,4)元素 0243b(2,2)元素 312v(2,3)元素 3243(2,4)元素 0114vbv(3,3)元素 23423(3,4)元素 4214vv(4,4)元素 43b由对于(1,1)、(2,2)、(3,3) 、(4,4)元素的等式和 , 有0ib4,321i, , ,012v2334vb由对于(1,3)、(2,4)
8、、(1,4) 元素的等式,有, ,1402413- 85 -由(1,2)、(2,3)、(3,4) 元素,有, ,241vb3v42vb因此,12341 1232,bv4v即, 121231234 00bV为对角线矩阵。因为 为半正定阵,所以要检查 在原点 以外的 是否满足系统状态方WTxWxx程。由于满足 的 同时满足 ,而 时,状态方程的解为0Tx044,所以满足 的状态方程的解只有 。 3214xTx0x由李雅普诺夫的稳定判据, 是渐近稳定的充要条件是对角矩阵 为正定阵。因V此 , , , 是求的充要条件。 01b203b45.9 下面的非线性微分方程式称为关于两种生物个体群的沃尔特纳(V
9、olterra)方程式1122dxaxt式中, 、 分别是生物个体数, 、 、 、 是不为零的实数。关于这个系统,1) 1x2试求平衡点;2) 在平衡点的附近线性化,试讨论平衡点的稳定性。解1) 由 , ,得01dtx2t0)(12122xx同时满足这二式的 、 有两组: 、 和 、 。即,系1x20/2x统的平衡点为:- 86 -平衡点(a) 、01x2平衡点(b) 、/x2) 分两种情况讨论。 平衡点(a)线性化的微分方程为 21210xxdt其特征方程式是 )(s、 时,平衡点(a)稳定,除此以外不稳定。0 平衡点(b)令 , ,得11/xx22/x 22/1/2 )()( xdt xx
10、 11/2/12 )()(txxx 因此,在平衡点(b)线性化的微分方程式是 21210/xxdt其特征方程式为 2s时,特征根是 ,为正、负实数,平衡点(b)不稳定。 时,特征根是0 0,为共轭纯虚数,平衡点(b)的稳定性在这样的线性化范围内不能决定。 j5.10 对于下面的非线性微分方程式试求平衡点;在各平衡点进行线性化,试判别平衡点是否稳定。 1212sinxx解 由 , ,知系统的平衡点是 ,02x0si21 ,201,x- 87 -。02x1) 在 , 处,将系统近似线性化得,421, 02x2121xdt其特征多项式是 。这是胡尔维茨多项式,因此这些平衡点渐近稳定。2s2)在 ,,
11、31x02x1122xdt特征多项式是 ,这不是胡尔维茨多项式。因此这些平衡点不稳定。 2s5.11 利用李雅普诺夫第二方法判断下列系统是否为大范围渐近稳定:123x解 令矩阵 12pP则由 TAPI得 12121033pp解上述矩阵方程,有 11222 174430 865pp即得 1274853pP- 88 -因为 121 757 1480 detdet034pP可知 是正定的。因此系统在原点处是大范围渐近稳定的。系统的李雅普诺夫函数及其沿轨迹的导数分别为 2212()(406)08TTVxPxxQA又因为 ,所以系统在原点处大范围渐近稳定。 lim()x5.12 给定连续时间的定常系统
12、122()xx试用李雅普诺夫第二方法判断其在平衡状态的稳定性。解 易知( )为其唯一的平衡状态。现取 ,且有:120, 21()Vx为正定2(i)Vx121212()i ()()xVx容易看出,除了两种情况(a) 任意,1x20(b) 任意, 1时 以外,均有 。所以, 为负半定。()0Vx()Vx()Vx(iii )检查 是否恒等于零。考虑到使得 的可能性只有上述两种0;,t()0x情况,所以问题归结为判断这两种情况是否为系统的受扰运动解。先考察情况(a): ,则由于 可导出 ,将此代入系统的方程可得:01(;,)(,Ttxt2()0xt2()t- 89 -12211()00()()(xtx
13、txtt这表明,除了点( )外, 不是系统的受扰运动解。再12, 0;,T考察情况(b): ,则由 可导出 ,将此代入系01(;)(,Ttxt2()xt2()0xt统的方程可得: 12211()0()()(xtxtxtt显然这是一个矛盾的结果,表明 也不是系统的受扰运动解。综上分0;,T析可知, 。0(;,)Vtx(iv)当 时,显然有 。212()Vx于是,可以断言,此系统的原点平衡状态是大范围渐近稳定的。 5.13 试用克拉索夫斯基定理判断下列系统是否是大范围渐近稳定的。1232xx解 显然 是系统的一个平衡点。021()3Fxx26()()T x由 和 知 。由克拉索夫斯基定理可知602
14、213406x()F系统在原点渐近稳定。又因为 23211lim()li()()Txxf xx所以原系统在原点处是大范围渐近稳定的。 5.14 试用克拉索夫斯基定理判断下列系统的稳定性。- 90 -2113233xx解 显然 是系统的一个平衡点。0x1236()9xF322340()()618T xxx由 , , ,知 。602480x3223236070xx()0Fx由克拉索夫斯基定理可知系统在原点渐近稳定。又因为 2222321313lim()li()()()Txxf xxx所以原系统在原点处是大范围渐近稳定的。 5.15 试用克拉索夫斯基定理确定使下列系统 1252xabx的原点为大范围
15、渐近稳定的参数 和 的取值范围。解 41()5aFxbx41()()25Tabx因为系统在原点渐近稳定,所以当 ,应有 , 又 时, 的0x()0Fx()0Fx充要条件为 。于是 应满足 。又因为系统大范围渐近稳20,a45ba1定,所以当 时,应有 。注意 , 时, 的充要条x()Fxx()x件为 ; 时, 的充要条件为 。综上, 的取值范围为:b100b,ab- 91 -,或 。 1,0ab1,0ab5.16 试用变量梯度法构造下列系统的李雅普诺夫函数 31122xx解 设 的梯度为V12ax于是 的导数为 121212()()()TVxaxax试取 则12,0,a221()xx当 时,
16、。注意到 满足旋度方程 ,所以可210xV12V122Vx知 122100xxdx由这个李亚普诺夫函数可看出,在 范围内,系统是渐近稳定的。 215.17 用变量梯度法求解下列系统的稳定性条件。 1212()()xattx解 设 的梯度为V12xa于是 的导数为 12()TVxx 2121()catxc- 92 -12112 2120()()catxxcat 显然,当 , 时, 。注意到 满足旋度方2()0cat12()0tV12cxV程 ,于是可知122Vx,12 22100()0xcdcx 12,0c由上式可看出,对于给定的 ,当 时,不难确定 使得 。从1)at212()cat而可得系统是渐近稳定的充分条件是 , 。 1(02()t