发明名称 一种位置传感器和姿态传感器的导航系统及其融合方法
摘要 本发明公开了一种位置传感器和姿态传感器的导航系统及其融合方法,将互补滤波与卡尔曼滤波相结合,简称互补卡尔曼滤波,将经过低通滤波后GPS/气压计测量的位置和速度信息,高通滤波后加速度估算的位置和速度信息作为卡尔曼滤波器的输入,通过卡尔曼滤波自动调节互补滤波的权值,有效的解决手动调节互补滤波的权值的难题,实现将GPS位置信息与姿态传感器的加速度信息融合输出横向位置与速度信息,将气压计高度信息与姿态传感器输出的z轴加速度进行融合并输出z轴位置和速度信息。融合后的输出的位置和速度信息和融合前普通GPS和气压计输出的位置信息相比更加精确,并且能提供较精确的速度信息。
申请公布号 CN105865453A 申请公布日期 2016.08.17
申请号 CN201610340135.1 申请日期 2016.05.20
申请人 南京航空航天大学 发明人 刘翎予;王新华;李腾
分类号 G01C21/20(2006.01)I 主分类号 G01C21/20(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 朱小兵
主权项 一种位置传感器和姿态传感器信息融合方法,其特征在于,将互补滤波与卡尔曼滤波相结合,将经过低通滤波后位置传感器测量的位置和速度信息、高通滤波后加速度估算的位置和速度信息作为卡尔曼滤波器的输入,位置传感器测量的位置信息和姿态传感器测量的加速度信息经过互补滤波与卡尔曼滤波后进行融合,得到融合后的位置信息和速度信息:1)融合后的速度信息通过以下方程计算得到:v(k|k‑1)=v(k‑1|k‑1)+a<sub>g</sub>(k)*dtP<sub>1</sub>(k|k‑1)=P<sub>1</sub>(k‑1|k‑1)+Q<maths num="0001"><math><![CDATA[<mrow><msub><mi>K</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>P</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>R</mi><mn>1</mn></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000995263840000011.GIF" wi="475" he="146" /></maths><maths num="0002"><math><![CDATA[<mrow><mi>v</mi><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mi>z</mi><mo>+</mo><msub><mi>K</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>*</mo><mi>v</mi><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mi>z</mi><mo>+</mo><msub><mi>K</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>*</mo><msub><mi>v</mi><mrow><mi>b</mi><mi>e</mi><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000995263840000012.GIF" wi="1133" he="150" /></maths>式中,v(k|k‑1)为根据加速度信息计算得到的k时刻速度信息,v(k‑1|k‑1)为k‑1时刻融合后的速度信息,a<sub>g</sub>(k)为k时刻姿态传感器测量的加速度信息,dt为k‑1时刻和k时刻之间的时间间隔,P<sub>1</sub>(k|k‑1)为k时刻根据加速度信息计算得到的速度信息的协方差,P<sub>1</sub>(k‑1|k‑1)为k‑1时刻根据加速度信息计算得到的速度信息的协方差,Q为加速度信息的协方差,K<sub>1</sub>(k)为第一卡尔曼滤波增益,v<sub>before</sub>(k)为k时刻位置传感器测量的速度信息,R<sub>1</sub>为v<sub>before</sub>(k)的协方差;v(k|k)为k时刻融合后的速度信息;2)融合后的位置信息通过以下方程计算得到:<maths num="0003"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>L</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>v</mi><mrow><mi>a</mi><mi>f</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>*</mo><mi>d</mi><mi>t</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>*</mo><msub><mi>a</mi><mi>g</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>*</mo><mi>d</mi><mi>t</mi><mo>*</mo><mi>d</mi><mi>t</mi></mrow>]]></math><img file="FDA0000995263840000013.GIF" wi="1196" he="116" /></maths>P<sub>2</sub>(k|k‑1)=P<sub>2</sub>(k‑1|k‑1)+Q<maths num="0004"><math><![CDATA[<mrow><msub><mi>K</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msub><mi>R</mi><mn>2</mn></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000995263840000014.GIF" wi="486" he="149" /></maths><maths num="0005"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><msub><mi>K</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mi>z</mi><mo>+</mo><msub><mi>K</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>*</mo><mi>L</mi><mrow><mo>(</mo><mi>k</mi><mo>|</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><mi>z</mi><mo>+</mo><msub><mi>K</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>*</mo><msub><mi>L</mi><mrow><mi>b</mi><mi>e</mi><mi>f</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000995263840000015.GIF" wi="1167" he="150" /></maths>式中,L(k|k‑1)为根据加速度信息计算得到的k时刻位置信息,L(k‑1|k‑1)为k‑1时刻融合后的位置信息,v<sub>after</sub>(k‑1)为k‑1时刻融合后的速度信息,P<sub>2</sub>(k|k‑1)为k时刻根据加速度信息计算得到的位置信息的协方差,P<sub>2</sub>(k‑1|k‑1)为k‑1时刻根据加速度信息计算得到的位置信息的协方差,K<sub>2</sub>(k)为第二卡尔曼滤波增益,L<sub>before</sub>(k)为k时刻位置传感器测量的位置信息,R<sub>2</sub>为L<sub>bef</sub>o<sub>re</sub>(k)的协方差;L(k|k)为k时刻融合后的位置信息。
地址 210016 江苏省南京市秦淮区御道街29号