发明名称 一种利用偏航机动的自主导航敏感器系统误差校正方法
摘要 一种利用偏航机动的自主导航敏感器系统误差校正方法,自主导航敏感器通过对地球进行观测得到卫星本体系地心方向矢量和地心距,自主导航敏感器通过对恒星进行观测得到惯性系到卫星本体系的姿态转换矩阵,二者结合得到惯性系下的地心方向和地心距,从而可以确定卫星轨道。本发明方法将自主导航敏感器系统误差扩充为状态变量,以地心方向和地心距作为测量值进行偏差自校准。借鉴惯导系统领域研究较多的动基座对准思想,利用偏航机动来提高导航系统对敏感器系统误差的可观性,从而可以更好地对敏感器系统误差进行估计和补偿。本发明方法操作简单,可以显著提高自主导航精度。
申请公布号 CN102519472B 申请公布日期 2014.07.02
申请号 CN201110409345.9 申请日期 2011.12.08
申请人 北京控制工程研究所 发明人 魏春岭;张斌;张青春;刘良栋;李果;亢淼
分类号 G01C21/24(2006.01)I 主分类号 G01C21/24(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 安丽
主权项 1.一种利用偏航机动的自主导航敏感器系统误差校正方法,其特征在于实现步骤如下:(1)通过自主导航敏感器对地球和恒星进行观测,得到惯性系下地心方向矢量和地心距自主导航敏感器对地球观测,得到卫星滚动角<img file="FSB0000124263910000011.GIF" wi="56" he="43" />和俯仰角θ<sub>m</sub>,导航敏感器测量存在系统误差分别为<img file="FSB0000124263910000012.GIF" wi="47" he="42" />和θ<sub>c</sub>,则卫星真实滚动角<img file="FSB0000124263910000013.GIF" wi="35" he="41" />和俯仰角θ分别为<img file="FSB0000124263910000014.GIF" wi="221" he="44" />θ=θ<sub>m</sub>-θ<sub>c</sub>计算得到卫星本体系下地心方向矢量<img file="FSB0000124263910000015.GIF" wi="47" he="49" /><img file="FSB0000124263910000016.GIF" wi="383" he="220" />自主导航敏感器对恒星观测,得到惯性系到卫星本体系的姿态转换矩阵<img file="FSB0000124263910000017.GIF" wi="79" he="77" />从而计算得到惯性系下地心方向矢量<img file="FSB0000124263910000018.GIF" wi="31" he="40" /><maths num="0001"><![CDATA[<math><mrow><mover><mi>u</mi><mo>&RightArrow;</mo></mover><mo>=</mo><msup><mrow><mo>(</mo><msubsup><mi>C</mi><mi>i</mi><mi>b</mi></msubsup><mo>)</mo></mrow><mi>T</mi></msup><msup><mover><mi>u</mi><mo>&RightArrow;</mo></mover><mi>b</mi></msup><mo>;</mo></mrow></math>]]></maths>由自主导航敏感器测量得到的地球视半径ρ,计算得到卫星到地心的距离r<maths num="0002"><![CDATA[<math><mrow><mi>r</mi><mo>=</mo><mfrac><mrow><msub><mi>R</mi><mi>e</mi></msub><mo>+</mo><msub><mi>h</mi><mrow><mi>co</mi><mn>2</mn></mrow></msub></mrow><mrow><msup><mi>sin</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>&rho;</mi></mrow></mfrac></mrow></math>]]></maths>其中R<sub>e</sub>为地球半径,h<sub>co2</sub>为二氧化碳层高度;(2)以卫星位置矢量、卫星速度矢量为状态变量,同时辅以导航敏感器测量存在的系统误差为状态变量,以地心方向矢量和地心距为测量值,进行导航滤波解算,最后得到位置估计值、速度估计值和系统误差估计值,具体如下:(a)状态变量<img file="FSB00001242639100000111.GIF" wi="408" he="76" />其中r=[x y z]为卫星的位置矢量,x、y、z分别为X轴、Y轴和Z轴的位置坐标值,<maths num="0003"><![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>为卫星的速度矢量,v<sub>x</sub>、v<sub>y</sub>、v<sub>z</sub>分别为X轴、Y轴和Z轴的速度坐标值,<img file="FSB0000124263910000021.GIF" wi="79" he="53" />θ<sub>c</sub>为导航敏感器系统误差,状态方程如下:<img file="FSB0000124263910000022.GIF" wi="704" he="963" />式中<img file="FSB0000124263910000023.GIF" wi="318" he="77" />μ为地球引力常数,Re为地球半径,J<sub>2</sub>为二阶带谐项系数,w<sub>x</sub>,x<sub>y</sub>,w<sub>z</sub>分别为X轴、Y轴和Z轴的系统噪声;(b)将惯性系下地心方向和地心距作为测量值,测量方程为:<maths num="0004"><![CDATA[<math><mrow><mi>Z</mi><mo>=</mo><mi>h</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>+</mo><mi>v</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mover><mi>u</mi><mo>&RightArrow;</mo></mover></mtd></mtr><mtr><mtd><mi>r</mi></mtd></mtr></mtable></mfenced><mo>+</mo><mi>v</mi></mrow></math>]]></maths>其中v为测量噪声,h(X)表示测量方程是状态变量的非线性函数,计算测量值对状态变量的偏导数矩阵为<img file="FSB0000124263910000025.GIF" wi="656" he="291" />其中<maths num="0005"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mover><mi>u</mi><mo>&RightArrow;</mo></mover></mrow><mrow><mo>&PartialD;</mo><mi>r</mi></mrow></mfrac><mfrac><mn>1</mn><mrow><mo>|</mo><mi>r</mi><mo>|</mo></mrow></mfrac><msub><mi>I</mi><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub><mo>-</mo><mi>r</mi><mfrac><msup><mi>r</mi><mi>T</mi></msup><msup><mrow><mo>|</mo><mi>r</mi><mo>|</mo></mrow><mn>3</mn></msup></mfrac></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mover><mi>u</mi><mo>&RightArrow;</mo></mover></mrow><mrow><mo>&PartialD;</mo><mover><mi>r</mi><mo>&CenterDot;</mo></mover></mrow></mfrac><mo>=</mo><msub><mn>0</mn><mrow><mn>3</mn><mo>&times;</mo><mn>3</mn></mrow></msub></mrow></math>]]></maths><img file="FSB0000124263910000031.GIF" wi="406" he="218" /><img file="FSB0000124263910000032.GIF" wi="381" he="219" /><maths num="0007"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>r</mi></mrow><mrow><mo>&PartialD;</mo><mi>r</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mrow><mo>|</mo><mi>r</mi><mo>|</mo></mrow></mfrac><msup><mi>r</mi><mi>T</mi></msup></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>r</mi></mrow><mrow><mo>&PartialD;</mo><mover><mi>r</mi><mo>&CenterDot;</mo></mover></mrow></mfrac><mo>=</mo><mn>0</mn></mrow></math>]]></maths><img file="FSB0000124263910000035.GIF" wi="160" he="128" /><maths num="0009"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>r</mi></mrow><mrow><mo>&PartialD;</mo><msub><mi>&theta;</mi><mi>c</mi></msub></mrow></mfrac><mo>=</mo><mn>0</mn></mrow></math>]]></maths>I<sub>3×3</sub>为3×3维的单位矩阵,0<sub>3×3</sub>为3×3维的零矩阵(c)导航滤波算法采用扩展卡尔曼滤波算法,即可实现对敏感器系统误差<img file="FSB0000124263910000037.GIF" wi="77" he="54" />θ<sub>c</sub>的实时估计和校正,从而得到位置估计值、速度估计值和系统误差估计值;(3)在导航滤波解算过程,通过地面指令使卫星进行90度偏航姿态机动,提高位置估计值、速度估计值和系统误差估计值的准确性,从而提高自主导航精度,实现利用偏航机动的自主导航敏感器系统误差校正。
地址 100190 北京市海淀区2729信箱