发明名称 无陀螺捷联惯导系统的初始姿态快速测量方法
摘要 本发明提供的是一种无陀螺捷联惯导系统初始姿态快速测量方法。将低成本无人飞行器上的无陀螺捷联惯性导航系统作为子惯导系统,低成本无人飞行器装载机构上的捷联惯性导航系统作为主惯导系统。利用主惯导输出的速度参考信息进行滤波,估测出主惯导系统与子惯导系统之间的水平方向安装偏差,再利用无陀螺捷联惯性导航系统中的加速度计输出信息对主惯导系统与子惯导系统之间的方位安装偏差进行匹配测量。本发明能快速测量低成本无人飞行器的初始姿态,从而提高低成本无人飞行器的快速反应能力,具有实用价值,本发明适用于低成本无人飞行器配备的中高精度无陀螺捷联惯导系统。
申请公布号 CN101694389A 申请公布日期 2010.04.14
申请号 CN200910073074.7 申请日期 2009.10.20
申请人 哈尔滨工程大学 发明人 于飞;奔粤阳;高伟;孙枫;周广涛;李倩
分类号 G01C21/16(2006.01)I 主分类号 G01C21/16(2006.01)I
代理机构 代理人
主权项 1.一种无陀螺捷联惯导系统的初始姿态快速测量方法,其特征是包括如下步骤:步骤1、将低成本无人飞行器上的无陀螺捷联惯性导航系统与高精度的低成本无人飞行器装载机构上的捷联惯性导航系统通过数据传输电缆相联通;步骤2、无陀螺捷联惯性导航系统进行预热后,采集加速度计输出的数据;步骤3、利用主惯导系统将初始速度参数和包括初始的经度、纬度的初始位置参数装订至无陀螺捷联惯性导航系统导航计算机中;步骤4、将主惯导系统输出的包括纵摇角、横摇角和航向角的姿态测量信息通过数据传输电缆传输至无陀螺捷联惯性导航系统,粗略测量出低成本无人飞行器的初始姿态,完成初始姿态的一次快速传递;步骤5、无陀螺捷联惯性导航系统由步骤3和步骤4得到的姿态初始值,利用无陀螺捷联惯性导航系统上加速度计输出,测量低成本无人飞行器的比力信息,再由比力解算出低成本无人飞行器的角速度信息;步骤6、在子惯导系统中,利用低成本无人飞行器的比力信息与角速度信息进行导航递推解算,测量出低成本无人飞行器的速度值;步骤7、将主惯导系统测量的速度与子惯导系统测量的速度作差,并将差值作为观测向量,采用卡尔曼滤波技术估测出主惯导系统与子惯导系统之间的水平方向安装偏差φ<sub>x</sub>、φ<sub>y</sub>;估测时间为120秒,其中前80秒为滤波估测稳定时间,由80秒到120秒应用平均滤波技术对于数据进行平滑处理;<maths num="0001"><![CDATA[<math><mrow><msub><mi>&phi;</mi><mi>x</mi></msub><mo>=</mo><mfrac><mrow><msub><mover><mi>&phi;</mi><mo>~</mo></mover><mi>x</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mover><mi>&phi;</mi><mo>~</mo></mover><mi>x</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mover><mi>&phi;</mi><mo>~</mo></mover><mi>x</mi></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow><mi>N</mi></mfrac></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>&phi;</mi><mi>y</mi></msub><mo>=</mo><mfrac><mrow><msub><mover><mi>&phi;</mi><mo>~</mo></mover><mi>y</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mover><mi>&phi;</mi><mo>~</mo></mover><mi>y</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mover><mi>&phi;</mi><mo>~</mo></mover><mi>y</mi></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow><mi>N</mi></mfrac></mrow></math>]]></maths>其中<img file="F2009100730747C0000013.GIF" wi="298" he="67" />表示由80秒到120秒数据平滑过程中,第k个滤波估测时间点的滤波估测值;在80秒到120秒数据平滑过程中,共有N个滤波估测时间点,即k=1~N;步骤8、利用步骤7估测出主惯导系统与子惯导系统之间的水平方向安装偏差φ<sub>x</sub>、φ<sub>y</sub>,以及子惯导系统测量比力值<img file="F2009100730747C0000021.GIF" wi="498" he="81" />主惯导上系统测量比力值<img file="F2009100730747C0000022.GIF" wi="474" he="82" />测量出主惯导系统与子惯导系统之间的方位安装偏差φ<sub>z</sub>;<maths num="0003"><![CDATA[<math><mrow><msub><mi>&phi;</mi><mi>z</mi></msub><mo>=</mo><mi>arcsin</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>k</mi><mn>1</mn></msub><msub><mi>k</mi><mn>4</mn></msub><mo>+</mo><msub><mi>k</mi><mn>3</mn></msub><msub><mi>k</mi><mn>2</mn></msub></mrow><mrow><msup><msub><mi>k</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>k</mi><mn>2</mn></msub><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>其中参数k<sub>1</sub>,k<sub>2</sub>,k<sub>3</sub>,k<sub>4</sub>为<maths num="0004"><![CDATA[<math><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>=</mo><msubsup><mi>f</mi><mi>ibx</mi><mi>m</mi></msubsup></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><msubsup><mi>f</mi><mi>iby</mi><mi>m</mi></msubsup></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>k</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><msubsup><mi>f</mi><mi>ibx</mi><mi>s</mi></msubsup><mo>+</mo><mi>sin</mi><msub><mi>&phi;</mi><mi>y</mi></msub><msub><mrow><mi>cos</mi><mi>&phi;</mi></mrow><mi>y</mi></msub><msubsup><mi>f</mi><mi>ibz</mi><mi>m</mi></msubsup></mrow><msub><mrow><mi>cos</mi><mi>&phi;</mi></mrow><mi>y</mi></msub></mfrac><mo>-</mo><mfrac><mrow><msub><mrow><mi>sin</mi><mi>&phi;</mi></mrow><mi>y</mi></msub><msub><mrow><mi>sin</mi><mi>&phi;</mi></mrow><mi>x</mi></msub><msubsup><mi>f</mi><mi>iby</mi><mi>s</mi></msubsup></mrow><mrow><msub><mrow><mi>cos</mi><mi>&phi;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>cos</mi><mi>&phi;</mi></mrow><mi>y</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msub><mrow><mi>sin</mi><mi>&phi;</mi></mrow><mi>y</mi></msub><msub><mrow><mi>sin</mi><mi>&phi;</mi></mrow><mi>x</mi></msub><msub><mrow><mi>sin</mi><mi>&phi;</mi></mrow><mi>x</mi></msub><msubsup><mi>f</mi><mi>ibz</mi><mi>m</mi></msubsup></mrow><mrow><msub><mrow><mi>cos</mi><mi>&phi;</mi></mrow><mi>x</mi></msub><mi>cos</mi><msub><mi>&phi;</mi><mi>y</mi></msub></mrow></mfrac><mo>;</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mi>k</mi><mn>4</mn></msub><mo>=</mo><mfrac><mrow><msub><mrow><mi>sin</mi><mi>&phi;</mi></mrow><mi>x</mi></msub><msubsup><mi>f</mi><mi>ibz</mi><mi>m</mi></msubsup><mo>-</mo><msubsup><mi>f</mi><mi>iby</mi><mi>s</mi></msubsup></mrow><msub><mrow><mi>cos</mi><mi>&phi;</mi></mrow><mi>x</mi></msub></mfrac></mrow></math>]]></maths>步骤9、利用步骤7估测出主惯导系统与子惯导系统之间的水平方向安装偏差φ<sub>x</sub>、φ<sub>y</sub>,以及步骤8测量出主惯导系统与子惯导系统之间的方位安装偏差φ<sub>z</sub>构造出安装偏差方向余弦矩阵C<sub>s</sub><sup>m</sup>;<maths num="0008"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>s</mi><mi>m</mi></msubsup><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mi>&phi;</mi><mi>z</mi></msub></mtd><mtd><mo>-</mo><msub><mi>&phi;</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mrow><mo>-</mo><mi>&phi;</mi></mrow><mi>z</mi></msub></mtd><mtd><mn>1</mn></mtd><mtd><msub><mi>&phi;</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&phi;</mi><mi>y</mi></msub></mtd><mtd><mo>-</mo><msub><mi>&phi;</mi><mi>x</mi></msub></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>步骤10、利用步骤8构造的主惯导系统与子惯导系统之间的安装偏差矩阵C<sub>s</sub><sup>m</sup>,以及主惯导输出的方向余弦矩阵C<sub>m</sub><sup>n</sup>构造出低成本无人飞行器的载体坐标系与导航坐标系之间的方向余弦矩阵C<sub>s</sub><sup>n</sup>;<maths num="0009"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>s</mi><mi>n</mi></msubsup><mo>=</mo><msubsup><mi>C</mi><mi>m</mi><mi>n</mi></msubsup><msubsup><mi>C</mi><mi>s</mi><mi>m</mi></msubsup><mo>;</mo></mrow></math>]]></maths>步骤11、由低成本无人飞行器的载体坐标系与导航坐标系之间的方向余弦矩阵C<sub>s</sub><sup>n</sup>测量出无人飞行器的初始姿态角。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号1号楼哈尔滨工程大学科技处知识产权办公室