线性系统的稳定性分析.doc

1、关于线性系统稳定性的进一步探究任何一个实际系统总是在各种偶然和持续的干扰下运动或工作的。显然，我们首先要考虑的问题是，当系统承受这种干扰之后，能否稳妥地保持预定的运动轨迹或者工作状态，这就是稳定性。此外，我们知道，描述系统的数学模型，绝大部分都是近似的，这或者是由于量测误差，或者是为使问题简化，而不得不忽略某些次要因素。近似的数学模型能否如实反映实际的运动，在某种意义上说，也是稳定性问题。系统的稳定性在控制中是一个很重要的问题。在学习完稳定性理论之后，对此有了更为深刻的理解，不单单停留在输出跟踪输入的浅显印象之上，获益匪浅。因此，本文根据黄琳院士较为精炼的数学讲解，描述了一些自己对该问题的直观

2、思考，并且结合线性系统和具体实例对稳定性作进一步分析，使内容不再过于抽象，更为深入地理解其应用价值。1 预备理论1.1 微分方程解的表示考虑微分方程 0(,)xft其解 是自变量 的函数，而 ， 变动时对应的解也随着变动，故它应()xttt该是自变量 与初值 、 的函数, 可记为 。0x0(;,)xt例如： 000(;,)()t tteex问题：当初值变动时，对应的解如何变动？在应用上的意义是：初值通常是用实验方法求得的，实验测得的数据不可能绝对准确，若微小的误差会引起对应解的巨大变动，那么所求的初值问题解的实用价值就很小。1.2 Lipschitz 条件012(,),)():xfttIxWR

3、的定义域记为 。若存在常数 L，使得对任何 都有(,)fxtI I,Wtxy(,),ftfytxy则称 在 上满足 Lipschitz 条件。这个定义可以推广到 W 为任意有限 n 维fWI空间的情形。注：满足 Lipschitz 条件可保证微分方程解的存在性和唯一性1.3 解的存在性、唯一性及对初值的连续依赖性定理 1-1 (存在性及唯一性定理)对于微分方程 (,)xft若 在 域内连续且满足 Lipschitz 条件，则对任意的初始条件(,)fxtWI总存在常数 ，使得有唯一解 ，在0,0a0(;,)xt上存在、对 t 连续 ，且满足初始条件 。ta 稳定性所要研究的是解的渐近性质，即当解

4、 在 时的性状。故总t假定在 上解是存在的。0,t定理 1-2 (解对初值的连续依赖性)在定理 1 的条件下，若 在域内连(,)fxt续且满足 Lipschitz 条件，则微分方程的解 作为 ， ， 的函数在它0(;,)xtt0的存在范围内是连续的，即， ， 00()tt，0(;,)(;,xtt0,abt以上定理说明：若在初始时刻 和 十分接近，则在定义域 内0)x0(t ,ab的解 和 也会十分接近。因此，1.1 中所提的问题也就迎刃而解了。()tt2 平衡状态的稳定性李雅普诺夫稳定性的概念是微分方程解对初值的连续依赖性这一概念在无穷时间区间上的推广和发展。因此下面讨论时均假定所研究方程的解

5、在无穷区间 满足存在和唯一性条件。0,t2.1 平衡状态考虑系统 (,)nxftR若随着时间 的变化，状态 保持不变（即恒为常数） ，则称这个状态te为系统的平衡状态。由于平衡状态也是系统的一个状态，故它是上述微分方程的一个解，即 的解。(,)0efxt例：微分方程 显然 是它的一个解并且是它的一个平衡状态。A0ex2.2 简化的平衡状态在初始时刻 时，干扰引起的状态向量 与平衡状态 之差 称0t 0xex0eyx为初始扰动向量。由 所决定的运动过程是 的解，成为被扰运动，记x(,)ft为 。0(),)xtt由于平衡状态和被扰动运动均为微分方程 的解。由此可导出扰动(,)xft向量 应服从微分

