自适应神经网络控制

自适应神经网络控制

基本思路

自适应控制率 $u^{*}$

$u^{*} =-\frac{1}{b( x)}[ a( x) +v] -\left[\frac{1}{\varepsilon b( x)} +\frac{1}{\varepsilon b^{2}( x)} -\frac{\dot{b( x)}}{2b^{2}( x)}\right] e_{s}$ ———— (1)

之所以这么设计的原因，一方面是基于反馈控制的基本思想，设计误差动态系统，另一方面是为了保证系统的稳定性，也就是为了使得基于误差系统的Lyapunov稳定性可以保证

误差 $e_{s}$ 动态方程

$\dot{e_{s}} =a( x) +b( x) u+v+d( t)$ ———— (2)
$v=-y^{( n)}_{d} +\left[ 0\ \Lambda ^{T}\right] e$

利用神经网络 $u^{*}_{1}( Z)$ 拟合控制 $u^{*}$

$u^{*} =u^{*}_{1}( Z) -\frac{1}{\varepsilon } e_{s}$ ———— (3)

HONN自适应控制

神经网络控制器$u_{1}(Z)$

$u_{1}(Z)=\hat{W} \cdot S( Z)$ ———— (4)

自适应律

$\mathcal{\dot{\hat{W}} =-} \Gamma \left[ S( Z) e_{s} +\sigma \hat{W}\right]$ ———— (5)

拟合误差

$u^{*}_{1}( Z) =W^{*} \cdot S( Z) +\mu _{l}$ ———— (6)
$W^{*} \cdot S( Z)$表示最优网络
$\mu _{l}$表示拟合误差

推导误差 $e_{s}$ 动态方程
1. 将控制(4)代入误差动态方程(2)

$\dot{e_{s}} =a( x) +b( x)\hat{W} \cdot S( Z) +v+d( t)$

1. 添加辅助项 $b( x) u^{*}$

   $\dot{e_{s}} =a( x) +b( x)\left[\hat{W} \cdot S( Z) - u^{*}\right] +b( x) u^{*} +v+d( t)$ ———— (7)

2. 将(6)代入

   $\dot{e_{s}} =a( x) +v+b( x)\left[\hat{W} \cdot S( Z) -W^{*} \cdot S( Z) -\mu _{l}\right] +b( x) u^{*} +d( t)$
   由于$\tilde{W} \cdot S( Z) =\hat{W} \cdot S( Z) -W^{*} \cdot S( Z)$，则
   $\dot{e_{s}} =a( x) +v+b( x)\left[\tilde{W} \cdot S( Z) -\mu _{l}\right] +b( x) u^{*} +d( t)$
   这里体现出线性神经网络的优点，便于化简

3. 将(1)代入

   由于$b( x) u^{*} =-[ a( x) +v] -b( x)\left[\frac{1}{\varepsilon b( x)} +\frac{1}{\varepsilon b^{2}( x)} -\frac{\dot{b( x)}}{2b^{2}( x)}\right] e_{s}$，则
   $\dot{e_{s}} =b( x)\left[\tilde{W} \cdot S( Z) -\mu _{l}\right] -b( x)\left[\frac{1}{\varepsilon b( x)} +\frac{1}{\varepsilon b^{2}( x)} -\frac{\dot{b( x)}}{2b^{2}( x)}\right]e_{s} +d( t)$
   $\dot{e_{s}} =b( x)\left[\frac{\dot{b( x)}}{2b^{2}( x)}e_{s} -\frac{1}{\varepsilon b( x)}e_{s} -\frac{1}{\varepsilon b^{2}( x)}e_{s} +\tilde{W} \cdot S( Z) -\mu _{l}\right] +d( t)$ ———— (8)

Lyapunov方程

$V_{s} =\frac{1}{2}\left[\frac{e^{2}_{s}}{b( x)} +\tilde{W}^{T} \Gamma \mathcal{^{-1}} \tilde{W}\right]$

稳定性证明

$\dot{V}_{s} =\frac{e_{s} \cdot \dot{e}_{s}}{b( x)} -\frac{\dot{b( x)}}{2b( x)} e^{2}_{s} +\tilde{W}^{T} \Gamma \mathcal{^{-1}\dot{\hat{W}}}$

将(5)(8)代入，得

$\dot{V}=e_{s}\left[\frac{\dot{b( x)}}{2b^{2}( x)}e_{s} -\frac{1}{\varepsilon b( x)}e_{s} -\frac{\dot{b( x)}}{2b( x)} e^{2}_{s} +\tilde{W} \cdot S( Z) -\mu _{l}\right] +e_{s} \cdot \dfrac{d( t)}{b( x)} -e^{2}_{s}\frac{\dot{b( x)}}{2b^{2}( x)} +$
$\ \ \ \ \ \ \ \ \tilde{W}^{T} \Gamma \mathcal{^{-1}\left[ - \Gamma \left[ S( z) e_{s} +\sigma \hat{W}\right]\right]}$
$\dot{V} =e^{2}_{s}\left( -\frac{1}{\varepsilon b( x)} -\frac{1}{\varepsilon b^{2}( x)}\right) -e_{s} \mu _{l} +e_{s} \cdot \dfrac{d( t)}{b( x)} -\tilde{W}^{T}\mathcal{\sigma \hat{W}}$ ———— (9)

