Skip to main content

李雅普诺夫直接法

李雅普诺夫直接法(Lyapunov's Direct Method,也称第二方法)是分析非线性系统稳定性的核心工具——无需求解方程,直接通过构造一个"能量函数"来判断稳定性。它特别适用于线性化方法失效的非双曲情形。

核心思想

为系统构造一个类似于"能量"的函数 V(X)V(X)。若 VV 沿系统轨线单调递减,则系统趋向不动点(稳定);若递增加,则远离(不稳定)。

稳定性定义

对自治系统 X=F(X)X' = F(X),设 X=0X^* = 0 是不动点(通过平移总能做到),考虑零解的稳定性:

类型定义
稳定(李雅普诺夫意义)ε>0,δ>0\forall \varepsilon > 0, \exists \delta > 0X(0)<δX(t)<ε,t0\|X(0)\| < \delta \Rightarrow \|X(t)\| < \varepsilon, \forall t \geq 0
渐近稳定稳定 + δ>0\exists \delta > 0X(0)<δlimtX(t)=0\|X(0)\| < \delta \Rightarrow \lim_{t \to \infty} X(t) = 0
指数稳定δ,α,C>0\exists \delta, \alpha, C > 0X(0)<δX(t)CeαtX(0)\|X(0)\| < \delta \Rightarrow \|X(t)\| \leq C e^{-\alpha t}\|X(0)\|
不稳定不满足稳定条件

李雅普诺夫函数

定义

V:URV: U \to \mathbb{R}C1C^1 函数,定义在原点邻域 UU 上。VV 沿系统 X=F(X)X' = F(X)全导数为:

V˙(X)=V(X)F(X)=i=1nVxifi(X)\dot{V}(X) = \nabla V(X) \cdot F(X) = \sum_{i=1}^{n} \frac{\partial V}{\partial x_i} f_i(X)

直观上,V˙(X)\dot{V}(X)VV 沿解轨线对时间的导数。

正定函数

V(X)V(X) 称为正定的,若 V(0)=0V(0) = 0V(X)>0V(X) > 0 对所有 X0X \neq 0(在原点邻域内)。

李雅普诺夫稳定性定理

定理一:稳定

若存在正定函数 V(X)V(X)C1C^1),使在原点某邻域内 V˙(X)0\dot{V}(X) \leq 0半负定),则零解是稳定的。

定理二:渐近稳定

若存在正定函数 V(X)V(X)C1C^1),使在原点某邻域内 V˙(X)<0\dot{V}(X) < 0 对所有 X0X \neq 0负定),则零解是渐近稳定的。

定理三:不稳定(Chetaev 定理)

若存在函数 V(X)V(X)V(0)=0V(0) = 0,且在原点任意小邻域内存在点使 V>0V > 0,且在这些点上 V˙>0\dot{V} > 0,则零解不稳定

判定表

V(X)V(X)V˙(X)\dot{V}(X)结论
正定负定渐近稳定
正定半负定稳定(需额外验证渐近性)
正定正定不稳定
正定不定不能直接判定

线性系统的李雅普诺夫方程

对于线性系统 X=AXX' = AX,可取二次型 V(X)=XPXV(X) = X^\top P XPP 对称正定)。则:

V˙(X)=X(AP+PA)X\dot{V}(X) = X^\top (A^\top P + PA) X

李雅普诺夫方程AP+PA=QA^\top P + PA = -Q。若对任意正定 QQ,存在正定解 PP,则 AA 的所有特征值实部为负(零解渐近稳定)。

逆命题也成立:若 AA 的特征值实部全负,则对任意正定 QQ,存在唯一正定解 PP

构造李雅普诺夫函数的常见技巧

1. 能量函数法

对力学系统,取总能量(动能 + 势能)作为 VV

2. 二次型试探

V(X)=XPXV(X) = X^\top P XPP 待定正定。

3. 变量梯度法

先假设 V=g(X)\nabla V = g(X)(待定),由可积性条件 gi/xj=gj/xi\partial g_i/\partial x_j = \partial g_j/\partial x_i 确定 gg,再积分得 VV

举例

分析 {x=yx3y=xy3\begin{cases} x' = -y - x^3 \\ y' = x - y^3 \end{cases} 零解的稳定性。

V(x,y)=x2+y2V(x, y) = x^2 + y^2(正定)。

V˙=2x(yx3)+2y(xy3)=2x42y4\dot{V} = 2x(-y - x^3) + 2y(x - y^3) = -2x^4 - 2y^4

V˙<0\dot{V} < 0 对所有 (x,y)(0,0)(x, y) \neq (0, 0)(负定)。

由李雅普诺夫定理,零解渐近稳定。注意:线性化系统的 Jacobi 矩阵特征值为纯虚数 ±i\pm i(非双曲中心),线性化方法失效——这正是李雅普诺夫直接法大显身手的场合。

不变原理(LaSalle 不变原理)

V˙0\dot{V} \leq 0(半负定)且不能证明 V˙<0\dot{V} < 0 时,可用 LaSalle 不变原理判断渐近稳定性:

LaSalle 不变原理:若 VV 正定、径向无界,V˙0\dot{V} \leq 0,且集合 {XV˙(X)=0}\{X \mid \dot{V}(X) = 0\} 中除原点外不含系统的不变集,则零解渐近稳定。

李雅普诺夫函数的意义

  • 无需求解方程即可判断稳定性
  • 特别适用于非双曲不动点(线性化失效)
  • 提供稳定性的充分条件(不是必要条件)
  • 本质是寻找系统的"单调递减的广义能量"