发明名称 一种激光惯组与里程计组合导航方法及系统
摘要 本发明公开了一种激光惯组与里程计组合导航方法及系统,方法包括步骤:S1、采集正常工作状态下N个采样周期内里程计脉冲数,求解惯组与里程计之间的俯仰安装偏角和航向安装偏角;S2、建立状态方程;S3、建立量测方程,其中量测量<img file="DDA0000758554310000011.GIF" wi="508" he="92" /><img file="DDA0000758554310000012.GIF" wi="129" he="84" />为惯组在第i个采样周期内导航系下的位移增量,<img file="DDA0000758554310000013.GIF" wi="118" he="83" />为第i个采样周期内里程计增量<img file="DDA0000758554310000014.GIF" wi="90" he="83" />在导航系下的位移投影;S4、利用状态方程和量测方程,进行卡尔曼滤波,以实时修正惯性导航系统参数误差和器件参数误差,完成组合导航。本发明还提供实现上述方法的系统。实施本发明能有效修正组合导航误差量,大幅度提高组合导航精度,尤其适用于惯组处于大角度斜置的情况。
申请公布号 CN104977004A 申请公布日期 2015.10.14
申请号 CN201510408912.7 申请日期 2015.07.13
申请人 湖北航天技术研究院总体设计所 发明人 李旦;刘明;牟方佳;杨道安;胡华峰
分类号 G01C21/20(2006.01)I;G01C25/00(2006.01)I 主分类号 G01C21/20(2006.01)I
代理机构 武汉东喻专利代理事务所(普通合伙) 42224 代理人 宋业斌
主权项 一种激光惯组与里程计组合导航方法,其特征在于,所述方法包括以下步骤:S1、采集正常工作状态下第i个采样周期内里程计脉冲数N<sub>i</sub>,求解惯组与里程计之间的俯仰安装偏角α<sub>θ</sub>和航向安装偏角α<sub>ψ</sub>,<img file="FDA0000758554280000011.GIF" wi="987" he="385" />所述<img file="FDA0000758554280000012.GIF" wi="874" he="102" />为矢量<img file="FDA0000758554280000013.GIF" wi="82" he="76" />第j个分量,j=1,2,3;<img file="FDA0000758554280000014.GIF" wi="133" he="84" />为第i个采样周期内惯组坐标系与导航坐标系之间的方向余弦矩阵,<img file="FDA0000758554280000015.GIF" wi="186" he="87" />为预设时间内惯组导航解算的位移增量,K<sub>od</sub>为里程计脉冲当量;S2、建立状态方程<img file="FDA0000758554280000016.GIF" wi="294" he="72" />其中,F为状态转移矩阵,w为系统噪声,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><msup><mfenced open = '[' close = ']'><mtable><mtr><mtd><msup><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><msup><mrow><mo>(</mo><msup><mi>&delta;v</mi><mi>n</mi></msup><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><msup><mrow><mo>(</mo><mi>&delta;</mi><mi>P</mi><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><msup><mrow><mo>(</mo><msubsup><mi>&epsiv;</mi><mi>b</mi><mi>b</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><msup><mrow><mo>(</mo><msubsup><mo>&dtri;</mo><mi>b</mi><mi>b</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><mrow><msup><msub><mi>&xi;</mi><mrow><mi>o</mi><mi>d</mi></mrow></msub><mi>T</mi></msup></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000758554280000017.GIF" wi="1087" he="119" /></maths>φ为惯导姿态误差,δv<sup>n</sup>为惯导速度误差,δP为惯导位置误差,<img file="FDA0000758554280000018.GIF" wi="62" he="84" />为陀螺漂移误差;<img file="FDA0000758554280000019.GIF" wi="77" he="82" />为加速度计偏置误差,里程计误差相关量ξ<sub>od</sub>=[δK<sub>od</sub>,δα<sub>θ</sub>,δα<sub>ψ</sub>],其中δK<sub>od</sub>为里程计刻度系数误差,δα<sub>θ</sub>为俯仰安装偏角误差,δα<sub>ψ</sub>为航向安装偏角误差;S3、建立量测方程z=HX+v,其中v为量测噪声;量测量<img file="FDA00007585542800000110.GIF" wi="506" he="107" />累积预设采样周期内惯导位移增量与里程计位移增量,并将二者之差作为量测量z,所述<img file="FDA00007585542800000111.GIF" wi="123" he="83" />为惯组在第i个采样周期内导航系下的位移增量,<img file="FDA00007585542800000112.GIF" wi="128" he="85" />为第i个采样周期内里程计增量<img file="FDA00007585542800000113.GIF" wi="91" he="81" />在导航系下的位移投影;量测矩阵<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mrow><mo>-</mo><mi>&Sigma;</mi><mrow><mo>(</mo><msubsup><mi>C</mi><mi>b</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><msubsup><mi>&Delta;S</mi><mi>i</mi><mi>b</mi></msubsup><mo>)</mo></mrow><mo></mo><mo>&times;</mo></mrow></mtd><mtd><msub><mi>I</mi><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mtd><mtd><mrow><mo>-</mo><msubsup><mi>&Sigma;C</mi><mi>b</mi><mi>n</mi></msubsup><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><msubsup><mi>M</mi><mi>i</mi><mi>b</mi></msubsup></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA00007585542800000114.GIF" wi="1161" he="96" /></maths>I为单位矩阵,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>M</mi><mi>i</mi><mi>b</mi></msubsup><mo>=</mo><msub><mi>N</mi><mi>i</mi></msub><msub><mi>K</mi><mrow><mi>o</mi><mi>d</mi></mrow></msub><mfenced open = '[' close = ']'><mtable><mtr><mtd><mrow><msub><mi>sin&alpha;</mi><mi>&psi;</mi></msub><msub><mi>cos&alpha;</mi><mi>&theta;</mi></msub></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>sin&alpha;</mi><mi>&psi;</mi></msub><msub><mi>sin&alpha;</mi><mi>&theta;</mi></msub></mrow></mtd><mtd><mrow><msub><mi>cos&alpha;</mi><mi>&psi;</mi></msub><msub><mi>cos&alpha;</mi><mi>&theta;</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>cos&alpha;</mi><mi>&psi;</mi></msub><msub><mi>cos&alpha;</mi><mi>&theta;</mi></msub></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>cos&alpha;</mi><mi>&psi;</mi></msub><msub><mi>sin&alpha;</mi><mi>&theta;</mi></msub></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>sin&alpha;</mi><mi>&psi;</mi></msub><msub><mi>cos&alpha;</mi><mi>&theta;</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>sin&alpha;</mi><mi>&theta;</mi></msub></mrow></mtd><mtd><mrow><msub><mi>cos&alpha;</mi><mi>&theta;</mi></msub></mrow></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000758554280000021.GIF" wi="1276" he="253" /></maths>S4、利用所述状态方程和量测方程,进行卡尔曼滤波,以实时修正惯性导航系统参数误差和器件参数误差,完成组合导航。
地址 430040 湖北省武汉市金山大道9号