发明名称 一种用于捷联惯性/多普勒组合导航系统的自适应滤波方法
摘要 本发明公开了一种用于捷联惯性/多普勒组合导航系统的自适应滤波方法,其目的是提高滤波器在动态条件下的响应速度和滤波精度,并提高捷联惯性/多普勒组合导航系统的定位精度。该方法通过引入关于新息协方差的限定窗口平滑器,利用该新息协方差平滑值直接修正滤波器中的增益矩阵,并引入调节因子对一步预测均方误差进行修正,进而达到改善滤波器动态响应速度以及提高滤波精度的目的。本发明所设计的自适应滤波方法用于捷联惯性/多普勒组合导航系统中,有效地提高了组合系统在动态条件下的导航定位精度。
申请公布号 CN103389095A 申请公布日期 2013.11.13
申请号 CN201310313649.4 申请日期 2013.07.24
申请人 哈尔滨工程大学 发明人 高伟;李敬春;陈嵩博;董会云;任晶;陈佳;陶晨斌;李敬国;刘适;尹冬寒
分类号 G01C21/20(2006.01)I;G06F19/00(2011.01)I 主分类号 G01C21/20(2006.01)I
代理机构 北京科亿知识产权代理事务所(普通合伙) 11350 代理人 汤东凤
主权项 1.一种用于捷联惯性/多普勒组合导航系统的自适应滤波方法,其特征在于,包括以下步骤:步骤一:利用捷联惯性导航系统中惯性测量单元测得沿载体系相对惯性空间的角速率和加速度分量信息,并与系统的初始位置、速度以及姿态信息一起进行捷联惯导解算,得到载体实时的位置、速度和姿态信息;步骤二:根据捷联惯导系统的误差方程选取状态变量,建立系统的状态方程,并选取经纬度误差、速度误差以及失准角作为状态变量;利用多普勒导航系统与捷联惯导系统所提供速度的差值作为量测变量,建立组合导航系统的量测方程;步骤三:经过连续系统离散化处理后,对组合导航系统进行滤波估计;首先更新状态一步预测值<img file="FDA00003561638300011.GIF" wi="126" he="89" />及其均方误差P<sub>k,k-1</sub>,并利用k时刻量测信息Z<sub>k</sub>和滤波器的一步预测值<img file="FDA00003561638300012.GIF" wi="152" he="88" />计算k时刻的新息序列ε<sub>k</sub>;再确定平滑窗口的宽度N,并借助量测新息序列ε<sub>k</sub>建立关于新息协方差的N步限定窗口平滑器<img file="FDA00003561638300013.GIF" wi="100" he="78" />所述的一步预测值:    <maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>X</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>&Phi;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mover><mi>X</mi><mo>^</mo></mover><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></math>]]></maths>所述的一步预测均方误差:    <maths num="0002"><![CDATA[<math><mrow><msub><mi>P</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>&Phi;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>P</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msubsup><mi>&Phi;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>T</mi></msubsup><mo>+</mo><msub><mi>&Gamma;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>Q</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msubsup><mi>&Gamma;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>T</mi></msubsup></mrow></math>]]></maths>所述的k时刻新息序列:<img file="FDA00003561638300016.GIF" wi="382" he="94" />所述的N步新息协方差平滑器:    <maths num="0003"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>k</mi><mi>N</mi></msubsup><mo>=</mo><msubsup><mi>C</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow><mi>N</mi></msubsup><mo>+</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><mo>[</mo><msub><mi>&epsiv;</mi><mi>k</mi></msub><msubsup><mi>&epsiv;</mi><mi>k</mi><mi>T</mi></msubsup><msub><mi>&epsiv;</mi><mrow><mi>k</mi><mo>-</mo><mi>N</mi></mrow></msub><msubsup><mi>&epsiv;</mi><mrow><mi>k</mi><mo>-</mo><mi>N</mi></mrow><mi>T</mi></msubsup><mo>]</mo></mrow></math>]]></maths>其中,P<sub>k-1</sub>表示k-1时刻的滤波均方误差,Q<sub>k-1</sub>表示k-1时刻系统噪声协方差,<img file="FDA00003561638300018.GIF" wi="94" he="87" />表示k-1时刻状态转移矩阵的转置,<img file="FDA00003561638300019.GIF" wi="85" he="88" />表示k-1时刻系统噪声系数矩阵的转置,ε<sub>k</sub>表示k时刻的新息序列,<img file="FDA000035616383000110.GIF" wi="76" he="79" />和<img file="FDA000035616383000111.GIF" wi="88" he="86" />分别是k和k-1时刻的新息协方差N步平滑值,N=10,且当k&lt;N时,<img file="FDA000035616383000112.GIF" wi="357" he="147" />步骤四:借助步骤三中得到的关于新息协方差的限定窗口平滑器<img file="FDA000035616383000113.GIF" wi="96" he="81" />对k时刻滤波器的增益矩阵K<sub>k</sub>进行修正,使得增益矩阵能够根据新息序列的变化直接做出调整;步骤五:利用步骤三中得到的限定窗口平滑器<img file="FDA000035616383000114.GIF" wi="100" he="86" />并引入调节因子ρ<sub>k</sub>,对滤波器的一步预测估计均方误差P<sub>k,k-1</sub>进行修正;步骤六:将步骤三至步骤五中得到的滤波参数带入到自适应滤波器中,并不断重复以上过程完成对导航参数误差的估计;利用输出校正实时修正捷联惯导系统输出的导航信息,从而得到更高精度的导航参数信息,即补偿后的位置、速度和姿态信息,直至捷联惯性/多普勒组合导航过程结束。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号
您可能感兴趣的专利