第五章李雅普诺夫稳定性分析

合集下载

第5章李雅普诺夫稳定性分析

第5章李雅普诺夫稳定性分析
3
第5章 李雅普诺夫稳定性分析
第五章 李雅普诺夫稳定性分析
5.1 李雅普诺夫意义下的稳定性 5.2 李雅普诺夫第一法(间接法) 5.3 李雅普诺夫第二法(直接法) 5.4 线性定常系统的李雅普诺夫稳定性分析
4
第5章 李雅普诺夫稳定性分析
5.1 李雅普诺夫意义下的稳定性
1.自治系统
没有外输入作用时的系统称为自治系统,可 用如下系统状态方程来描述:
如果时变函数V(x,t)有一个正定函数作为下限, 也就是说,存在一个正定函数W(x) ,使得
V ( x ,t) W ( x), V (0,t) 0, t t0
则称时变函数V(x,t)在域S(域S包含状态空间的 原点)内是正定的。
24
第5章 李雅普诺夫稳定性分析
3. 负定函数:如果-V(x)是正定函数,则标量函数 V(x)为负定函数。
则称平衡状态xe在李雅普诺夫意义下是稳定的。
在上述稳定的定义中,实数δ通常与ε和初始时
刻t0都有关,如果δ只依赖于ε ,而和t0的选取无关,
则称平衡状态是一致稳定的。
9
第5章 李雅普诺夫稳定性分析
5. 渐近稳定性
若系统的平衡状态xe不仅具有李雅普诺夫意 义下的稳定性,且有
lim
t
||
x(t;
x0 ,
(s)
则 m(s) 为矩阵A的最小多项式。
注:换言之,矩阵A的最小多项式就是(sI-A)-1
中所有元素的最小公分母。
17
第5章 李雅普诺夫稳定性分析
例5-1(补充):判断下述线性定常系统的稳定性
0 0 0
x 0 0
0
x
0 0 1
解:1)系统矩阵A为奇异矩阵,故系统存在无穷

5李雅普诺夫稳定性分析.ppt

5李雅普诺夫稳定性分析.ppt
本章主要讨论李雅普诺夫稳定性分析,首先介绍了稳定性的定义,即系统受到外界干扰后仍能自动恢复平衡态的能力。进而阐述了李雅普诺夫稳定性的基本定理,为分析系统状态稳定性提供了理论基础。针对线性系统,虽然经典控制理论中的劳斯-赫尔维茨判据和奈奎斯特判据等关系,未涉及系统内部状态的稳定性。因此,李雅普诺夫理论在解决复杂系统,如非线性或时变系统的稳定性问题时具有更广泛的应用。其中,李雅普诺夫第二法通过定义李雅普诺夫函数来分析判别稳定性,无需解方程即可直接判断,为稳定性分析提供了便捷途径。此外,本章还介绍了Matlab在李亚普诺夫稳定性问题中的计算与程序设计,为实际应用提供了有力支持。

第五章 李雅普诺夫稳定性分析4.8

第五章 李雅普诺夫稳定性分析4.8

二、非线性系统的稳定性
1、非线性系统线性化 设系统的状态方程为:x f ( x, t )
xe 为平衡状态;f ( x, t ) 为与 x 同维的矢量函数,并且对 x
具有连续的偏导数。
将非线性矢量函数 f ( x, t ) 在 xe 邻域内展开成泰勒级数:
f x ( x xe ) R x x
塞尔维斯特(Sylvester)定理: V x xT Px
为正定的充要条件是P的所有顺序主子行列式都
是正的。如果P的所有主子行列式为非负的(其 中有的为零),那么V(x)为半正定的。 如果V(x)是正定的(半正定的),则-V(x)将是负定 的 (半负定的)。
例5.2.3
证明下列二次型函数是正定的。
图5.1(a)、(b)、(c)分别表示平衡状态为稳定、 渐近稳定和不稳定时初始扰动所引起的典型轨迹。
5.2 李雅普诺夫稳定性理论
5.2.1 李雅普诺夫第一方法
(间接法,通过系统状态方程的解来判定系统的稳定性。)
一、线性系统的稳定性
内部稳定 (平衡状态xe=0渐进稳定) BIBO稳定 系统矩阵A的所有特征值 均具有负实部 传递函数的所有极点均位 于s的左半平面
一个因果系统,如果对于任意一个有界输入
u (t ) 1 , t (t0 , )
对应的输出均有界
y (t ) 2 , t (t0 , )
则称该系统为外部稳定。 线性定常连续系统,BIBO稳定的充分必要条件为 其传递函数矩阵G(s)的所有极点都具有负实部。
5.1 几个稳定性概念
p11 p 21 P p n1
于是有:
P 为权矩阵(常取对称矩阵)。 式中,x T 为 x 的转置,

第5章 “控制系统的李雅普诺夫稳定性分析”练习题

第5章 “控制系统的李雅普诺夫稳定性分析”练习题

