李雅普诺夫稳定性分析
一些概念:
- 能量函数:是一个大于等于0的量,在数学上可以用一个二次型函数来描述,形如 V ( x ) = x T P x V(x)={{x}^{T}}Px V(x)=xTPx,其中 x x x为一系列状态量或误差量
- 系统稳定:可以定义对于一个非线性系统 x ˙ ( t ) = f ( x ( t ) , t ) \dot{x}(t)=f(x(t),t) x˙(t)=f(x(t),t),若存在能量函数 V ( x ) V(x) V(x)满足条件 V ( x ) V(x) V(x)是正定的且 V ˙ ( x ) \dot{V}(x) V˙(x)是半负定的,则系统在原点处的平衡状态是李雅普诺夫意义下稳定的
解题步骤
- 求解平衡点
- 凑出对称矩阵P中各元素关系
- 凑出能量函数形式
- 反向求解证明能量函数导数半负定,以证明系统稳定性
例题
1.考虑系统
{
x
˙
1
=
x
2
x
˙
2
=
−
x
1
−
(
1
+
x
2
2
)
3
x
2
5
\left.\left\{\begin{aligned}\dot{x}_1&=x_2\\\dot{x}_2&=-x_1-(1+x_2^2)^3x_2^5\end{aligned}\right.\right.
{x˙1x˙2=x2=−x1−(1+x22)3x25
试采用李雅普诺夫第二法判断系统关于平衡点
(
x
1
=
0
,
x
2
=
0
)
(x_1=0,x_2=0)
(x1=0,x2=0)的稳定性。
解:
草稿运算
设李雅普诺夫函数 V ( x ) = x T P x V(x)={{x}^{T}}Px V(x)=xTPx,求导则有 V ˙ ( x ) = 2 P 11 x ˙ 1 x 1 + 2 P 22 x ˙ 2 x 2 + ( P 12 + P 21 ) x ˙ 1 x 2 + + ( P 12 + P 21 ) x ˙ 2 x 1 \dot{V}(x)=2P_{11}\dot{x}_1x_1+2P_{22}\dot{x}_2x_2+(P_{12}+P_{21})\dot{x}_1x_2++(P_{12}+P_{21})\dot{x}_2x_1 V˙(x)=2P11x˙1x1+2P22x˙2x2+(P12+P21)x˙1x2++(P12+P21)x˙2x1
带入题设条件,则有
V ˙ ( x ) = 2 P 11 x 2 x 1 + 2 P 22 ( − x 1 − ( 1 + x 2 2 ) 3 x 2 5 ) x 2 + ( P 12 + P 21 ) x 2 x 2 + ( P 12 + P 21 ) ( − x 1 − ( 1 + x 2 2 ) 3 x 2 5 ) x 1 \dot{V}(x)=2P_{11}x_2x_1+2P_{22}(-x_1-(1+x_2^2)^3x_2^5)x_2+(P_{12}+P_{21})x_2x_2+(P_{12}+P_{21})(-x_1-(1+x_2^2)^3x_2^5)x_1 V˙(x)=2P11x2x1+2P22(−x1−(1+x22)3x25)x2+(P12+P21)x2x2+(P12+P21)(−x1−(1+x22)3x25)x1
化简有
V ˙ ( x ) = 2 ( P 11 − P 22 ) x 1 x 2 − 2 P 22 ( 1 + x 2 2 ) 3 x 2 6 + ( P 12 + P 21 ) ( − x 1 − ( 1 + x 2 2 ) 3 x 2 5 ) x 1 \dot{V}(x)=2(P_{11}-P_{22})x_1x_2-2P_{22}(1+x_2^2)^3x_2^6+(P_{12}+P_{21})(-x_1-(1+x_2^2)^3x_2^5)x_1 V˙(x)=2(P11−P22)x1x2−2P22(1+x22)3x26+(P12+P21)(−x1−(1+x22)3x25)x1
则当 P 11 = P 22 > 0 P_{11}=P_{22}>0 P11=P22>0且 P 12 = P 21 = 0 P_{12}=P_{21}=0 P12=P21=0时
V ˙ ( x ) = − 2 P 22 ( 1 + x 2 2 ) 3 x 2 6 < = 0 \dot{V}(x)=-2P_{22}(1+x_2^2)^3x_2^6<=0 V˙(x)=−2P22(1+x22)3x26<=0
设李雅普诺夫函数
V
(
x
)
=
x
T
P
x
V(x)={{x}^{T}}Px
V(x)=xTPx,其中
P
=
[
1
2
0
0
1
2
]
P=\begin{bmatrix}\dfrac{1}{2}&0\\0&\dfrac{1}{2}\end{bmatrix}
P=
210021
,对
V
V
V求导,则有
V
˙
(
x
)
=
−
P
22
(
1
+
x
2
2
)
3
x
2
6
<
=
0
\dot{V}(x)=-P_{22}(1+x_2^2)^3x_2^6<=0
V˙(x)=−P22(1+x22)3x26<=0
又P为正定矩阵,则有
V
(
x
)
>
0
V(x)>0
V(x)>0
故系统在原点附近稳定
2.考虑系统
m
x
¨
+
c
x
˙
+
sin
(
x
)
+
x
=
0
m\ddot{x}+c\dot{x}+\sin(x)+x=0
mx¨+cx˙+sin(x)+x=0
其中
m
,
c
m,c
m,c均为正常数。
(1)求系统平衡点;
(2)试应用李雅普诺夫第二法分析系统关于平衡点的稳定性
解:
(1) 选取
x
1
=
x
,
x
2
=
x
1
˙
=
x
˙
x_1=x,x_2=\dot{x_1}=\dot{x}
x1=x,x2=x1˙=x˙
则有状态方程
{
x
˙
1
=
x
2
x
˙
2
=
−
1
m
(
c
x
2
+
s
i
n
x
1
+
x
1
)
\left.\left\{\begin{aligned}\dot{x}_1&=x_2\\\dot{x}_2&=-\frac1{m}(cx_2+sinx_1+x_1)\end{aligned}\right.\right.
⎩
⎨
⎧x˙1x˙2=x2=−m1(cx2+sinx1+x1)
令
x
˙
1
=
x
˙
2
=
0
\dot x_1=\dot x_2=0
x˙1=x˙2=0,则有
x
1
=
x
2
=
0
x_1=x_2=0
x1=x2=0
故平衡点为原点
(2)
草稿运算
设李雅普诺夫函数 V ( x ) = x T P x V(x)={{x}^{T}}Px V(x)=xTPx,求导则有 V ˙ ( x ) = 2 P 11 x ˙ 1 x 1 + 2 P 22 x ˙ 2 x 2 + ( P 12 + P 21 ) x ˙ 1 x 2 + + ( P 12 + P 21 ) x ˙ 2 x 1 \dot{V}(x)=2P_{11}\dot{x}_1x_1+2P_{22}\dot{x}_2x_2+(P_{12}+P_{21})\dot{x}_1x_2++(P_{12}+P_{21})\dot{x}_2x_1 V˙(x)=2P11x˙1x1+2P22x˙2x2+(P12+P21)x˙1x2++(P12+P21)x˙2x1
带入题设条件,则有
V ˙ ( x ) = 2 P 11 x 2 x 1 + 2 P 22 ( − 1 m ( c x 2 + s i n x 1 + x 1 ) ) x 2 + ( P 12 + P 21 ) x 2 x 2 + ( P 12 + P 21 ) ( − 1 m ( c x 2 + s i n x 1 + x 1 ) ) x 1 \dot{V}(x)=2P_{11}x_2x_1+2P_{22}(-\frac1{m}(cx_2+sinx_1+x_1))x_2+(P_{12}+P_{21})x_2x_2+(P_{12}+P_{21})(-\frac1{m}(cx_2+sinx_1+x_1))x_1 V˙(x)=2P11x2x1+2P22(−m1(cx2+sinx1+x1))x2+(P12+P21)x2x2+(P12+P21)(−m1(cx2+sinx1+x1))x1
P 12 = P 21 = 0 P_{12}=P_{21}=0 P12=P21=0时,化简有
V ˙ ( x ) = 2 ( P 11 − 1 m P 22 ) x 1 x 2 − 2 P 22 ( 1 m c x 2 2 ) − 2 P 22 1 m x 2 s i n x 1 \dot{V}(x)=2(P_{11}-\frac1{m}P_{22})x_1x_2-2P_{22}(\frac1{m}cx_2^2)-2P_{22}\frac1{m}x_2sinx_1 V˙(x)=2(P11−m1P22)x1x2−2P22(m1cx22)−2P22m1x2sinx1
当 P 11 = 1 m P 22 > 0 P_{11}=\frac1{m}P_{22}>0 P11=m1P22>0时
V ˙ ( x ) = − 2 P 22 ( 1 m c x 2 2 ) − 2 P 22 1 m x 2 s i n x 1 = − 2 P 22 ( 1 m c x 2 2 ) + 2 P 22 1 m c o s , ( x 1 ) \dot{V}(x)=-2P_{22}(\frac1{m}cx_2^2)-2P_{22}\frac1{m}x_2sinx_1=-2P_{22}(\frac1{m}cx_2^2)+2P_{22}\frac1{m}cos^,(x_1) V˙(x)=−2P22(m1cx22)−2P22m1x2sinx1=−2P22(m1cx22)+2P22m1cos,(x1)
可以看到多了一项+2P_{22}\frac1{m}cos^,(x_1),故需要凑李雅普诺夫函数来消掉该项
设李雅普诺夫函数
V
(
x
)
=
x
T
P
x
−
c
o
s
x
1
+
1
V(x)={{x}^{T}}Px-cosx_1+1
V(x)=xTPx−cosx1+1,其中
P
=
[
1
2
0
0
1
2
m
]
P=\begin{bmatrix}\dfrac{1}{2}&0\\0&\dfrac{1}{2}m\end{bmatrix}
P=
210021m
,
对
V
V
V求导,则有
V
˙
(
x
)
=
−
2
P
22
(
1
m
c
x
2
2
)
+
2
P
22
1
m
c
o
s
,
(
x
1
)
−
x
2
s
i
n
x
1
=
−
c
x
2
2
<
=
0
\dot{V}(x)=-2P_{22}(\frac1{m}cx_2^2)+2P_{22}\frac1{m}cos^,(x_1)-x_2sinx_1=-cx_2^2<=0
V˙(x)=−2P22(m1cx22)+2P22m1cos,(x1)−x2sinx1=−cx22<=0
又
V
(
x
)
=
x
T
P
x
−
c
o
s
x
1
+
1
>
=
0
V(x)={{x}^{T}}Px-cosx_1+1>=0
V(x)=xTPx−cosx1+1>=0
故系统在平衡点附近稳定
习题
习题1.考虑系统
x
¨
=
−
x
˙
(
c
o
s
x
)
2
−
k
x
\ddot{x}=-\dot{x}(cosx)^2-kx
x¨=−x˙(cosx)2−kx
其中
k
>
0
k>0
k>0为常数。试采用李雅普诺夫第二法判断系统关于平衡点
x
=
x
˙
=
0
x=\dot x=0
x=x˙=0的稳定性。
习题2.考虑系统
x
¨
+
2
x
˙
+
x
1
3
=
0
\ddot{x}+2\dot x+x^{\frac1{3}}=0
x¨+2x˙+x31=0
(1)求系统平衡点;
(2)试应用李雅普诺夫第二法分析系统关于平衡点的稳定性
习题3.考虑系统
x
¨
=
−
x
˙
−
c
(
x
−
1
)
1
3
\ddot{x}=-\dot{x}-c(x-1)^{\frac1{3}}
x¨=−x˙−c(x−1)31,其中c为正常数
(1)求系统平衡点;
(2)试应用李雅普诺夫第二法分析系统关于平衡点的稳定性
习题4.考虑系统
m
x
¨
+
x
˙
(
c
o
s
x
)
2
+
k
(
s
i
n
x
+
x
)
=
0
m\ddot{x}+\dot{x}(cosx)^2+k(sinx+x)=0
mx¨+x˙(cosx)2+k(sinx+x)=0,其中m,k为正常数
(1)求系统平衡点;
(2)试应用李雅普诺夫第二法分析系统关于平衡点的稳定性
欢迎讨论,如有错误敬请指出