6、方程():eytx()(,)(,):(,)e edyxftfytGytt称为关于平衡状态 的扰动方程，即e (,)yt其中， 满足 。这是因为(,)Gyt(0,)t,(,)0eGtfxt因此，在下面考虑一般的时变、非线性、多变量系统时，我们总假定它的微分方程(2-1)F(,)xt满足(2-2)(0,)t其中 为 n 维向量， 为 n 维的函数向量。这时方程（2-1）有解x(,)Fxt(满足 ),称为（ 2-1）的显然解或零解。00t在以下讨论平衡状态的稳定性时，只需要讨论零解这个平衡状态的稳定性就可以了。2.3 李雅普诺夫稳定性定义设有一个初始扰动，使系统的状态偏离了平衡状态 。若初始扰动为0

7、x，显然在这个初始扰动作用下，方程（2-1 ）所决定的运动是下列初值0xt问题 0(,)xFt的解。将这个解表示为 。(),)xtt例：考虑微分方程 4x显然， 是它的一个平衡状态。现若有初始扰动 则其解为0ex x(0)=.1，。可见，即使初始值微小地偏离了平衡状态，且在任意有限的4.1,tx时间内其解有界，但最终将发散。例：考虑微分方程 4x显然 是它的一个平衡状态，先若有初始扰动0ex x(0)=1则其解为 。事实上无论初始扰动多么大，最终将收敛到平衡状41,t态。以上两个例子是熟悉的线性系统的稳定和不稳定的例子。从第一个例子还可以看到，尽管在任意有限的时间内解是有界的，但若讨论时间趋于