第5章 “控制系统的李雅普诺夫稳定性分析”练习题及答案5.1 判断下列函数的正定性1) 2221231213()2322V x x x x x x x =++-+x 2) 222123121323()82822V x x x x x x x x x =++-+-x 3) 22131223()2V x x x x x x =+-+x解1) T T 211()130101V A -⎡⎤⎢⎥==-⎢⎥⎢⎥⎣⎦x x x x x , 因为顺序主子式2120,50,13->=>- 2111302011--=> 所以0>A ,()V x 为正定函数。

2) T T 841()421111V -⎡⎤⎢⎥==--⎢⎥⎢⎥-⎣⎦x x Ax x x , 因为主子式8481218,2,10,0,70,10,421111-->==>=>--841421164421680111---=++---<- 所以A 不定,()V x 为不定函数。

3) T T 1212110()1001V -⎡⎤⎢⎥==-⎢⎥⎢⎥⎣⎦x x Ax x x , 因为顺序主子式1110,10,1->=-<- 121211011001041--=--<所以A 为不定矩阵,()V x 为不定函数。

5.2 用李雅普诺夫第一方法判定下列系统在平衡状态的稳定性。

2211211222212212()()x x x x x x x x x x x x =-+++=--++解解方程组 22121122212212()0()0x x x x x x x x x x ⎧-+++=⎨--++=⎩只有一个实孤立平衡点(0,0)。

在(0,0)处将系统近似线性化,得**1111x x -⎡⎤=⎢⎥--⎣⎦,由于原系统为定常系统,且矩阵1111-⎡⎤⎢⎥--⎣⎦的特征根1s i =-±均具有负实部,于是根定理5.3可知系统在原点(0,0)附近一致渐近稳定。

第五章李雅普诺夫稳定性分析

第五章李雅普诺夫稳定性分析
即 x e = f (xe , t) = 0 。
从定义可知,平衡状态的各分量相对于时间不再发生变化。
线性定常系统:x = Ax
A非奇异:Axe = 0 xe = 0 是唯一零解 A奇异:Axe = 0 xe 有无穷多个解
非线性系统:x = f (x,t)
x = f (xe , t) = 0 xe 可能有一个也可能有多个平衡状态
5-2 李雅普诺夫稳定性的基本概念
一、 平衡状态
系统x = f (x,t) ,X为n 维状态向量,且显含时间变量t,x = f (x,t)为线性或
非线性、定常或时变的n
维向量函数,假定方程的解为
x(t;
x
0
,
t 0
)
,式中
x
0
和 t0 分别为初始状态和初始时刻。
定义:系统 x = f (x,t) 的平衡状态是使x = 0的那一类状态,并用 xe 表示,
1 2
Mx22

若用标量函数 V (x) 表示系统的能量。则
V
(x)
=
1 2
Kx12
+
1 2
Mx22
V (x) = Kx1x1 + Mx2x2
=
Kx1x2
+ Mx2 (−
K M
x1

f M
x2 )
= − fx22 0
结论:坐标原点处的平衡状态是渐近稳定的。
一、标量函数及其定号性
1.标量函数 V (x) 的符号和性质
+ ... +
a1
+
a0
=
0
如何判断系统的渐近稳定性?
5-4 李雅普诺夫第二方法
李雅普诺夫第二方法,建立在用能量观点分析稳定性的基础上: 若系统的某个平衡状态是渐近稳定的,则系统储存的能量将随时

第五章 控制系统的李雅普诺夫稳定性分析汇总

第五章 控制系统的李雅普诺夫稳定性分析汇总
则状态方程的解为: x(t ) e At x(0) ( R1e1t ... Rnent ) x(0)
Re(i ) 0, (i 1, 2,..., n) lim x(t ) 0, 系统渐近稳定。
t
如果只有一个(或一对)特征值的实部等于0,其余特征值实 部均小于0,则系统仅仅可能是李亚普诺夫意义下的稳定性。
线性定常系统的特征值判据: 系统 x Ax 渐近稳定的充要条件是A的特征值均具有负实 部,即:Re( i ) 0 (i 1,2,, n) 证明:假定A有相异特征值 1 ,..., n 根据凯莱哈密顿定理:矩阵指数eAt为 e1t ,..., ent的线性组合
e At R1e1t ... Rn ent
x xe ( x1 xe1 ) 2 ... ( xn xen ) 2
2
2
2
由范数的定义可知,向量 ( x xe ) 的范数可写成
通常又将 x xe 称为 围之内时,则记为
x 与 xe 的距离。当向量 ( x xe ) 的范数限定在某一范
x xe
0
xe
与经典控制理论的区别: 1. 2. 3. 4. 5. 6. 平衡点/BIBO; 状态稳定/输出稳定; 经典控制的稳定大致对应于现代控制的渐进稳定; 即便输出稳定,状态可能不稳定; 李雅普诺夫意义下的稳定在经典中是不稳定的; 经典控制不需要一致性、全局性概念。
5.2 李雅普诺夫稳定性理论 一、李雅普诺夫第一方法 李雅普诺夫第一法的基本思想是利用状态方程解的性质来 判断系统的稳定性。通常又称为间接法。它适用于线性定常系 统以及线性时变系统和非线性系统可以线性化的情况。
意义:当系统运动到xe点时,系统状态各分量将维持平衡, 不再随时间变化。 平衡点:由系统状态在状态空间中所确定的点 求法:1、线性定常系统

