推迟参数设计的自适应反步控制和自适应神经网络的反步控制设计

推迟参数设计的自适应反步控制和自适应神经网络的反步控制设计

前言

在上一篇文章中，我们从反步法出发，分析基于反步法的跟踪控制设计，以及含有未知参数的自适应反步控制的控制器设计。

上一篇文章链接：Adaptive Backstepping Besign自适应反步控制设计_ADi_hhh的博客

但是上一篇文章中的系统满足匹配条件，比较容易设计。本篇文章针对不满足匹配条件的系统进行自适应反步设计，并针对系统中含有未知参数的系统进行自适应神经网络的的反步控制设计。

匹配与非匹配

在控制理论中的不确定性可以分为两类。当控制输入可以克服不确定性时，系统就是在匹配条件下的，或者看网上有人说不确定性是在控制矩阵B所张成的空间中；而当控制输入不在控制输入下时，系统就是在不匹配条件下的，不能通过控制量的设计来抵消其影响。

举例如下：

上一篇文章中：
$$\begin{array}{rl}{\dot{x}}_1& =f_1(x_1)+x_2\\ {\dot{x}}_2& =f_2(x_1,x_2)\theta +bu\\ y& =x_1\end{array}$$
这样的系统就是匹配条件下的系统，因而我们能通过自适应反步控制设计，设计出控制器和参数更新率。

而形如
$$\begin{array}{rl}{\dot{x}}_1& =f_1(x_1)\theta +x_2\\ {\dot{x}}_2& =u\\ y& =x_1\end{array}(1)$$
则是非匹配的，控制输入无法直接抵消不确定量$\theta$所带来的影响。需要利用反步设计、鲁棒控制等其他方法来设计控制器。

1. 基于自适应反步控制的非匹配条件下的系统控制器设计

问题描述

针对形如公式(1)的系统进行控制器设计，其中$f_1(x)$是已知的函数，$\theta$为未知参数，设计控制器使得$y\to y_d$，并且$y_d$的导数可知。

我们首先按照前面的自适应反步控制方法来设计控制器。

控制器设计

1. 定义系统误差：$z_1=y-y_d=x_1-y_d$(为什么不用e来表示误差了？别问，问就是显得专业)

2. 系统误差的动态方程：
   $$\begin{array}{rl}{\dot{z}}_1& ={\dot{x}}_1-{\dot{y}}_d\\ & =f_1(x_1)\theta +x_2-{\dot{y}}_d\end{array}$$

3. 考虑如下候选李亚普诺夫函数：$V_1=\frac{1}{2}{z}_1^2$

   对其求导得：
   $$\begin{array}{rl}{\dot{V}}_1& =z_1{\dot{z}}_1\\ & =z_1(f_1(x_1)\theta +x_2-{\dot{y}}_d)\end{array}$$
   这时候我们可以把$x_2$当作一个输入量，那么假如
   $$x_2=-k_1{z}_1-f_1(x_1)\hat{\theta }+{\dot{y}}_d$$
   其中$\hat{\theta }$是对$\theta$的估计，估计误差为$\tilde{\theta }=\hat{\theta }-\theta$，引入一虚拟控制量${\alpha }_1$，即令$\alpha_1=-k_1z_1-f_1(x_1) \hat\theta+\dot y_d $

   此时可以定义出第二个误差为$z_2=x_2-{\alpha }_1$

   因为引入了新的估计误差$\tilde{\theta }$，因此考虑如下候选李雅普诺夫函数：
   $$V_2=V_1+\frac{1}{2\gamma }{\tilde{\theta }}^2$$
   对其求导得：
   $$\begin{array}{rl}{\dot{V}}_2& ={\dot{V}}_1+\frac{1}{\gamma }\tilde{\theta }\dot{\hat{\theta }}\\ & =z_1(f_1(x_1)\theta +x_2-{\dot{y}}_d)+\frac{1}{\gamma }\tilde{\theta }\dot{\hat{\theta }}\\ & =z_1(f_1(x_1)\theta +z_2+{\alpha }_1-{\dot{y}}_d)+\frac{1}{\gamma }\tilde{\theta }\dot{\hat{\theta }}\\ & =z_1(f_1(x_1)\theta +z_2+(-k_1{z}_1-f_1(x_1)\hat{\theta }+{\dot{y}}_d)-{\dot{y}}_d)+\frac{1}{\gamma }\tilde{\theta }\dot{\hat{\theta }}\\ & =-k_1{z}_1^2+z_1{z}_2+\tilde{\theta }(-z_1{f}_1(x_1)+\frac{1}{\gamma }\dot{\hat{\theta }})\end{array}(\ast )$$
   则可以设计出参数更新率
   $$\dot{\hat{\theta }}=\gamma z_1{f}_1(x_1)(2)$$
   此时${\dot{V}}_2=-k_1{z}_1^2+z_1{z}_2$，$z_1、z_2$为互连项。