补充说明1

Lyapunov稳定性证明可以分为两种思路
第一直接设计的控制率使得 $\dot{V} < 0$
第二构造形如 $\dot{V} < V+\beta ( x)$

补充说明2

控制率和自适应律的设计也分为两部分
第一部分为了消除Lyapunov中的某些项
第二部分是为了满足稳定性条件，此部分有两种思路，第一种为了满足 $\dot{V} < 0$ ，第二种是为了构造 ${V}$ ，使得 $\dot{V} < V+\beta ( x)$ 的条件满足

从化简的过程可以出来
(5)自适应律的设计分为两部分，第一部分 $S( z) e_{s}$ 是为了抵消Lyapunov $\dot{V}$的非线性项，第二部分 $\sigma \hat{W}$ 是为了构造 $V$

再对(9)继续化简之前，需要考虑以下不等式

$2ab\leq a^{2} +b^{2}$
$2\dfrac{a}{\varepsilon }( \varepsilon b) \leq \left(\dfrac{a}{\varepsilon }\right)^{2} +( \varepsilon b)^{2}$

---

Lyapunov方程中的项采用如下不等式进行化简

$2\tilde{W}^{T}\mathcal{\hat{W} =} \left\Vert \tilde{W}\right\Vert ^{2} +\left\Vert \mathcal{\hat{W}} \right\Vert ^{2} -\left\Vert W^{*}\right\Vert ^{2} \geq \left\Vert \tilde{W}\right\Vert ^{2} -\left\Vert W^{*}\right\Vert ^{2}$
$\dfrac{d( t)}{b( x)} e_{s} =\dfrac{e_{s}}{b( x)}\left(\frac{1}{\sqrt{\dfrac{\varepsilon }{2}}}\right) \cdot \left(\sqrt{\dfrac{\varepsilon }{2}} d( t)\right) \leq \dfrac{1}{2}\left[\dfrac{e_{s}}{b( x)}\left(\frac{1}{\sqrt{\dfrac{\varepsilon }{2}}}\right)\right]^{2} +\dfrac{1}{2}\left(\sqrt{\dfrac{\varepsilon }{2}} d( t)\right)^{2}$
$\ \ \ \ \ \ \ \ \ \ \ \ =\dfrac{e^{2}_{s}}{\varepsilon b^{2}( x)} +\dfrac{\varepsilon }{4} d^{2}( t)$
$| \mu _{l} e_{s}| \leq \dfrac{e^{2}_{s}}{2\varepsilon b( x)} +\dfrac{\varepsilon }{2} \mu ^{2}_{l} b( x)$
$| \mu _{l}| \leq \mu _{0}$
$|d_{t}| \leq d_{0} \ $
$0< b( x) < b_{1}$

因此Lyapunov方程可以化简为

$\dot{V} \leq -\frac{e^{2}_{s}}{2\varepsilon b^{2}( x)} +\dfrac{\varepsilon }{2} \mu ^{2}_{0} b_{1} +\dfrac{\varepsilon }{4} d^{2}_{0} -\dfrac{\sigma }{2}\left\Vert \tilde{W}\right\Vert ^{2} +\dfrac{\sigma }{2}\left\Vert W^{*}\right\Vert ^{2}$

对于该Lyapunov方程的处理我们采用构造 $\dot{V} < V+\beta ( x)$ 的方式进行稳定性证明
再进一步处理之前，我们需要利用如下不等式对Lyapunov函数 ${V}_{s}$ 进行变换

$\tilde{W}^{T} \Gamma \mathcal{^{-1}} \tilde{W} \leq \lambda _{max}\left\Vert \tilde{W}\right\Vert ^{2}$
$-\dfrac{\sigma }{2}\left\Vert \tilde{W}\right\Vert ^{2} \leq -\dfrac{\sigma }{2}\dfrac{\tilde{W}^{T} \Gamma \mathcal{^{-1}} \tilde{W}}{\lambda _{max}}$

则 $\dot{V}$ 可以转化为

$\dot{V} \leq -\frac{e^{2}_{s}}{2\varepsilon b^{2}( x)} -\dfrac{\sigma }{2}\dfrac{\tilde{W}^{T} \Gamma \mathcal{^{-1}} \tilde{W}}{\lambda _{max}} +\dfrac{\varepsilon }{2} \mu ^{2}_{0} b_{1} +\dfrac{\varepsilon }{4} d^{2}_{0} +\dfrac{\sigma }{2}\left\Vert W^{*}\right\Vert ^{2}$

假设

$\alpha _{0} =max\left\{\varepsilon ,\dfrac{\lambda _{max}}{\sigma }\right\}$
$\sigma _{0} =\dfrac{\varepsilon }{2} \mu ^{2}_{0} b_{1} +\dfrac{\varepsilon }{4} d^{2}_{0} +\dfrac{\sigma }{2}\left\Vert W^{*}\right\Vert ^{2}$

则

$\dot{V} \leq -\dfrac{1}{\alpha _{0}} V+\dfrac{1}{2\sigma _{0}}$