8、无穷(或在工程上，当时间“ 很长” ）时系统的行为，则这种发散的特性就是完全不能接受的了。Lyapunov 稳定性就是要研究微分方程的解在 上的有界性。0,t根据微分方程解对初值的连续依赖性质，可知只要 充分小，对于x之间的任一时刻， 偏离 （平衡状态）也可以任意小。现在0,tT0(;,)xtx要研究这一性质是否对 均成立。定义 2-1 对于任意的 0 都存在 ，使得当 时有0(,)t00(),)xtt(;,xt成立。则称系统关于平衡状态（或原点） 是（李雅普诺夫意义下）稳定的。x定义 2-2 若定义 2-1 中的 ，即 与 无关(关于 一致)，则称所()0t0t定义的稳定为一致稳定。定义 2

9、-1（李雅普诺夫意义下稳定）的图示： 0()xt0t t图 1 李雅普诺夫意义下的稳定(1) 此处 随着 、 而变化；0t(2) ， ，初值变化充分小时，解的变化 可任意小0(;,)xt0t(不是无变化) ；(3) 显然， 。0(,)t李雅普诺夫意义下稳定也可表示为 0x()xt0(,)t图 2 李雅普诺夫意义下的稳定例如：讨论下列系统的稳定性和一致稳定性： 10xx其解为 0 0()()11220,t teext任给 ，取 (与 无关)，则只要00t1201020()()xtxtt，就有 0 0()()1212202()()t ttteextxt，故系统是(李氏) 稳定的。又由于 ，即与 无

10、关，系统还是一致稳定的。0t定义 2-3 （不稳定的定义）若对任意给定的 ，无论 多么小，总可以找到满足 的某一初值 ，使得从它出发的运动轨线 在某一0()xt0x 0(;,)xt时刻 ，有 ，则称系统(2-1)的零解是不稳定的。1t0;,0()xt0t t1t图 3 不稳定定义 2-4 （渐进稳定定义）若 (a) 是稳定的；(b)存在 ，使得0x0()t对任意的 ，存在 ，当 ， 时，有00(,)Ttx()tt0,tTx，则称 为渐近稳定。0(;,)xt(a) 是稳定的， 在 的行为已决定；0t(b) 是 充分大时的性质。t0()t0t t00(,)tTtx图 4 渐进稳定(1) 此处 是固

11、定的一个范围(称为吸引区，不是任意小的 )；0)t(2) ，;,x00(,)tTtx讨论：(1) 定义 2-4 的第二部分(b)又称为关于零解是吸引的。它反映的是解的渐近性质。可以将(b)改成：存在 ，使得 蕴涵0()t00()xtt0lim,tt(2) 稳定和吸引（即(a)和(b)）是相互独立的概念，对于一般的系统，它们之间不存在蕴涵关系。苏联人给出了一个著名的反例 (参见黄琳“稳定性理论” ，1992，p.7 ） ，表明一个微分方程的解是吸引的但却不是稳定的。(3) 正数 称为系统渐近稳定的吸引区。若吸引区是整个空间，称系统0()t是关于原点全局渐近稳定的。定义 2-5 （一致渐进稳定定义

12、）若(a) 是一致稳定的。(b)存在 ，0x0使得对任意的 ，存在 ，当 ，当 时有0()xt()T0()tt0()tT，则称 为一致渐近稳定，即0(;,)xt00(,)txxt 关 于 、 均 一 致这里，一致性在于： 不依赖于 、且 T 仅依赖于 ，不依赖于 、 。0tx定义 2-6 （指数渐进稳定定义）若存在 ，对任意的 ，存在 ，()使得当 ，就有0()xt0()00(;,)txtet成立。则称 是按指数渐近稳定的。显然，以上定义关于 、 是一致的。0t这里所定义的稳定、一致稳定、渐近稳定、一致渐近稳定和按指数渐近稳定都是局部的概念，即定义中的条件只要在 的附近成立即可。但在工程技0x

13、术上，特别是在控制系统中，所发生的初始偏差并非任意的小，而是有限的或是任意大的。幸好，就我们所讨论的线性系统而言，全局和局部是一致的。指数渐近稳定一致稳定稳 定指数渐近稳定渐近稳定图 5 各种稳定性之间的关系3 运动的稳定性前一节讨论了动态系统的一种特殊的运动平衡状态的稳定性，现在来讨论系统(3-1)F(,)xt任一运动的稳定性问题。我们已经知道，每一个初始状态 确定唯一的解 。一个系0()t 0(,)xt统随着初始条件不同可以有很多不同的运动。现在，设我们关心(2-1)的某一个运动： 00(),()xftxft我们欲研究这个运动的稳定性。我们称这个运动为给定运动，或未被扰运动。进而，设于初始

14、时刻 ，系统受到干扰，状态由 变成 。从这一初0t 00xy始状态出发的运动，即初值问题 00()F,()xttxfty的解，称为被扰运动。类比于平衡状态的稳定性（李氏稳定、一致稳定、渐近稳定等等） ，我们也可以相应地定义相对于给定运动的稳定性（李氏稳定、一致稳定、渐近稳定等等） 。定义 3-1 对于任意的 都存在 ，使得当00(,)t时有00(,)(,)xtft000(;,)(;,),xtftxt成立。则称系统关于给定运动 是（李雅普诺夫意义下）稳定的。但需要指出，关于给定运动的稳定性可以变换成关于零解的稳定性问题，故上述定义事实上是不必要的。为此，考虑变换 ，则扰动方程定义为：yxf(3-

15、2)()()F,(),:ytxftxtftGy则显然 (0,)t结论：这说明，通过上述变换可以将给定运动（或称为未被扰运动）的稳定性问题化为(3-2)的零解稳定性问题。也就是说，今后讨论运动的稳定性时，可先列出其扰动方程，然后讨论扰动方程(3-2)零解的稳定性就可以了，而没有必要再给出运动稳定性的其它定义。4 线性系统稳定性的特点考虑(4-1)A()Bxttu其对应的齐次方程为(4-2)()t（4-1）式比一般的方程（3-1）式的结构要简单，因此它在稳定性方面有更多的简单特性。定理 4-1 对于方程(4-1) 所表示的线性系统，若有一个运动稳定，则其所有运动稳定。因此，对线性系统而言，今后可笼

16、统地说“系统是稳定的” ，而一般的非线性系统并不具备这一特性。讨论线性系统 A()B(.1)xttuB在任意输入 u 作用下任一实际运动的稳定性，等价于讨论其所对应的齐次方程 () (.2)t关于零解的稳定性且(B.1)具有什么性质的稳定性等价于(B.2)具有同一种性质的稳定性。例：讨论如下系统的稳定性： 05,()xtx根据上面的分析，只需要讨论所对应的齐次方程的零解稳定性即可。齐次方程渐近稳定，故原系统渐近稳定。此外，注意到在这个例子中，系统的响应是无界的。这是由于输入信号是无界的。这和系统的稳定性不是同一个概念。t未被扰运动 被扰运动图 6 无界响应5 线性系统的稳定性判据由于线性动态方

17、程的稳定性等价于其对应的齐次方程的零解的稳定性，故这里只讨论齐次方程(5-1)A()xt对于（5-1 ）零解的稳定性问题。由于 不是常量矩阵，因此一般不能用特征值来讨论系统运动的性质，而应该用与系统运动关系密切的状态转移矩阵 。0(,)t定理 5-1 设 是连续（或分段连续）的函数矩阵，则有以下充分必要条()At件成立：(1) （5-1）稳定 存在某常数 ，使得对于任意的 和()xt0()Nt 0t有000(,)tt(2) （5-1 ）一致稳定 (1)中的 与 无关(3) （5-1)渐近稳定 0lim(,)tt(4) （5-1 ）一致渐近稳定 存在 N、 C 0，使得对于任意的 和 有0t0t