4. 考虑系统第二个误差的动态方程：
   $$\begin{array}{rl}{\dot{z}}_2& ={\dot{x}}_2-{\dot{\alpha }}_1\\ & =u-{\dot{\alpha }}_1\end{array}$$
   值得注意的是此时${\dot{\alpha }}_1$是一个关于$x_1,y_d,{\dot{y}}_d,\hat{\theta }$的复合函数，并且其中含有未知项，这也是上一篇文章所提到的问题。
   $$\begin{array}{rl}{\dot{\alpha }}_1& =\frac{\partial \alpha }{\partial x_1}{\dot{x}}_1+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}\\ & =\frac{\partial \alpha }{\partial x_1}(f_1(x_1)\theta +x_2)+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}\\ & =\frac{\partial \alpha }{\partial x_1}{f}_1(x_1)\theta +\frac{\partial \alpha }{\partial x_1}{x}_2+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}\\ & =w\theta +\varphi \end{array}$$
   其中令$\frac{\partial \alpha }{\partial x_1}{f}_1(x_1)=w$,$\frac{\partial \alpha }{\partial x_1}{x}_2+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}=\varphi$，这两部分都是已知的。

   则：
   $$\begin{array}{rl}{\dot{z}}_2& =u-{\dot{\alpha }}_1\\ & =u-w\theta -\varphi \end{array}$$

5. 考虑如下候选李雅普诺夫函数$V_3=V_2+\frac{1}{2}{z}_2^2$

   对其求导得：
   $$\begin{array}{rl}{\dot{V}}_3& ={\dot{V}}_2+z_2{\dot{z}}_2\\ & =-k_1{z}_1^2+z_1{z}_2+z_2(u-w\theta -\varphi )\end{array}(3)$$
   不知道看到公式(3)，你有什么想法？我的想法则是估计过的参数怎么又出现了，这该怎么处理那？甚至有点后悔自己为什么那么早处理那个未知的参数，但是做到这个地方了，只能再次估计这个未知参数了。采用${\hat{\theta }}_2$来估计$\theta$。则引入新的估计误差${\tilde{\theta }}_2={\hat{\theta }}_2-\theta$.

6. 考虑新的候选李雅普诺夫函数$V_4=V_3+\frac{1}{2{\gamma }_2}{\tilde{\theta }}_2^2$

   对其求导得：
   $$\begin{array}{rl}{\dot{V}}_4& ={\dot{V}}_3+\frac{1}{{\gamma }_2}{\tilde{\theta }}_2\dot{{\hat{\theta }}_2}\\ & =-k_1{z}_1^2+z_1{z}_2+z_2(u-w\theta -\varphi )+\frac{1}{{\gamma }_2}{\tilde{\theta }}_2\dot{{\hat{\theta }}_2}\\ & =-k_1{z}_1^2+z_2(z_1+u-w\theta -\varphi )+\frac{1}{{\gamma }_2}{\tilde{\theta }}_2\dot{{\hat{\theta }}_2}\end{array}(4)$$