第5章“控制系统的李雅普诺夫稳定性分析”练习题及答案.doc

第5章“控制系统的李雅普诺夫稳定性分析”练习题及答案.doc

第5章 “控制系统的李雅普诺夫稳定性分析”练习题及答案5.1 判断下列函数的正定性1) 2221231213()2322V x x x x x x x =++-+x 2) 222123121323()82822V x x x x x x x x x =++-+-x 3) 22131223()2V x x x x x x =+-+x解1) T T 211()130101V A -⎡⎤⎢⎥==-⎢⎥⎢⎥⎣⎦x x x x x , 因为顺序主子式2120,50,13->=>- 2111302011--=> 所以0>A ,()V x 为正定函数。

2) T T 841()421111V -⎡⎤⎢⎥==--⎢⎥⎢⎥-⎣⎦x x Ax x x , 因为主子式8481218,2,10,0,70,10,421111-->==>=>--841421164421680111---=++---<- 所以A 不定,()V x 为不定函数。

3) T T 1212110()1001V -⎡⎤⎢⎥==-⎢⎥⎢⎥⎣⎦x x Ax x x , 因为顺序主子式1110,10,1->=-<- 121211011001041--=--<所以A 为不定矩阵,()V x 为不定函数。

5.2 用李雅普诺夫第一方法判定下列系统在平衡状态的稳定性。

2211211222212212()()x x x x x x x x x x x x =-+++=--++解解方程组 22121122212212()0()0x x x x x x x x x x ⎧-+++=⎨--++=⎩只有一个实孤立平衡点(0,0)。

在(0,0)处将系统近似线性化,得**1111x x -⎡⎤=⎢⎥--⎣⎦,由于原系统为定常系统,且矩阵1111-⎡⎤⎢⎥--⎣⎦的特征根1s i =-±均具有负实部,于是根定理5.3可知系统在原点(0,0)附近一致渐近稳定。

《自动控制原理》第五章:系统稳定性

《自动控制原理》第五章:系统稳定性

