发明名称 一种基于共基座安装的传递对准精度定量评估方法
摘要 本发明属于惯导系统性能评估领域,具体涉及一种针对捷联惯导系统传递对准精度的基于共基座安装的传递对准精度定量评估方法。本发明包括:将子、主惯导共基座安装于同一铝板,构建子、主惯导间方位安装失准角;将安装有子、主惯导的铝板置于高精度转台,利用光学瞄准镜使主惯导的y轴指向正北,利用激光跟踪仪测定子、主惯导y轴指向间的方位安装误差角,作为方位安装误差角基准值;主惯导开机并完成自对准。该方法获取对准精度评估的基准在传递对准之前,通过基准值和估计值的对比即可实现对准精度的定量评估。避免了对准完成后子惯导进行导航解算和平滑解算,大大降低了计算复杂度和数据存储需求;可以实现短时间内多次的重复试验。
申请公布号 CN105973268A 申请公布日期 2016.09.28
申请号 CN201610293287.0 申请日期 2016.05.06
申请人 哈尔滨工程大学 发明人 程建华;陈岱岱;费再慧;于天琦;董楠楠;蒋国桉;董铭涛
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 代理人
主权项 一种基于共基座安装的传递对准精度定量评估方法,其特征在于,包括如下步骤:(1)将子、主惯导共基座安装于同一铝板,构建子、主惯导间方位安装失准角;(2)将安装有子、主惯导的铝板置于高精度转台,利用光学瞄准镜使主惯导的y轴指向正北,利用激光跟踪仪测定子、主惯导y轴指向间的方位安装误差角,作为方位安装误差角基准值;(3)主惯导开机并完成自对准,进入导航状态;(4)子惯导开机;(5)主惯导向子惯导传递位置、速度、姿态矩阵;(6)子惯导进行导航解算,同步采集子、主惯导的速度、姿态信息;(7)利用卡尔曼滤波进行“速度+姿态”匹配快速传递对准解算,估计出子惯导系统误差状态;所涉及的子惯导系统误差状态向量为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>&phi;</mi><mi>mx</mi></msub></mtd><mtd><msub><mi>&phi;</mi><mi>my</mi></msub></mtd><mtd><msub><mi>&phi;</mi><mi>mz</mi></msub></mtd><mtd><msub><mi>&delta;V</mi><mi>x</mi></msub></mtd><mtd><msub><mi>&delta;V</mi><mi>y</mi></msub></mtd><mtd><msub><mi>&phi;</mi><mi>ax</mi></msub></mtd><mtd><msub><mi>&phi;</mi><mi>ay</mi></msub></mtd><mtd><msub><mi>&phi;</mi><mi>az</mi></msub></mtd><mtd><msub><mo>&dtri;</mo><mi>x</mi></msub></mtd><mtd><msub><mo>&dtri;</mo><mi>y</mi></msub></mtd><mtd><msub><mo>&dtri;</mo><mi>z</mi></msub></mtd><mtd><msub><mi>&epsiv;</mi><mi>x</mi></msub></mtd><mtd><msub><mi>&epsiv;</mi><mi>y</mi></msub></mtd><mtd><msub><mi>&epsiv;</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow>]]></math><img file="FDA0000982816220000018.GIF" wi="1580" he="90" /></maths>式中,φ<sub>mx</sub>,φ<sub>my</sub>,φ<sub>mz</sub>为主惯导系到计算得到的子惯导坐标系之间的相对姿态误差;δV<sub>x</sub>,δV<sub>y</sub>为子惯导速度误差;φ<sub>ax</sub>,φ<sub>ay</sub>,φ<sub>az</sub>为主惯导载体坐标系到实际的子惯导载体坐标系的相对姿态误差;<img file="FDA0000982816220000019.GIF" wi="211" he="63" />为子惯导三个轴向的加速度计常值零偏;ε<sub>x</sub>,ε<sub>y</sub>,ε<sub>z</sub>为子惯导三个轴向的陀螺常值漂移;所涉及的快速传递对准系统误差方程为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mover><mi>&phi;</mi><mo>&CenterDot;</mo></mover><mi>m</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mi>&phi;</mi><mi>m</mi></msub><mo>-</mo><msub><mi>&phi;</mi><mi>a</mi></msub><mo>)</mo></mrow><mo>&times;</mo><msubsup><mover><mi>&omega;</mi><mo>&OverBar;</mo></mover><mrow><mi>n</mi><mi>m</mi></mrow><mi>s</mi></msubsup><mo>+</mo><msup><mi>&epsiv;</mi><mi>s</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>&delta;</mi><mover><mi>V</mi><mo>&CenterDot;</mo></mover><mo>=</mo><msubsup><mi>C</mi><mover><mi>s</mi><mo>^</mo></mover><mi>n</mi></msubsup><mrow><mo>(</mo><mo>&lsqb;</mo><msub><mi>&phi;</mi><mi>m</mi></msub><mo>&times;</mo><mo>&rsqb;</mo><mo>-</mo><mo>&lsqb;</mo><msub><mi>&phi;</mi><mi>a</mi></msub><mo>&times;</mo><mo>&rsqb;</mo><mo>)</mo></mrow><msubsup><mover><mi>f</mi><mo>&OverBar;</mo></mover><mi>m</mi><mi>s</mi></msubsup><mo>-</mo><mrow><mo>(</mo><mn>2</mn><msubsup><mover><mi>&omega;</mi><mo>&OverBar;</mo></mover><mrow><mi>i</mi><mi>e</mi></mrow><mi>n</mi></msubsup><mo>+</mo><msubsup><mover><mi>&omega;</mi><mo>&OverBar;</mo></mover><mrow><mi>e</mi><mi>n</mi></mrow><mi>n</mi></msubsup><mo>)</mo></mrow><mo>&times;</mo><mi>&delta;</mi><mi>V</mi><mo>+</mo><msubsup><mi>C</mi><mover><mi>s</mi><mo>^</mo></mover><mi>n</mi></msubsup><msubsup><mo>&dtri;</mo><mi>s</mi><mi>s</mi></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><msub><mover><mi>&phi;</mi><mo>&CenterDot;</mo></mover><mi>a</mi></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mover><mo>&dtri;</mo><mo>&CenterDot;</mo></mover><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><mover><mi>&epsiv;</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000982816220000011.GIF" wi="1061" he="413" /></maths>式中,φ<sub>m</sub>=[φ<sub>mx</sub>,φ<sub>my</sub>,φ<sub>mz</sub>]<sup>T</sup>,φ<sub>a</sub>=[φ<sub>ax</sub>,φ<sub>ay</sub>,φ<sub>az</sub>]<sup>T</sup>,δV=[δV<sub>x</sub>,δV<sub>y</sub>]<sup>T</sup>,ε<sup>s</sup>=[ε<sub>x</sub>,ε<sub>y</sub>,ε<sub>z</sub>]<sup>T</sup>,<img file="FDA0000982816220000012.GIF" wi="70" he="55" />为主惯导坐标系m相对东北天地理坐标系n的旋转角速度在子惯导坐标系s下的投影;<img file="FDA0000982816220000013.GIF" wi="54" he="62" />为子惯导计算坐标系<img file="FDA0000982816220000014.GIF" wi="27" he="47" />转换至东北天地理坐标系n的转换矩阵;<img file="FDA0000982816220000015.GIF" wi="54" he="69" />为主惯导比力在子惯导系的投影;<img file="FDA0000982816220000016.GIF" wi="58" he="55" />为地球自转角速度;<img file="FDA0000982816220000017.GIF" wi="62" he="63" />为东北天地理坐标系相对于地球系的旋转角速度;(8)滤波结束后,存储主惯导载体坐标系到实际的子惯导载体坐标系的相对姿态误差φ<sub>ax</sub>,φ<sub>ay</sub>,φ<sub>az</sub>,其中φ<sub>az</sub>作为方位安装误差角估计值;(9)通过对比步骤(2)中测得的方位安装误差角基准值与步骤(8)中得到的方位安装误差角估计值,实现传递对准定量精度评估。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室