发明名称 一种自主导航系统误差校正方法
摘要 一种自主导航系统误差校正方法,基于紫外地球敏感器和星敏感器的自主导航系统中,采用紫外敏感器测量得到卫星本体系中地心矢量,采用星敏感器测量得到惯性系到卫星本体系的姿态转换矩阵。由于敏感器不可避免存在相对安装误差,滤波器无法消除该项系统误差,导致导航精度变差。本发明方法对紫外敏感器和星敏感器的相对安装误差进行了建模,将安装误差扩充为状态变量进行导航滤波。假定卫星一段时间内能通过GPS接收机或地面定轨获得高精度位置测量信息,则通过滤波可以对系统误差进行实时估计和校正。本发明方法操作简单,可以显著提高导航精度。
申请公布号 CN102175260B 申请公布日期 2012.11.14
申请号 CN201010623851.3 申请日期 2010.12.31
申请人 北京控制工程研究所 发明人 魏春岭;张斌;黄翔宇;王大轶
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 安丽
主权项 1.一种自主导航系统误差校正方法,其特征在于步骤如下:(1)基于紫外地球敏感器和星敏感器的自主导航系统中,紫外地球敏感器对地球进行观测,得到卫星本体系下的地心矢量值<img file="FSB00000852408400011.GIF" wi="71" he="54" />星敏感器对恒星观测,通过星图识别,得到惯性系到卫星本体系的姿态转换矩阵<img file="FSB00000852408400012.GIF" wi="82" he="68" />其中惯性系原点在地球中心,X轴正方向沿地球赤道面和黄道面的交线指向春分点,Z轴为赤道平面的法线,正方向指向北极方向,Y轴与X、Z轴成右手系;卫星本体系的原点在卫星质心,对于对地定向三轴稳定卫星Z轴指向地心,X轴与Z轴垂直并指向卫星速度方向,Y轴与X、Z轴成右手系;(2)设地球敏感器和星敏感器的安装误差角为Ф=[φ<sub>x</sub> φ<sub>y</sub> φ<sub>z</sub>],则惯性系到卫星本体系的姿态转换矩阵测量值<img file="FSB00000852408400013.GIF" wi="54" he="68" />为<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>i</mi><msup><mi>b</mi><mo>&prime;</mo></msup></msubsup><mo>=</mo><mo>{</mo><msub><mi>I</mi><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub><mo>-</mo><mo>[</mo><mi>&Phi;</mi><mo>&times;</mo><mo>]</mo><mo>}</mo><msubsup><mi>C</mi><mi>i</mi><mi>b</mi></msubsup><mo>,</mo></mrow></math>]]></maths>式中<img file="FSB00000852408400015.GIF" wi="50" he="60" />为惯性系到卫星本体系的真实的姿态转换矩阵,[Ф×]为向量Ф构成的反对称矩阵,<maths num="0002"><![CDATA[<math><mrow><mo>[</mo><mi>&Phi;</mi><mo>&times;</mo><mo>]</mo><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><msub><mrow><mo>-</mo><mi>&phi;</mi></mrow><mi>z</mi></msub></mtd><mtd><msub><mi>&phi;</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&phi;</mi><mi>z</mi></msub></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>&phi;</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>&phi;</mi><mi>y</mi></msub></mtd><mtd><msub><mi>&phi;</mi><mi>x</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>(3)根据步骤(1)和步骤(2)的结果,得到惯性系下地心方向矢量<img file="FSB00000852408400017.GIF" wi="38" he="46" />的测量方程为<maths num="0003"><![CDATA[<math><mrow><msup><mover><mi>u</mi><mo>^</mo></mover><mi>I</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mi>i</mi><msup><mi>b</mi><mo>&prime;</mo></msup></msubsup><mo>)</mo></mrow><mi>T</mi></msup><msup><mover><mi>u</mi><mo>^</mo></mover><mi>B</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mi>i</mi><mi>b</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><mo>{</mo><msub><mi>I</mi><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub><mo>+</mo><mo>[</mo><mi>&Phi;</mi><mo>&times;</mo><mo>]</mo><mo>}</mo><msup><mover><mi>u</mi><mo>^</mo></mover><mi>B</mi></msup></mrow></math>]]></maths>(4)通过GPS接收机或地面定轨获得高精度位置测量信息,将安装误差角Ф的三个分量扩充为滤波状态变量<img file="FSB00000852408400019.GIF" wi="363" he="57" />其中r=[x y z]为卫星的位置矢量,<maths num="0004"><![CDATA[<math><mrow><mover><mi>r</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>v</mi><mi>x</mi></msub></mtd><mtd><msub><mi>v</mi><mi>y</mi></msub></mtd><mtd><msub><mi>v</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>为卫星的速度矢量,将三个安装误差角与卫星位置矢量、卫星速度矢量一起作为状态变量,进行UKF导航滤波解算,对安装误差进行实时估计和校正。
地址 100080 北京市2729信箱