发明名称 一种模糊自适应变分贝叶斯无迹卡尔曼滤波方法
摘要 一种模糊自适应变分贝叶斯无迹卡尔曼滤波方法。本发明包括估计一步预测目标状态<img file="2014103650444100004dest_path_image002.GIF" wi="72" he="24" />及其协方差阵<img file="dest_path_image004.GIF" wi="72" he="22" />;迭代估计测量噪声的方差<img file="dest_path_image006.GIF" wi="37" he="24" />;计算当前时刻的残差方差阵的真实值<img file="dest_path_image008.GIF" wi="64" he="25" />,估计值<img file="dest_path_image010.GIF" wi="37" he="26" />,匹配程度指标<img file="dest_path_image012.GIF" wi="48" he="25" />,调整量<img file="dest_path_image014.GIF" wi="32" he="25" />以及调整后的测量噪声方差<img file="dest_path_image016.GIF" wi="22" he="28" />;计算目标状态的估计<img file="dest_path_image018.GIF" wi="54" he="24" />及其误差协方差<img file="dest_path_image020.GIF" wi="56" he="22" />。本发明提出的方法具备在线估计测量噪声的统计方差能力,降低了由于噪声统计特性未知时变而导致的滤波误差,提高了非线性滤波估计精度。与此同时,利用基于新息的协方差匹配技术的模糊逻辑方法,实时调整变分贝叶斯方法估计的量测噪声方差,抑制滤波器的发散,增强了滤波方法的鲁棒性。
申请公布号 CN104112079A 申请公布日期 2014.10.22
申请号 CN201410365044.4 申请日期 2014.07.29
申请人 洛阳理工学院 发明人 王国勇;王剑;李冠峰;李明照;崔文;孙昭峰;王帆;张红霞
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 洛阳公信知识产权事务所(普通合伙) 41120 代理人 苗强
主权项 一种模糊自适应变分贝叶斯无迹卡尔曼滤波方法,其特征在于:步骤1 设置滤波初始条件,具体包括:(1.1) 初始状态及其协方差阵<img file="2014103650444100001dest_path_image001.GIF" wi="98" he="24" />;(1.2) 模糊逻辑方法中的移动窗口大小W;(1.3) VB迭代次数N,初始化参数<img file="2014103650444100001dest_path_image002.GIF" wi="14" he="22" />,<img file="2014103650444100001dest_path_image003.GIF" wi="31" he="28" />,<img file="2014103650444100001dest_path_image004.GIF" wi="32" he="25" />;步骤2 进行滑动窗口内的变分贝叶斯无迹卡尔曼滤波(VB‑UKF)计算,具体包括:(2.1) 设置循环控制变量<img file="2014103650444100001dest_path_image005.GIF" wi="9" he="18" />的初值,令<img file="2014103650444100001dest_path_image006.GIF" wi="112" he="19" />,迭代循环开始;(2.2) 时间更新,估计一步预测目标状态<img file="2014103650444100001dest_path_image007.GIF" wi="61" he="21" />及其协方差阵<img file="2014103650444100001dest_path_image008.GIF" wi="65" he="22" />,其中,<img file="dest_path_image009.GIF" wi="42" he="20" />表示用<img file="2014103650444100001dest_path_image010.GIF" wi="30" he="18" />时刻的目标信息估计第<img file="16092dest_path_image005.GIF" wi="9" he="18" />时刻的目标信息;(2.3) 量测更新,具体包括:(2.3.1) 计算测量值的预测估计值<img file="dest_path_image011.GIF" wi="60" he="26" />;(2.3.2) 计算状态和测量值的互协方差矩阵<img file="dest_path_image012.GIF" wi="73" he="24" />;(2.3.3) 采用变分贝叶斯方法迭代计算测量噪声方差阵<img file="dest_path_image013.GIF" wi="29" he="26" />,均方根新息协方差阵<img file="dest_path_image014.GIF" wi="74" he="20" />,增益阵<img file="dest_path_image015.GIF" wi="29" he="19" />、最优估计<img file="dest_path_image016.GIF" wi="44" he="20" />及其误差协方差<img file="dest_path_image017.GIF" wi="45" he="19" />,迭代过程如下:(I) 设置循环控制变量<i>t</i>的初值,令<i>t</i>=1,根据迭代次数N的值,迭代循环开始;(II) 计算测量噪声方差阵<img file="dest_path_image019.GIF" wi="39" he="26" />,其中上标<i>t</i>表示第<i>t</i>次迭代时的值;(III) 计算均方根新息协方差阵<img file="26380dest_path_image021.GIF" wi="73" he="24" />和增益阵<img file="dest_path_image023.GIF" wi="38" he="23" />;(IV) 计算第t次迭代的最优估计<img file="dest_path_image025.GIF" wi="49" he="21" />及其误差协方差<img file="dest_path_image027.GIF" wi="49" he="24" />;(V) 如果<img file="dest_path_image029.GIF" wi="38" he="21" />,令<i>t</i>=<i>t+</i>1,然后返回(II),否则执行(VI);(VI) 结束VB迭代过程,求得滑动窗口内第<img file="143985dest_path_image005.GIF" wi="9" he="18" />次VB‑UKF滤波估计结果:<img file="dest_path_image031.GIF" wi="89" he="23" />,<img file="dest_path_image033.GIF" wi="111" he="23" />,<img file="dest_path_image035.GIF" wi="118" he="25" />;(2.4) 如果<img file="dest_path_image037.GIF" wi="33" he="16" />,令<img file="dest_path_image039.GIF" wi="53" he="18" />,然后返回步骤2,否则结束滑动窗口内VB‑UKF计算,执行步骤3;步骤3利用模糊逻辑方法动态调整<img file="dest_path_image041.GIF" wi="33" he="23" />,具体过程包括:(3.1) 计算当前时刻的残差方差阵的真实值<img file="dest_path_image043.GIF" wi="62" he="24" />,估计值<img file="598362dest_path_image045.GIF" wi="33" he="23" />以及匹配程度指标<img file="457734dest_path_image047.GIF" wi="43" he="23" />;(3.2) 以<img file="66569dest_path_image047.GIF" wi="43" he="23" />为模糊推理系统(FIS)的输入参数,计算测量噪声方差阵的调整量<img file="50181dest_path_image049.GIF" wi="29" he="23" />,并以此对VB方法估计的测量噪声方差<img file="123180dest_path_image051.GIF" wi="29" he="22" />进行调整,获得当前时刻测量噪声方差<img file="91136dest_path_image053.GIF" wi="19" he="24" />;步骤4 代回调整后的测量噪声方差<img file="317dest_path_image053.GIF" wi="19" he="24" />,在当前时刻执行一次UKF算法,获得目标状态的估计<img file="dest_path_image055.GIF" wi="52" he="23" />及其误差协方差<img file="dest_path_image057.GIF" wi="52" he="20" />。
地址 471000 河南省洛阳市高新技术产业开发区芳泽路