李雅普诺夫第二法求解思路及例题

本文详细介绍了李雅普诺夫稳定性分析的概念,通过例题展示了如何使用李雅普诺夫第二法判断非线性系统在特定平衡点的稳定性,包括求解能量函数、导数的半负定性等步骤。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

李雅普诺夫稳定性分析

一些概念:

  • 能量函数:是一个大于等于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)是半负定的,则系统在原点处的平衡状态是李雅普诺夫意义下稳定的

解题步骤

  1. 求解平衡点
  2. 凑出对称矩阵P中各元素关系
  3. 凑出能量函数形式
  4. 反向求解证明能量函数导数半负定,以证明系统稳定性

例题

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(P11P22)x1x22P22(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(P11m1P22)x1x22P22(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)=xTPxcosx1+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)=xTPxcosx1+1>=0
故系统在平衡点附近稳定

习题

习题1.考虑系统

x ¨ = − x ˙ ( c o s x ) 2 − k x \ddot{x}=-\dot{x}(cosx)^2-kx x¨=x˙(cosx)2kx
其中 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(x1)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)试应用李雅普诺夫第二法分析系统关于平衡点的稳定性

欢迎讨论,如有错误敬请指出

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值