现代控制理论 ——05 控制系统的李雅普诺夫稳定性分析

# 控制系统的李雅普诺夫稳定性分析

# 李氏稳定性的定义

何为平衡状态?对于一个系统 x˙=f(x,t)\dot{x}=f(x,t),若果存在状态 xex_e 满足 x˙e0\dot{x}_e\equiv 0,那么该状态即为平衡状态。

  1. 稳定:对于任意实数 ε>0\varepsilon>0,都存在一个实数 δ>0\delta>0 满足 x0xeδ||x_0-x_e||\leq\delta,从任意 x0x_0 触发的解都能够满足 xtxeε||x_t-x_e||\leq\varepsilon,则称 xex_e 在李雅普诺夫意义下是稳定的。
  2. 渐进稳定:当上述解能够满足 xtxeμ||x_t-x_e||\leq\mu,也就是能够收敛到 xex_e 时,则称系统渐进稳定
  3. 不稳定:无论 δ\delta 有多小,都会使得 xtxe>ε||x_t-x_e||>\varepsilon,则称系统不稳定
  1. 大范围渐进稳定:从状态空间中所有初始点出发的轨迹都具有渐进稳定性,那么状态 xex_e大范围渐进稳定

  2. 正定函数:对于函数 V(x)V(x),在区域 SS 内的所有 xx 都有:① V(x)V(x) 中的各分量的偏导均存在;②V(0)=0V(0)=0 ;③当 x0x\neq0 时, V(x)>0(V(x)0)V(x)>0 (V(x)\geq0) 。则称该函数是正定 (半正定) 的。

# 李雅普诺夫第一方法

  1. 线性系统的稳定性判据

李雅普诺夫稳定的充要条件:系统矩阵 AA 的全部特征值实部大于 0,即位于复平面左半部。

  1. 非线性系统的稳定性判据

对于非线性系统 x˙=f(x)\dot{x}=f(x),讨论其在可能平衡状态 xex_e 的稳定性。引入新向量 y=xxey=x-x_e,那么系统的状态方程转换为 y˙=Ay+G(y)y\dot{y}=Ay+G(y)y,其中 AA 为雅克比矩阵。

......

# 李雅普诺夫第二方法

对于状态方程为 x˙=f(x,t),f(0,t)=0\dot{x}=f(x,t),f(0,t)=0 的系统,存在一个具有连续偏导的标量函数 V(x,t)V(x,t),满足

  1. V(x,t)V(x,t) 正定,V˙(x,t)\dot{V}(x,t) 半正定,则系统在原点一致稳定;在此基础上,若对于任意 t0t_0x00x_0\neq0,在 tt0t\geq t_0 时不恒等于 0,则系统在原点渐进稳定;在此基础上,若随着 x||x||\to \inftyV(x,t)V(x,t)\to \infty,则系统在原点大范围渐进稳定

  2. V(x,t)V(x,t) 正定,V˙(x,t)\dot{V}(x,t) 正定,则系统在原点不稳定

例题

用李雅普诺夫第二方法判断以下系统的稳定性。

x˙1=(x1+x2)x22x˙2=(x1+x2)+x1x2\begin{aligned}&\dot{x}_1=-(x_1+x_2)-x_2^2\\&\dot{x}_2=-(x_1+x_2)+x_1x_2\end{aligned}

系统存在的唯一可能平衡状态为 x1=0,x2=0x_1=0,x_2=0,取标量函数 V(x)=x12+x22V(x)=x_1^2+x_2^2,显然 V(x)V(x) 正定,求导有

V˙(x)=2x1x˙1+2x2x˙2=2(x1+x2)2\dot{V}(x)=2x_1\dot{x}_1+2x_2\dot{x}_2=-2(x_1+x_2)^2

负定。除原点外有 x1=x2x_1=-x_2 使得 V˙(x)=0\dot{V}(x)=0,但是系统状态仍在转移中,故 V˙(x)\dot{V}(x) 不会恒定等于 0。且随着 x||x||\to \inftyV(x)V(x)\to \infty,故系统在原点大范围渐进稳定。

:::