5.2 稳定的条件
当σi和λi均为负数,即特征根的 σi和λi均为负数, 均为负数 实部为负数,系统是稳定的; 实部为负数,系统是稳定的; 或极点均在左平面。 或极点均在左平面。
5.3 代数稳定性判据
定常线性系统稳定的充要条件 定常线性系统稳定的充要条件是特征方程的根具有负 充要条件是特征方程的根具有负 实部。因此,判别其稳定性,要解系统特征方程的根。为 实部。因此,判别其稳定性,要解系统特征方程的根。 避开对特征方程的直接求解,可讨论特征根的分布, 避开对特征方程的直接求解,可讨论特征根的分布,看其 是否全部具有负实部,并以此来判别系统的稳定性,这样 是否全部具有负实部,并以此来判别系统的稳定性, 也就产生了一系列稳定性判据。 也就产生了一系列稳定性判据。 其中最主要是E.J.Routh(1877 )h和Hurwitz( 其中最主要是E.J.Routh(1877年)h和Hurwitz(1895 E.J.Routh(1877年 年)分别提出的代数判据。 分别提出的代数判据 代数判据。
习题讲解: 习题讲解:
µ
G1
Q21
G1
h2
k1 k1 G1 ( s ) = , G1 ( s ) = (T1s + 1) (T1s + 1) k1k 2 G0 ( s ) = (T1s + 1)(T2 s + 1)
kp
G0 ( s ) G(s) = 1 + G0 ( s ) K p
5.4 Nyquist稳定性判据 Nyquist稳定性判据
系统稳定的条件? 系统稳定的条件?
5.2 稳定的条件
d n y (t ) d ( n −1) y (t ) dy (t ) 线性系统微分方程: 线性系统微分方程: n a + an −1 + L + a1 + a0 y (t ) n ( n −1) dt dt dt d m x(t ) d ( m −1) x(t ) dx(t ) = bm + bm−1 + L + b1 + b0 x(t ) m ( m −1) dt dt dt d n y (t ) d ( n −1) y (t ) dy (t ) + a( n −1) + L + a1 + a0 y (t ) = 0 齐次微分方程: 齐次微分方程: an n ( n −1) dt dt dt an s n + an −1s n −1 + L + a1s + a0 = 0 设系统k 设系统k个实根

第五章李雅普诺夫稳定性分析讲诉

第五章李雅普诺夫稳定性分析讲诉

第六章 李雅普诺夫稳定性分析在反馈控制系统的分析设计中,系统的稳定性是首先需要考虑的问题之一。

因为它关系到系统是否能正常工作。

经典控制理论中已经建立了劳斯判据、Huiwitz 稳定判据、Nquist 判据、对数判据、根轨迹判据等来判断线性定常系统的稳定性,但不适用于非线性和时变系统。

分析非线性系统稳定性及自振的描述函数法,则要求系统的线性部分具有良好的滤除谐波的性能;而相平面法则只适合于一阶、二阶非线性系统。

1892年俄国学者李雅普诺夫(Lyapunov )提出的稳定性理论是确定系统稳定性的更一般的理论,它采用状态向量来描述,不仅适用于单变量、线性、定常系统,还适用于多变量、非线性、时变系统。

§6-1 外部稳定性和内部稳定性系统的数学模型有输入输出描述(即外部描述)和状态空间描述(即内部描述),相应的稳定性便分为外部稳定性和内部稳定性。

一、外部稳定性1、定义(外部稳定性):若系统对所有有界输入引起的零状态响应的输出是有界的,则称该系统是外部稳定的。

(外部稳定性也称为BIBO (Bounded Input Bounded Output )稳定性) 说明:(1)所谓有界是指如果一个函数)(t h ,在时间区间],0[∞中,它的幅值不会增至无穷,即存在一个实常数k ,使得对于所有的[]∞∈0t ,恒有∞<≤k t h )(成立。

(2)所谓零状态响应,是指零初始状态时非零输入引起的响应。

2、系统外部稳定性判据线性定常连续系统∑),,(C B A 的传递函数矩阵为Cxy Bu Ax x=+=BUA sI X BU X A sI CX Y BU AX sX 1)()(--==-=+=B A sIC s G 1)()(--=当且仅当)(s G 极点都在s 的左半平面内时,系统才是外部稳定(或BIBO 稳定)的。

【例6.1.1】已知受控系统状态空间表达式为u x x ⎥⎦⎤⎢⎣⎡-+⎥⎦⎤⎢⎣⎡-=121160 , []x y 10= 试分析系统的外部稳定性。

第5章现代控制理论之系统运动的稳定性分析

第5章现代控制理论之系统运动的稳定性分析
当然,对于线性系统, 从不稳定平衡状态出发的轨 迹,理论上趋于无穷远。
由稳定性定义知,球域S(δ) 限制着初始状态x0的取值,球域
S(ε)规定了系统自由运动响应 xt xt; x0的, t0边 界。
简单地说:1.如果 x t; x0, t0 有界,则称 xe 稳定;
2.如果 x t; x0, t0 不仅有界,而且当t→∞时收敛于原点,则
5.1.1 平衡状态
李雅普诺夫关于稳定性的研究均针对平衡状态而言。
1. 平衡状态的定义
设系统状态方程为: x f x,t , x Rn
若对所有t ,状态 x 满足 x 0 ,则称该状态x为平衡状
态,记为xe。故有下式成立:f xe ,t 0
由平衡状态在状态空间中所确定的点,称为平衡点。
2.平衡状态的求法
由定义,平衡状态将包含在 f x,t 这样0 一个代数方程组
中。
对于线性定常系统 x A,x其平衡状态为 xe 应满足代数
方程 。Ax 0
只有坐标原点处是线性系统的平衡状态点。
对于非线性系统,方程 方程而定。
如:
x1 x2
x1 x1
x2
x
3 2
f x的,t 解 可0 能有多个,视系统
稳定性是系统的重要特性,是系统正常工作的必要条件。
稳定性是指系统在平衡状态下受到扰动后,系统自由运动 的性质。因此,系统的稳定性是相对于系统的平衡状态而 言的。它描述初始条件下系统方程是否具有收敛性,而不 考虑输入作用。
1. 线性系统的稳定性只取决于系统的结构和参数,与系统 初始条件及外作用无关; 2. 非线性系统的稳定性既取决于系统的结构和参数,也与 系统初始条件及外作用有关;
当稳定性与 t0 的选择无关时,称一致全局渐近稳定。

第五章 李雅普诺夫稳定性理论

第五章 李雅普诺夫稳定性理论

非线性系统的稳定性是相对系统的平衡态而言的, 非线性系统的稳定性是相对系统的平衡态而言的,很难 笼统地讨论非线性系统在整个状态空间的稳定性。 笼统地讨论非线性系统在整个状态空间的稳定性。 对于非线性系统, 对于非线性系统,其不同的平衡态有着不同的稳定 故只能分别讨论各平衡态附近的稳定性。 性,故只能分别讨论各平衡态附近的稳定性。 对于稳定的线性系统,由于只存在唯一的孤立平衡 对于稳定的线性系统, 态,所以只有对线性系统才能笼统提系统的稳定性 问题。 问题。 李雅普诺夫稳定性理论讨论的是动态系统各平衡态 附近的局部稳定性问题。 附近的局部稳定性问题。 它是一种具有普遍性的稳定性理论, 它是一种具有普遍性的稳定性理论,不仅适用于线 性定常系统,而且也适用于非线性系统、时变系统、 性定常系统,而且也适用于非线性系统、时变系统、 分布参数系统。 分布参数系统。
5.1 动态系统的外部稳定性
控制系统的外部稳定性,常称为有界输入有界输出稳定性。 控制系统的外部稳定性,常称为有界输入有界输出稳定性。 在讨论系统的外部稳定性时,一般只适用于线性动态系统, 在讨论系统的外部稳定性时,一般只适用于线性动态系统, 而且必须假定系统的初始条件为零。 而且必须假定系统的初始条件为零。 外部稳定性的定义是,初始条件为零的线性系统, 外部稳定性的定义是,初始条件为零的线性系统,在任何一 个有界的输入作用下,若系统所产生的输出也是有界的, 个有界的输入作用下,若系统所产生的输出也是有界的,就 称该动态系统是外部稳定的,又称为BIBO稳定。 稳定。 称该动态系统是外部稳定的,又称为 稳定 对于单输入单输出线性定常系统而言, 对于单输入单输出线性定常系统而言,在经典控制理论中定 义的传递函数正是表征了系统在零初始条件下, 义的传递函数正是表征了系统在零初始条件下,输出量与输 入量两者间的关系。因此,对线性定常系统, 入量两者间的关系。因此,对线性定常系统,具有外部稳定 性的充要条件等价于其传递函数的所有极点都位于s平面的 性的充要条件等价于其传递函数的所有极点都位于 平面的 左半边。 左半边。

李雅普诺夫稳定性分析的方法

李雅普诺夫稳定性分析的方法
f ( x, t ), 且f (0, t ) 0 • 定理一.设系统的状态方程: x (坐标原点为平衡状态)如果上述给定系统 存在一个有连续偏导数的标量函数V(x)并 满足下列条件:
1).对所有 x 0 时V(x)>0 ( x) 0 ,则平衡点x=0是渐 2).对所有 x 0 时V 近稳定的. 3).除满足1),2)外,如果 x ,V ( x) 则x=0是大范围渐近稳定的.
(3)用李氏方法分析的必要性 • 以一个例子说明:用特征值来判断线性时变 系统一般稳定性是会失效的.
1 e 2 t x x 0 1
• 其中特征值为 -1,-1.
• 但由于其解为
et x(t ) 0 (e t e t ) / 2 x(0) t e
• 当 x(0) 0时,若 t 则必有 x • 故平衡状态是不稳定的,即系统的实际表现 是不收敛和发散的.从而采用特征值判断失 效.
一.系统运动稳定性的性质.
• 运动稳定性的实质,归结为系统平衡状态的 稳定性. • 平衡状态的稳定性问题实际就是:偏离平衡 状态的受扰运动能否只依靠系统内部的结 构因素,或者使之限制在平衡状态的有限临 域内,或者使之同时返回平衡状态.
1.预备知识
1).标量函数V(x)性质意义: 令V(x)是向量x的标量函数,Ω是x空间包含 原点的封闭有限区域. (1).如果对所有区域Ω中的非零向量x,有 V(x)>0,且在x=0处有V(x)=0则在域Ω内称 V(x)为正定.
(2).如果V(x)除原点以及某些状态等于零 外,在域Ω内其余状态处都是正的,则V(x)称
• 定理二
前提如定理一. 1).对所有 x 0 时V(x)>0
( x) 0 ,但不恒等于零,则 2).对所有 x 0 时 V

