-
摘要: 针对比例—积分—微分(Proportional-integral-differential, PID)控制器的整定问题, 提出了自耦PID (Self-coupling PID, SC-PID)控制方法.该方法将系统动态和内外不确定性定义为总和扰动, 从而将非线性不确定系统变换为线性不确定系统, 进而构建了总和扰动反相激励下的误差动态系统, 据此设计了SC-PID控制律模型和整定规则, 进而设计了自适应速度因子(Adaptive speed factor, ASF)模型.数值仿真结果表明, SC-PID具有快的响应速度、高的控制精度、良好的抗总和扰动鲁棒性等诸多优点. SC-PID整定规则为现有PID整定结果的技术评估与技术升级提供了科学的理论依据, 在国防和工业控制领域具有广泛的应用价值.Abstract: Aiming at tuning problem for proportional-integral-differential (PID) controller, the self-coupling PID (SC-PID) control method was proposed. Dynamic and uncertainty both inside and outside system was defined as a total disturbance in the proposed method, so as to transform nonlinear uncertain system to linear uncertain system, and then error dynamic system (EDS) inversely driven by the total disturbance was built. Control law and its tuning rules for SC-PID were designed on the basis of EDS, and then adaptive speed factor (ASF) model was designed. Numerical simulation results showed that SC-PID has fast response, high control accuracy, good global stability robustness and good total disturbance rejection robustness etc. many advantages. The SC-PID tuning rules will provide a scientific theoretical basis for the technical evaluation and technical upgrading of the existing PID tuning results, and has extensive application value in national defense and industrial control.
-
Key words:
- Nonlinear uncertain system /
- adaptive speed factor /
- SC-PID control /
- tuning rule /
- global stability
-
近年来, 多传感器信息融合技术得到广泛关注.经典Kalman滤波需要已知模型参数和噪声统计特性[1].而在实际应用中, 系统可能会出现模型参数或噪声统计特性未知情况.文献[2-4]研究了带未知噪声方差系统的自校正估计问题.文献[5]和[6]研究了模型参数和噪声方差未知系统的自校正融合估计问题.文献[7]针对带未知模型参数和噪声统计特性且噪声相关的多传感器系统, 分别应用RELS算法、Gevers-Wouters算法对未知模型参数和噪声统计特性进行辨识, 并提出了自校正融合估值器.文献[5-7]中在对未知模型参数进行融合处理时, 采用加权平均方法获得最终融合辨识器.该方法没有考虑不同传感器的局部参数辨识的差异, 不能保证融合后的参数估计精度都优于每个局部参数估计.上述文献所提出的辨识和估计算法都是基于完整的传感器观测数据, 而没有考虑数据的不完整现象.
在实际的网络化系统或传感器网络中, 由于传感器老化或故障、以及传感器的观测数据在通信传输过程中由于带宽有限, 可能会出现数据丢失、衰减、延迟等问题, 使得估值器所收到的传感器数据具有不完整现象[8].文献[9]研究了带多丢包系统的最优估计问题.文献[10]将文献[9]的单传感器系统推广到了多传感器系统, 对带有不同丢包率的多传感器系统, 提出了集中式和分布式融合估值器.文献[11-14]考虑数据在传输过程中存在的丢失和延迟现象, 设计了相应的最优估值器.其中文献[13]在文献[12]的基础上研究了带随机乘性噪声、多丢包和滞后系统的最优估计问题.文献[14]考虑了过程噪声和观测噪声具有一步自相关和互相关的情况.文献[15]和[16]对带有数据包丢失的多传感器系统, 应用协方差信息方法设计了分布式融合估值器.在上述文献[9-17]中数据丢失现象均由一组满足伯努利分布的随机变量所描述, 这种数据丢失可以看作是观测数据发生衰减现象的一种特殊情况.文献[18]和[19]考虑了带随机参数矩阵、相关噪声和衰减观测系统的最优和次优估计问题.文献[20]中研究了带衰减观测系统的Kalman滤波估计问题, 同时分析了误差协方差的有界性和稳态特性.文献[21-22]对带衰减观测的随机不确定系统, 分别设计了多传感器分布式和序贯融合估值器.上述文献都是在假设观测丢失率或衰减率已知的情况下得到的最优估计结果.而没有考虑观测丢失率或衰减率未知的自校正估计问题.目前, 有关同时考虑带有未知模型参数和未知衰减观测率多传感器系统的自校正融合估计问题的报导甚少.
基于以上文献分析, 本文将对带有未知模型参数和衰减观测率的多传感器随机系统, 应用相关函数和递推增广最小二乘算法, 分别在线辨识衰减观测的数学期望、方差和模型参数.应用线性无偏最小方差估计准则, 分别提出了分布式模型参数融合辨识器和自校正融合状态滤波器, 并分析算法的收敛性.
1. 问题阐述
考虑带衰减观测的多传感器随机系统
$$ \begin{equation} {\boldsymbol x(t+1)} = {\Phi}{\boldsymbol x(t)}+{\Gamma}{\boldsymbol w}(t) \end{equation} $$ (1) $$ \begin{equation} {{y}_{i}}(t) = {\mu_{i}(t)}{ h_{i}}{\boldsymbol x}(t)+{v_{i}(t)}, \ i = 1, 2, \cdots, L \end{equation} $$ (2) 其中状态$ \boldsymbol x(t)\in {\bf R}^n $, 观测$ {{y}_{i}}(t)\in {{\bf R}} $, 过程噪声$ \boldsymbol w(t)\in {\bf R}^r $, 观测噪声$ {{v}_{i}}(t)\in {{\bf R}} $, 下标$ i $表示第$ i $个传感器, $ L $表示传感器的个数. $ \left\{ {{\mu _i}(t)} \right\} $为一组在${[0, 1]}$区间取值的用来描述第$ i $个传感器衰减观测现象的标量随机变量.其中$ {\rm E}{[\mu_{i}(t)] = \alpha_{i}} $, $ {\rm Cov}[\mu_{i}(t)] = \sigma_{i}^2 $, $ {\rm E} $为数学期望符号, $ {\rm Cov} $为协方差符号. $ \left\{ {{\mu _i}(t)} \right\} $与其他随机变量不相关. $ \Phi $, $ \Gamma $, $ { h_{i}} $是适当维数矩阵.
假设1. $ {\boldsymbol w}(t) $和$ {v_{i}(t)} $为零均值、方差阵分别为$ Q_{{\boldsymbol w}} $和$ Q_{v_{i}} $的不相关白噪声.
假设2. 初值$ {\boldsymbol x(0)} $不相关于$ {\boldsymbol w}(t) $和$ {v_{i}(t)} $, 且$ {\rm{E}}\left\{ {\boldsymbol x(0)} \right\} = {\boldsymbol u_0}, {\rm{E}}\left\{ {\left[{\boldsymbol x(0) - {\boldsymbol u_0}} \right]{{\left[{\boldsymbol x(0) - {\boldsymbol u_0}} \right]}^{\rm{T}}}} \right\} = {P_0} $.其中$ {\rm T} $为转置号.
假设3. $ \Phi $为稳定矩阵, ($ \Phi, { h_{i}} $)是完全可观对, ($ \Phi, \Gamma $)是完全可控对.
假设4. $ \Phi $中部分参数未知, $ \left\{ {{\mu _i}(t)} \right\} $的数学期望$ {\alpha_i} $和方差$ \sigma_{i}^2 $未知.
问题是基于观测$ ({y}_{i}(1), \cdots, {y}_{i}(t)) $, $ i = 1, 2, $ $ \cdots, L $, 辨识$ \Phi $中未知参数、$ \left\{ {{\mu _i}(t)} \right\} $的数学期望$ {\alpha_i} $和方差$ \sigma_{i}^2 $, 并求$ \Phi $中未知参数的融合辨识器$ {{\rm{\hat \Phi }}_o}(t) $和状态$ {\boldsymbol x(t)} $的自校正融合滤波器$ {\hat{\boldsymbol {x}}_s(t|t)} $.
注1. 为了阅读方便, 这里对后文中经常遇到一些术语给予解释.局部滤波器, 即基于单个传感器的观测数据所获得的滤波器; 最优滤波器, 即系统的模型参数以及衰减观测期望和方差已知情况下, 获得的线性最小方差意义下的滤波器; 自校正滤波器, 即系统模型含有未知参数, 以及衰减观测期望和方差未知情况下, 通过辨识这些未知参数, 然后代入最优滤波算法中获得的滤波器.分布式融合滤波器, 即基于各个传感器的局部滤波器, 应用线性无偏最小方差矩阵加权融合估计算法[23]获得的融合滤波器.
2. 分布式最优融合滤波
当系统模型参数、$ \left\{ {{\mu _i}(t)} \right\} $的数学期望$ {\alpha_i} $和方差$ \sigma_{i}^2 $已知时, 应用线性无偏最小方差意义下的矩阵加权融合估计算法[23]可获得分布式最优融合滤波器.下面给出实现过程.
由式(2)可得
$$ \begin{equation} {{y}_{i}}(t) = {\alpha_i}{ h_{i}}{\boldsymbol x(t)}+{V_i(t)} \end{equation} $$ (3) 其中
$$ \begin{equation} {V_i(t)} = (\mu_{i}(t)-{\alpha_i}){h_{i}}{\boldsymbol x(t)}+v_i(t) \end{equation} $$ (4) 可计算其方差阵为
$$ \begin{equation} Q_{{V_i(t)}} = {\rm E}[V_i^2(t)] = {{{\sigma}}_{i}^2}{ h_{i}}X(t){ h_{i}^{\rm T}}+Q_{v_i} \end{equation} $$ (5) 状态二阶矩$ X(t) = {\rm E}[{\boldsymbol x}(t){\boldsymbol x}^{\rm T}(t)] $可递推计算如下:
$$ \begin{equation} X(t + 1) = {\rm{\Phi }}X(t){{\rm{\Phi }}^{\rm{T}}} + {\rm{\Gamma }}{Q_{\boldsymbol w}}{{\rm{\Gamma }}^{\rm{T}}} \end{equation} $$ (6) 初值为$ X(0) = {\boldsymbol \mu_0}{\boldsymbol \mu}_0^{\rm{T}}+ {P_0} $.由假设3可知$ X(t) $有界.
下面引理1给出了最优局部滤波器算法; 引理2给出了互协方差计算公式; 引理3给出了分布式最优加权融合滤波算法.
引理1[24]. 在假设1$ \sim $3下, 随机系统(1)和(3)基于每个传感器的观测有最优局部滤波器
$$ \begin{equation} {\hat {\boldsymbol x}_i}(t + 1|t + 1) = {\Psi _{fi}}(t + 1){\hat {\boldsymbol x}_i}(t|t)+{{K}_i}(t + 1){{ y}_i}(t + 1) \end{equation} $$ (7) $$ \begin{equation} {{ K}_i}(t + 1) = {\Sigma _i}(t + 1|t) F_i^{\rm{T}}Q_{{{{C}}_i}}^{ - 1}(t + 1) \end{equation} $$ (8) $$ \begin{equation} {\Sigma _i}(t + 1|t) = \Phi {P_i}(t|t){\Phi ^{\rm{T}}} + \Gamma {Q_{\boldsymbol w}}{\Gamma ^{\rm{T}}} \end{equation} $$ (9) $$ \begin{equation} {Q_{{{{C}}_i}}}(t + 1) = { F_i}{\Sigma _i}(t + 1|t)F_i^{\rm{T}} + {Q_{{V_i}}}(t + 1) \end{equation} $$ (10) $$ \begin{equation} {P_i}(t + 1|t + 1) = [{I_n} - {{ K}_i}(t + 1){ F_i}]{\Sigma _i}(t + 1|t) \end{equation} $$ (11) 其中$ { F_i} = {\alpha _i}{ h_i}, {\Psi _{fi}}(t + 1) = {\rm{[}}{I_n}-{{ K}_i}(t + 1){F_i}]\Phi $. $ {\hat {\boldsymbol x}_i}(t|t) $为第$ i $个传感器的局部滤波器, $ {K_i}(t + 1) $为相应的滤波增益, $ {P_i}(t|t) $为滤波误差方差阵.初值为$ {\hat {\boldsymbol x}_i}(0|0) = {\boldsymbol u_0}, {P_i}(0|0) = {P_0} $.
引理2[24]. 任意两个局部滤波误差之间的互协方差阵$ {P_{ij}}(t|t) = {\rm E}[{\tilde {\boldsymbol x}_i}(t|t)\tilde {\boldsymbol x}_j^{\rm T}(t|t)] $ (其中滤波误差$ {\tilde {\boldsymbol x}_i}(t|t) = {\boldsymbol x}(t)-{\hat {\boldsymbol x}_i}(t|t) $)可递推计算如下:
$$ \begin{align} {P_{ij}}&(t + 1|t + 1) = [{I_n} - { K_i}(t + 1){ F_i}]\times \\ &[\Phi{P_{ij}}(t|t){\Phi ^{\rm{T}}} + \Gamma {Q_{\boldsymbol w}}{\Gamma ^{\rm{T}}}]{[{I_n} - { K_j}(t + 1){ F_j}]^{\rm{T}}} \end{align} $$ (12) 初值为$ {P_{ij}}(0|0) = {P_0} $.
引理3[23]. 基于引理1的各局部滤波器和引理2的任意两个局部滤波误差之间的互协方差阵, 分布式最优矩阵加权融合滤波器可计算如下:
$$ \begin{equation} {\hat {\boldsymbol x}_o}(t|t) = \sum\limits_{i = 1}^L {W_i^{}(t){{\hat {\boldsymbol x}}_i}(t|t)} \end{equation} $$ (13) 加权矩阵计算为
$$ \begin{equation} [W_1^{}(t), \cdots , W_L^{}(t)]{\kern 1pt} {\kern 1pt} {\kern 1pt} = {(e_{}^{\rm{T}}P_{}^{ - 1}(t|t)e)^{ - 1}}e_{}^{\rm{T}}P_{}^{ - 1}(t|t) \end{equation} $$ (14) 其中$ e = [I_{n}, \cdots, I_{n}]^{\mathrm{T}}, P(t|t) = {\left[{{P_{ij}}(t|t)} \right]_{nL \times nL}} $是以$ {P_{ij}}(t|t) $为第$ (i, j) $元素的分块矩阵.融合滤波器的估计误差方差阵计算为
$$ \begin{equation} {P_o}(t|t) = {(e_{}^{\rm T}P_{}^{ - 1}(t|t)e)^{ - 1}} \end{equation} $$ (15) 且有$ {P_o}(t|t) \le {P_i}(t|t), i = 1, \cdots, L. $
3. 未知模型参数融合辨识器
上一节我们针对系统模型精确已知时给出了分布式最优融合估计算法.而在实际应用中, 系统模型可能含有未知参数.当$ \Phi $中含有未知参数时, 本节采用RELS算法辨识未知模型参数, 并对辨识得到的$ L $组参数估值进行加权融合, 获得模型参数的分布式融合辨识器.下面给出具体实现过程.
由式(1)可得
$$ \begin{equation} {\boldsymbol x}(t) = {({I_n} - {q^{ - 1}}\Phi )^{ - 1}}{q^{ - 1}}\Gamma {\boldsymbol w}(t) \end{equation} $$ (16) 式中$ q^{-1} $为单位滞后算子, 即$ {q^{ - 1}}{\boldsymbol x}(t) = {\boldsymbol x}(t - 1) $.将式(16)代入式(3)得
$$ \begin{equation} {y_i}(t) = {\alpha _i}{ h_i}{({I_n} - {q^{ - 1}}\Phi )^{ - 1}}{q^{ - 1}}\Gamma {\boldsymbol w}(t) + {V_i}(t) \end{equation} $$ (17) 将式(17)进一步化简得
$$ \begin{equation} A({q^{ - 1}}){ y_i}(t) = {\alpha _i}{B_i}({q^{-1}}){\boldsymbol w}(t)+ A({q^{ - 1}}){V_i}(t) \end{equation} $$ (18) 其中$ A({q^{-1}}) = \det ({I_n}-{q^{ - 1}}\Phi), {B_i}({q^{ - 1}}) = { h_i}{\rm{adj}}({I_n}-{q^{ - 1}}\Phi){q^{ - 1}}\Gamma $, 式中det和adj分别表示矩阵行列式和伴随矩阵. $ A(q^{-1}) $和$ B_i(q^{-1}) $具有如下多项式形式:
$$ \begin{equation} A({q^{ - 1}}) = 1 + {a_1}{q^{ - 1}} + \cdots + {a_{{n_a}}}{q^{ - {n_a}}} \end{equation} $$ (19) $$ \begin{equation} {B_i}({q^{ - 1}}) = {B_{i1}}{q^{ - 1}} + \cdots + {B_{i{n_{{B_i}}}}}{q^{ - {n_{{B_i}}}}} \end{equation} $$ (20) 其中$ {a_k}, k = {\rm{1}}, \cdots, {n_a} $和$ {B_{ik}}, k = {\rm{1}}, \cdots, {n_{{B_i}}} $是多项式系数. $ n_A, n_{B_i} $分别为$ A(q^{-1}) $和$ B_i(q^{-1}) $的阶次.式(18)等号右侧两个滑动平均过程可以等价为一个稳定的滑动平均过程$ {D_i}({q^{ - 1}}){\varepsilon _i}(t) $[24], 即
$$ \begin{equation} {D_i}({q^{ - 1}}){\varepsilon _i}(t) = {\alpha _i}{B_i}({q^{ - 1}}){\boldsymbol w}(t) + A({q^{ - 1}}){V_i}(t) \end{equation} $$ (21) 其中$ {\varepsilon _i}(t) $是零均值且带有未知噪声方差$ \sigma _{{\varepsilon _i}}^2 $的白噪声, $ {D_i}({q^{ - 1}}) $具有如下多项式形式:
$$ \begin{equation} {D_i}({q^{ - 1}}) = 1 + {d_{i1}}{q^{ - 1}} + \cdots + {d_{i{n_{{D_i}}}}}{q^{ - {n_{{D_i}}}}} \end{equation} $$ (22) 其中$ {d_{ik}}, k = 1, \cdots, {n_{{D_i}}} $是多项式$ {D_i}({q^{ - 1}}) $的系数, $ {n_{{D_i}}} $是$ {D_i}({q^{ - 1}}) $的阶次.
将式(18)重写为
$$ \begin{equation} A({q^{ - 1}}){y_i}(t) = {D_i}({q^{ - 1}}){\varepsilon _i}(t) \end{equation} $$ (23) 令$ \varphi _i^{\rm{T}}(t) = [-{y_i}(t- 1), \cdots, -{y_i}(t- {n_a}), {\hat \varepsilon _i}(t - 1), \cdots, {\hat \varepsilon _i}(t -{n_{{D_i}}})], {\vartheta _i} = [{a_1}, \cdots, {a_{{n_a}}}, {d_{i1}}, \cdots, {d_{i{n_{{D_i}}}}}{{\rm{]}}^{\rm{T}}} $, 则式(23)可表示为
$$ \begin{equation} {y_i}(t) = \varphi _i^{\rm{T}}(t){\vartheta _i} + {\varepsilon _i}(t) \end{equation} $$ (24) 参数$ {a_k}, k = 1, \cdots, {n_a}; {d_{ik}}, k = 1, \cdots, {n_{{D_i}}} $未知.
基于每个单传感器的观测数据, 应用RELS算法[24]可得到局部参数估计为
$$ \begin{equation} {\hat \vartheta _i}(t + 1) = {\hat \vartheta _i}(t) + {M_i}(t + 1){\hat \varepsilon _i}(t{\rm{ + }}1) \end{equation} $$ (25) $$ \begin{equation} {\hat \varepsilon _i}(t + 1) = {y_i}(t + 1) - \varphi _i^{\rm{T}}(t + 1){\hat \vartheta _i}(t) \end{equation} $$ (26) $$ \begin{equation} {M_i}(t + 1) = \frac{{{Z_i}(t){\varphi _i}(t + 1)}}{{1 + \varphi _i^{\rm{T}}(t + 1){Z_i}(t){\varphi _i}(t + 1)}} \end{equation} $$ (27) $$ \begin{equation} {Z_i}(t + 1) = [{I_{{n_a} + {n_{{D_i}}}}} - {M_i}(t + 1)\varphi _i^{\rm{T}}(t + 1)]{Z_i}(t) \end{equation} $$ (28) 代初值$ {\hat \vartheta _i}{\rm{(}}0) = 0, \; {Z_i}(0) = {\beta _i}I, {\beta _i} $为充分大的正数, 且规定$ {\hat \varepsilon _i}(j) = 0, \; {y_i}(j) = 0\; (j \le 0) $.
由文献[24]可知, RELS算法参数估计是一致的, 即$ {\hat \vartheta _i}(t) \to {\vartheta _i}, t \to \infty, w.p.1 $.符号"$ w.p.1. $"表示"以概率1".
由式(25)可得基于单传感器的局部参数估计误差$ {\tilde \vartheta _i}(t) = {\vartheta _i} - {\hat \vartheta _i}(t) $满足如下方程:
$$ \begin{align} {\tilde \vartheta _i}(t + 1) = \, & [{I_{{n_A} + {n_{{D_i}}}}} - {M_i}(t + 1)\varphi _i^{\rm{T}}(t + 1)]{\tilde \vartheta _i}(t) -\\ & {M_i}(t + 1){\varepsilon _i}(t + 1) \end{align} $$ (29) 于是, 任意两个局部参数估计之间的估计误差协方差阵$ {P_{{\vartheta _{ij}}}}(t) = {\rm{E}}[{\tilde \vartheta _i}(t)\tilde \vartheta _j^{\rm{T}}(t)] $可计算如下:
$$ \begin{align} {P_{{\vartheta _{ij}}}}(t + 1) = \, &[{I_{{n_A} + {n_{{D_i}}}}} - {M_i}(t + 1)\varphi _i^{\rm{T}}(t + 1)]{P_{{\vartheta _{ij}}}}(t) \\ &{[{I_{{n_A} + {n_{{D_j}}}}} - {M_j}(t + 1)\varphi _j^{\rm{T}}(t + 1)]^{\rm{T}}} +\\ & {M_i}(t + 1)\hat \sigma _{{\varepsilon _{ij}}}^2(t + 1)M_j^{\rm{T}}(t + 1) \end{align} $$ (30) 当$ i = j $时, $ {P_{{\vartheta _{ii}}}}(t) $即为局部参数估计误差方差阵$ {P_{{\vartheta _i}}}(t) $. $ {\varepsilon _i}(t) $与$ {\varepsilon _j}(t) $之间的互协方差$ \sigma _{{\varepsilon _{ij}}}^2 $可近似计算如下:
$$ \begin{equation} \hat \sigma _{{\varepsilon _{ij}}}^2(t) = \frac{1}{t}\sum\limits_{k = 1}^t {\hat \varepsilon _i^{}(k)\hat \varepsilon _j^{}(k)} \end{equation} $$ (31) 它可递推地计算为
$$ \begin{equation} \hat \sigma _{{\varepsilon _{ij}}}^2(t) = \hat \sigma _{{\varepsilon _{ij}}}^2(t - 1) + \frac{1}{t}[\hat \varepsilon _i^{}(t)\hat \varepsilon _j^{}(t) - \hat \sigma _{{\varepsilon _{ij}}}^2(t - 1)] \end{equation} $$ (32) 初值为$ \hat \sigma _{{\varepsilon _{ij}}}^2(0) = {y_i}(0){y_j}(0), i, j = 1, \cdots, L $.
令$ {\vartheta _A} = {[{a_1}, \cdots, {a_{{n_a}}}]^{\rm{T}}} $, 则有$ {\vartheta _{A_i}} = [{I_{{n_a}}}, 0]{\vartheta _i} $.于是, 我们有参数$ {\vartheta _A} $基于传感器$ i $的局部估计和估计误差协方差阵如下:
$$ \begin{equation} {\hat \vartheta _{Ai}}(t) = [{I_{{n_a}}}, 0]{\hat \vartheta _i}(t) \end{equation} $$ (33) $$ \begin{equation} {P_{{\vartheta _{Aij}}}}(t) = [{I_{{n_a}}}, 0]{P_{{\vartheta _{ij}}}}(t){[{I_{{n_a}}}, 0]^{\rm{T}}} \end{equation} $$ (34) 当$ i = j $时, $ {P_{{\vartheta _{Aii}}}}(t) $即为局部参数$ {\vartheta _A} $的估计误差方差阵$ {P_{{\vartheta _{Ai}}}}(t) $.
由式(18)可知, 参数$ {\vartheta _A} = {[{a_1}, \cdots, {a_{{n_a}}}]^{\rm{T}}} $是$ \Phi $中未知参数的函数.假设$ \Phi $中未知模型参数组成的列向量为$ {\Lambda ^{[{\rm{\Phi }}]}} \in {{\boldsymbol{\rm R}}^{{n_\Phi }}}, {n_\Phi } \le {n_a} $, 且可由$ {\vartheta _A} $唯一确定.设${\Lambda ^{[{\rm{\Phi }}]}}$与${\vartheta _A}$之间满足如下关系:
$$ \begin{equation} {\Lambda ^{[{\rm{\Phi }}]}} = f({\vartheta _A}) \end{equation} $$ (35) 其中$ f({\vartheta _A}) $为关于$ {\vartheta _A} $的线性或非线性函数.
1) 如果$ f({\vartheta _A}) $是线性函数, 我们将式(35)重写为
$$ \begin{equation} {\Lambda ^{{\rm{[\Phi }}]}} = S{\vartheta _A} + \gamma \end{equation} $$ (36) 式中$ S, \gamma $为适当维数的系数阵.
那么, 基于传感器$ i $的数据获得的$ \Phi $中未知模型参数在$ t $时刻的局部估计为
$$ \begin{equation} \hat \Lambda _i^{[{\rm{\Phi }}]}(t) = S{\hat \vartheta _{Ai}}(t) + \gamma \end{equation} $$ (37) 定义局部估计误差方差$ {P_{\Lambda _i^{[{\rm{\Phi }}]}}}(t) = {\rm{E}}[\tilde \Lambda _i^{[{\rm{\Phi }}]}(t){(\tilde \Lambda _i^{[{\rm{\Phi }}]}(t))^{\rm{T}}}] $, 任意两个局部参数估值器之间的估计误差互协方差为$ {P_{\Lambda _{ij}^{[{\rm{\Phi }}]}}}(t) = {\rm{E}}[\tilde \Lambda _i^{{\rm{[\Phi }}]}(t){(\tilde \Lambda _j^{[{\rm{\Phi]}}}(t))^{\rm{T}}}] $, 其中估计误差$ \tilde \Lambda _i^{[\Phi]}(t) = \Lambda _{}^{[\Phi]} - \hat \Lambda _i^{[\Phi]}(t) = S{\tilde \vartheta _{Ai}}(t) $, 则可获得
$$ \begin{equation} {P_{\Lambda _i^{[{\rm{\Phi }}]}}}(t) = S{P_{{\vartheta _{Ai}}}}(t){S^{\rm{T}}}, {P_{\Lambda _{ij}^{[{\rm{\Phi }}]}}}(t) = S{P_{{\vartheta _{Aij}}}}(t){S^{\rm{T}}} \end{equation} $$ (38) 2) 如果$ f({\vartheta _A}) $是非线性函数, 我们将$ f({\vartheta _A}) $在点$ {\hat \vartheta _{Ai}}(t - 1) $处进行线性化, 有
$$ \begin{equation} \Lambda _i^{[\Phi ]}(t) \approx S({\hat \vartheta _{Ai}}(t - 1)){\vartheta _A} + \gamma ({\hat \vartheta _{Ai}}(t - 1)) \end{equation} $$ (39) 那么, 基于传感器$ i $的数据获得的$ \Phi $中未知模型参数在$ t $时刻的局部估计为
$$ \begin{equation} \hat \Lambda _i^{[\Phi ]}(t) \approx S({\hat \vartheta _{Ai}}(t - 1)){\hat \vartheta _{Ai}}(t) + \gamma ({\hat \vartheta _{Ai}}(t - 1)) \end{equation} $$ (40) 类似地, 可以得到相应的估计误差方差$ {P_{\Lambda _i^{[{\rm{\Phi }}]}}}(t) $和互协方差$ P_{ij}^{{\Lambda ^{[{\rm{\Phi }}]}}}(t) $分别为
$$ \begin{align} &{P_{\Lambda _i^{[{\rm{\Phi }}]}}}(t) = S({\hat \vartheta _{Ai}}(t - 1)){{{P}}_{{\vartheta _{Ai}}}}(t){S^{\rm{T}}}({\hat \vartheta _{Ai}}(t - 1)) \\ &{P_{\Lambda _{ij}^{[{\rm{\Phi }}]}}}(t) = S({\hat \vartheta _{Ai}}(t - 1)){P_{{\vartheta _{Aij}}}}(t){S^{\rm{T}}}({\hat \vartheta _{Aj}}(t - 1)) \end{align} $$ (41) 其中$ i, j = 1, \cdots, L $.通过以上算法, 基于$ L $个传感器的数据可获得$ \Phi $中未知参数在时刻$ t $处的局部估值$ \hat \Lambda _i^{[\Phi]}(t), i, j = 1, \cdots, L $.由于对$ \Phi $中未知参数估计了$ L $次, 因此我们可应用线性无偏最小方差加权融合估计算法[23]将它们进行融合处理.下面定理1给出了未知参数分布式融合估计的结果.
定理1. 基于局部参数估计$ \hat \Lambda _i^{[{\rm{\Phi }}]}(t) $、局部参数估计误差方差阵$ {P_{\Lambda _i^{[{\rm{\Phi }}]}}}(t) $, 以及参数估计误差互协方差$ {P_{\Lambda _{ij}^{[{\rm{\Phi }}]}}}(t) $, 可得在线性无偏最小方差意义下的矩阵加权参数融合辨识器如下:
$$ \begin{equation} \hat \Lambda _o^{{\rm{[\Phi ]}}}(t) = \sum\limits_{i = 1}^L {{W_{\Lambda _i^{{\rm{[\Phi }}]}}}(t)\hat \Lambda _i^{[{\rm{\Phi }}]}(t)} \end{equation} $$ (42) 参数融合加权矩阵计算为
$$ \begin{align} [{W_{\Lambda _i^{[{\rm{\Phi }}]}}}(t), &\cdots , {W_{\Lambda _L^{[{\rm{\Phi }}]}}}(t)] = \\ &{(e_{{\Lambda ^{[\Phi ]}}}^{\rm{T}}P_{{\Lambda ^{[\Phi ]}}}^{ - 1}(t){e_{{\Lambda ^{[\Phi ]}}}})^{ - 1}}e_{{\Lambda ^{[\Phi ]}}}^{\rm{T}}P_{{\Lambda ^{[\Phi ]}}}^{ - 1}(t) \end{align} $$ (43) 其中$ {e_{{\Lambda ^{[\Phi]}}}} = {[{\begin{array}{*{20}{c}} {{I_{{n_\Phi }}}} & \cdots & {{I_{{n_\Phi }}}} \\ \end{array}}]^{\rm T}}_{{n_\Phi }L \times {n_\Phi }}, {P_{{\Lambda ^{[\Phi]}}}}(t) = {[{{P_{\Lambda _{ij}^{[{\rm{\Phi }}]}}}(t|t)}]_{{n_\Phi }L \times {n_\Phi }L}} $是以$ {P_{\Lambda _{ij}^{[{\rm{\Phi }}]}}}(t|t) $为第$ (i, j) $元素的分块矩阵.
参数融合辨识器的估计误差方差阵计算为
$$ \begin{equation} {P_{\Lambda _o^{[{\rm{\Phi }}]}}}(t) = {(e_{{\Lambda ^{[\Phi ]}}}^{\rm{T}}P_{{\Lambda ^{[\Phi ]}}}^{ - 1}(t){e_{{\Lambda ^{[\Phi ]}}}})^{ - 1}} \end{equation} $$ (44) 且$ {P_{\Lambda _o^{{\rm{[\Phi]}}}}}(t) \le {P_{\Lambda _i^{[{\rm{\Phi }}]}}}(t), i = 1, \cdots, L. $
基于以上算法可以获得未知模型参数融合辨识器$ \Lambda _o^{[{\rm{\Phi]}}}(t) $, 进而可获得$ \Phi $的融合估值$ {{\rm{\hat \Phi }}_o}(t) $.由前面的分析, 有$ {\hat \Phi _o}(t) \to \Phi, t \to \infty, w.p.1 $.
注2. 文献[5-7, 24]在进行模型参数融合辨识时, 将各传感器辨识得到的模型参数采用加权平均方法进行融合处理.该方法不能保证所获得的参数融合辨识器的估计精度不低于所有的局部参数估计.而本文采用线性无偏最小方差分布式矩阵加权融合算法[23]对各传感器辨识得到的参数进行融合处理.所获得的参数融合辨识器的估计精度不低于所有的局部参数估计.因此, 本文的矩阵加权融合的参数估计精度高于加权平均融合的参数估计精度.这在后面的仿真研究中也能看到.
4. $ \boldsymbol {\mu _i(t)} $的数学期望与方差辨识
当各传感器的衰减观测率未知时, 为了能应用第2节中的算法获得状态估计, 我们需要辨识描述衰减观测的随机变量$ \{ {\mu _i}(t)\} $的均值和方差.下面我们采用相关函数来辨识它们.
将模型参数融合辨识器$ {\hat \Phi _o}(t) $代入式(6)中, 得到
$$ \begin{equation} \hat X(t + 1) = {\hat \Phi _o}(t) \hat X(t)\hat \Phi _o^{\rm{T}}(t) + {\rm{\Gamma }}{Q_{\boldsymbol w}}{{\rm{\Gamma }}^{\rm{T}}} \end{equation} $$ (45) 由式(2)可计算零步相关函数$ R_i(t, 0) = {\rm E}[y_i^2(t)] $为
$$ \begin{equation} {R_i}(t, 0) = (\alpha _i^2 + \sigma _i^2){h_{i}}{{\hat X}}(t)h_{i}^{\rm{T}} + {{{Q}}_{{v_i}}} \end{equation} $$ (46) 一步相关函数$ {{{R}}_i}(t, {\rm{1}}) = {\rm{E}}[{y_i}(t){y_i}(t- {\rm{1}})] $为
$$ \begin{equation} {R_i}(t, 1) = \alpha _i^2{ h_{i}}{\hat \Phi _o}(t - 1)\hat X(t - 1) h_{i}^{\rm{T}} \end{equation} $$ (47) 零阶和一阶相关函数阵$ {{{R}}_i}(t, r), r = 0, {\rm{1}} $可通过如下采样相关函数$ {{{\hat R}}_i}(t, r) $逼近:
$$ \begin{align} {\hat R_i}(t, r) \approx &\frac{1}{t}\sum\limits_{s = 1}^t {{{ y}_i}(s){{y}_i}(s - r)} = {\hat R_i}(t - 1, r) + \\ & \frac{1}{t}[{y_i}(t){y_i}(t - r) - {\hat R_i}(t - 1, r)] \end{align} $$ (48) 初值为$ {{{\hat R}}_i}(0, 0) = {\rm E}[y_i^2{\rm{(0)]}}, \; {{{\hat R}}_i}(0, 1){\rm{ = 0}} $.
最后, 由式(47)可求出随机变量$ {\rm{\{ }}{\mu _i}{\rm{(}}t{\rm{)\} }} $的数学期望为
$$ \begin{equation} {\hat \alpha _i}(t) = \sqrt {\frac{{{{\hat R}_i}(t, 1)}}{{{{ h}_i}{{\hat \Phi }_o}(t - 1)\hat X(t - 1){ h}_i^{\rm{T}}}}} \end{equation} $$ (49) 将式(49)代入式(46)可得随机变量$ {\rm{\{ }}{\mu _i}{\rm{(}}t{\rm{)\} }} $的方差为
$$ \begin{equation} \hat \sigma _i^2(t) = \frac{{{{\hat R}_i}(t, 0) - {{{Q}}_{{v_i}}}}}{{ {h_{i}}{{\hat X}}(t)h_{i}^{\rm{T}}}} - \hat \alpha _i^2(t) \end{equation} $$ (50) 通过式(49)可获得$ {\alpha _i} $在$ t $时刻的估值.由第3节可知$ {\hat \Phi _o}(t) \to \Phi, t \to \infty $.由$ \Phi $的稳定性和式(45)可得$ \hat X(t) \to X(t), t \to \infty $.又由随机过程的遍历性可得$ {{{\hat R}}_i}(t, r) \to {{{R}}_i}(t, r), t \to \infty $.因而, 有估值$ {\hat \alpha _i}(t) $, $ \hat \sigma _i^2(t) $具有一致性, 即
$$ \begin{equation} {\hat \alpha _i}(t) \to {\alpha _i}, \hat \sigma _i^2(t) \to \sigma _i^2, t \to \infty , w.p.1 \end{equation} $$ (51) 进一步, 有$ {\hat Q_{{V_i}}}(t) = \hat \sigma _i^2(t){{{{ h}}}_i}\hat X(t){{ h}_i^{\rm{T}}} + {Q_{{v_i}}} \to {Q_{{V_i}}}(t) = \sigma _i^2{{{ h}}_i}X(t){{{ h}}_i^{\rm{T}}} + {Q_{{v_i}}}, t \to \infty, w.p.1 $.
将每时刻辨识的$ {\hat \Phi _o}(t) $、$ {\hat \alpha _i}(t) $和$ {\hat Q_{{V_i}}}(t) $替代第2节中的局部最优滤波器、互协方差阵和分布式融合滤波算法中的$ \Phi $、$ {\alpha _i} $和$ {Q_{{V_i}}}(t) $, 可得到相应的自校正状态估计算法.为了后文阐述方便, 记相应的自校正局部滤波器、预报误差方差、滤波误差协方差、增益、融合滤波器分别为$ {\hat {\boldsymbol x}_{si}}(t|t) $、$ {\hat \Sigma _{si}}(t|t - 1) $、$ {\hat P_{sij}}(t|t) $、$ {\hat K_{si}}(t) $和$ {\hat {\boldsymbol x}_s}(t|t) $.
5. 自校正滤波器的收敛性分析
下面的引理4给出了DESA方法, 应用DESA方法可证明自校正融合状态滤波器$ {\hat {\boldsymbol x}_s}(t|t) $的收敛性.
引理4[24]. 考虑动态误差系统
$$ \begin{equation} \delta (t) = T(t)\delta (t - 1) + u(t) \end{equation} $$ (52) 其中$ t \ge 0 $, 输出$ \delta (t) \in{{\bf R}}^n $, 输入$ u(t) \in {{\bf R}}^n $, 并且矩阵$ T(t) \in {{{\bf R}}^{n \times n}} $是一致渐近稳定的.若$ u(t) $是有界的, 则$ \delta {\rm{(}}t{\rm{)}} $是有界的.当$ t \to \infty $时, 若$ u(t) \to 0 $, 则$ \delta (t) \to 0 $.
记带已知参数$ \Phi $、$ {\alpha _i} $和$ \sigma _i^2 $的系统(1)和(3)为$ ({\rm{\Phi }}, {\rm{\Gamma}}, {F_i}, {Q_{\boldsymbol w}}, {Q_{{V_i}}}(t)) $, 相应的带时变参数$ {{\rm{\hat \Phi }}_o}(t) $、$ {\hat \alpha _i}(t) $和$ \hat \sigma_i^2(t) $的系统(1)和(3)为$ ({{\rm{\hat \Phi }}_o}(t), {\rm{\Gamma }}, {\hat F_i}(t), {Q_{\boldsymbol w}}, {\hat Q_{{V_i}}}(t)) $.
引理5[24]. 在假设1$ \, \sim\, $4下, 系统$ ({\rm{\Phi }}, {\rm{\Gamma }}, $ $ {F_i}, $ $ {Q_{\boldsymbol w}}, $ $ {Q_{{V_i}}}(t)) $和$ ({{\rm{\hat \Phi }}_o}(t), \Gamma, {\hat F_i}(t), {Q_{\boldsymbol w}}, {\hat Q_{{V_i}}}(t)) $是一致完全可观和一致完全可控的; 系统$ ({\rm{\Phi }}, \Gamma, {F_i}, {Q_{\boldsymbol w}}, {Q_{{V_i}}}(t)) $的最优局部滤波器的状态转移阵$ {\Psi _{fi}}(t) = [{I_n}-{K_i}(t){F_i}]{\rm{\Phi }} $是一致渐近稳定的, 且$ {K_i}(t) $是有界的; 系统$ ({{\rm{\hat \Phi }}_o}(t), \Gamma, {\hat F_i}(t), {Q_{\boldsymbol w}}, {\hat Q_{{V_i}}}(t)) $的自校正局部滤波器的状态转移阵$ {\hat \Psi _{sfi}}(t) = [{I_n}-{\hat K_{si}}(t){\hat F_i}{\rm{(}}t{\rm{)}}]{{\rm{\hat \Phi }}_o}(t) $是一致渐近稳定的, 且$ {\hat K_{si}}(t) $是有界的.而且$ \Delta {\hat K_i}(t) = {\hat K_{si}}(t)-{K_i}(t)\to 0 $和$ \Delta {\hat \Psi _{fi}}(t) = {\hat \Psi _{sfi}}(t)-{\Psi _{fi}}(t) \to 0 $, $ t \to\infty $.自校正局部滤波误差方差阵$ {\hat P_{si}}(t|t) $收敛于最优局部滤波误差方差阵$ {P_i}(t|t) $, 自校正滤波误差互协方差阵$ {\hat P_{sij}}(t|t) $收敛于最优滤波误差互协方差阵$ {P_{ij}}(t|t) $.
定理2. 在假设1$ \, \sim\, $4下, 自校正局部滤波器$ {\hat {\boldsymbol x}_{si}}(t|t) $收敛于最优局部滤波器$ {\hat {\boldsymbol x}_i}(t|t) $, 即
$$ \begin{equation} [{\hat {\boldsymbol x}_{si}}(t|t) - {\hat {\boldsymbol x}_i}(t|t)] \to 0, t \to \infty \end{equation} $$ (53) 证明. 由式(7)可得自校正滤波器为
$$ \begin{equation} {\hat {\boldsymbol x}_{si}}(t|t) = {\hat \Psi _{sfi}}(t){\hat {\boldsymbol x}_{si}}(t - 1|t - 1)+{\hat K_{si}}(t){y_i}(t) \end{equation} $$ (54) 注意到$ {\hat K_{si}}(t) $和$ {y_i}(t) $是有界的, $ {\hat \Psi _{sfi}}(t) $是一致渐近稳定的, 应用引理4得到$ {\hat {\boldsymbol x}_{si}}(t|t) $是有界的.令$ {\delta _i}(t) = {\hat {\boldsymbol x}_{si}}(t|t) - {\hat {\boldsymbol x}_i}(t|t) $, 式(54)减式(7)有动态误差方程
$$ \begin{equation} {\delta _i}(t) = {\Psi _{fi}}(t){\delta _i}(t - 1) + {u_i}(t) \end{equation} $$ (55) 其中$ {u_i}(t) = \Delta {\hat \Psi _{fi}}(t){\hat {\boldsymbol x}_{si}}(t - 1|t - 1) + \Delta {\hat K_i}(t){y_i}(t) $.根据$ {\hat {\boldsymbol x}_{si}}(t|t) $和$ { y_i}(t) $的有界性, 由$ \Delta {\hat K_i}(t)\to 0 $和$ \Delta{\hat \Psi _{fi}}(t){\to0} $, 有$ {u_i}(t)\to0 $.对式(55)应用引理4, 当$ t \to \infty $时, 有$ {\delta _i}(t) \to 0 $, 即式(53)成立.
定理3. 在假设1$ \, \sim\, $4下, 自校正加权融合滤波器$ {\hat {\boldsymbol x}_s}(t|t) $收敛于最优加权融合滤波器$ {\hat {\boldsymbol x}_o}(t|t) $, 即
$$ \begin{equation} [{\hat {\boldsymbol x}_s}(t|t) - {\hat {\boldsymbol x}_o}(t|t)] \to 0, t \to \infty \end{equation} $$ (56) 证明. 由引理3和引理5有$ {W_i}(t) $有界且$ \Delta {\hat W_i}(t) = [{\hat W_{si}}(t) -{W_i}(t)] \to 0 $.由式(53)以及$ {\hat {\boldsymbol x}_{si}}(t|t) $的有界性可得
$$ \begin{align} {\hat {\boldsymbol x}_s}(t|t) & - {\hat {\boldsymbol x}_o}(t|t) = \sum\limits_{i = 1}^L {{W_i}(t)[{{\hat {\boldsymbol x}}_{si}}(t|t) - {{\hat {\boldsymbol x}}_i}(t|t)} ]+ \\ &\sum\limits_{i = 1}^L {\Delta {{\hat W}_i}(t){{\hat {\boldsymbol x}}_{si}}(t|t)} \to 0 \end{align} $$ (57) 即式(56)成立.
6. 仿真例子
考虑带三传感器系统(1)和(2), 其中系数阵为$ \Phi = \left[{\begin{array}{*{20}{c}} {{a_{11}}} & {{a_{12}}} \\ {0.4} & { - 0.8} \\ \end{array}} \right], \Gamma = \left[{\begin{array}{*{20}{c}} {0.5} \\ {0.6} \\ \end{array}} \right], {{h}_1} = \left[{\begin{array}{*{20}{c}} {0.5} & {1.2} \\ \end{array}} \right] $, $ {{ h}_2} = \left[{\begin{array}{*{20}{c}} {0.6} & {1.9} \\ \end{array}}\right], {{ h}_3} = \left[{\begin{array}{*{20}{c}} {1.4} & 2 \\ \end{array}} \right], $噪声方差为$ {Q_{\boldsymbol w}} = 3, {Q_{{v_1}}} = 2, {Q_{{v_2}}} = 0.4, {Q_{{v_3}}} = 1, \left\{ {{\mu _i}(t)} \right\} $, $ i = 1, 2, 3 $为在$ [0, 1] $区间取值的标量随机变量.取初值$ {\hat{\boldsymbol x}_i}(0|0) = 0, {P_{ij}}(0|0) = 0.1{I_2}, i, j = 1, 2, 3 $.假设$ {a_{11}}, {a_{12}}, \{ {\mu _i}(t)\} $的数学期望$ {\alpha _i} $与方差$ \sigma _i^2, i = 1, 2, 3 $, 未知.目的是辨识未知模型参数$ {a_{11}} $和$ {a_{12}} $、期望$ {\alpha _i} $和方差$ \sigma _i^2, i = 1, 2, 3 $并求自校正融合状态滤波器.
在仿真中假设未知模型参数$ {a_{11}} = 0.6 $, $ {a_{12}} = -0.2 $, $ {\mu _i}(t), i = 1, 2, 3 $, 的概率分布分别为$ P{\rm{\{ }}{\mu _1}(t){\rm{ = 0.3\} = 0.3}} $, $ P{\rm{\{ }}{\mu _1}(t){\rm{ = 0.5\} = 0.2}} $, $ P{\rm{\{ }}{\mu _1}(t){\rm{ = 1\} = 0.5}} $, $ P{\rm{\{ }}{\mu_2}(t){\rm{ = 0.4\} = 0.4}} $, $ P{\rm{\{ }}{\mu_2}(t) {\rm{ = 0.7\} = 0.3}} $, $ P{\rm{\{ }}{\mu_2}(t){\rm{ = 0.9\} = 0.3}} $, $ P{\rm{\{ }}{\mu_3}(t){\rm{ = 0.1\}}} $ $ {{ = 0.2}} $, $ P{\rm{\{ }}{\mu_3}(t){\rm{ = 0.6\} = 0.6}} $, $ P{\rm{\{ }}{\mu_3}(t){\rm{ = 0.9\} = 0.2}} $.我们可以计算$ {\mu _i}(t) $, $ i = 1, 2, 3 $的数学期望和方差分别为$ {\alpha _1} = 0.69 $, $ {\alpha _2}{\rm{ = }}0.64 $, $ {\alpha_3} = 0.56 $, $ \sigma _1^2 = 0.1009 $, $ \sigma_2^2 = 0.0444 $, $ \sigma_3^2 = 0.0664 $.
为了与文献[5-7, 24]中的参数加权平均融合算法相比较.如下给出了模型参数$ {a_{11}} $和$ {a_{12}} $的局部估计误差方差、加权平均估计误差方差以及分布式加权融合估计误差方差算法.
a) 局部参数估计和分布式加权融合估计:
根据前面的第3节, 可知未知模型参数与$ {\vartheta _A} = {[{a_1}, {a_2}]^{\rm{T}}} $有如下线性关系
$$ \begin{equation} {\Lambda ^{[\Phi ]}} = S{\vartheta _A} + \gamma \end{equation} $$ (58) 其中$ {\Lambda ^{[\Phi]}} = \left[{\begin{array}{*{20}{c}} {{a_{11}}} \\ {{a_{12}}} \\ \end{array}} \right] $, $ S = \left[{\begin{array}{*{20}{c}} { - 1} & 0 \\ {{ {{ - {a_{22}}} \over {{a_{21}}}}}} & {{ {{ - 1} \over {{a_{21}}}}}} \\ \end{array}} \right] $, $ \gamma = \left[{\begin{array}{*{20}{c}} {{a_{22}}} \\ {{ {{a_{_{22}}^2} \over {{a_{21}}}}}} \\ \end{array}} \right]. $
${\Lambda ^{[{\rm{\Phi }}]}}$的局部估值为
$$ \begin{equation} \hat \Lambda _i^{[\Phi ]}(t) = S{\hat \vartheta _{Ai}}(t) + \gamma , i = 1, 2, 3 \end{equation} $$ (59) $ {\hat \vartheta _{Ai}}(t) $可由最小二乘算法辨识, 将$ S, {\hat \vartheta _{Ai}}(t), \gamma $代入式(38)中, 可分别获得局部估值$ \hat \Lambda _i^{[\Phi]}(t) $和局部估计误差方差$ {P_{\Lambda _i^{[\Phi]}}}(t) $以及误差互协方差$ {P_{\Lambda _{ij}^{[\Phi]}}}(t), i, j = 1, 2, 3 $.再根据定理1中参数融合估计算法, 可得参数融合估计误差方差$ {P_{\Lambda _0^{[{\rm{\Phi }}]}}}(t) $.
b) 参数的加权平均融合估计:
${\Lambda ^{[{\rm{\Phi }}]}}$的加权平均估值为
$$ \begin{equation} \bar \Lambda _{}^{[\Phi ]}(t) = { {1 \over 3}}\sum\limits_{i = 1}^3 {\hat \Lambda _i^{[\Phi ]}(t)} \end{equation} $$ (60) 定义加权平均估计误差方差为$ {P_{\bar \Lambda _{}^{[\Phi]}}}(t) = {\rm{E}}[\tilde {\bar \Lambda}_{}^{[\Phi]}(t)(\tilde {\bar \Lambda} _{}^{[\Phi]}(t)){^{\rm{T}}}] $, 其中$ \tilde {\bar \Lambda} _{}^{[\Phi]}(t) = {\Lambda ^{[\Phi]}} -\bar \Lambda _{}^{[\Phi]}(t) = { {1 \over 3}}\sum_{i = 1}^3 {(S{{\tilde \vartheta }_{Ai}}(t))} $.于是, 加权平均估计误差方差可计算为
$$ \begin{equation} {P_{\bar \Lambda _{}^{[\Phi ]}}}(t) = { {1 \over 9}}S\left(\sum\limits_{i = 1}^3 {\sum\limits_{j = 1}^3 ({P_{{\vartheta _{Aij}}}}}(t))\right){S^{\rm{T}}} \end{equation} $$ (61) 图 1给出了应用部分3参数融合辨识算法获得的未知模型参数融合辨识结果.由图可知随着时间的增长, 辨识结果收敛于真值.
图 2和图 3分别给出了未知模型参数$ {a_{11}} $和$ {a_{12}} $的局部估计误差方差、加权平均估计误差方差和分布式融合估计误差方差的比较结果.由图可知分布式加权融合辨识误差方差小于各局部辨识误差方差和加权平均辨识误差方差.图中$ {S_i}, i = 1, 2, 3 $表示第$ i $个传感器的局部辨识的误差方差, DWF表示分布式加权融合辨识的误差方差, WAEV表示加权平均融合辨识的误差方差.
图 4和图 5给出了应用部分4中辨识算法分别对不同传感器的随机变量$ \{ {\mu _i}(t)\}, i = 1, 2, 3 $的数学期望和方差进行辨识的结果.曲线表示辨识结果, 直线表示相应的真值.由图可知随着时间的增长, 辨识结果收敛于真值.图 6和图 7给出了自校正融合状态滤波器, 可见自校正融合估计具有有效性.
图 4 $\mu_{i}(t)$的数学期望辨识Fig. 4 Identification of Mathematical expectation of $\mu_{i}(t)$图 8和图 9给出了局部和融合的最优与自校正状态估计误差方差图.由图可见, 各局部自校正误差方差收敛于局部最优误差方差, 自校正融合误差方差收敛于最优融合误差方差, 即自校正滤波器具有渐近最优性.而且自校正融合滤波器比各局部自校正滤波器具有更高精度.图中$ {S_i}, i = 1, 2, 3 $表示第$ i $个传感器的局部自校正估计误差方差, SF表示自校正融合估计误差方差, 直线表示相应的最优方差.
图 8 局部、融合最优与自校正状态分量1的滤波误差方差Fig. 8 Variance of the first state component of local, fusion optimal and self-tuning filters
图 9 局部、融合最优与自校正状态分量2的滤波误差方差Fig. 9 Variance of the second state component of local, fusion optimal and self-tuning filters目前参考文献[5-7, 24]中的自校正滤波算法大都没有考虑传感器的衰减观测现象.图 10给出了传感器存在衰减观测而没有给予考虑的自校正融合滤波器与本文考虑衰减观测的自校正融合滤波器在30次蒙特卡洛实验下均方误差迹的比较.可见, 在传感器存在衰减观测时, 本文考虑衰减观测的自校正融合滤波器具有更高的精度.
图 10 考虑衰减观测与没有考虑衰减观测自校正融合滤波器的均方误差的迹Fig. 10 Trace of mean square error of the self-tuning fusion filters with/without considering fading measurements7. 结论
对带未知模型参数和衰减观测率的多传感器随机系统, 应用RELS算法和相关函数分别对未知模型参数、描述衰减观测现象的随机变量的数学期望和方差进行在线实时辨识, 提出了线性无偏最小方差矩阵加权融合模型参数辨识器.与已有文献的加权平均融合模型参数辨识算法相比, 本文所提出的线性无偏最小方差矩阵加权融合参数辨识算法具有更高的估计精度.将实时辨识的模型参数、数学期望和方差代入到最优局部和融合状态估计算法中获得了相应的自校正状态滤波算法.利用DESA方法证明了自校正状态滤波器收敛于最优状态滤波器.与现有文献的带未知模型参数的自校正估计算法相比, 本文还考虑了传感器的衰减观测现象, 并给出了采用相关函数辨识衰减观测的数学期望和方差的算法.
-
[1] 柴天佑.自动化科学与技术发展方向.自动化学报, 2018, 44(11): 1923-1930 doi: 10.16383/j.aas.2018.c180252Chai Tian-You. Development directions of automation science and technology. Acta Automatica Sinica, 2018, 44(11): 1923-1930 doi: 10.16383/j.aas.2018.c180252 [2] 王维洲, 吴志伟, 柴天佑.电熔镁砂熔炼过程带输出补偿的PID控制.自动化学报, 2018, 44(7): 1282-1292 doi: 10.16383/j.aas.2018.c170620Wang Wei-Zhou, Wu Zhi-Wei, Chai Tian-You. PID control with output compensation for the fused magnesia smelting process. Acta Automatica Sinica, 2018, 44(7): 1282-1292 doi: 10.16383/j.aas.2018.c170620 [3] 王兰豪, 贾瑶, 柴天佑.再磨过程的泵池液位和给矿压力双速率区间控制.自动化学报, 2017, 43(6): 993-1006 doi: 10.16383/j.aas.2017.c170134Wang Lan-Hao, Jia Yao, Chai Tian-You. Dual-rate interval control of pump pool level and feeding pressure during regrinding. Acta Automatica Sinica, 2017, 43(6): 993-1006 doi: 10.16383/j.aas.2017.c170134 [4] 杨辉, 郝丽娜, 陈洋, 薛帮灿.针对气动肌肉仿生肘关节抖振现象的Kalman-PID控制.控制理论与应用, 2017, 34(4): 477-482 https://www.cnki.com.cn/Article/CJFDTOTAL-KZLY201704008.htmYang Hui, Hao Li-Na, Chen Yang, Xue Bang-Can. Kalman-PID control for chattering phenomena of bionic elbow joint actuated by pneumatic artificial muscles. Control Theory & Applications, 2017, 34(4): 477-482 https://www.cnki.com.cn/Article/CJFDTOTAL-KZLY201704008.htm [5] Wang Z D, Xu Y, Liu G, Wang J X, Jiang W, Wang Y, et al. Simulation and experiment of PID applied to the automatic voltage control of gyrotron traveling wave tubes. IEEE Transactions on Plasma Science, 2018, 46(7): 2446-2451 doi: 10.1109/TPS.2017.2723245 [6] Zhang Y, Jia Y, Chai T Y, Wang D H, Dai W, Fu J. Data-driven PID controller and its application to pulp neutralization process. IEEE Transactions on Control Systems Technology, 2018, 26(3): 828-841 doi: 10.1109/TCST.2017.2695981 [7] 韩京清.自抗扰控制器及其应用.控制与决策, 1998, 13(1): 19-23 doi: 10.3321/j.issn:1001-0920.1998.01.005Han Jing-Qing. Active disturbance rejection controller and applications. Control and Decision, 1998, 13(1): 19-23 doi: 10.3321/j.issn:1001-0920.1998.01.005 [8] 韩京清.控制理论—模型论还是控制论.系统科学与数学, 1989, 9(4): 328-335 https://www.cnki.com.cn/Article/CJFDTOTAL-STYS198904005.htmHan Jing-Qing. Control theory, is it a model analysis approach or a direct control approach. Journal of Systems Science and Mathematical Sciences, 1989, 9(4): 328-335 https://www.cnki.com.cn/Article/CJFDTOTAL-STYS198904005.htm [9] 韩京清.自抗扰控制技术.北京:国防工业出版社, 2009, XIHan Jing-Qing. Active Disturbance Rejection Control Technology, Beijing: National Defense Industry Press, 2009, XI [10] 魏伟, 梅生伟, 张雪敏.先进控制理论在电力系统中的应用综述及展望.电力系统保护与控制, 2013, 41(12): 143-153 doi: 10.7667/j.issn.1674-3415.2013.12.023Wei Wei, Mei Sheng-Wei, Zhang Xue-Min. Review of advanced control theory and application in power system. Power System Protection and Control, 2013, 41(12): 143-153 doi: 10.7667/j.issn.1674-3415.2013.12.023 [11] 周延延, 吴晓燕, 李刚.基于BP神经网络的PID控制器研究.空军工程大学学报, 2007, 8(4): 45-48 doi: 10.3969/j.issn.1009-3516.2007.04.013Zhou Yan-Yan, Wu Xiao-Yan, Li Gang. Study on PID controller based on BP neural network. Journal of Air Force Engineering University, 2007, 8(4): 45-48 doi: 10.3969/j.issn.1009-3516.2007.04.013 [12] 王明涛, 张百浩.基于专家PID的燃气机转速控制试验.热能动力工程, 2015, 30(6): 848-852 https://www.cnki.com.cn/Article/CJFDTOTAL-RNWS201506011.htmWang Ming-Tao, Zhang Bai-Hao. Experiment of the rotating speed control over a gas engine based on an expert PID (Proportional, Integral and Differential) control method. Journal of Engineering for Thermal energy and power, 2015, 30(6): 848-852 https://www.cnki.com.cn/Article/CJFDTOTAL-RNWS201506011.htm [13] 李明, 封航, 张延顺.基于UMAC的RBF神经网络PID控制.北京航空航天大学学报, 2018, 44(10): 2063-2070 https://www.cnki.com.cn/Article/CJFDTOTAL-BJHK201810006.htmLi Ming, Feng Hang, Zhang Yan-Shu. RBF neural network tuning PID control based on UMAC. Journal of Beijing University of Aeronautics and Astronautics, 2018, 44(10): 2063-2070 https://www.cnki.com.cn/Article/CJFDTOTAL-BJHK201810006.htm [14] 王青山, 梁得亮, 杜锦华.交流稳压电源的改进神经网络PID控制.电机与控制学报, 2017, 21(2): 1-9 https://www.cnki.com.cn/Article/CJFDTOTAL-DJKZ201702001.htmWang Qing-Shan, Liang De-Liang, Du Jin-Hua. Improved neural network PID controller for regulated power supply. Electric Machines and Control, 2017, 21(2): 1-9 https://www.cnki.com.cn/Article/CJFDTOTAL-DJKZ201702001.htm [15] 邱占芝, 李世峰.基于神经网络的PID网络化控制系统建模与仿真.系统仿真学报, 2018, 30(4): 1423-1432 https://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ201804027.htmQiu Zhan-Zhi, Li Shi-Feng. Modeling and simulation of PID networked control systems based on neural network. Journal of System Simulation, 2018, 30(4): 1423-1432 https://www.cnki.com.cn/Article/CJFDTOTAL-XTFZ201804027.htm [16] 曾喆昭, 肖亚芬, 郝逸轩.非线性类PID神经元网络控制器.北京科技大学学报, 2012, 34(1): 12-15 https://www.cnki.com.cn/Article/CJFDTOTAL-BJKD201201002.htmZeng Zhe-Zhao, Xiao Ya-Fen, Hao Yi-Xuan. Nonlinear PID-like neuron network controller. Journal of Beijing University of Science and Technology, 2012, 34(1): 12-15 https://www.cnki.com.cn/Article/CJFDTOTAL-BJKD201201002.htm [17] 李桂梅, 曾喆昭.一种基于神经网络算法的非线性PID控制器.中南大学学报(自然科学版), 2010, 41(5): 1865-1870 https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201005038.htmLi Gui-Mei, Zeng Zhe-Zhao. A nonlinear PID controller based on neural network algorithm. Journal of Central South University (Science and Technology), 2010, 41(5): 1865-1870 https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201005038.htm [18] Shen J, Xin B. Control of single-axis rotation INS by tracking differentiator based fuzzy PID. IEEE Transactions on Aerospace and Electronic Systems, 2017, 53(6): 2976-2986 doi: 10.1109/TAES.2017.2722558 [19] 潘泽跃, 程健, 吴嘉珉, 程景行.基于FPGA广义预测比例—积分—微分控制在特种电源中的应用.电工技术学报, 2018, 33(10): 2376-2382 https://www.cnki.com.cn/Article/CJFDTOTAL-DGJS201810023.htmPan Ze-Yue, Cheng Jian, Wu Jia-Min, Cheng Jing-Xing. The application of generalized predictive proportion-integration-differentiation control in special power supply based on FPGA. Transactions of China Electrotechnical Society, 2018, 33(10): 2376-2382 https://www.cnki.com.cn/Article/CJFDTOTAL-DGJS201810023.htm [20] 谢宏, 杨鹏, 陈海滨, 张小刚, 陈俊辉, 谭阳红.遗传优化模糊PID融合算法的5自由度机械手控制.电子测量与仪器学报, 2015, 29(1): 21-30 https://www.cnki.com.cn/Article/CJFDTOTAL-DZIY201501003.htmXie Hong, Yang Peng, Chen Hai-Bin, Zhang Xiao-Gang, Chen Jun-Hui, Tan Yang-Hong. Fuzzy PID control system optimized by genetic algorithm for 5-freedom robot arm. Journal of Electronic Measurement and Instrumentation, 2015, 29(1): 21-30 https://www.cnki.com.cn/Article/CJFDTOTAL-DZIY201501003.htm [21] 王镇道, 张乐, 彭子舜.基于PSO优化算法的模糊PID励磁控制器设计.湖南大学学报(自然科学版), 2017, 44(8): 106-111 https://www.cnki.com.cn/Article/CJFDTOTAL-HNDX201708016.htmWang Zhen-Dao, Zhang Le, Peng Zi-Shun. Design of fuzzy PID excitation control based on PSO optimization algorithm. Journal of Hunan University (Natural Science), 2017, 44(8): 106-111 https://www.cnki.com.cn/Article/CJFDTOTAL-HNDX201708016.htm [22] Rabiatuladawiah A H, Siti F T, Salmiah A, Mohd K H. Swarm-intelligence tuned current reduction for power-assisted steering control in electric vehicles. IEEE Transactions on Industrial Electronics, 2018, 65(9): 7202-7210 doi: 10.1109/TIE.2017.2784344 [23] Marín A, Hernández R J A, Jiménez J A. Tuning multivariable optimal PID controller for a continuous stirred tank reactor using an evolutionary algorithm. IEEE Latin America Transactions, 2018, 16(2): 422-427 doi: 10.1109/TLA.2018.8327395 [24] 曾喆昭, 贺莹, 张畅, 李霖.非线性PID自学习控制方法研究.计算机工程, 2014, 40(10): 224-227 doi: 10.3969/j.issn.1000-3428.2014.10.042Zeng Zhe-Zhao, He Ying, Zhang Chang, Li Lin. Research on nonlinear PID control method for self-learning. Computer Engineering, 2014, 40(10): 224-227 doi: 10.3969/j.issn.1000-3428.2014.10.042 [25] 周勇, 曾喆昭.自学习非线性PID抗扰控制原理研究.控制工程, 2017, 24(6): 1180-1185 https://www.cnki.com.cn/Article/CJFDTOTAL-JZDF201706014.htmZhou Yong, Zeng Zhe-Zhao. Research on self-learning nonlinear PID disturbance rejection control principle. Control Engineering of China, 2017, 24(6): 1180-1185 https://www.cnki.com.cn/Article/CJFDTOTAL-JZDF201706014.htm [26] 曾喆昭, 吴亮东, 杨振源, 唐欢.非仿射系统的自学习滑模抗扰控制.控制理论与应用, 2016, 33(7): 980-987 https://www.cnki.com.cn/Article/CJFDTOTAL-KZLY201607017.htmZeng Zhe-Zhao, Wu Liang-Dong, Yang Zhen-Yuan, Tang Huan. Self-learning sliding-mode disturbance rejection control for non-affine systems. Control Theory & Applications, 2016, 33(7): 980-987 https://www.cnki.com.cn/Article/CJFDTOTAL-KZLY201607017.htm [27] Cheng Chun-Hua, Hu Yun-An, Wu Jin-Hua, Li Jing. Auto disturbance rejection controller for non-affine nonlinear systems with adaptive observers. Control Theory & Applications, 2014, 31(2): 148-158 -
下载:
下载:
计量
- 文章访问数: 2857
- HTML全文浏览量: 519
- PDF下载量: 617
- 被引次数: 0