7. 则依据公式(4)，可以设计出控制率$u$
   $$u=-k_2{z}_2-z_1+\varphi +w{\hat{\theta }}_2(5)$$
   代入到(4)中，
   $$\begin{array}{rl}{\dot{V}}_4& =-k_1{z}_1^2+z_2(z_1+u-w\theta -\varphi )+\frac{1}{{\gamma }_2}{\tilde{\theta }}_2\dot{{\hat{\theta }}_2}\\ & =-k_1{z}_1^2-k_2{z}_2^2+{\tilde{\theta }}_2(z_{2}w+\frac{1}{\gamma }\dot{{\hat{\theta }}_2})\end{array}(6)$$

8. 依据(6)可以设计出参数更新率为：
   $$\dot{{\hat{\theta }}_2}=-{\gamma }_2{z}_{2}w(7)$$

则公式(2)、(5)、(7)共同构成了系统的控制率和参数更新率。

小结

显然再不匹配条件下很离谱啊，一个未知参数按照在匹配条件下的做法设计，居然需要两个参数来估计，这显然是不合理的。假如这个系统是n阶的，那么就会需要n个估计参数来估计这一个不确定项，这就不可接受了。我们也自然而然地想到能不能将未知参数放到最后来一起处理那？有的，这叫做推迟参数设计。

接下来我们将通过推迟参数设计来重新设计上面这个系统的控制器和参数更新率。

2. 基于自适应反步控制和推迟参数设计的非匹配条件下的系统控制器设计

问题描述

仍旧针对公式(1)所代表的不匹配系统进行控制器设计，其中$f_1(x)$是已知的函数，$\theta$为未知参数，设计控制器使得$y\to y_d$，并且$y_d$的导数可知。

控制器设计

我们设计的步骤前面与“1. 基于自适应反步控制的非匹配条件下的系统控制器设计”的一致，一直到公式(*)，即我们得到了
$$\begin{array}{rl}{\dot{V}}_2& =-k_1{z}_1^2+z_1{z}_2+\tilde{\theta }(-z_1{f}_1(x_1)+\frac{1}{\gamma }\dot{\hat{\theta }})\end{array}(\ast \ast )$$
但是在这个时候不设计参数更新率，直接考虑第二个误差的动态方程。

1. 考虑系统第二个误差的动态方程：
   $$\begin{array}{rl}{\dot{z}}_2& ={\dot{x}}_2-{\dot{\alpha }}_1\\ & =u-{\dot{\alpha }}_1\end{array}$$
   值得注意的是此时${\dot{\alpha }}_1$是一个关于$x_1,y_d,{\dot{y}}_d,\hat{\theta }$的复合函数，并且其中含有未知项，这也是上一篇文章所提到的问题。
   $$\begin{array}{rl}{\dot{\alpha }}_1& =\frac{\partial \alpha }{\partial x_1}{\dot{x}}_1+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}\\ & =\frac{\partial \alpha }{\partial x_1}(f_1(x_1)\theta +x_2)+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}\\ & =\frac{\partial \alpha }{\partial x_1}{f}_1(x_1)\theta +\frac{\partial \alpha }{\partial x_1}{x}_2+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}\\ & =w\theta +\psi +\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}\end{array}$$
   其中令$\frac{\partial \alpha }{\partial x_1}{f}_1(x_1)=w$,$\frac{\partial \alpha }{\partial x_1}{x}_2+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d=\psi$，这两部分都是已知的。

   则：
   $$\begin{array}{rl}{\dot{z}}_2& =u-{\dot{\alpha }}_1\\ & =u-w\theta -\psi -\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}\end{array}$$