现代控制理论5.4 非线性系统的李雅普诺夫稳定性分析

现代控制理论5.4 非线性系统的李雅普诺夫稳定性分析

克拉索夫斯基法 (1/7)
5.4.1 克拉索夫斯基法
� 设非线性定常连续系统的状态方程为
̇ (t ) = f ( x ) x
� 对该系统有如下假设: 1) 所讨论的平衡态xe=0; 2) f(x)对状态变量x是连续可微的,即存在雅可比矩阵
J ( x ) = ∂f ( x ) / ∂xτ
� 对上述非线性系统 ,有如下判别渐近稳定性的克拉索夫斯 基定理。
0
1
x1
x2
0
(x1 , x2 ,0,⋯ ,0)
dx2 + ⋯ + ∫ ∇Vn (x , x ,⋯, x ) dxn
0
1 2
xn
n
变量梯度法 (5/10)
� 按变量梯度法构造李雅普诺夫函数方法的步骤如下。 1) 将李雅普诺夫函数V(x)的梯度假设为
⎡ a11 x1 + a12 x2 + ⋯ + a1n xn ⎤ ⎢a x + a x + ⋯ + a x ⎥ 22 2 2n n ⎥ grad V = ⎢ 21 1 ⎢ ⎥ ⋮ ⎢ ⎥ a x + a x + ⋯ + a x ⎣ n1 1 2n 2 nn n ⎦
非线性系统的李雅普诺夫稳定性分析(2/4)
� 本节主要研究Lyapunov方法在非线性系统中的应用。 � 由于非线性系统千差万别,没有统一的描述,目前也不存在 统一的动力学分析方法,因此对其进行稳定性分析是困难 的。 � 对于非线性系统,李雅普诺夫第二法虽然可应用于非线性 系统的稳定性判定,但其只是一个充分条件,并没有给出建 立李雅普诺夫函数的一般方法。 � 而只能针对具体的非线性系统进行具体分析。
̇1 = x2 ⎧x ⎨ ̇2 = − x2 − x13 ⎩x

第五章 控制系统的李雅普诺夫稳定性分析汇总

