【MFAC】基于全格式动态线性化的无模型自适应控制

来源：侯忠生教授的《无模型自适应控制：理论与应用》（2013年科学出版社）。

👉对应书本 3.4 单输入单输出系统(SISO)全格式动态线性化(FFDL) 和 4.4 单输入单输出系统(SISO)全格式动态线性化(FFDL)的无模型自适应控制(MFAC)

PFDL仅考虑了在下一时刻的系统输出变化量与当前时刻的一个固定长度滑动时间窗口内的输入变化量之间的动态时变关系，而FFDL还考虑了当前时刻的具有某个长度的滑动时间窗口内的所有控制输入变化量和系统输出变化量对下一时刻的系统输出变化量的影响。

FFDL

全格式动态线性化
(full form dynamic linearization)

数据模型：
$\Delta y\left(k+1\right)={\varphi }_{f,L_y,L_u}^T(k)\Delta H_{L_y,L_u}(k)$

伪偏导(PPD) ${\varphi }_{f,L_y,L_u}(k)$ 的下标 f 表示full
$L_y$ 和 $L_u$ 为伪阶数，$0\le L_y\le n_y$，$1\le L_u\le n_u$.
当 $L_y=0$ 和 $L_u=L$ 时，FFDL数据模型可变为PFDL数据模型；
当 $L_y=0$ 和 $L_u=1$ 时，FFDL数据模型可变为CFDL数据模型。

${\varphi }_{f,L_y,L_u}(k)={\left[{\varphi }_1(k),{\varphi }_2(k),\dots {,\varphi }_{L_y}(k),{\varphi }_{L_y+1}(k),\dots ,{\varphi }_{L_y+L_u}(k)\right]}^T$
$\Delta H_{L_y,L_u}(k)={\left[\Delta y(k),\dots ,\Delta y\left(k-L_y+1\right),\Delta u(k),\dots ,\Delta u\left(k-L_u+1\right)\right]}^T$

SISO-FFDL-MFAC

学习控制律

$u(k)=u\left(k-1\right)+\frac{{\rho }_{L_y+1}{\varphi }_{L_y+1}(k)\left[y^{\ast }\left(k+1\right)-y(k)\right]}{{|{\varphi }_{L_y+1}(k)|}^2+\lambda }-\frac{{\varphi }_{L_y+1}(k)\sum _{i=1}^{L_y}{\rho }_i{\varphi }_i(k)\Delta u(k-i+1)}{{|{\varphi }_{L_y+1}(k)|}^2+\lambda }-\frac{{\varphi }_{L_y+1}(k)\sum _{i=L_y+2}^{L_u+L_y}{\rho }_i{\varphi }_i(k)\Delta u\left(k-L_y-i+1\right)}{{|{\varphi }_{L_y+1}(k)|}^2+\lambda }$

步长因子 $\rho \in (0,1](i=1,2,.\dots ,L_y+L_u)$

PPD参数估计算法

{ ϕ ^ f , L y , L u ( k ) = ϕ ^ f , L y , L u ( k − 1 ) + η Δ H L y , L u ( k − 1 ) μ + ∣ | Δ H L y , L u ( k − 1 ) | ∣ 2 [ Δ y ( k ) − ϕ ^ f , L y , L u T ( k − 1 ) Δ H L y , L u ( k − 1 ) ] \{\hat{\phi}_{f,L_{y},L_{u}}(k) = {\hat{\phi}}_{f,L_{y},L_{u}}\left( {k - 1} \right) + \frac{\eta\Delta H_{L_{y},L_{u}}\left( {k - 1} \right)}{\mu + \left| \middle| \Delta H_{L_{y},L_{u}}\left( {k - 1} \right) \middle| \right|^{2}}\left\lbrack \Delta y(k) - {\hat{\phi}}_{f,L_{y},L_{u}}^{T}(k - 1)\Delta H_{L_{y},L_{u}}\left( {k - 1} \right) \right\rbrack {ϕ^​f,Ly​,Lu​​(k)=ϕ^​f,Ly​,Lu​​(k−1)+μ+∣∣ΔHLy​,Lu​​(k−1)∣∣2ηΔHLy​,Lu​​(k−1)​[Δy(k)−ϕ^​f,Ly​,Lu​T​(k−1)ΔHLy​,Lu​​(k−1)]

其中，μ>0为权重因子，为了让控制算法更具一般性，引入步长因子 η∈(0,2] .

PPD参数重置算法

如果
$\Vert \widehat{{\varphi }_{f,L_y,L_u}}(k)\Vert \le \epsilon$
或
$\Vert \Delta H_{L_y,L_u}\Vert \le \epsilon$
或
$sign\left(\widehat{{\varphi }_{f,L_y,L_u}}(k)\right)\ne sign\left(\widehat{{\varphi }_{f,L_y,L_u}}(1)\right)$
则
$\widehat{{\varphi }_{f,L_y,L_u}}(k)=\widehat{{\varphi }_{f,L_y,L_u}}(1)$

算法重置机制的引入是为了使PPD估计算法具有更强的对时变参数的跟踪能力。

仿真实验

【MFAC】基于全格式动态线性化的无模型自适应控制（Matlab代码）