18、0()0(,)tte结论：对线性系统(a) 李氏稳定等价于状态转移矩阵范数的有界性；(b) 一致稳定等价于状态转移矩阵范数的一致有界性；(c) 渐近稳定等价于状态转移矩阵范数趋向于零；(d) 一致渐近稳定等价于状态转移矩阵按指数规律稳定。讨论：(1) 定理 5-1 所给出的线性系统的重要性质，完全是由 中，00(,)(,)xtt对 的线性关系所致。状态转移矩阵 决定了一切性质。0(,)xt0x一般地，对于非线性系统，定理 5-1 的结论均不成立。(2) 线性系统的稳定性具有全局性质。定义 系统 的零解称为是全局（一致）渐近稳定的，若其零解是F(,)xt（一致）渐近稳定的且无论初始扰动多大，均有

19、 lim()0tx定义 对任意的 x(0) , 均有 x(t)有界,则称 =A(t)x 的零解是李雅普诺夫意义下稳定的。定理 5-1 之(3)、(4) 清楚地表明，对于线性系统 =A(t)x 而言，若其零解是x（一致) 渐近稳定的，那么由状态空间任一点为起点的运动轨线都要收敛到原点，即原点的渐定稳定的吸引区遍及整个状态空间，这就是上面定义所述的全局（一致）渐近稳定或大范围（一致）渐近稳定的概念。(3) 对于线性系统而言，零解的吸引性蕴涵其稳定性，而一般的非线性系统则不具备这一性质（此性质将进一步讨论） 。(4) 再回到方程 A()Bxttu已经证明，其扰动运动的稳定性等价于对应的齐次方程零解的

20、稳定性。注：对于线性系统，零解的吸引性蕴涵稳定性。参考文献1 高为炳编著：运动稳定性基础，高等教育出版社， 1987 年 5 月.2 黄琳：稳定性理论，北京大学出版社，1992 年 7 月.3 秦元勋、王慕秋、王联：运动稳定性理论与应用, 科学出版社，1980 年.4 王柔怀、伍卓群编：常微分方程讲义, 人民教育出版社, 1978 年 5 月.5 黄琳：稳定性与鲁棒性的理论基础,科学出版社，2003 年 2 月.6 LaSalle, J. P., Stability by Lyapunov direct method, New York: Academic Press, 1961.7 Hahn, W., Stability of motion, New York, Springer-Verlag, 1967.8 Desoer, C.A. and Vidyasagar, M., Feedback systems: Input-output properties, New York: Academic Press, 1975.