第五章 控制系统的李雅普诺夫稳定性分析汇总
几何意义:从 S(d ) 发出的轨迹, 在t
t 0 的任何时刻
总不会超出 S ( )
定义中对 、d 的大小没有具体要求,只要是有限的实数就可 以。因此,若状态解是等幅振荡的自由运动,在经典理论中 是不稳定的,而在李雅普诺夫的稳定性定义中是稳定的。
1、一致稳定性 如果平衡状态是稳定的,且 d 与 t 0 无关,则称该平衡状态是 一致稳定的。 因此,若定常系统的平衡状态是稳定的,则一定是一致稳定的。
x xe ( x1 xe1 ) 2 ... ( xn xen ) 2
2
2
2
由范数的定义可知,向量 ( x xe ) 的范数可写成
通常又将 x xe 称为 围之内时,则记为
x 与 xe 的距离。当向量 ( x xe ) 的范数限定在某一范
x xe
0
0 xe2 1 0 xe3 1 0 xe1 0
因此该系统有3个平衡状态:
二、范数的概念 范数:衡量(度量)状态空间距离的大小向量x的长度称 为向量x的范数: n维状态空间中,向量x的长度称为向量x的范数,用 x 表示,则 1
x x1 x2 ... xn ( xT x) 2
xe
与经典控制理论的区别: 1. 2. 3. 4. 5. 6. 平衡点/BIBO; 状态稳定/输出稳定; 经典控制的稳定大致对应于现代控制的渐进稳定; 即便输出稳定,状态可能不稳定; 李雅普诺夫意义下的稳定在经典中是不稳定的; 经典控制不需要一致性、全局性概念。
5.2 李雅普诺夫稳定性理论 一、李雅普诺夫第一方法 李雅普诺夫第一法的基本思想是利用状态方程解的性质来 判断系统的稳定性。通常又称为间接法。它适用于线性定常系 统以及线性时变系统和非线性系统可以线性化的情况。

《现代控制理论》李雅普诺夫稳定性分析

《现代控制理论》李雅普诺夫稳定性分析
向量和矩阵的范数
1、向量空间上的欧几里德范数(即向量长度)
其欧几里德范数定义为:
一般
一、向量和矩阵的范数
预备知识
矩阵范数
矩阵 的范数定义为:
【例】
Hale Waihona Puke , 则即:矩阵每个元素平方和开根号
预备知识
2、矩阵范数
1.二次型函数:由n个变量
组成的二次齐次多项式,称(n元)二次型函数
2.二次型函数的矩阵表示
则系统在原点处的平衡状态是不稳定的。
为唯一的平衡状态。
定理4:设系统状态方程为
李雅普诺夫主要的稳定性定理
例题
[例] 设系统状态方程为
试确定系统的稳定性。
解 xe=0
,
是该系统惟一的平衡状态。
由于当

,所以系统在原点处的平衡状态是
大范围渐近稳定的。
选取
李雅普诺夫主要的稳定性定理
例题
[例] 已知定常系统状态方程为
定义:若所有有界输入引起的零状态响应输出有界,则称系统为有界输入输出稳定。
李雅普诺夫第一方法—间接法
定理3:连续定常系统 传递函数为: 系统 BIBO 稳定的充要条件为:传递函数的所有极点均位于S左半平面。
【例】试分析系统渐近稳定和BIBO稳定。
李雅普诺夫主要的稳定性定理
讨论续
这是一个矛盾的结果,表明
也不是系统的
受扰运动解。综合以上分析可知,

时,显然有
根据定理9-12可判定系统的原点平衡状态是大范围渐近稳定的。
李雅普诺夫主要的稳定性定理
线性系统稳定性分析
一.线性定常系统李雅普诺夫稳定性分析
线性定常连续系统
系统状态方程为

李雅普诺夫稳定性分析

李雅普诺夫稳定性分析

李雅普诺夫稳定性分析常微分⼤作业--李雅普诺夫稳定性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 xe e (2)的状态x e 称为平衡状态(⼜称为平衡点)。

平衡状态的各分量不再随时间变化。

若已知状态⽅程,令0=x所求得的解x ,便是平衡状态。

对于线性定常系统Ax x= ,其平衡状态满⾜0=e Ax ,如果A ⾮奇异,系统只有惟⼀的零解,即存在⼀个位于状态空间原点的平衡状态。

⾄于⾮线性系统,0),(=t x f e 的解可能有多个,由系统状态⽅程决定。

李雅普诺夫稳定性分析

李雅普诺夫稳定性分析
第四章
李雅普诺夫稳定性分析
李雅普诺夫稳定性理论
李雅普诺夫理论在建立一系列关于稳定性概念的基础上,提出了判断 系统稳定性的两种方法: 间接法:利用线性系统微分方程的解来判断系统稳定性,又称之为李 雅普诺夫第一法; 直接法:首先利用经验和技巧来构造李雅普诺夫函数,进而利用李雅 普诺夫函数来判断系统稳定性,又称为李雅普诺夫第二法。
这表明, 当且仅当‖eAt‖≤ k <∞ 时,对任给的一个实数ε > 0,都对应存在和初始时 刻无关的一个实数 δ(ε)= ε /k,使得由满足不等式 ||x0 — xe|| ≤ δ(ε) (4-391) 的任一初态x0出发的受扰运动都满足不等式 xt; x0 ,0 xe e At x0 xe k , t 0 (4 392)
2)
证明 1) 设 xe 为线性定常系统(4-388+)的平衡状态,则由性质 e 0 和 Axe 0 x 可知,对于所有 t≥0 均有(可通过等式两边求微分证明下式)
xe e At xe (4 389) (4 390)
于是,考虑到 x(t; x0, 0) = eAtx0,有
x(t; x0 ,0) xe e At ( x0 xe ), t 0
2 李雅普诺夫意义下的稳定性
设系统初始状态位于以平衡状态xe为球心、δ为半径的闭球域S(δ)内,即 ||x0 - xe|| ≤ δ, t =t0 (4-385) 若能使系统方程的解x(t;x0,t0)在t→∞的过程中,都位于以xe为球心、任意规 定的半径为ε的闭球域S(ε)内,即 ||x(t;x0,t0)-xe|| ≤ ε,t≥t0 (4-386) 则称系统的平衡状态xe在李雅普诺夫意义下是稳定的。式中||· ||为欧几里德范 数,其几何意义是空间距离的尺度。 例如: ||x0 - xe||表示状态空间中, x0 点至 xe 点之间距离的尺度,数学表达式 为: ||x0 - xe|| = [(x10 – x1e)2+ (x20 – x2e)2+… +(xn0 – xne)2]1/2 (4-385)