2. 考虑如下候选李雅普诺夫函数$V_3=V_2+\frac{1}{2}{z}_2^2$

   对其求导得：
   $$\begin{array}{rl}{\dot{V}}_3& ={\dot{V}}_2+z_2{\dot{z}}_2\\ & =-k_1{z}_1^2+z_1{z}_2+\tilde{\theta }(-z_1{f}_1(x_1)+\frac{1}{\gamma }\dot{\hat{\theta }})+z_2(u-w\theta -\psi -\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }})\\ & =-k_1{z}_1^2+\tilde{\theta }(-z_1{f}_1(x_1)+\frac{1}{\gamma }\dot{\hat{\theta }})+z_2(z_1+u-w\theta -\psi -\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }})\end{array}(8)$$

3. 则依据公式(8)，可以设计出控制率$u$
   $$u=-k_2{z}_2-z_1+\psi +w\hat{\theta }+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }}(9)$$
   代入到(8)中，
   $$\begin{array}{rl}{\dot{V}}_3& =-k_1{z}_1^2+\tilde{\theta }(-z_1{f}_1(x_1)+\frac{1}{\gamma }\dot{\hat{\theta }})+z_2(z_1+u-w\theta -\psi -\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }})\\ & =-k_1{z}_1^2+\tilde{\theta }(-z_1{f}_1(x_1)+\frac{1}{\gamma }\dot{\hat{\theta }})+z_2(z_1+(-k_2{z}_2-z_1+\psi +w\hat{\theta }+\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }})-w\theta -\psi -\frac{\partial \alpha }{\partial \hat{\theta }}\dot{\hat{\theta }})\\ & =-k_1{z}_1^2-k_2{z}_2^2+\tilde{\theta }(-z_1{f}_1(x_1)+\frac{1}{\gamma }\dot{\hat{\theta }})+z_{2}w\tilde{\theta }\\ & =-k_1{z}_1^2-k_2{z}_2^2+\tilde{\theta }(-z_1{f}_1(x_1)+\frac{1}{\gamma }\dot{\hat{\theta }}+z_{2}w)\end{array}(10)$$

4. 依据(10)可以设计出参数更新率为：
   $$\begin{array}{rl}\dot{\hat{\theta }}& =\gamma (z_1{f}_1(x_1)-z_{2}w)\\ & =\gamma (z_1{f}_1(x_1)-z_2\frac{\partial \alpha }{\partial x_1}{f}_1(x_1))\end{array}(11)$$

则公式(9)、(11)共同构成了系统的控制率和参数更新率。

小结

实际上推迟参数设计只是讲未知参数放到最后再进行处理，防止后续过程中再次出现未知参数。上述增加了推迟参数设计的过程中可以看到只用一个参数估计就可以解决非匹配条件下的控制系统，从而避免过参问题。

3. 基于自适应神经网络+反步控制+推迟参数设计的非匹配条件下的系统控制器设计

问题描述

$$\begin{array}{rl}{\dot{x}}_1& =f(x_1)+x_2\\ {\dot{x}}_2& =u\\ y& =x_1\end{array}$$

$f_1(x_1)$ 是未知的光滑函数 ，设计设计控制器使得$y\to y_d$，并且$y_d$的导数可知。

控制器设计

由前面设计自适应神经网络控制的经验可知，我们可以采用RBF神经网络估计$f(x)$
$$f(x_1)=W^{\ast T}S(x_1)+ϵ(x_1)$$
其中$W^{\ast }$是理想的连接权值，$S(x)$是径向基函数，$ϵ(x)$是逼近误差。

显然理想的连接权值我们也是无法直接获取的，因此也需要估计。设计$\hat{W}$为$W^{\ast }$的估计值，其误差表示为$\tilde{W}=\hat{W}-W^{\ast }$

那么接下来我们就能来继续设计控制器了。

1. 定义系统误差：$z_1=y-y_d=x_1-y_d$

