非线性干扰观测器简介
最近要用到非线性干扰观测器,搜了搜网上的资料,关于干扰观测器的资料不少,但是关于 非线性干扰观测器的不多,原因在于非线性干扰观测器较难设计。本文在此简单介绍一下非线性干扰观测器的基本原理及其稳定性证明。
1. 系统方程
这里依然采用系统的标准柯西形式:
{
x
˙
1
=
x
2
x
˙
2
=
f
+
g
U
+
d
(1)
\begin{cases} \dot x_1 = x_2 \\ \dot x_2 = f + gU + d \end{cases} \tag{1}
{x˙1=x2x˙2=f+gU+d(1)其中
U
U
U为控制量,
d
d
d为干扰。
对于不同的实际系统,表达式(1)可能有所不同,但最终都可以通过数学化简得到(1)式的结果。
现不考虑干扰
d
d
d,可以得到系统的理想微分方程组:
{
x
˙
1
=
x
2
x
˙
2
=
f
+
g
U
(2)
\begin{cases} \dot x_1 = x_2 \\ \dot x_2 = f + gU \end{cases} \tag{2}
{x˙1=x2x˙2=f+gU(2)
2. 非线性干扰观测器的设计
设计干扰观测器如下:
z
˙
=
−
L
z
−
L
(
L
x
2
+
f
+
g
U
)
d
^
=
z
+
L
x
2
(3)
\begin{aligned} \dot z &= -L z - L \left(L x_2 + f + gU \right) \\ \hat d &= z + L x_2 \tag{3} \end{aligned}
z˙d^=−Lz−L(Lx2+f+gU)=z+Lx2(3)利用上式得出干扰的估计值
d
^
\hat d
d^,将其引入系统中加以消除,即可以将系统补偿为趋近于理想系统(2)的情况。
3. 非线性干扰观测器稳定性证明
不出意外地,依然使用李雅普诺夫稳定性理论进行证明。
对
d
^
\hat d
d^进行求导:
d
^
˙
=
z
˙
+
L
x
˙
2
=
−
L
z
−
L
(
L
x
2
+
f
+
g
U
)
+
L
(
f
+
g
U
+
d
)
=
−
L
z
−
L
2
x
2
−
L
(
f
+
g
U
)
+
L
(
f
+
g
U
)
+
L
d
=
−
L
(
z
+
L
x
2
)
+
L
d
=
−
L
d
^
+
L
d
=
L
(
d
−
d
^
)
=
L
d
~
(4)
\begin{aligned} \dot{\hat d} &= \dot z + L \dot x_2 \\ &= -Lz - L \left(L x_2 + f + gU \right) + L \left( f + gU + d \right) \\ &= - Lz - L^2 x_2 - L \left( f + gU \right) + L \left( f + gU \right) + Ld \\ &= - L \left( z + L x_2 \right) + Ld \\ &= - L \hat d + Ld \\ &= L \left( d - \hat d \right) \\ &= L \tilde d \end{aligned} \tag{4}
d^˙=z˙+Lx˙2=−Lz−L(Lx2+f+gU)+L(f+gU+d)=−Lz−L2x2−L(f+gU)+L(f+gU)+Ld=−L(z+Lx2)+Ld=−Ld^+Ld=L(d−d^)=Ld~(4)设李雅普诺夫函数
V
=
1
2
d
~
2
V = \frac{1}{2} \tilde d^2
V=21d~2对其求导:
V
˙
=
d
~
d
~
˙
=
d
~
(
d
˙
−
d
^
˙
)
\dot V = \tilde d \dot{\tilde d} = \tilde d \left( \dot d - \dot{\hat d} \right)
V˙=d~d~˙=d~(d˙−d^˙)代入式(4):
V
˙
=
d
~
(
d
˙
−
d
^
˙
)
=
d
~
(
d
˙
−
L
d
~
)
=
−
L
d
~
2
+
d
~
d
˙
(5)
\dot V = \tilde d \left( \dot d - \dot{\hat d} \right) = \tilde d \left( \dot d - L \tilde d \right) = - L \tilde d^2 + \tilde d \dot d \tag{5}
V˙=d~(d˙−d^˙)=d~(d˙−Ld~)=−Ld~2+d~d˙(5)根据基本不等式
a
+
b
≥
2
a
b
a +b \geq 2 \sqrt{ab}
a+b≥2ab对式(5)有
V
˙
≤
−
L
d
~
2
+
d
~
2
+
d
˙
2
2
=
−
(
L
−
1
2
)
d
~
2
+
1
2
d
˙
2
(6)
\dot V \leq - L \tilde d^2 + \frac{\tilde d^2 + \dot d^2}{2} = - \left( L - \frac{1}{2} \right) \tilde d^2 + \frac{1}{2} \dot d^2 \tag{6}
V˙≤−Ld~2+2d~2+d˙2=−(L−21)d~2+21d˙2(6)又有以下定理:
定理1:若李雅普诺夫函数
V
V
V满足
V
˙
≤
−
κ
V
λ
+
θ
,
κ
,
λ
,
θ
>
0
,
且
θ
有界
\dot V \leq - \kappa V^{\lambda} + \theta, \qquad \kappa, \lambda, \theta > 0,\quad 且\theta有界
V˙≤−κVλ+θ,κ,λ,θ>0,且θ有界则系统为稳定的。
则(6)式为稳定的,其中 κ = L − 1 2 , λ = 1 , θ = 1 2 d ˙ 2 \kappa = L - \frac{1}{2}, \lambda = 1, \theta = \frac{1}{2} \dot d^2 κ=L−21,λ=1,θ=21d˙2。
4. 总结
非线性干扰观测器的方程核心即为(3)式,通过计算(3)式可以得到干扰的估计值 d ^ \hat d d^,再引入系统中与干扰 d d d相减,即可得到近似理想的无干扰系统。但使用该干扰观测器仍需注意:
- 干扰必须为有界的,否则定理1将不成立。
- 由于式(3)为迭代计算过程,在实际中需要注意代数环的产生,必要时可使用延迟模块delay进行规避。