李雅普诺夫稳定性分析

李雅普诺夫稳定性分析

第5章李雅普诺夫稳定性分析本章讨论李雅普诺夫稳定性分析。

主要介绍李雅普诺夫稳定性的定义以及分析系统状态稳定性的李雅普诺夫理论和方法;着重讨论李雅普诺夫第二法及其在线性系统和3类非线性系统的应用、李雅普诺夫函数的构造、李亚普诺夫代数(或微分)方程的求解等。

最后介绍李亚普诺夫稳定性问题的Matlab 计算与程序设计。

一个自动控制系统要能正常工作,必须首先是一个稳定的系统,即当系统受到外界干扰时它的平衡被破坏,但在外界干扰去掉以后,它仍有能力自动地恢复在平衡态下继续工作。

系统的这种性能,叫做稳定性。

例如,电压自动调解系统中保持电机电压为恒定的能力、电机自动调速系统中保持电机转速为一定的能力以及火箭飞行中保持航向为一定的能力等。

具有稳定性的系统称为稳定系统,不具有稳定性的系统称为不稳定系统。

也可以说,系统的稳定性就是系统在受到外界干扰后,系统状态变量或输出变量的偏差量(被调量偏离平衡位置的数值)过渡过程的收敛性,用数学方法表示就是ε≤Δ∞→)(Lim t x t 式中,)(t x Δ为系统被调量偏离其平衡位置的变化量;ε为任意小的规定量。

如果系统在受到外扰后偏差量越来越大,显然它不可能是一个稳定系统。

在经典控制理论中,借助于常微分方程稳定性理论,产生了许多线性定常系统的稳定性判据,如劳斯-胡尔维茨(Routh-Hurwitz)判据和奈奎斯特判据等,都给出了既实用又方便的稳定性判别及设计方法。

但这些稳定性判据仅限于讨论SISO 线性定常系统输入输出间动态关系,讨论的是有界输入有界输出(BIBO)稳定性,未研究系统的内部状态变化的稳定性。

再则,对于非线性或时变系统,虽然通过一些系统转化方法,上述稳定判据尚能在某些特定系统和范围内应用,但是难以胜任一般系统。

现代控制系统的结构比较复杂,大都存在非线性或时变因素,即使是系统结构本身, 往往也需要根据性能指标的要求而加以改变,才能适应新的情况,保证系统的正常或最佳运行状态。

第5章系统的稳定性

第5章系统的稳定性

经典控制论中,系统稳定性判据

代数判据


Routh(劳斯)判据 Hurwitz(古尔维茨)判据 Nyquist判据 Bode判据

几何判据

5.2 Routh(劳斯)稳定判据
Routh稳定判据
不求解特征方程的根,直接根据特征方程的系 数,判断系统的稳定性,回避了求解高次方程根 的困难。

【例】D(s) s 4 3s3 4s 2 12s 16
【解】:Routh表为: s4 s3 s2 s1 s0 1 3 4 16 12
12 48 48
0( ) 16 12 48 0

很小时,

12

0
16
【结论】:系统不稳定,并有两个正实部根。
【情况2】:
n n n an1 an2 an3 si, si s j, an an an i 1 i j i j k i 1, j 2
(1)
n
s
i 1
n
i
si s j sk,
i 1, j 2, k 3
n a0 n , (1) si an i 1
系统稳定的必要条件:特征方程中所有项的系数均大 于0,只要有一项等于或小于0,则为不稳定系统。 充分必要条件:Routh表第一列元素均大于0。

必要条件证明
D(s) an s n an1s n1 an1 n1 得:s s an 再展开,得
n
a1s a0 0,两端同除以an,并分解因式, (s sn )
其中N+为:正穿越与半次正穿越次数的和。 其中N-为:负穿越与半次负穿越次数的和。

李雅普洛夫稳定性分析

李雅普洛夫稳定性分析
G ( s ) C ( s A I ) 1 B 0 1 s 1 s 6 1 1 1 2 ( s ( s 2 ) s 2 ) 3 ( ) ( s 1 3 )
极点位于s左半平面,s=2的极点被对消掉了。 系统是有界输入有界输出稳定的。
(2)
求系统的特征方程:
4)(半)正定和(半)负定间的关系
V(x)为正定,则-V(x)为负定; V(x)为半正定,则-V(x)为半负定;
5)不定性 如果无论取多么小的零点的某个邻域,V(X)可为正 值也可为负值,则V(x)为不定。
4.3 二次型标量函数XTPX 1、二次型函数V(x):
p11 p12 p1nx1
V(x)[x1,x2, ,xn]p21
如果 1 p 1 10 , 2p p 1 21 1p p 1 2 2 20 , ,n P 0
则P为正定,即V(x)正定。
2)二次型 V(x)xTPx为负定,或实对称阵P为负定的
充要条件是P的主子行列式满足 i 0(i为奇数;)i 0
( i为偶数)i=1,2,3,…,n。
4.4 稳定性判据
判据1:设系统的状态方程为 xf(x)
dI e A t ) ( 1 6 1 ( 2 )( 3 ) 0
求 得 12 , : 2 3
系统不是渐近稳定的。
3.2 非线性系统的李亚普洛夫第一法
对非线性系统 Xf(X,t)
当f (X,t)为与X 同维的矢量函数,且对X 具有连续偏导数,则可将
其作线性化处理,在Xe 领域内展开成泰勒级数 X X fT(XX e)R (X ) 。
A为非奇异阵时,x=0是其唯一的平衡状态。 A为奇异阵时,系统有无穷多个平衡状态。 2、对于非线性系统,有一个或多个平衡状态。 3、对任意 xe 0 ,总可经过一定的坐标变换,把它化到坐标原 点(即零状态)。一般将平衡状态取为状态空间原点。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第六章 李雅普诺夫稳定性分析在反馈控制系统的分析设计中,系统的稳定性是首先需要考虑的问题之一。

因为它关系到系统是否能正常工作。

经典控制理论中已经建立了劳斯判据、Huiwitz 稳定判据、Nquist 判据、对数判据、根轨迹判据等来判断线性定常系统的稳定性,但不适用于非线性和时变系统。

分析非线性系统稳定性及自振的描述函数法,则要求系统的线性部分具有良好的滤除谐波的性能;而相平面法则只适合于一阶、二阶非线性系统。

1892年俄国学者李雅普诺夫(Lyapunov )提出的稳定性理论是确定系统稳定性的更一般的理论,它采用状态向量来描述,不仅适用于单变量、线性、定常系统,还适用于多变量、非线性、时变系统。

§6-1 外部稳定性和内部稳定性系统的数学模型有输入输出描述(即外部描述)和状态空间描述(即内部描述),相应的稳定性便分为外部稳定性和内部稳定性。

一、外部稳定性1、定义(外部稳定性):若系统对所有有界输入引起的零状态响应的输出是有界的,则称该系统是外部稳定的。

(外部稳定性也称为BIBO (Bounded Input Bounded Output )稳定性) 说明:(1)所谓有界是指如果一个函数)(t h ,在时间区间],0[∞中,它的幅值不会增至无穷,即存在一个实常数k ,使得对于所有的[]∞∈0t ,恒有∞<≤k t h )(成立。

(2)所谓零状态响应,是指零初始状态时非零输入引起的响应。

2、系统外部稳定性判据线性定常连续系统∑),,(C B A 的传递函数矩阵为Cxy Bu Ax x=+=BUA sI X BU X A sI CX Y BU AX sX 1)()(--==-=+=B A sIC s G 1)()(--=当且仅当)(s G 极点都在s 的左半平面内时,系统才是外部稳定(或BIBO 稳定)的。

【例6.1.1】已知受控系统状态空间表达式为u x x ⎥⎦⎤⎢⎣⎡-+⎥⎦⎤⎢⎣⎡-=121160 , []x y 10= 试分析系统的外部稳定性。

解:系统为SISO 系统,传递函数为B A sIC s G 1)()(--=[]⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡+--=-12116101s s )3)(2(2+--=s s s31+=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.2】已知受控系统状态空间表达式为u x x ⎥⎦⎤⎢⎣⎡-+⎥⎦⎤⎢⎣⎡-=121160 , []x y 10= 试分析系统的内部稳定性。

解:该系统为线性定常系统,其特征方程为:0)3)(2(6)1(=+-=-+=-λλλλλA I于是系统的特征值为21=λ,32-=λ,故系统不是内部稳定(渐近稳定)的。

三、内部稳定性与外部稳定性的关系1、若系统是内部稳定(渐近稳定)的,则一定是外部稳定(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 xe 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 点之间的距离的尺度,其数学表达式为2021100)()(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),;(lim00→-∞→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 内正定。

如2221)(x x x V +=是正定的。

0x - 初始状态图6-3 二维空间不稳定的几何解释示意图(2)负定性标量函数)(x V 在域s 中对所有非零状态(0≠x )有0)(<x V 且0)0(=V ,则称)(x V 在域s 内负定。

如)()(2221x x x V +-=是正定的。

相关文档
最新文档