发明名称 一种捷联惯性导航系统双轴旋转调制方法
摘要 本发明属于捷联惯导系统调制技术领域,涉及一种捷联惯导系统双轴旋转调制方法,目的是是降低复杂性,实现无需高角加速度换向机构的速率偏频。它包括设定坐标系、安装陀螺和加速度计、安装惯性测量单元、双轴旋转调制、惯性测量单元导航计算和捷联惯导系统导航数据解调计算六个步骤。本发明采用双轴旋转调制,实现了对捷联惯导系统内惯性器件误差的调制,提高了捷联惯导系统的导航精度。对于使用激光陀螺的捷联惯导系统,本发明不需要高角加速度的换向机构,降低了系统实现的复杂度,进而降低了系统实现成本。
申请公布号 CN101900559B 申请公布日期 2013.07.03
申请号 CN200910210714.4 申请日期 2009.11.06
申请人 北京自动化控制设备研究所 发明人 邱宏波;练涛;崔鹏程;周章华;黄妍妍;张良军;陈明刚
分类号 G01C21/16(2006.01)I;G01C21/20(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅
主权项 1.一种捷联惯导系统双轴旋转调制方法,包括如下步骤:第一步:设定坐标系;第二步:安装陀螺和加速度计;将陀螺和加速度计安装在惯性测量单元上,使加速度计坐标系X<sup>a</sup>Y<sup>a</sup>Z<sup>a</sup>与惯性测量单元的坐标系X<sup>m</sup>Y<sup>m</sup>Z<sup>m</sup>重合,3个加速度计相互正交,X轴向加速度计(7),Y轴向加速度计(8)、Z轴向加速度计(9)分别安装在与加速度计坐标系X<sup>a</sup>Y<sup>a</sup>Z<sup>a</sup>的X<sup>a</sup>、Y<sup>a</sup>、Z<sup>a</sup>的3个轴平行的位置上;陀螺坐标系X<sup>g</sup>Y<sup>g</sup>Z<sup>g</sup>的X<sup>g</sup>、Y<sup>g</sup>、Z<sup>g</sup>的3个轴与惯性测量单元的坐标系X<sup>m</sup>Y<sup>m</sup>Z<sup>m</sup>的Y<sup>m</sup>轴的夹角都是54.74°;X轴向陀螺(10),Y轴向陀螺(11)、Z轴向陀螺(12)相互正交,X轴向陀螺(10),Y轴向陀螺(11)、Z轴向陀螺(12)分别安装在与陀螺坐标系X<sup>g</sup>Y<sup>g</sup>Z<sup>g</sup>的X<sup>g</sup>、Y<sup>g</sup>、Z<sup>g</sup>的3个轴平行的位置上;第三步:安装惯性测量单元将惯性测量单元安装在双轴旋转机构上,使得旋转机构的内轴(3)与惯性测量单元的坐标系X<sup>m</sup>Y<sup>m</sup>Z<sup>m</sup>的Y<sup>m</sup>轴平行;第四步:双轴旋转调制控制双轴旋转机构内轴(3)以恒定角速度连续旋转,角速率的范围为10°/s~1000°/s;控制外轴(4)进行周期性的往复运动,往复运动1圈到100圈,往复运动的周期为36s~36000s;第五步:惯性测量单元导航计算计算惯性测量单元在地理系的北向速度、东向速度、天向速度、经度、纬度、高度、姿态矩阵<img file="FDA00002963065100021.GIF" wi="112" he="76" />第六步:捷联惯导系统导航数据解调计算对捷联惯导系统导航数据进行解调计算,得到捷联惯导系统的航向角、俯仰角、滚动角、北向速度、东向速度、天向速度、经度、纬度和高度;包括如下步骤:(1)确定捷联惯导系统的北向速度、东向速度、天向速度、经度、纬度和高度;(2)确定捷联惯导系统的航向角、俯仰角和滚动角包括如下步骤:(a)计算捷联惯导系统在载体坐标系下的姿态矩阵捷联惯导系统载体坐标系的姿态矩阵为:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>n</mi><mi>b</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>cos</mi><msub><mi>&phi;</mi><mi>x</mi></msub></mtd><mtd><mn>0</mn></mtd><mtd><mi>sin</mi><msub><mi>&phi;</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><mi>sin</mi><msub><mi>&phi;</mi><mi>y</mi></msub><mi>sin</mi><msub><mi>&phi;</mi><mi>x</mi></msub></mtd><mtd><mi>cos</mi><msub><mi>&phi;</mi><mi>y</mi></msub></mtd><mtd><mi>sin</mi><msub><mi>&phi;</mi><mi>y</mi></msub><mi>cos</mi><msub><mi>&phi;</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><mi>cos</mi><msub><mi>&phi;</mi><mi>y</mi></msub><mi>sin</mi><msub><mi>&phi;</mi><mi>x</mi></msub></mtd><mtd><mo>-</mo><mi>sin</mi><msub><mi>&phi;</mi><mi>y</mi></msub></mtd><mtd><mi>cos</mi><msub><mi>&phi;</mi><mi>y</mi></msub><mi>cos</mi><msub><mi>&phi;</mi><mi>x</mi></msub></mtd></mtr></mtable></mfenced><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup></mrow></math>]]></maths>其中,<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>1,1</mn><mo>)</mo></mrow></mtd><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>1,2</mn><mo>)</mo></mrow></mtd><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>1,3</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>2,1</mn><mo>)</mo></mrow></mtd><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>2,2</mn><mo>)</mo></mrow></mtd><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>2,3</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>3,1</mn><mo>)</mo></mrow></mtd><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>3,2</mn><mo>)</mo></mrow></mtd><mtd><msubsup><mi>C</mi><mi>n</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mn>3,3</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>φ<sub>x</sub>双轴旋转机构的内轴转角;φ<sub>y</sub>为双轴旋转机构的外轴转角;<img file="FDA00002963065100024.GIF" wi="77" he="71" />为惯性测量单元的姿态矩阵;(b)计算捷联惯导系统的航向角、俯仰角和滚动角捷联惯导系统的航向角φ、俯仰角θ和滚动角γ如下:<maths num="0003"><![CDATA[<math><mrow><mi>&phi;</mi><mo>=</mo><msup><mi>tan</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msubsup><mi>C</mi><mi>n</mi><mi>b</mi></msubsup><mrow><mo>(</mo><mn>1,3</mn><mo>)</mo></mrow></mrow><mrow><msubsup><mi>C</mi><mi>n</mi><mi>b</mi></msubsup><mrow><mo>(</mo><mn>1,1</mn><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>&theta;</mi><mo>=</mo><msup><mi>sin</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msubsup><mi>C</mi><mi>n</mi><mi>b</mi></msubsup><mrow><mo>(</mo><mn>1,2</mn><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mi>&gamma;</mi><mo>=</mo><msup><mi>tan</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msubsup><mi>C</mi><mi>n</mi><mi>b</mi></msubsup><mrow><mo>(</mo><mn>3,2</mn><mo>)</mo></mrow></mrow><mrow><msubsup><mi>C</mi><mi>n</mi><mi>b</mi></msubsup><mrow><mo>(</mo><mn>2,2</mn><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>.</mo></mrow></math>]]></maths>
地址 100074 北京市丰台区云岗北区西里1号院