发明名称 一种基于惯性和地磁传感器的目标姿态跟踪方法
摘要 本发明公开了一种基于惯性和地磁传感器的目标姿态跟踪方法,无需任何参考标志物和特定跟踪环境即可实现姿态跟踪,方法简单易行;采用陀螺仪、加速度计以及地磁传感器分别采集目标当前姿态所对应的角速度、加速度和磁强度数据在三个敏感轴上的分量,利用加速度计和磁传感器校正陀螺仪的姿态跟踪结果,消除了漂移误差,从而提高跟踪结果的精度;本发明使用高效离线运行的卡尔曼滤波算法,融合了多传感器跟踪结果,实现了实时在线姿态跟踪。
申请公布号 CN104764451A 申请公布日期 2015.07.08
申请号 CN201510196626.9 申请日期 2015.04.23
申请人 北京理工大学 发明人 刘越;贺长宇;闫达远;常军;翁冬冬
分类号 G01C21/16(2006.01)I;G01C21/08(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 李微微;仇蕾安
主权项 一种基于惯性和地磁传感器的目标姿态跟踪方法,其特征在于,采用陀螺仪、加速度计以及地磁传感器分别采集目标当前姿态所对应的角速度、加速度和磁强度数据在三个敏感轴上的分量,并采用卡尔曼滤波器对上述三个传感器数据进行融合并计算目标的姿态信息,具体方法如下:步骤1、通过卡尔曼滤波器对目标姿态信息进行时间更新,具体为:S10、在时间更新方程中对目标的状态矢量<img file="FDA0000704552360000011.GIF" wi="58" he="75" />进行估计:<img file="FDA0000704552360000012.GIF" wi="218" he="73" />其中,A为所述陀螺仪测得的目标在三个敏感轴上角速度测量值构成的滤波器增益矩阵;<img file="FDA0000704552360000013.GIF" wi="63" he="75" />为当前时刻目标状态矢量估计值,<img file="FDA0000704552360000014.GIF" wi="90" he="76" />为卡尔曼滤波器上一时刻输出的目标状态矢量;S11、对卡尔曼滤波器的传递函数的估计值<img file="FDA0000704552360000015.GIF" wi="60" he="67" />进行更新:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>p</mi><mi>k</mi><mo>-</mo></msubsup><mo>=</mo><msub><mi>Ap</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msup><mi>A</mi><mi>T</mi></msup><mo>+</mo><mi>Q</mi></mrow>]]></math><img file="FDA0000704552360000016.GIF" wi="323" he="70" /></maths>其中,Q值为系统噪声协方差矩阵;A<sup>T</sup>为增益矩阵A的转置矩阵;P<sub>k‑1</sub>为卡尔曼滤波器上一时刻输出的传递函数;步骤2、通过卡尔曼滤波器对目标姿态信息进行测量更新,具体为:S20、基于步骤1中更新得到的传递函数的估计值<img file="FDA0000704552360000017.GIF" wi="92" he="76" />获得卡尔曼滤波器的测量置信参数K<sub>k</sub>:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>K</mi><mi>k</mi></msub><mo>=</mo><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup><msup><mi>H</mi><mi>T</mi></msup><msup><mrow><mo>(</mo><msubsup><mi>HR</mi><mi>k</mi><mo>-</mo></msubsup><msup><mi>H</mi><mi>T</mi></msup><mo>+</mo><mi>R</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup></mrow>]]></math><img file="FDA0000704552360000018.GIF" wi="453" he="71" /></maths>其中,H为测量模型与估计模型间转换矩阵:<img file="FDA0000704552360000019.GIF" wi="500" he="170" />[q<sub>1</sub> q<sub>2</sub> q<sub>3</sub> q<sub>4</sub>]表示组成目标的状态矢量估计值<img file="FDA00007045523600000110.GIF" wi="61" he="76" />的旋转四元数,即<img file="FDA00007045523600000111.GIF" wi="387" he="75" />R为测量误差模型:<img file="FDA00007045523600000112.GIF" wi="607" he="77" />其中ACC<sub>x</sub>、ACC<sub>y</sub>和ACC<sub>z</sub>分别为加速度传感器测得的目标在三个敏感轴上的加速度分量,g为重力加速度;S21、对目标的状态矢量的测量结果<img file="FDA0000704552360000021.GIF" wi="62" he="83" />进行更新:<img file="FDA0000704552360000022.GIF" wi="445" he="77" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Z</mi><mi>k</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&theta;</mi></mtd></mtr><mtr><mtd><mi>&gamma;</mi></mtd></mtr><mtr><mtd><mi>&psi;</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000704552360000023.GIF" wi="161" he="195" /></maths>表示目标的姿态,其中,θ=arcsin(ACC<sub>x</sub>)表示目标的俯仰角,<img file="FDA0000704552360000024.GIF" wi="326" he="116" />表示目标的横滚角,<img file="FDA0000704552360000025.GIF" wi="240" he="115" />表示目标的航向角,A=MAG<sub>x</sub>cos(θ)+MAG<sub>y</sub>sin(θ)+MAG<sub>z</sub>cos(θ)sin(γ),B=MAG<sub>y</sub>cos(γ)+MAG<sub>z</sub>sin(γ),MAG<sub>x</sub>、MAG<sub>y</sub>和MAG<sub>z</sub>分别表示磁场传感器获得的目标所在位置的地球磁场在体坐标系的三个敏感轴上的分量;S22、对传递函数的测量结果P<sub>k</sub>进行更新:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>P</mi><mi>k</mi></msub><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>K</mi><mi>l</mi></msub><mi>H</mi><mo>)</mo></mrow><msubsup><mi>P</mi><mi>k</mi><mo>-</mo></msubsup></mrow>]]></math><img file="FDA0000704552360000026.GIF" wi="349" he="77" /></maths>步骤3、目标状态矢量的测量结果<img file="FDA0000704552360000027.GIF" wi="62" he="85" />即为当前目标姿态的跟踪结果,将当前时刻测量结果<img file="FDA0000704552360000028.GIF" wi="70" he="81" />和传递函数的测量结果P<sub>k</sub>输出并作为下一时刻滤波器的输入结果,返回步骤1,进行下一时刻的目标姿态跟踪。
地址 100081 北京市海淀区中关村南大街5号