# 李雅普诺夫方程判断线性系统的稳定性

  1. 在连续系统 x˙=Ax\dot{x}=Ax 中,在平衡状态 x=0x=0 处是大范围渐进稳定的充要条件:对于给定的正定对称实矩阵 QQ,存在一个正定实对称矩阵 PP ,满足 ATP+PA=QA^TP+PA=-Q。(其中 xTPxx^TPx 就是李雅普诺夫函数 )

  2. 在离散系统 x(k+1)=Gx(k)x(k+1)=Gx(k),在平衡状态 x=0x=0 处是渐进稳定的充要条件:对于给定的正定对称实矩阵 QQ,存在一个正定实对称矩阵 PP ,满足 GTPGP=QG^TPG-P=-Q。(其中 xTPxx^TPx 就是李雅普诺夫函数 )

例题

以下系统的平衡状态在坐标原点,判断其渐进稳定性。

x˙=[0111]x\dot{x}=\begin{bmatrix}0&1\\-1&-1\end{bmatrix}x

P=[p11p12p12p22]P=\begin{bmatrix}p_{11}&p_{12}\\p_{12}&p_{22}\end{bmatrix},由 ATP+PA=IA^TP+PA=-I 有,[0111][p11p12p12p22]+[p11p12p12p22][0111]=[1001]\begin{bmatrix}0&-1\\1&-1\end{bmatrix}\begin{bmatrix}p_{11}&p_{12}\\p_{12}&p_{22}\end{bmatrix}+\begin{bmatrix}p_{11}&p_{12}\\p_{12}&p_{22}\end{bmatrix}\begin{bmatrix}0&1\\-1&-1\end{bmatrix}=\begin{bmatrix}-1&0\\0&-1\end{bmatrix},则 {2p12=1p11p12p22=02p122p22=1{p11=32p12=1p22=12\left\{\begin{aligned}&-2p_{12}=-1\\&p_{11}-p_{12}-p_{22}=0\\&2p_{12}-2p_{22}=-1\end{aligned} \right.\to \left\{\begin{aligned}&p_{11}=\frac{3}{2}\\&p_{12}=1\\&p_{22}=\frac{1}{2}\end{aligned} \right.,得到 P=[3212121]P=\begin{bmatrix}\frac{3}{2}&\frac{1}{2}\\\frac{1}{2}&1\end{bmatrix},验证各阶主子行列式是否大于 0:P11=32>0P_{11}=\frac{3}{2}>0P22=det[p11p12p12p22]=54>0P_{22}=det\begin{bmatrix}p_{11}&p_{12}\\p_{12}&p_{22}\end{bmatrix}=\frac{5}{4}>0,故矩阵 PP 正定。故系统在原点大范围渐进稳定。

以下系统的平衡状态在坐标原点,判断其渐进稳定性。

[x1(k+1)x2(k+1)]=[00.50.51][x1(k)x2(k)]\begin{bmatrix}x_1(k+1)\\x_2(k+1)\end{bmatrix}=\begin{bmatrix}0&0.5\\-0.5&-1\end{bmatrix}\begin{bmatrix}x_1(k)\\x_2(k)\end{bmatrix}

P=[p11p12p12p22]P=\begin{bmatrix}p_{11}&p_{12}\\p_{12}&p_{22}\end{bmatrix},由 GTPGP=IG^TPG-P=-I 有,[00.50.51][p11p12p12p22][00.50.51][p11p12p12p22]=[1001]\begin{bmatrix}0&-0.5\\0.5&-1\end{bmatrix}\begin{bmatrix}p_{11}&p_{12}\\p_{12}&p_{22}\end{bmatrix}\begin{bmatrix}0&0.5\\-0.5&-1\end{bmatrix}-\begin{bmatrix}p_{11}&p_{12}\\p_{12}&p_{22}\end{bmatrix}=\begin{bmatrix}-1&0\\0&-1\end{bmatrix},得到 P=[52274027402710027]P=\begin{bmatrix}\frac{52}{27}&\frac{40}{27}\\\frac{40}{27}&\frac{100}{27}\end{bmatrix},验证各阶主子行列式均大于 0,故矩阵 PP 正定。故系统在原点大范围渐进稳定。

更新于 阅读次数