留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于X射线脉冲星的环火探测器同步定位和授时仿真分析

魏二虎 王凌轩 张帅 刘经南

魏二虎, 王凌轩, 张帅, 刘经南. 基于X射线脉冲星的环火探测器同步定位和授时仿真分析[J]. 武汉大学学报 ● 信息科学版, 2017, 42(8): 1035-1039. doi: 10.13203/j.whugis20150437
引用本文: 魏二虎, 王凌轩, 张帅, 刘经南. 基于X射线脉冲星的环火探测器同步定位和授时仿真分析[J]. 武汉大学学报 ● 信息科学版, 2017, 42(8): 1035-1039. doi: 10.13203/j.whugis20150437
WEI Erhu, WANG Lingxuan, ZHANG Shuai, LIU Jingnan. Simulation Analysis of Mars Probe Autonomous Positioning and Timing Based on X-Ray Pulsars[J]. Geomatics and Information Science of Wuhan University, 2017, 42(8): 1035-1039. doi: 10.13203/j.whugis20150437
Citation: WEI Erhu, WANG Lingxuan, ZHANG Shuai, LIU Jingnan. Simulation Analysis of Mars Probe Autonomous Positioning and Timing Based on X-Ray Pulsars[J]. Geomatics and Information Science of Wuhan University, 2017, 42(8): 1035-1039. doi: 10.13203/j.whugis20150437

基于X射线脉冲星的环火探测器同步定位和授时仿真分析

doi: 10.13203/j.whugis20150437
基金项目: 

国家自然科学基金 41374012

详细信息
    作者简介:

    魏二虎, 博士, 教授, 现主要从事空间大地测量和地球动力学研究。ehwei@sgg.whu.edu.cn

  • 中图分类号: P228

Simulation Analysis of Mars Probe Autonomous Positioning and Timing Based on X-Ray Pulsars

Funds: 

The National Natural Science Foundation of China 41374012

More Information
    Author Bio:

    WEI Erhu, PhD, professor, specializes in space geodesy and geodynamics. E-mail: ehwei@sgg.whu.edu.cn

图(4) / 表(1)
计量
  • 文章访问数:  910
  • HTML全文浏览量:  43
  • PDF下载量:  335
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-03-29
  • 刊出日期:  2017-08-05

基于X射线脉冲星的环火探测器同步定位和授时仿真分析

doi: 10.13203/j.whugis20150437
    基金项目:

    国家自然科学基金 41374012

    作者简介:

    魏二虎, 博士, 教授, 现主要从事空间大地测量和地球动力学研究。ehwei@sgg.whu.edu.cn

  • 中图分类号: P228

摘要: 利用X射线脉冲星观测量不仅能够为深空探测器提供稳定可靠的位置、速度和姿态信息,还能够对探测器上的原子钟进行长期误差校正。火星探测任务中,存在探测器飞行时间长,钟差随时间累积严重降低脉冲星自主定位系统的精度等问题。针对此,提出了一种环火探测器同步定位和授时方法,对钟差进行建模,将钟差模型参数纳入到滤波系统中,利用自适应扩展卡尔曼滤波器同时计算探测器的位置、速度和钟差模型参数。仿真实验表明,经滤波校正后的钟差保持在300 ns以内,探测器的位置和速度精度分别优于200 m和0.03 m/s,精度得到明显改善。

English Abstract

