主权项 |
一种位置传感器和姿态传感器信息融合方法,其特征在于,将互补滤波与卡尔曼滤波相结合,将经过低通滤波后位置传感器测量的位置和速度信息、高通滤波后加速度估算的位置和速度信息作为卡尔曼滤波器的输入,位置传感器测量的位置信息和姿态传感器测量的加速度信息经过互补滤波与卡尔曼滤波后进行融合,得到融合后的位置信息和速度信息: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时刻融合后的位置信息。 |