第五章李雅普诺夫稳定性分析
第六章 李雅普诺夫稳定性分析
在反馈控制系统的分析设计中，系统的稳定性是首先需要考虑的问题之一。因为它关系到系统是否能正常工作。
经典控制理论中已经建立了劳斯判据、Huiwitz 稳定判据、Nquist 判据、对数判据、根轨迹判据等来判断线性定常系统的稳定性，但不适用于非线性和时变系统。分析非线性系统稳定性及自振的描述函数法，则要求系统的线性部分具有良好的滤除谐波的性能；而相平面法则只适合于一阶、二阶非线性系统。
1892年俄国学者李雅普诺夫（Lyapunov ）提出的稳定性理论是确定系统稳定性的更一般的理论，它采用状态向量来描述，不仅适用于单变量、线性、定常系统，还适用于多变量、非线性、时变系统。
§6-1 外部稳定性和内部稳定性
系统的数学模型有输入输出描述（即外部描述）和状态空间描述（即内部描述），相应的稳定性便分为外部稳定性和内部稳定性。
一、外部稳定性
1、定义（外部稳定性）：
若系统对所有有界输入引起的零状态响应的输出是有界的，则称该系统是外部稳定的。
(外部稳定性也称为BIBO （Bounded Input Bounded Output ）稳定性) 说明：
（1）所谓有界是指如果一个函数)(t h ，在时间区间],0[∞中，它的幅值不会增至无穷，即存在一个实
常数k ，使得对于所有的[]∞∈0
t ，恒有∞<≤k t h )(成立。
（2）所谓零状态响应，是指零初始状态时非零输入引起的响应。
2、系统外部稳定性判据
线性定常连续系统
∑),,(C B A 的传递函数矩阵为
Cx
y Bu Ax x =+=
BU
A sI X BU X A sI CX Y BU AX sX 1)()(--==-=+=
B A sI
C s G 1
)()(--=
当且仅当)(s G 极点都在s 的左半平面内时，系统才是外部稳定（或BIBO 稳定）的。
【例6.1.1】已知受控系统状态空间表达式为
u x x ??
????-+??????-=121160 ， []x y 10= 试分析系统的外部稳定性。
解：系统为SISO 系统，传递函数为
B A sI
C s G 1
)()(--=[]??
????-??????+--=-12116101
s s
)
3)(2(2
+--=
s s s
3
1+=
s 由于传递函数的极点位于s 左平面，故系统是外部稳定的。
二、内部稳定性
对于线性定常系统 Bu Ax x
+= ， 00)(x t x =
Cx y =
如果外部输入0)(=t u ，初始条件0x 为任意，且由0x 引起的零输入响应为
00),()(x t t t x φ=
满足
0),(lim 00=∞
→x t t t φ
则称系统是内部稳定的，或称为系统是渐近稳定的。
说 明：线性定常系统的渐近稳定与经典控制理论中的稳定性一致。
【例6.1.２】已知受控系统状态空间表达式为
u x x ??
????-+??????-=121160 ， []x y 10=
试分析系统的内部稳定性。
解：该系统为线性定常系统，其特征方程为：
0)3)(2(6)1(=+-=-+=-λλλλλA I
于是系统的特征值为21=λ，32-=λ，故系统不是内部稳定（渐近稳定）的。
三、内部稳定性与外部稳定性的关系
１、若系统是内部稳定（渐近稳定）的，则一定是外部稳定（BIBO 稳定）的。
2、若系统是外部稳定（BIBO 稳定）的，且又是可控可观测的，则系统是内部稳定（渐近稳定）的。此时内部稳定和外部稳定是等价的。
§6-2 李雅普诺夫稳定性的基本概念
一、自治系统
没有外界输入作用的系统叫自治系统。
自治系统可用如下的显含时间t 的状态方程来描述
),(t x f x = ， 00)(x t x =，0t t ≥………………………… (6-1)
其中x 为n 维状态向量。),(t x f 为线性或非线性、定常或时变的n 维向量函数。假定方程的解为
),;(00t x t x ，式中0x 和0t 分别为初始状态向量和初始时刻，那么初始条件0x 必满足0000),;(x t x t x =。
如果系统为线性系统，则(6-1)方程中的),(t x f 为x 的线性向量函数，或按习惯表示为：
x t A x
)(= ， 00)(x t x =，0t t ≥………………………… (6-2)
二、平衡状态
设控制系统的齐次状态方程为：
),(t x f x = ， 00)(x t x =，0t t ≥
对于所有t ，如果存在某个状态e x ，满足：
0),(==t x f x
e e
则称e x 为系统的一个平衡点或平衡状态。
平衡状态的各分量相对时间不再发生变化。若已知系统状态方程，令0=x
所求得的解x ，便是
平衡状态。
在大多数情况下，0=e x （状态空间原点）为系统的一个平衡状态。当然，系统也可以有非零平衡状态。如果系统的平衡状态在状态空间中表现为彼此分隔的孤立点，则称其为孤立平衡状态。对于
孤立平衡状态，总是可以通过移动坐标系而将其转换为状态空间的原点，所以在下面的讨论中，假定原点即0=e x 为平衡状态。
所谓系统运动的稳定性，就是研究其平衡状态的稳定性，也即偏离平衡状态的受扰运动，能否只依靠系统内部的结构因素而返回到平衡状态，或者限制在平衡状态的附近。
线性定常系统Ax x
= ，其平衡状态满足0=e Ax ，只要A 非奇异，系统只有唯一的零解，即存
在一个位于状态空间原点的平衡状态；当A 为奇异矩阵时，0=e Ax 有无数解，也就是系统有无数个
平衡状态。对于非线性系统，0),(=t x f e 的解可能有多个，由系统状态方程决定。
三、李雅普诺夫意义下稳定
设系统初始状态0x 位于以平衡状态e x 为球心、半径为δ的闭球域)(δS 内，即 ),(00t x x e εδ≤- 0t t =
若能使系统方程的解),;(00t x t x 在∞→t 的过程中，都位于以e x 为球心、任意规定的半径为ε的闭球域)(εS 内，即
ε≤-e x t x t x ),;(00 0t t ≥
则称该e x 是稳定的，通常称e x 为李雅普诺夫意义下稳定的平衡状态。以二维系统为例，上述定义的平面几何表示如图6-1所示。
式中?称为向量的范数，其几何意义是空间距离的尺度。如e x x -0表示状态空间中0x 至e x 点之间的距离的尺度，其数学表达式为
2
021100)()(ne n e e x x x x x x -++-=-
在上述稳定性的定义中，如果δ只依赖于ε而和初始时刻0t 的选取无关，则称平衡状态e x 是一致稳定的。对于定常系统，e x 的稳定等价于一致稳定。但对于时变系统，e x 的稳定并不意味着其为一致稳定。
要注意到，按李雅普诺夫意义下的稳定性定义，当系统作不衰减的振荡运动时，将在平面描绘出一条封闭曲线，但只要不超过)(εS ，则认为稳定，这同经典控制理论中线性定常系统稳定性的定义是有差异的。
四、渐近稳定
0x - 初始状态
e x - 平衡状态
图6-1 二维空间李雅普诺夫意义下稳定性的几何解释示意图
设e x 是系统),(t x f x
= ， 00)(x t x =，0t t ≥的一个孤立平衡状态，如果 （1）e x 是李雅普诺夫意义下稳定的； （2）0),;(lim
00→-∞
→e t x t x t x
则称此平衡状态是渐近稳定的。
实际上，渐近稳定即为工程意义下的稳定，也就是经典控制理论中所讨论的稳定性。当δ与0t 无关时，称平衡状态e x 是一致渐近稳定的。
五、大范围（全局）渐近稳定
当初始条件扩展到整个状态空间，且具有渐近稳定性时，称此平衡状态e x 是大范围渐近稳定的。 对于严格线性系统，如果它是渐近稳定的，必具有大范围渐近稳定性，这是因为线性系统稳定性
与初始条件的大小无关。一般非线性系统的稳定性与初始条件的大小密切相关，其δ总是有限的，故通常只能在小范围内渐近稳定。当δ与0t 无关时，称平衡状态e x 是大范围一致渐近稳定。
六、不稳定
不管把域)(δS 取得多么小，也不管把域)(εS 取得如何的大，只要在)(δS 内存在一个非零初始状态0x ，使得有0x 出发的运动轨迹超出域)(εS 以外，则称平衡状态e x 是不稳定的。
线性系统的平衡状态不稳定，表征系统不稳定。非线性系统的平衡状态不稳定，只说明存在局部发散的轨迹，至于是否趋于无穷远，要看)(εS 域外是否存在其它平衡状态，若存在，如有极限环，则系统仍是李雅普诺夫意义下稳定的。
下面介绍李雅普诺夫理论中判断系统稳定性的方法。
§6-3 李雅普诺夫稳定性判别方法
0x - 初始状态
e x - 平衡状态
图6-2 二维空间渐近稳定性的几何解释示意图
一、李雅普诺夫第一法（间接法）
这是利用状态方程解的特性来判别系统稳定性的方法，它适用于线性定常、线性时变以及非线性函数可线性化的情况。由于本章主要研究线性定常系统，所以在此仅介绍线性定常系统的特征值判据。
线性定常系统的特征值判据：
对于线性定常系统Ax x
= ，0)0(x x =，0≥t 有
（1）系统的平衡状态是在李雅普诺夫意义下稳定的充分必要条件是，
A 的所有特征值均具有非正
（负或零）实部，且具有零实部的特征值为A 的最小多项式的单根。
（2）系统的惟一平衡状态0=e x 是渐近稳定的充分必要条件是，A 的所有特征值均具有负实部。
二、李雅普诺夫第二法（直接法）
根据古典力学中的振动现象，若系统能量（含动能与位能）随时间推移而衰减，系统迟早回到达平衡状态，但要找到实际系统的能量函数表达式并非易事。李雅普诺夫提出，可虚构一个能量函数（后来被称为李雅普诺夫函数），一般它与n x x x ,,,21 及t 有关，记为),(t x V 。若不显含t ，则记为)(x V 。
它是一个标量函数，考虑到能量函数总是大于零，故为正定函数。能量衰减特性用),(t x V
或)(x V 表示。李雅普诺夫第二法利用V 及V
的符号特征，直接对平衡状态稳定性作出判断，无需求出系统状态方程的解，故称直接法。用此方法解决了一些用其它稳定性判据难以解决的非线性系统的稳定性问题，遗憾的是对一般非线性系统仍未形成构造李雅普诺夫函数的通用方法。对于线性系统，通常用二次型函数Px x
T
作为李雅普诺夫函数。
1、标量函数)(x V 符号性质的几个定义
（1）正定性
标量函数)(x V 在域s 中对所有非零状态（0≠x ）有0)(>x V 且0)0(=V ，则称)(x V 在域s 内
正定。如2
221)
(x x x V +=是正定的。
0x - 初始状态
图6-3 二维空间不稳定的几何解释示意图
（2）负定性
标量函数)(x V 在域s 中对所有非零状态（0≠x ）有0)( 22 1x x x V +-=是正定的。 （3）正半定性 0)0(=V ，且标量函数)(x V 在域s 内某些非零状态处有0)(=x V ，而在其它非零状态处有 0)(>x V 且，则称)(x V 在域s 内正半定。如221)2()(x x x V +=，当212x x -=时有0)(=x V ；当 212x x -≠时有0)(>x V ，故)(x V 为正半定。 （4）负半定性 0)0(=V ，且标量函数)(x V 在域s 内某些状态处有0)(=x V ，而在其它状态处有0)( 21)2()(x x x V +-=是负半定的。 （5）不定性 标量函数)(x V 在域s 内可正可负，则称)(x V 不定。如21)(x x x V =是不定的。 2、标量函数)(x V 取二次型时的符号 Px x x V T =)([]? ????? ??????????????????=n nn n n n n n x x x p p p p p p p p p x x x 212 1 22221 112112 1 式中P 为对称矩阵，有ji ij p p =。显然满足0)0(=V 。当P 阵的每一个元都为实数时，称作实二次型。 实二次型)(x V 是正定的充要条件是矩阵P 的各顺序主子行列式均大于零（赛尔维斯特准则），即 011>p ， 022 21 1211>p p p p ，…，01111>nn n n p p p p 则)(x V 正定，且称P 为正定矩阵。 当矩阵P 的各顺序主子行列式负、正相间时，即 011 022 21 1211>p p p p ，…，0)1(1111>-nn n n n p p p p 则)(x V 负定，且称P 为负定矩阵。 若矩阵P 的各顺序主子行列式含有等于零的情况，则)(x V 为正半定或负半定。不属于以上所有情 况的，)(x V 为不定。 【例6.3.1】证明下列的二次型是正定的 3132212 32221422410) (x x x x x x x x x x V --+++= 证明：上式用矩阵形式表示为 Px x x V T =)([]???? ? ???????????????----=32132 1 1121412110x x x x x x 由于01011>=p 0394 1 11022 21 1211>== p p p p 0171 1 214 1 2 11033 32 31 232221 131211 >=----=p p p p p p p p p 所以，)(x V 是正定的。 3、李雅普诺夫第二法 定理1 设系统的状态方程为 )(x f x = 0)0(=f 如果存在一个标量函数)(x V ，它有连续的一阶偏导数，而且满足：（1）)(x V 是正定的；（2）)(x V 是负定的，则系统在原点处的平衡状态是渐近稳定的。 【例6.3.2】设系统的状态方程是 )(2 221121x x x x x +-= )(2 22 1212x x x x x +--= 试分析系统在原点处的平衡状态是否为渐近稳定的。 解：令01=x 及02=x ，解得01=x ，02=x ，故原点为平衡状态，且只有一个平衡状态。 设2 22 1)(x x x V +=，则221122)(x x x x x V +=2 2221)(2x x +-= 显然，对于0≠x 存在0)( 以及0)0(=V ，故)(x V 是负定的。又由于)(x V 是正定的，所以该系统在原点处的平衡状态是渐进稳定的。因为只有一个平衡状态，故该系统是渐进稳定的。 【例6.3.4】设系统的状态方程是 ???????? ????--=??????21211110x x x x 试分析系统在原点处的平衡状态是否为渐近稳定的。 解：这个系统的平衡点是021==x x 。 （1）先猜想一个能代表系统能量的正定函数 22 212)(x x x V += 则：221124)(x x x x x V += 由原状态方程知：21x x = ，212x x x --= ，代入上式得 22 2122)(x x x x V -= 它是不定的，因而无法判断系统的稳定性。 （2）重新构造正定函数)(x V ： 2 221)(x x x V += 因此：2 22211222)(x x x x x x V -=+= 它是负半定的。但是因为)(x V 不恒为零，所以可知该系统是渐近稳定的。)(x V 是李亚普诺夫函数。 （3）再取函数 [] 2 2 212212)(2 1)(x x x x x V +++= 因此：)()(2 2 2 1x x x V +-= 是负定的。 且： ∞→x ，∞→)(x V ，故是大范围渐近稳定。 从上述例子可以看出，找到一个系统的李亚普诺夫函数证明系统是稳定时，则系统必定是稳定的；反之，如果用李亚普诺夫第二法不能判断系统是稳定时，并不能得出系统是不稳定的。 4.3 线性定常系统的稳定性分析 线性定常系统可以应用各种方法，诸如劳斯—霍尔维茨，奈奎斯特法来判断系统的稳定性。李亚普诺夫直接法也提供了一种稳定性判据的方法。在上一节的例4-3中，我们用猜想和试探、验证的方法来找李亚普诺夫函数，这是很费时的，对于高阶的系统尤其费时，有时甚至是困难的。在本节中，采用二次型函数来判定系统稳定性，这种方法较为简捷，并且还可以利用这种方法作为基础来解参数最优问题以及某些控制系统的设计问题。 设线性定常系统为 Ax x = (4-18) 式中x 是n 维状态向量。 如果采用二次型函数作为李亚普诺夫函数，即 ()Px x x V T = （4-19） 式中 P 是正定的实对称矩阵。那么对上式求导得 ()Px x x P x x V T T += （4-20） 将式（4-18）代入（4-20）式得 ()Px A x PAx x x V T T T += () x P A PA x T T += （4-21） 令 P A PA Q T +=- （4-22） 代入式（4-21）得 ()Qx x x V T -= （4-23） 从上式可以看出，为了判定(x)V 是否负定，只需要判定Q 是否正定。 因此，判定线性定常系统的稳定性可如下进行：选定一个正定实对称矩阵P ，按式（4-22）计算Q ， 如果Q 是正定的，则式（4-19）所表达的李亚普诺夫函数证明系统是稳定的。 实际应用这个方法时，常不是先假定P ，而是先假设一个正定矩阵Q,例如取Q=I,再按式（4-22） 计算P ，应用赛尔维斯特法则检查P 是否正定。如果P 是正定的，则系统是渐进稳定的。这样做法，有 时要方便得多。如果 ，()Qx x x V T -= 沿任意一条轨迹不恒等于零，那么Q 可取正半定的， [例4-4] 设系统的状态方程为 ????????????????????---=??????????321321x x x 10k 120010x x x 求使系统稳定的K 值。 [解] 假设 ?? ??? ?????=100000000Q 由式（4-23）得 ()[]23321321x x x x 100000000x ,x ,x x V -=?? ??? ???????????????-= 上式中()x V 不恒为零，所以这样选取Q 是合适的。 按（4-22）式计算P: ????? ?????---+??????????---??????????110021K 1010K 120010p p p p p p p p p 333231232221131211 ???? ? ?????-=??????????100000000p p p p p p p p p 3332 31 232221131211 解出上式得 ?????????? ?? ?? ? ?-------+=2K 1262K 12K 02K 12K 2K 123K 2K 126K 2K 126K 2K 1212K K P 2 根据赛尔维斯法则，如果P 是正定的，需使 0K 02K -12>> 所以 6K 0<< 作业：p82: 4-2 课堂练习： 4-1判定下列二次型函数的正定性 （1）3132212 32221242113x x x x x x x x x Q --+---= 123112321233()(32)(211)Q x x x x x x x x x x x x =-+-+--+---， ???? ??????-------=112123111 1P 常微分大作业--李雅普诺夫稳定性 11091059洪一洲 从19世纪末以来，李雅普诺夫稳定性理论一直指导着关于稳定性的研究和应用。不少学者遵循李雅普诺夫所开辟的研究路线对第二方法作了一些新的发展。一方面，李雅普诺夫第二方法被推广到研究一般系统的稳定性。例如，1957年，В．И．祖博夫将李雅普诺夫方法用于研究度量空间中不变集合的稳定性。随后，J.P.拉萨尔等又对各种形式抽象系统的李雅普诺夫稳定性进行了研究。在这些研究中，系统的描述不限于微分方程或差分方程，运动平衡状态已采用不变集合表示，李雅普诺夫函数是在更一般意义下定义的。1967年，D.布肖对表征在集合与映射水平上的系统建立了李雅普诺夫第二方法。这时，李雅普诺夫函数已不在实数域上取值，而是在有序定义的半格上取值。另一方面，李雅普诺夫第二方法被用于研究大系统或多级系统的稳定性。此时，李雅普诺夫函数被推广为向量形式，称为向量李雅普诺夫函数。用这种方法可建立大系统稳定性的充分条件。 1.李雅普诺夫稳定性概念 忽略输入后，非线性时变系统的状态方程如下 ),(t x f x = （1） 式中，x 为n 维状态向量；t 为时间变量；),(t x f 为n 维函数，其展开式为 12(,,,,)i i n x f x x x t = n i ,,1 = 假定方程的解为 ),;(00t x t x ，x 0和t 0 分别为初始状态向量和初始时刻，0000),;(x t x t x =。 平衡状态 如果对于所有t ，满足 0),(==t x f x e e （2） 的状态x e 称为平衡状态（又称为平衡点）。平衡状态的各分量不再随时间变化。若已知状态 方程，令0=x 所求得的解x ，便是平衡状态。 对于线性定常系统Ax x = ，其平衡状态满足0=e Ax ，如果A 非奇异，系统只有惟一的零解，即存在一个位于状态空间原点的平衡状态。至于非线性系统，0),(=t x f e 的解可能有多个，由系统状态方程决定。 控制系统李雅普诺夫意义下的稳定性是关于平衡状态的稳定性，反映了系统在平衡状态附近的动态行为。鉴于实际线性系统只有一个平衡状态，平衡状态的稳定性能够表征整个系统的稳定性。对于具有多个平衡状态的非线性系统来说，由于各平衡状态的稳定性一般并不相同，故需逐个加以考虑，还需结合具体初始条件下的系统运动轨迹来考虑。 第六章 李雅普诺夫稳定性分析 在反馈控制系统的分析设计中，系统的稳定性是首先需要考虑的问题之一。因为它关系到系统是否能正常工作。 经典控制理论中已经建立了劳斯判据、Huiwitz 稳定判据、Nquist 判据、对数判据、根轨迹判据等来判断线性定常系统的稳定性，但不适用于非线性和时变系统。分析非线性系统稳定性及自振的描述函数法，则要求系统的线性部分具有良好的滤除谐波的性能；而相平面法则只适合于一阶、二阶非线性系统。 1892年俄国学者李雅普诺夫（Lyapunov ）提出的稳定性理论是确定系统稳定性的更一般的理论，它采用状态向量来描述，不仅适用于单变量、线性、定常系统，还适用于多变量、非线性、时变系统。 §6-1 外部稳定性和内部稳定性 系统的数学模型有输入输出描述（即外部描述）和状态空间描述（即内部描述），相应的稳定性便分为外部稳定性和内部稳定性。 一、外部稳定性 1、定义（外部稳定性）： 若系统对所有有界输入引起的零状态响应的输出是有界的，则称该系统是外部稳定的。 (外部稳定性也称为BIBO （Bounded Input Bounded Output ）稳定性) 说明： （1）所谓有界是指如果一个函数)(t h ，在时间区间],0[∞中，它的幅值不会增至无穷，即存在一个实 常数k ，使得对于所有的[]∞∈0 t ，恒有∞<≤k t h )(成立。 （2）所谓零状态响应，是指零初始状态时非零输入引起的响应。 2、系统外部稳定性判据 线性定常连续系统 ∑),,(C B A 的传递函数矩阵为 Cx y Bu Ax x =+= BU A sI X BU X A sI CX Y BU AX sX 1)()(--==-=+= B A sI C s G 1 )()(--= 当且仅当)(s G 极点都在s 的左半平面内时，系统才是外部稳定（或BIBO 稳定）的。 李雅普诺夫稳定性方法 李雅普诺夫第一方法又称间接法，它是通过系统状态方程的解来判断系统的稳定性。如果其解随时间而收敛，则系统稳定；如果其解随时间而发散，则系统不稳定。 李雅普诺夫第二方法又称直接法，它不通过系统状态方程的解来判断系统的稳定性，而是借助李雅普诺夫函数对稳定性作出判断，是从广义能量的观点进行稳定性分析的。例如有阻尼的振动系统能量连续减小（总能量对时间的导数是负定的），系统会逐渐停止在平衡状态，系统是稳定的。由于李雅普诺夫第一方法求解通常很烦琐，因此李雅普诺夫第二方法获得更广泛的应用。李雅普诺夫第二方法的难点在于寻找李雅普诺夫函数。迄今为止，尚没有通用于一切系统的构造李雅普诺夫函数的方法。 对于系统[]t ,f x x = ，平衡状态为,0e =x 满足()0f e =x 。如果存在一个标量函数()x V ,它满足()x V 对所有x 都具有连续的 一阶偏导数；同时满足()x V 是正定的；则 （1）若()x V 沿状态轨迹方向计算的时间导数()dt /)(dV V x x = 为半负定，则平衡状态稳定； （2） 若()x V 为负定，或虽然()x V 为半负定，但对任意初始状 态不恒为零，则平衡状态渐近稳定。进而当∞→∞→)(V x x 时，，则系统大范围渐近稳定； （3） 若()x V 为正定，则平衡状态不稳定。 判断二次型 x x x P )(V τ=的正定性可由赛尔维斯特 （Sylvester ）准则来确定，即正定（记作V(x)>0）的充要条件为P 的所有主子行列式为正。如果P 的所有主子行列式为非负，为正半定（记作V(x)≥0）；如果-V(x)为正定,则V(x)为负定（记作V(x)<0）；如果-V(x)为正半定,则V(x)为负半定（记作V(x)≤0）。 例： []正定。 则)(V 0 1121412110 ,0411 10,010x x x 1121412110x x x )(V 321321x x >---->>----=??? ????????????? 例： )x x (x x x ) x x (x x x 2 2212122221121+--=+-= (0,0）是唯一的平衡状态。设正定的标量函数为 ∞→∞→<+-=+--++-=+=??+??=+=)V(,且当0 )x 2(x )]x (x x x [2x )]x (x x [x 2x x 2x x 2x dt dx x V dt dx x V )(V x x )V(2222122212122221121221122112 2 21x x x x 故系统在坐标原点处为大范围渐近稳定。 基于正定二次型的 李雅普诺夫稳定性分析 张俊超 （控制科学与工程、控制理论与控制工程、2010010215） 摘要：李雅普诺夫稳定性理论以状态向量描述为基础，不仅适用于单变量、 线性、定常系统，而且适用于多变量、非线性、时变系统。但要应用李氏判据判断系统稳定性，就要涉及到系统矩阵A特征值的求解以及根据系统状态方程构造正定二次型的李雅普诺夫函数来判断系统稳定性。 1.问题的提出 我们在处理实际工程问题时，经常需要判断系统稳定性，一般稳定性判据都有一定局限性，李雅普诺夫稳定性理论是确定系统稳定性的一般的理论，不仅适用于单变量、线性、定常系统，而且适用于多变量、非线性、时变系统，它以状态向量描述为基础，结合正定二次型的相关知识对系统稳定性进行判断。 2.问题的求解 李雅普诺夫稳定性理论分析系统稳定性的两种方法： （1）利用线性系统微分方程的解来判断系统的稳定性 ——李雅普诺夫第一法（间接法） 李雅普诺夫第一法的主要内容 1）用一次近似式表示状态方程，即：X=AX+B(x) 如果A的全部特征值都具有负实部，则系统在平衡点xe=0处是稳定的， 且系统的稳定性与高阶项B(x)无关。 2）如果X=AX+B(x)的A的特征值中至少有一个具有正实部，则无论B(x)如何，系统在平衡点xe=0处为不稳定的。 3）如果X=AX+B(x)的A的含有等于零的特征值，则系统的稳定性由B(x)决定。李雅普诺夫第一法是根据系统矩阵A的特征值来判断系统的稳定性的。 （2）构造李雅普诺夫函数，利用构造的李氏函数判断系统稳定性 ——李雅普诺夫第二法（直接法） 观察振动现象，若系统能量(含动能和位能)随时间推移而衰减，则系统迟早会达到平衡状态。基本思想：若系统内部能量随时间↑而↓,最终到达静止状态，系统稳定。虚构一个能量函数（李雅普诺夫函数) V(x，t)=f（x 1，x 2 ， (x) n ，t） V(x)=f（x 1，x 2 ， (x) n ） V(x,t)或V(x)是一个标量函数。能量总大于零，故为正定函数。能量随随时间增加而衰减，即：V(x，t)或V(x)的导数小于零。 控制系统的李雅普诺夫稳定性分析 内容提要 稳定性是系统的又一重要特性。所谓系统的稳定性，就是系统在受到小的外界扰动后，被调量与规定量之间的偏差值的过渡过程的收敛性。显然，稳定性是系统的一个动态属性。在控制理论和控制工程中，无论是调节器理论、观测器理论还是滤波预测、自适应理论，都不可避免的要遇到系统稳定性问题。稳定性问题一直是一个最基本的和最重要的问题。 随着控制理论与控制工程所涉及的领域由线性时不变系统扩展为时变系统和非线性系统，稳定性分析的复杂程度也在急剧的增长。直到目前，虽然有许多判据可应用于线性时不变系统或其它各自相应类型的问题中，以判断系统稳定情况，但能同时有效地适用于线性、非线性、定常、时变等各类系统的方法，则是俄国数学家李雅普诺夫(Lyaponov)在19世纪所提出的方法。这就是控制系统稳定性分析的李雅普诺夫方法。李雅普诺夫稳定性理论是稳定性分析、应用与研究的最重要基础。 习题与解答 5.1 判断下列函数的正定性 1)222 1231213()2322V x x x x x x x x =++-+ 2)222 123121323()82822V x x x x x x x x x x =++-+- 3)22 131223()2V x x x x x x x =+-+ 4)222 123122313()104224V x x x x x x x x x x =+++-- 5)222 123122313()311242V x x x x x x x x x x =++-++ 解 1) 210()131011T T V x x Ax x x -?? ??==-?????? , 因为顺序主子式 2120, 50,1 3 ->=>- 2101 11300 1 1 --=> 所以0A >，()V x 为正定函数。 李雅普诺夫稳定性定理的应用—— 设计模型参考自适应律 2010.04.14 理论依据 李雅普李雅普诺夫直接法一致渐近稳定的条件：接致渐稳定条件V (x , t 正定V (x , t 负定? á假设可调系统与参考模型在数学模型的结构上完全相同，该设计要求设计可调参数的变化规律（自适应律），以使得可调系统的外特性能够完全趋于参考模型的外特性。 例题 试用李雅普诺夫稳定性理论设计参数可调试用李雅普诺夫稳定性理论设计参数调的模型参考自适应律，其中参考模型和可， 调系统的传递函数分别是： k ?(s =g 参考模型：s +a k v ?v ( s =可调系统：g s +a v 解：给予参考模型和可调系统以相同的输入u ，假设它们的输出分别是y 和y v ，当然它们都是可以直接量测的所要求的模型参考自是可以直接量测的。所要求的模型参考自适应律就是当 a v =a v (t , u , y , y v 及k v =k v (t , u , y , y v 时可调系统实现对参考模型的自适应时，可调系统实现对参考模型的自适应，即： =k ?k →0?k v ? =a ?a v →0?a ?e =y ?y → 0v ? 将参考模型和可调系统都写成微分方程的形式： y (t + a ? y (t = k ? u (t yv (t + av yv (t = kvu (t 于是：e (t = y (t ? y v (t = ku (t ? a ? y (t ? k v u (t + a v y v (t = k u (t ? k v u (t ? a y (t + a y v (t ? a y v (t + a v y v (t = ( k ? k v u (t ? a[ y (t ?y v (t ] ? ( a ? a v y v (t ~ ~ = ? a ? e (t + k ? u (t ? a ? y (t v 设系统的广义状态变量是 ~ ~ x (t = [ e(t k (t a (t ]T 则前述自适应的目标就是广义系统渐近稳定。为此取李雅普诺夫函数 ~2 ~ v( x = P e (t + P2 k (t + P3 a 2 (t 1 2 Pi > 0 显然是正定泛函，另一方面观察 ~~ ~~ v( x = 2 P e e + 2 P2 k k + 2 P3 a a 1 ~ ~~ ~ y + 2P k k + 2P a a ~~ = 2 P e (?a e + k u ? a v 1 2 3 ~ ~~ 2 ~ ~~ = ?2 P a e + 2 P e k u + 2 P2 k k ? 2 P e a yv + 2 P3 a a 1 1 1 ~ ~ ~ ~ = ?2 P a e 2 + 2 k ( P e u + P2 k ? 2a ( P e yv ? P3 a 1 1 1 显然，只要保证a > 0, ~ P e u + P2 k = 0, 1 ~ P e yv ? P3 a = 0 1 就能确保 v( x < 0 ，即为负定泛函。即为负定泛函即可求出~ P k = ? 1 e u, P2 ~ = P ey 1 a v P3 最 第5章 李雅普诺夫稳定性分析 本章讨论李雅普诺夫稳定性分析。主要介绍李雅普诺夫稳定性的定义以及分析系统状态稳定性的李雅普诺夫理论和方法;着重讨论李雅普诺夫第二法及其在线性系统和3类非线性系统的应用、李雅普诺夫函数的构造、李亚普诺夫代数(或微分)方程的求解等。最后介绍李亚普诺夫稳定性问题的Matlab 计算与程序设计。 一个自动控制系统要能正常工作,必须首先是一个稳定的系统,即当系统受到外界干扰时它的平衡被破坏,但在外界干扰去掉以后,它仍有能力自动地恢复在平衡态下继续工作。系统的这种性能,叫做稳定性。例如,电压自动调解系统中保持电机电压为恒定的能力、电机自动调速系统中保持电机转速为一定的能力以及火箭飞行中保持航向为一定的能力等。具有稳定性的系统称为稳定系统，不具有稳定性的系统称为不稳定系统。也可以说,系统的稳定性就是系统在受到外界干扰后,系统状态变量或输出变量的偏差量(被调量偏离平衡位置的数值)过渡过程的收敛性,用数学方法表示就是 ε≤Δ∞→)(Lim t x t 式中，)(t x Δ为系统被调量偏离其平衡位置的变 化量;ε为任意小的规定量。 如果系统在受到外扰后偏差量越来越大,显然它不可能是一个稳定系统。在经典控制理论中,借助于常微分方程稳定性理论,产生了许多线性定常系统的稳定性判据,如劳斯-胡尔维茨(Routh-Hurwitz)判据和奈奎斯特判据等，都给出了既实用又方便的稳定性判别及设计方法。但这些稳定性判据仅限于讨论SISO 线性定常系统输入输出间动态关系,讨论的是有界输入有界输出(BIBO)稳定性,未研究系统的内部状态变化的稳定性。再则，对于非线性或时变系统,虽然通过一些系统转化方法,上述稳定判据尚能在某些特定系统和范围内应用,但是难以胜任一般系统。现代控制系统的结构比较复杂,大都存在非线性或时变因素,即使是系统结构本身, 往往也需要根据性能指标的要求而加以改变,才能适应新的情况,保证系统的正常或最佳运行状态。在解决这类复杂系统的稳定性问题时,最通常的方法是基于李雅普诺夫第二法而得到的一些稳定性理论。 早在1892年,俄国学者李雅普诺夫就发表了题为“运动稳定性一般问题”的著名文献,建立了关于运动稳定性研究的一般理论。 李雅普诺夫把分析系统稳定性的方法归纳为两类，分别称为李雅普诺夫第一法和李雅普诺夫第二法。李雅普诺夫第一法(亦称间接法)是解描述系统动力学的微分方程式,然后根据解的性质来判断系统的稳定性的方法。对于线性定常系统,主要是根据系统极点的分布来判断系统的稳定性,即为经典控制理论的稳定性判 第六章李雅普诺夫稳定性分析 在反馈控制系统的分析设计中，系统的稳定性是首先需要考虑的问题之一。因为它关系到系统是否能正常工作。 经典控制理论中已经建立了劳斯判据、Huiwitz 稳定判据、Nquist 判据、对数判据、根轨迹判据等来判断线性定常系统的稳定性，但不适用于非线性和时变系统。分析非线性系统稳定性及自振的描述函数法，则要求系统的线性部分具有良好的滤除谐波的性能；而相平面法则只适合于一阶、二阶非线性系统。 1892 年俄国学者李雅普诺夫( Lyapunov )提出的稳定性理论是确定系统稳定性的更一般的理论，它采用状态向量来描述，不仅适用于单变量、线性、定常系统，还适用于多变量、非线性、时变系统。 §6-1 外部稳定性和内部稳定性系统的数学模型有输入输出描述(即外部描述)和状态空间描述(即内部描述) ，相应的稳定性便分为外部稳定性和内部稳定性。 一、外部稳定性 1、定义(外部稳定性) ：若系统对所有有界输入引起的零状态响应的输出是有界的，则称该系统是外部稳定 的。 ( 外部稳定性也称为BIBO( Bounded Input Bounded Output )稳定性) 说明： (1)所谓有界是指如果一个函数h(t) ，在时间区间[0, ]中，它的幅值不会增至无穷，即存在一个实 常数k ，使得对于所有的t 0 ，恒有h(t) k 成立。 (2)所谓零状态响应，是指零初始状态时非零输入引起的响应。 2、系统外部稳定性判据 线性定常连续系统(A,B,C) 的传递函数矩阵为 x Ax Bu y Cx sX AX BU Y CX (sI A)X BU X (sI A) 1BU G(s) C(sI A) 1 B 当且仅当G(s) 极点都在s的左半平面内时，系统才是外部稳定(或BIBO稳定)的。 例6.1.1 】已知受控系统状态空间表达式为 李雅普诺夫稳定性定理的应用 ——设计模型参考自适应律 2010.04.14 理论依据 李雅普接致渐稳定条件李雅普诺夫直接法一致渐近稳定的条件：(),V x t 正定 (),V x t ?负定 假设可调系统与参考模型在数学模型的结á 构上完全相同，该设计要求设计可调参数的变化规律（自适应律），以使得可调系统的外特性能够完全趋于参考模型的外特性。 例题 试用李雅普诺夫稳定性理论设计参数调试用李雅普诺夫稳定性理论设计参数可调的模型参考自适应律，其中参考模型和可， 调系统的传递函数分别是： k ?参考模型：a s s g +=)(可调系统：v v v a s k s g +=)( ? 解：给予参考模型和可调系统以相同的输入u ，假设它们的输出分别是y 和，当然它们都是可以直接量测的所要求的模型参考自v y 是可以直接量测的。所要求的模型参考自适应律就是当 及时可调系统实现对参考模型的自适应(,,,)v v v a a t u y y =) ,,,(v v v y y u t k k =时，可调系统实现对参考模型的自适应，即： ?00v k k k =?→? 0v a a a e =?→??=?→v y y ? 将参考模型和可调系统都写成微分方程的形式： ?? )()()()()()(t u k t y a t y t u k t y a t y v v v v =+=+ 于是： ) ()()()()()()(t y a t u k t y a t ku t y t y t e v v v v +???=?=)()()()()()(t y a t y a t y a t y a t u k t u k v v v v v ?????=+?+??=)(~)(~)() ()()]()([)()(t y a t u k t e a t y a a t y t y a t u k k v v v v v ???+?? =李雅普诺夫稳定性分析
第五章李雅普诺夫稳定性分析
李雅普诺夫稳定性方法
李雅普诺夫稳定性分析
李雅普诺夫稳定性分析报告
李雅普诺夫稳定性定理的应用汇总
李雅普诺夫稳定性分析
李雅普诺夫稳定性分析
李雅普诺夫稳定性定理的应用