发明名称 基于匹配技术的地磁辅助惯性的载体姿态测量方法
摘要 一种基于匹配技术的地磁辅助惯性的载体姿态测量方法,包括:1)从地磁数据库中读出载体所在位置的地磁矢量信息,建立空间直角坐标系,计算出每个姿态下磁传感器的理论值;2)在一定时间内,从惯性器件得到运动载体的N个测量姿态,从磁强计得到N个磁场测量值;3)根据步骤2得到的载体N个测量姿态,在步骤1所建立的姿态角坐标系中,找到对应的N个点,结合磁传感器N个磁场测量值,采用匹配算法进行匹配,得到载体姿态的匹配结果;4)把匹配结果作为载体的真实姿态来修正惯性器件的输出结果;5)在下一测量时刻,重复执行步骤2-4,实现连续测姿。本发明能够有效的矫正惯性器件的积累误差,实现长时间、高精度的载体姿态测量。
申请公布号 CN103196445B 申请公布日期 2015.12.02
申请号 CN201310063251.X 申请日期 2013.02.07
申请人 哈尔滨工业大学 发明人 解伟男;屈桢深;夏红伟;白俊林
分类号 G01C21/16(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 代理人
主权项 一种基于匹配技术的地磁辅助惯性的载体姿态测量方法,其特征在于包括如下步骤:步骤1:从地磁数据库中读出载体所在位置的地磁矢量信息,以载体的滚转角、偏航角和俯仰角为坐标轴建立空间直角坐标系,并计算出每个姿态下磁强计的理论值;按照如下步骤建立以载体的滚转角、偏航角和俯仰角为坐标轴的空间直角坐标系,首先从地磁数据库中读出载体所在地理位置的地磁矢量信息M,由(1)式所示,下标xyz分别指向北东地方向:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>M</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>M</mi><mi>X</mi></msub></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>Y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>Z</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000140301180000011.GIF" wi="1103" he="202" /></maths>以载体的滚转角γ、俯仰角θ和偏航角<img file="FSB0000140301180000012.GIF" wi="42" he="46" />为坐标轴建立空间直角坐标系,由(2)式计算载体在任意姿态时,与载体捷联的磁强计的理论测量值B:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><mi>L</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mi>L</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>M</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>M</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000140301180000013.GIF" wi="1159" he="213" /></maths>其中,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000140301180000014.GIF" wi="1143" he="198" /></maths><img file="FSB0000140301180000015.GIF" wi="1468" he="202" />L(0)把地磁数据库中磁场值从北东地坐标系转换到东北天坐标系,即地理坐标系;L(1)是方向余弦矩阵,用于把地理坐标系下的地磁分量转换成载体坐标系下的磁场分量;L(2)为三维行向量,表示磁强计在载体坐标系下的方位角,用于把地理坐标系下的地磁分量转换到磁强计的轴线上,其具体数值在磁强计装配好后,对装配方向进行标定得到;步骤2:在t‑N+1到t的N个时刻,从惯性器件得到运动载体的N个测量姿态A,从磁强计得到N个磁场测量值M<sub>b</sub>,如下所示:<img file="FSB0000140301180000016.GIF" wi="1378" he="228" />M<sub>b</sub>=[m<sub>1</sub> m<sub>2</sub> … m<sub>N</sub>]                 (6)步骤3:根据步骤2得到的载体N个测量姿态,在步骤1所建立的空间直角坐标系中,找到对应的N个点,并将该N个点作为参考点集,结合磁强计的N个磁场测量值,采用匹配算法进行匹配,得到载体姿态的匹配结果;步骤4:把匹配结果作为载体的真实姿态来修正惯性器件的输出结果;步骤5:在下一测量时刻t=t+1,重复执行步骤2~4,从而实现连续测姿。
地址 150001 黑龙江省哈尔滨市南岗区一匡街2号哈尔滨工业大学科学园E2栋空间控制与惯性技术研究中心