2. 系统误差的动态方程：
   $$\begin{array}{c}\begin{array}{rl}{\dot{z}}_1& ={\dot{x}}_1-{\dot{y}}_d\\ & =f_1(x_1)\theta +x_2-{\dot{y}}_d\\ & =W^{\ast T}S(x_1)+ϵ(x_1)+x_2-{\dot{y}}_d\end{array}\end{array}$$

3. 考虑如下候选李亚普诺夫函数：$V_1=\frac{1}{2}{z}_1^2+\frac{1}{2}{\tilde{W}}^T{\Gamma }^{-1}\tilde{W}$

   对$V_1$进行求导：
   $$\begin{array}{rl}{\dot{V}}_1& =z_1{\dot{z}}_1+{\Gamma }^{-1}{\tilde{W}}^T\dot{\hat{W}}\\ & =z_1(W^{\ast T}S(x_1)+ϵ(x_1)+x_2-{\dot{y}}_d)+{\Gamma }^{-1}{\tilde{W}}^T\dot{\hat{W}}\end{array}(12)$$
   引入虚拟控制量${\alpha }_1$,(大多数论文中都是应该直接引入虚拟控制量，而不会告诉你为什么这么做，这里也不再赘述了，上面也解释过了)

   令
   $${\alpha }_1=-k_1{z}_1-{\hat{W}}^{T}S(x_1)+{\dot{y}}_d$$
   定义误差为：$z_2=x_2-{\alpha }_1$，代入公式(12)，得：
   $$\begin{array}{rl}{\dot{V}}_1& =z_1(W^{\ast T}S(x_1)+ϵ(x_1)+x_2-{\dot{y}}_d)+{\Gamma }^{-1}{\tilde{W}}^T\dot{\hat{W}}\\ & =z_1(W^{\ast T}S(x_1)+ϵ(x_1)+z_2+{\alpha }_1-{\dot{y}}_d)+{\Gamma }^{-1}{\tilde{W}}^T\dot{\hat{W}}\\ & =z_1(W^{\ast T}S(x_1)+ϵ(x_1)+z_2+(-k_1{z}_1-{\hat{W}}^{T}S(x_1)+{\dot{y}}_d)-{\dot{y}}_d)+{\Gamma }^{-1}{\tilde{W}}^T\dot{\hat{W}}\\ & =-k_1{z}_1^2+z_1{z}_2-z_1{\tilde{W}}^{T}S(x_1)+z_1ϵ(x_1)+{\Gamma }^{-1}{\tilde{W}}^T\dot{\hat{W}}\\ & =-k_1{z}_1^2+z_1{z}_2+z_1ϵ(x_1)+{\tilde{W}}^T({\Gamma }^{-1}\dot{\hat{W}}-z_{1}S(x_1))\end{array}(13)$$
   考虑推迟参数设计，暂时不对这里得参数进行参数更新率的设计。

   考虑第二个误差的动态方程：
   $$\begin{array}{rl}{\dot{z}}_2& ={\dot{x}}_2-{\dot{\alpha }}_1\\ & =u-{\dot{\alpha }}_1\end{array}$$
   同样值得注意的是此时${\dot{\alpha }}_1$是一个关于$x_1,y_d,{\dot{y}}_d,\hat{W}$的复合函数，并且其中含有未知项。
   $$\begin{array}{rl}{\dot{\alpha }}_1& =\frac{\partial \alpha }{\partial x_1}{\dot{x}}_1+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}}\\ & =\frac{\partial \alpha }{\partial x_1}(W^{\ast T}S(x_1)+ϵ(x_1)+x_2)+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}}\\ & =\frac{\partial \alpha }{\partial x_1}(W^{\ast T}S(x_1)+ϵ(x_1))+\frac{\partial \alpha }{\partial x_1}{x}_2+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d+\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}}\\ & =\frac{\partial \alpha }{\partial x_1}{W}^{\ast T}S(x_1)+\frac{\partial \alpha }{\partial x_1}ϵ(x_1)+\psi +\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}}\end{array}$$
   其中：$\frac{\partial \alpha }{\partial x_1}{x}_2+\frac{\partial \alpha }{\partial y_d}{\dot{y}}_d+\frac{\partial \alpha }{\partial {\dot{y}}_d}{\ddot{y}}_d=\psi$是已知的。

   则对应的${\dot{z}}_2=u-\frac{\partial \alpha }{\partial x_1}{W}^{\ast T}S(x_1)-\frac{\partial \alpha }{\partial x_1}ϵ(x_1)-\psi -\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}}$。