魏二虎, 王凌轩, 张帅, 刘经南. 基于X射线脉冲星的环火探测器同步定位和授时仿真分析[J]. 武汉大学学报 ● 信息科学版, 2017, 42(8): 1035-1039. doi: 10.13203/j.whugis20150437
引用本文: 魏二虎, 王凌轩, 张帅, 刘经南. 基于X射线脉冲星的环火探测器同步定位和授时仿真分析[J]. 武汉大学学报 ● 信息科学版, 2017, 42(8): 1035-1039. doi: 10.13203/j.whugis20150437
WEI Erhu, WANG Lingxuan, ZHANG Shuai, LIU Jingnan. Simulation Analysis of Mars Probe Autonomous Positioning and Timing Based on X-Ray Pulsars[J]. Geomatics and Information Science of Wuhan University, 2017, 42(8): 1035-1039. doi: 10.13203/j.whugis20150437
Citation: WEI Erhu, WANG Lingxuan, ZHANG Shuai, LIU Jingnan. Simulation Analysis of Mars Probe Autonomous Positioning and Timing Based on X-Ray Pulsars[J]. Geomatics and Information Science of Wuhan University, 2017, 42(8): 1035-1039. doi: 10.13203/j.whugis20150437
  • 火星探测任务中航天器的精确导航定位至关重要,而火星和地球之间的遥远距离和探测器所处的复杂太空环境对基于地面测控通讯站的传统导航方法而言是一个巨大挑战,航天器的自主导航将是突破该瓶颈的关键技术之一[1]

    X射线脉冲星自主导航技术具有可靠性高、自主性强、精确性良好和使用范围广等优点,在火星探测器自主导航方面极具应用潜力。国内外一些学者在利用X射线脉冲星观测量为执行火星探测任务的飞行器提供自主导航服务方面开展了相关研究,包括相关时空基准的定义和转换、适用于火星探测器的脉冲星测量模型的建立、火星探测器轨道动力学研究和导航滤波算法等内容[2, 3]。本文在此基础上,针对初始状态误差较大会导致普通的滤波算法收敛较慢甚至发散,以及火星探测任务周期长,探测器上搭载的原子钟的钟差随时间累积,严重损害了自主定位系统的精度等问题,提出了利用X射线脉冲星观测量同步滤波估计探测器的位置、速度和钟差模型参数的同步定位和授时方法,并利用设计的自适应扩展卡尔曼滤波器进行了仿真实验。

    • 由X射线脉冲星自主导航原理可知,该方法是一种测定探测器相对于某一特定参考系的相对测量方法。针对火星轨道探测器,该方法选择火星质心作为基准点,建立J2000.0火星质心天球坐标系。此时,X射线脉冲星观测量即脉冲信号到达探测器和火星质心的时间延迟τ可以表示为探测器在J2000.0火星质心天球坐标系中的位置参数r和钟差ΔtS的函数。应用于工程计算的简化脉冲星轨道测量模型[1]如式(1) 所示:

      $$ \tau = \frac{1}{c}\mathit{\boldsymbol{n}} \cdot \mathit{\boldsymbol{r + }}\Delta {t_S} + v $$ (1)

      式中,n为脉冲星的方向矢量;v为测量噪声。若同时获得多于3颗脉冲星的测量信息,则可计算出探测器的位置参数r和钟差ΔtS

      在航天器导航中,通常结合探测器的轨道动力学信息和观测信息,利用相应的卡尔曼滤波算法进行状态参数的预测。

    • 在导航滤波算法中,探测器的轨道动力学模型主要是在获取新的观测量之前,对探测器轨道进行预报。火星轨道探测器受到的主要作用力是火星的中心引力,对于中高轨道探测器而言,研究表明火星非球型摄动中的J2项为主要摄动项[2, 4]

      在J2000.0火星质心天球坐标系中,火星轨道探测器的动力学模型可表示为:

      $$ \mathit{\boldsymbol{\ddot r = }} - \frac{{{\mu _M}\mathit{\boldsymbol{r}}}}{{{r^3}}} + {\mathit{\boldsymbol{a}}_{{\mathit{\boldsymbol{J}}_2}}} + \mathit{\boldsymbol{\omega }} $$ (2)

      式中,r=[x y z]T为探测器位置矢量;$\mathit{\boldsymbol{\ddot r}}$为其二阶导数,r=||r||; μM是火星引力常数;w为满足高斯分布的随机噪声;aJ2代表探测器受到的J2项摄动力加速度。

      $$ {\mathit{\boldsymbol{a}}_{{\mathit{\boldsymbol{J}}_2}}} = \left[ {\begin{array}{*{20}{c}} { - \frac{{{\mu _M}x}}{{{\mathit{\boldsymbol{r}}^3}}}{\mathit{\boldsymbol{J}}_2}{{\left( {\frac{{{R_M}}}{\mathit{\boldsymbol{r}}}} \right)}^2}\frac{3}{2}\left( {1 - 5\frac{{{z^2}}}{{{\mathit{\boldsymbol{r}}^2}}}} \right)}\\ { - \frac{{{\mu _M}y}}{{{\mathit{\boldsymbol{r}}^3}}}{\mathit{\boldsymbol{J}}_2}{{\left( {\frac{{{R_M}}}{\mathit{\boldsymbol{r}}}} \right)}^2}\frac{3}{2}\left( {1 - 5\frac{{{z^2}}}{{{\mathit{\boldsymbol{r}}^2}}}} \right)}\\ { - \frac{{{\mu _M}z}}{{{\mathit{\boldsymbol{r}}^3}}}{\mathit{\boldsymbol{J}}_2}{{\left( {\frac{{{R_M}}}{\mathit{\boldsymbol{r}}}} \right)}^2}\frac{3}{2}\left( {1 - 5\frac{{{z^2}}}{{{\mathit{\boldsymbol{r}}^2}}}} \right)} \end{array}} \right] $$ (3)

      式中,RM为火星平均半径。

      在给定探测器初始轨道参数的情况下,即可利用轨道动力学模型进行数值积分计算下一时刻的探测器状态。

    • 一般情况下原子钟的钟差可以用形如式(4) 的二阶多项式表示[5]

      $$ \begin{array}{l} {x_1}\left( {{t_{k + 1}}} \right) = {x_1}\left( {{t_k}} \right) + \tau \cdot {x_2}\left( {{t_k}} \right) + \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\frac{{{\tau ^2}}}{2}{x_3}\left( {{t_k}} \right) + \omega \left( k \right) \end{array} $$ (4)

      式中,x1(tk)、x2(tk)和x3(tk)分别表示tk时刻原子钟的钟差、频率漂移和频率漂移变化率,它们组成了钟差参数;ω(k)为相应的原子钟噪声。

      以时钟钟差、时钟频率漂移和频率漂移变化率为状态量的钟差状态方程为:

      $$ {\left[ {\begin{array}{*{20}{c}} {{x_1}}\\ {{x_2}}\\ {{x_3}} \end{array}} \right]_{k + 1}} = {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_{k + 1,k}}{\left[ {\begin{array}{*{20}{c}} {{x_1}}\\ {{x_2}}\\ {{x_3}} \end{array}} \right]_k} = {\mathit{\boldsymbol{\omega }}_k} $$ (5)

      式中,ωk为钟差噪声向量;状态转移矩阵Φk+1, k为:

      $$ {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_{k + 1,k}} = \left[ {\begin{array}{*{20}{c}} 1&\varepsilon &{{\varepsilon ^2}/2}\\ 0&1&\varepsilon \\ 0&0&1 \end{array}} \right] $$ (6)

      式中,ε为时间间隔;ωk满足以下统计学模型:

      $$ \begin{array}{*{20}{c}} {{Q_k}\left( \varepsilon \right) = E\left[ {{\mathit{\boldsymbol{\omega }}_k}\mathit{\boldsymbol{\omega }}_k^{\rm{T}}} \right] = }\\ {\left[ {\begin{array}{*{20}{c}} {{q_1}\varepsilon + \frac{1}{3}{q_2}{\varepsilon ^3} + \frac{1}{{20}}{q_3}{\varepsilon ^5}}&{\frac{1}{2}{q_2}{\varepsilon ^2} + \frac{1}{8}{q_3}{\varepsilon ^4}}&{\frac{1}{6}{q_3}{\varepsilon ^3}}\\ {\frac{1}{2}{q_2}{\varepsilon ^2} + \frac{1}{8}{q_3}{\varepsilon ^4}}&{{q_2}\varepsilon + \frac{1}{3}{q_3}{\varepsilon ^3}}&{\frac{1}{2}{q_3}{\varepsilon ^2}}\\ {\frac{1}{6}{q_3}{\varepsilon ^3}}&{\frac{1}{2}{q_3}{\varepsilon ^2}}&{{q_3}\varepsilon } \end{array}} \right]} \end{array} $$ (7)

      式(7) 中,q1q2q3为连续过程噪声功率谱密度。

    • 在J2000.0火星质心天球坐标系中,设探测器位置、速度以及3个钟差参数组成火星轨道探测器同步定位和授时系统状态向量为:

      $$ \mathit{\boldsymbol{x}}\left( t \right) = {\left[ {x\;y\;z\;{v_x}\;{v_y}\;{v_z}\;{x_1}\;{x_2}\;{x_3}} \right]^{\rm{T}}} $$ (8)

      结合火星轨道探测器动力学模型(2) 和钟差模型(4) 可得同步定位和授时系统状态方程为:

      $$ \mathit{\boldsymbol{\dot x}}\left( t \right) = \mathit{\boldsymbol{F}}\left( t \right) \cdot \mathit{\boldsymbol{x}}\left( t \right) + \mathit{\boldsymbol{W}}\left( t \right) $$ (9)

      式中,W(t)为系统状态噪声矩阵,状态方程系数矩阵F(t)为:

      $$ \begin{array}{*{20}{c}} {\mathit{\boldsymbol{F}}\left( t \right) = }\\ {\left[ {\begin{array}{*{20}{c}} 0&0&0&1&0&0&0&0&0\\ 0&0&0&0&1&0&0&0&0\\ 0&0&0&0&0&1&0&0&0\\ {\frac{{\partial {a_x}}}{{\partial x}}}&{\frac{{\partial {a_x}}}{{\partial y}}}&{\frac{{\partial {a_x}}}{{\partial z}}}&0&0&0&0&0&0\\ {\frac{{\partial {a_x}}}{{\partial x}}}&{\frac{{\partial {a_y}}}{{\partial y}}}&{\frac{{\partial {a_z}}}{{\partial z}}}&0&0&0&0&0&0\\ {\frac{{\partial {a_z}}}{{\partial x}}}&{\frac{{\partial {a_z}}}{{\partial y}}}&{\frac{{\partial {a_z}}}{{\partial z}}}&0&0&0&0&0&0\\ 0&0&0&0&0&0&0&1&0\\ 0&0&0&0&0&0&0&0&1\\ 0&0&0&0&0&0&0&0&0 \end{array}} \right]} \end{array} $$ (10)
    • 若同一时段观测了p颗脉冲星,参考脉冲星测量模型(1),则环火探测器自主定位和授时系统的测量方程为:

      $$ \mathit{\boldsymbol{z = }}{\left[ {\Delta {t_1} \cdots \cdots \Delta {t_p}} \right]^{\rm{T}}} = \mathit{\boldsymbol{H}} \cdot \mathit{\boldsymbol{x + \eta }} $$ (11)

      式中,z表示脉冲星观测量;方程矩阵H为:

      $$ \mathit{\boldsymbol{H = }}\left[ \begin{array}{l} \begin{array}{*{20}{c}} {\frac{{{\mathit{\boldsymbol{n}}_1}}}{c}}&0&0&0&1&0&0 \end{array}\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; \vdots \\ \begin{array}{*{20}{c}} {\frac{{{\mathit{\boldsymbol{n}}_p}}}{c}}&0&0&0&1&0&0 \end{array} \end{array} \right] $$ (12)

      完成上述同步定位和授时系统状态方程以及观测方程的建立后,便可以利用相应的自适应扩展卡尔曼滤波算法进行状态参数的解算[6, 7]

    • 在选择导航脉冲星时,除了应选择经过可见性条件检查的脉冲星,还要考虑脉冲星的测量误差及脉冲星组合的位置结构对定位精度的影响[8, 9]。实验中,利用脉冲星优选因子(POSI)来定量计算脉冲星测量误差和位置分布对定位精度的影响值:

      $$ \begin{array}{*{20}{c}} {{\rm{Cov}}\left( {\mathit{\boldsymbol{\hat x}}} \right) = E\left( {\mathit{\boldsymbol{\hat x}}{{\mathit{\boldsymbol{\hat x}}}^{\rm{T}}}} \right) = }\\ {{{\left( {{\mathit{\boldsymbol{H}}^{\rm{T}}}\mathit{\boldsymbol{H}}} \right)}^{ - 1}}{\mathit{\boldsymbol{H}}^{\rm{T}}}E\left( {\Delta \mathit{\boldsymbol{\rho }}\Delta {\mathit{\boldsymbol{\rho }}^{\rm{T}}}} \right)\mathit{\boldsymbol{H}}{{\left( {{\mathit{\boldsymbol{H}}^{\rm{T}}}\mathit{\boldsymbol{H}}} \right)}^{ - 1}}}\\ {{{\left( {{\mathit{\boldsymbol{H}}^{\rm{T}}}\mathit{\boldsymbol{H}}} \right)}^{ - 1}}{\mathit{\boldsymbol{H}}^{\rm{T}}}{\rm{diag}}\left( {\sigma _1^2,\sigma _2^2, \cdots \sigma _p^2} \right)\mathit{\boldsymbol{H}} \cdot }\\ {{{\left( {{\mathit{\boldsymbol{H}}^{\rm{T}}}\mathit{\boldsymbol{H}}} \right)}^{ - 1}} = {D_p}} \end{array} $$ (13)

      式中,$\mathit{\boldsymbol{\hat x}}$ 为状态向量x的修正量;Δρ为相应的脉冲星测量误差,其方差矩阵为diag(σ12, σ22, …σp2);H为简化的脉冲星观测方程矩阵:

      $$ \mathit{\boldsymbol{H = }}{\left[ {{\mathit{\boldsymbol{n}}_1}\;{\mathit{\boldsymbol{n}}_2}\; \cdots \;{\mathit{\boldsymbol{n}}_p}} \right]^{\rm{T}}} $$ (14)

      式中,n表示相应脉冲星的方向矢量。

      此次环火探测器同步定位和授时实验中,我们从澳大利亚脉冲星数据库中挑选了14颗信噪比高、轮廓清晰的脉冲星,利用式(13) 进行了定量优选,最终确定B1937+21、B1821-24、B0531+21、B0540-69、B1957+20、B0614+091作为导航脉冲星,相应参数列于表 1中。

      表 1  脉冲星测距精度表

      Table 1.  Pulsar Ranging Accuracy

      脉冲星编号 测距精度σ/m
      B1937+21 689.5
      B1821-24 495.8
      B0531+21 88.3
      B0540-69 1 493.0
      B1957+20 932.0
      B0614+091 282.0
    • 利用Satellite Tool Kit(STK)软件[10]仿真得到火星轨道探测器的“标称轨道”,采样间隔为3 600 s,总采样数为2 400个,初始轨道参数为:a=15 000 km, e=0.005, i=30°,Ω=30°,ω=30°,M0=0°。

      在观测量模拟过程中,原子钟的钟差由钟差模型式(4) 模拟得到,初始钟差设为2×10-6 s,时钟频率漂移设为4.0×10-11 s/s, 时钟漂移变化率设为6.0×10-18 s/s2图 1为模拟的钟差累积示意图。从图 1中可知,若不进行定期校正,钟差将严重损害定位精度。

      图  1  模拟钟差累积示意图

      Figure 1.  Schematic of Simulate Clock Correction Accumulation

      其次,在观测量模拟过程中加入了-2σ~2σ的随机噪声,σ为相应脉冲星测距精度中误差。

      导航滤波的初始误差设置为(10 km, 10 km, 10 km)、(2 m/s, 2 m/s, 2 m/s)和(2×10-7s, 0, 0),时钟频率漂移及其变化率同上。

    • 完成上述模型建立及观测量模拟后便可以利用相应的自适应扩展卡尔曼滤波(AEKF)算法[11-13]进行导航滤波计算。

      图 2为同步定位和授时算法滤波修正后的钟差示意图。从图 2中可以看出,经过15个历元左右滤波误差迅速收敛至400 ns以内,最终稳定在300 ns。

      图  2  滤波修正后的钟差示意图

      Figure 2.  Clock Correction After Filtering

      图 3为滤波修正后的探测器位置误差示意图,图 4为相应的速度误差示意图。从图 3~4中可知,探测器在3个坐标轴方向上的位置、速度误差经几个历元后便可以迅速收敛到1 000 m、0.1 m/s以内,经过十几个历元以后稳定在200 m、0.03 m/s以内。

      图  3  滤波修正后的位置误差示意图

      Figure 3.  Position Error After Filtering

      图  4  滤波修正后的速度误差示意图

      Figure 4.  Velocity Error After Filtering

      图 3图 4可知,3个坐标轴方向上的滤波位置、速度精度水平相当,没有出现自主定位系统中X轴方向滤波精度差于Y、Z轴方向的精度情况[5]

      统计结果表明,基于X射线脉冲星的环火探测器同步定位和授时滤波算法能够有效地克服初始状态误差的影响,获得较高精度的探测器位置、速度和钟差参数滤波值,验证了该算法在自主定位和授时方面的高效性。

    • 本文研究了基于X射线脉冲星的火星轨道探测器同步定位和授时原理,建立了相应的同步定位和授时模型,利用POSI进行导航脉冲星的优选,并开展了基于自适应扩展卡尔曼滤波的仿真实验。实验表明, 同步定位授时算法能够有效地校正原子钟钟差,削弱钟差对脉冲星观测量的影响,改善探测器位置和速度参数的精度。

参考文献 (13)

目录

    /

    返回文章
    返回