线性定常系统的稳定性分析:研究线性定常系统的渐近稳定性的方法有很多。例如,对于用如下方程描述的连续时间系统:
x
˙
=
A
x
{\dot{x}}{=}Ax
x˙=Ax
在原点处渐近稳定性的充分必要条件是:
A
A
A的所有特征值都有负实部,或者说特征多项式:
∣
s
??
I
?
????
A
??
∣
=
??
s
n
??
+
??
a
1
s
n
?
1
??
+
??
?
??
+
??
a
n
?
1
s
??
+
??
a
n
\mid s\!I-\!\!A\!\mid=\!s^{n}\!+\!a_{1}s^{n-1}\!+\!\cdots\!+\!a_{n-1}s\!+\!a_{n}
∣sI?A∣=sn+a1?sn?1+?+an?1?s+an?
的所有零点都有负实部。
同样,对于用如下方程描述的离散时间系统:
x
(
?
k
??
+
??
1
)
=
G
x
(
?
k
?
)
{{x}}(\,k\!+\!1)={G}{{x}}(\,k\,)
x(k+1)=Gx(k)
在原点处渐近稳定性的充分必要条件是:
G
G
G的所有特征值的模都小于
1
1
1,或者说特征多项式
∣
?
z
I
?
G
∣
=
z
n
+
a
1
z
n
?
1
+
?
+
a
n
?
1
z
+
a
n
|\,z I-G|=z^{n}+a_{1}z^{n-1}+\cdots+a_{n-1}z+a_{n}
∣zI?G∣=zn+a1?zn?1+?+an?1?z+an?
的零点都在平面上以原点为中心的单位圆内。
然而,在高阶系统中或特征多项式的系数不是数值的情况下,特征根的求解会比较困难。在这种情况下,应采用朱利稳定性判据或者劳斯-赫维茨稳定性判据。
对于离散时间系统,使用前向差分 V ( ? x ( ? k ?? + ?? 1 ) ? T ) ?? ? ?? V ( ? x ( ? k ??? T ) ? ) V(\,{x}(\,k\!+\!1)\,T)\!-\!V(\,{x}(\,k\!\,T)\,) V(x(k+1)T)?V(x(kT))代替 V ˙ ( x ) \dot V(x) V˙(x)。
定理:
离散时间系统表示如下:
x
(
?
(
k
+
1
)
?
T
)
=
f
(
x
(
k
T
)
?
)
{{x}}(\,({{k}}+1)\,T)={{f}}({{x}}({{k}}T)\,)
x((k+1)T)=f(x(kT))
其中 x x x: n n n维向量
f ( x ) f(x) f(x): n n n维向量,且 f ( 0 ) = 0 f(0)=0 f(0)=0。
T T T:采样周期。
假定存在一个在 x x x处连续的标量函数 V ( x ) V(x) V(x),使得
1. V ( x ) > 0 ( x ≠ 0 ) V(x)>0(x \neq 0) V(x)>0(x=0)。
2. Δ V ( x ) < 0 ( x ≠ 0 ) \Delta V(x)<0(x \neq 0) ΔV(x)<0(x=0),其中 Δ V ( x ( k T ) ) ? = ?? V ( x ( k ?? + ?? 1 ) T ) ?? ? ?? V ( x ( k T ) ) ?? = ???? V ( f ( x ( k T ) ) ) ?? ? ?? V ( x ( k T ) ) \Delta V({{x}}(k T))\,=\!V({{x}}(k\!+\!1)T)\!-\!V({{x}}(k T))\!=\!\!V(f({{x}}(k T)))\!-\!V({{x}}(k T)) ΔV(x(kT))=V(x(k+1)T)?V(x(kT))=V(f(x(kT)))?V(x(kT))
3.$ V(0)=0$。
4.当 ∣ x ∣ → ∞ |x| \rightarrow \infty ∣x∣→∞时, V ( x ) → ∞ V(x)\rightarrow \infty V(x)→∞。那么平衡状态 x = 0 x=0 x=0是全局渐近稳定的,且 V ( x ) V(x) V(x)是一个李雅普诺夫函数。
注意,此定理中,条件2可以替换为:
2’.对于任意 x x x, Δ V ( x ) ≤ 0 \Delta V(x) \leq 0 ΔV(x)≤0,且对于满足系统的任意解序列 x ( k T ) {x(kT)} x(kT), Δ V ( x ) \Delta V(x) ΔV(x)都不恒等于零。
这意味着,如果对差分方程的任意解序列, Δ V ( x ) \Delta V(x) ΔV(x)都不恒等于零,那么就不必要求 Δ V ( x ) \Delta V(x) ΔV(x)一定是负定的。
离散时间系统描述如下:
x
(
k
+
1
)
=
G
x
(
k
)
x(k+1)=Gx(k)
x(k+1)=Gx(k)
其中,
x
x
x是状态向量(
n
n
n维向量),
G
G
G是
n
×
n
n×n
n×n非奇异常系数矩阵。原点
x
=
0
x=0
x=0是平衡状态。下面将采用李雅普诺夫第二方法来研究该状态的稳定性。选取一个可能的李雅普诺夫函数为
V
(
x
(
k
)
)
=
x
?
(
k
)
P
x
(
k
)
V(x(k))=x^{*}(k)Px(k)
V(x(k))=x?(k)Px(k)
其中,
P
P
P是一个正定的埃尔米特(或正定实对称)矩阵。那么
Δ
V
(
x
(
k
)
?
)
=
????
V
(
x
(
k
??
+
??
1
)
?
)
??
?
??
V
(
x
(
k
)
?
)
=
????
x
?
?
(
k
??
+
??
1
)
P
x
(
k
??
+
??
1
)
??
?
??
x
?
?
(
k
)
P
x
(
x
)
=
????
[
G
x
(
k
)
?
]
?
?
P
[
G
x
(
k
)
?
]
??
?
??
x
?
?
(
k
)
P
x
(
k
)
=
????
x
?
?
(
k
)
G
?
?
P
G
x
(
k
)
??
?
??
x
?
?
(
k
)
P
x
(
k
)
=
????
x
?
?
(
k
)
?
(
G
?
?
P
G
??
?
??
P
)
?
x
(
k
)
=
?
x
?
?
(
k
)
Q
x
(
k
)
\begin{array}{l l}{\Delta V({{x}}({{k}})\,)=\!\!V({{x}}({{k}}\!+\!1)\,)\!-\!V({{x}}({{k}})\,)}\\ {=\!\!{{x}}^{*}\ ({{k}}\!+\!1){{P}}{{x}}({{k}}\!+\!1)\!-\!{{x}}^{*}\ ({{k}}){{P}}{{x}}({{x}})}\\ {=\!\![{{G}}{{x}}({{k}})\,]^{*}\,{{P}}[{{G}}{{x}}({{k}})\,]\!-\!{{x}}^{*}\ ({{k}}){{P}}{{x}}({{k}})}\\ {=\!\!{{x}}^{*}\ ({{k}}){{G}}^{*}\,{{P}}{{G}}{{x}}({{k}})\!-\!{{x}}^{*}\ ({{k}}){{P}}{{x}}({{k}})}\\ {=\!\!{{x}}^{*}\ ({{k}})\,({{G}}^{*}\,{{P}}{{G}}\!-\!{{P}})\,{{x}}({{k}})}\\ =-{{x}}^{*}\ ({{k}})Q{{x}}({{k}}) \end{array}
ΔV(x(k))=V(x(k+1))?V(x(k))=x??(k+1)Px(k+1)?x??(k)Px(x)=[Gx(k)]?P[Gx(k)]?x??(k)Px(k)=x??(k)G?PGx(k)?x??(k)Px(k)=x??(k)(G?PG?P)x(k)=?x??(k)Qx(k)?
其中
Q
=
?
(
G
?
?
P
G
??
?
??
P
)
Q=-({{G}}^{*}\,{{P}}{{G}}\!-\!{{P}})
Q=?(G?PG?P),只要¥正定,系统就是渐进稳定的。
定理:
离散时间系统定义如下:
x
(
?
k
??
+
??
1
?
)
=
??
G
x
(
k
)
{{x}}(\,k\!+\!1\,)=\!G{{x}}({{k}})
x(k+1)=Gx(k)
其中,
x
x
x是状态向量(
n
n
n维向量),
G
G
G是
n
×
n
n×n
n×n非奇异常系数矩阵。
平衡状态
x
=
0
x=0
x=0全局渐近稳定的充分必要条件是:给定任意的正定埃尔米特(或任意正定实对称)矩阵
Q
Q
Q,存在一个正定的埃尔米特(或正定实对称)矩阵
P
P
P,使得
G
?
?
P
G
?
P
=
?
Q
G^{\ast}\,P G{-}P{=}-Q{}
G?PG?P=?Q
标量函数
x
?
?
P
x
x^{\ast}\,P x
x?Px是该系统的一个李雅普诺夫函数。
如果 : Δ V ( x ( k ) ) = ? x ? ? ( k ) Q x ( k ) :\Delta V(x(k))=-x^{*}\ (k)Q x(k) :ΔV(x(k))=?x??(k)Qx(k)沿其任一解序列都不恒等于零,那么 Q Q Q可以选为半正定矩阵。
如果用状态空间方程来描述系统,则通过把连续时间系统离散化得到的离散时间系统,其平衡状态的渐近稳定性等同于原连续时间系统平衡状态的渐近稳定性。
连续时间系统表示如下:
x
˙
=
A
x
\dot {x}=Ax
x˙=Ax
其对应的离散时间系统为
x
(
(
k
+
1
)
T
)
=
G
x
(
k
T
)
x((k+1)T)=Gx(kT)
x((k+1)T)=Gx(kT)
其中
G
=
e
A
T
G=\mathrm{e}^{AT}
G=eAT
如果上述连续时间系统是渐近稳定的,即如果矩阵
A
A
A的所有特征值都有负实部,则IG"→0,当n→0时
且离散化后得到的系统也是渐近稳定的。
这是因为,如果 λ i \lambda_i λi?都是 A A A的特征值,那么 e λ i T \mathrm{e}^{\lambda_iT} eλi?T就都是 G G G的特征值(如果 λ i T \lambda_i T λi?T是负的,那么 ∣ e λ i T ∣ < 1 |\mathrm{e}^{\lambda_i T}|<1 ∣eλi?T∣<1)。
对于
x
≠
0
x≠0
x=0的某些集合和某些范数,如果函数
f
(
x
)
f(x)
f(x)满足
f
(
0
)
=
0
f(0)=0
f(0)=0,且
∣
∣
f
(
x
)
∣
∣
<
∣
∣
x
∣
∣
||f(x)|| <||x||
∣∣f(x)∣∣<∣∣x∣∣
那么函数
f
(
x
)
f(x)
f(x)称为压缩映射。
对于离散时间系统,范数 ∣ ∣ x ∣ ∣ ||x|| ∣∣x∣∣可用做李雅普诺夫函数。
离散时间系统表示如下: x ( k + 1 ) = f ( x ( k ) ) , f ( 0 ) = 0 x(k+1)=f(x(k)),f(0)=0 x(k+1)=f(x(k)),f(0)=0其中, x x x和 f ( x ) f(x) f(x)都是 n n n维向量。假设对所有的 x x x和某些范数, f ( x ) f(x) f(x)是压缩映射。那么上述系统的原点是广义渐近稳定的,且它的一个李雅普诺夫函数为 V ( x ) = ∣ ∣ x ∣ ∣ V(x)=||x|| V(x)=∣∣x∣∣这点可证明如下:首先V(x)=||x||是正定的。又因为 f ( x ) f(x) f(x)是对所有的 x x x和某些范数的压缩,有 Δ V ( x ( k ) ? ) ?? = ???? V ( f ( x ( k ) ? ) ? ) ?? ? ?? V ( x ( k ) ? ) ?? = ?? ∥ f ( x ) ∥ ???? ? ∥ x ∥ \Delta V({{x}}({{k}})\,)\!=\!\!V(f({{x}}({{k}})\,)\,)\!-\!V({{x}}({{k}})\,)\!=\!\parallel f({{x}})\parallel\!\!-\parallel{{x}}\parallel ΔV(x(k))=V(f(x(k)))?V(x(k))=∥f(x)∥?∥x∥是负定的。因此平衡点是渐进稳定的。