4. 考虑新的李雅普诺夫函数：$V_2=V_1+\frac{1}{2}{z}_2^2$

   对其求导得：
   $$\begin{array}{rl}{\dot{V}}_2& ={\dot{V}}_1+z_2{\dot{z}}_2\\ & =-k_1{z}_1^2+z_1{z}_2+z_1ϵ(x_1)+{\tilde{W}}^T({\Gamma }^{-1}\dot{\hat{W}}-z_{1}S(x_1))+z_2(u-\frac{\partial \alpha }{\partial x_1}{W}^{\ast T}S(x_1)-\frac{\partial \alpha }{\partial x_1}ϵ(x_1)-\psi -\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}})\\ & =-k_1{z}_1^2+z_1ϵ(x_1)+{\tilde{W}}^T({\Gamma }^{-1}\dot{\hat{W}}-z_{1}S(x_1))+z_2(z_1+u-\frac{\partial \alpha }{\partial x_1}{W}^{\ast T}S(x_1)-\frac{\partial \alpha }{\partial x_1}ϵ(x_1)-\psi -\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}})\end{array}(14)$$
   针对公式(14)可以设计控制器$u$

$$u=-k_2{z}_2-z_1+\frac{\partial \alpha }{\partial x_1}{\hat{W}}^{T}S(x_1)+\psi +\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}}(15)$$

​ 将公式(15)代入到公式(14)中去：

​
$$\begin{array}{rl}{\dot{V}}_2& =-k_1{z}_1^2+z_1ϵ(x_1)+{\tilde{W}}^T({\Gamma }^{-1}\dot{\hat{W}}-z_{1}S(x_1))+z_2(z_1+u-\frac{\partial \alpha }{\partial x_1}{W}^{\ast T}S(x_1)-\frac{\partial \alpha }{\partial x_1}ϵ(x_1)-\psi -\frac{\partial \alpha }{\partial \hat{W}}\dot{\hat{W}})\\ & =-k_1{z}_1^2-k_2{z}_2^2+z_1ϵ(x_1)+{\tilde{W}}^T({\Gamma }^{-1}\dot{\hat{W}}-z_{1}S(x_1))+z_2(\frac{\partial \alpha }{\partial x_1}{\tilde{W}}^{T}S(x_1)-\frac{\partial \alpha }{\partial x_1}ϵ(x_1))\\ & =-k_1{z}_1^2-k_2{z}_2^2+{\tilde{W}}^T({\Gamma }^{-1}\dot{\hat{W}}-z_{1}S(x_1)+z_2\frac{\partial \alpha }{\partial x_1}S(x_1))+ϵ(x_1)(z_1-z_2\frac{\partial \alpha }{\partial x_1})\end{array}(16)$$

5. 根据公式(16)可以设计出参数更新率为：
   $$\dot{\hat{W}}=\Gamma (z_{1}S(x_1)-z_2\frac{\partial \alpha }{\partial x_1}S(x_1))(17)$$
   则对应的${\dot{V}}_2$为：
   $$\begin{array}{rl}{\dot{V}}_2& =-k_1{z}_1^2-k_2{z}_2^2+ϵ(x_1)(z_1-z_2\frac{\partial \alpha }{\partial x_1})\end{array}(18)$$
   之后再对$ϵ(x_1)(z_1-z_2$