留言板

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

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

行星际深空探测中双程测速的高精度计算方法

金炜桐 李斐 杨轩 鄢建国 叶茂 郝卫峰 曲春凯

金炜桐, 李斐, 杨轩, 鄢建国, 叶茂, 郝卫峰, 曲春凯. 行星际深空探测中双程测速的高精度计算方法[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1483-1489. doi: 10.13203/j.whugis20170020
引用本文: 金炜桐, 李斐, 杨轩, 鄢建国, 叶茂, 郝卫峰, 曲春凯. 行星际深空探测中双程测速的高精度计算方法[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1483-1489. doi: 10.13203/j.whugis20170020
JIN Weitong, LI Fei, YANG Xuan, YAN Jianguo, YE Mao, HAO Weifeng, QU Chunkai. High Precision Computational Method of Two-Way Range-Rate in Long-Distance Deep Space Exploration[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1483-1489. doi: 10.13203/j.whugis20170020
Citation: JIN Weitong, LI Fei, YANG Xuan, YAN Jianguo, YE Mao, HAO Weifeng, QU Chunkai. High Precision Computational Method of Two-Way Range-Rate in Long-Distance Deep Space Exploration[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1483-1489. doi: 10.13203/j.whugis20170020

行星际深空探测中双程测速的高精度计算方法

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

国家自然科学基金 41374024

国家自然科学基金 41604004

国家自然科学基金 41174019

湖北省自然科学基金重点项目创新群体类项目 2015CFA011

宇航动力学国家重点实验室开放基金 2016ADL-DW0103

中国博士后科学基金 2016M602360

详细信息
    作者简介:

    金炜桐, 硕士生, 主要从事行星探测器精密定轨方面的研究。jwt_whdx@whu.edu.cn

    通讯作者: 叶茂, 博士。mye@whu.edu.cn
  • 中图分类号: P228

High Precision Computational Method of Two-Way Range-Rate in Long-Distance Deep Space Exploration

Funds: 

The National Natural Science Foundation of China 41374024

The National Natural Science Foundation of China 41604004

The National Natural Science Foundation of China 41174019

the Hubei Province Natural Science Foundation Innovation Group Project 2015CFA011

the Open Research Fund of State Key Laboratory of Astronautic Dynamics 2016ADL-DW0103

the China Postdoctoral Science Foundation 2016M602360

More Information
    Author Bio:

    JIN Weitong, postgraduate, specializes in planetary spacecraft precision orbit determination. E-mail: jwt_whdx@whu.edu.cn

    Corresponding author: YE Mao, PhD. E-mail: mye@whu.edu.cn
  • 摘要: 对于行星际深空探测(距地球1亿km以上)任务,由于受到计算机字长的限制,传统双程测速模型的计算精度无法满足高精度定轨的需要,其最大误差源于多普勒频移周计数终点和始点上行几何距离之间和下行几何距离之间差分值的计算过程。对此建立行星际双程测速模型,高精度地计算了两个差分值,推导模型的计算公式并给出详细步骤,同时给出计算过程中需要的切比雪夫差分多项式递推公式的形式。将该模型在深空探测器精密定轨与重力场解算软件系统(Wuhan University deep-space orbit determination and gravity recovery system,WUDOGS)中进行了实现,并以欧空局火星快车号(Mars express,MEX)探测任务为背景,利用该软件进行仿真测试,从计算精度和定轨结果两个方面验证该模型的优越性。结果表明,该模型将双程测速的计算值在计算机中表达的精度提高2个数量级,同时避免了定轨过程中引入额外的数值误差,可以为后续高精度的行星际深空探测任务的定轨提供参考。
  • 图  1  星际模型详细计算步骤

    Figure  1.  Detailed Calculation Steps of Interplanetary Model

    图  2  两种模型计算值残差图

    Figure  2.  Residual Plot of Calculation Values Using Two Models

    图  3  两种模型模拟定轨最后一次迭代残差图

    Figure  3.  Post-fit Residuals of Orbit Determination Using Two Models

    表  1  计算值残差统计信息

    Table  1.   Statistical Information of Calculation Residuals

    模型 计算值残差/mm·s-1
    最大值 均值 RMS
    传统双程测速 0.111 22 -0.007 005 0.040 993
    星际双程测速 5.000 0×10-3 -3.429 86×10-3 9.419 85×10-4
    下载: 导出CSV

    表  2  两种模型模拟定轨的初轨改正数

    Table  2.   Initial Orbit Corrections of POD Using Two Models, Respectively

    模型 ΔPx/m ΔPy/m ΔPz/m ΔVx/m·s-1 ΔVy/m·s-1 ΔVz/m·s-1
    传统模型 100.097 80 -100.104 90 -99.798 70 2.67×10-5 -6.32×10-7 -1.29×10-5
    星际模型 99.999 98 -99.999 98 -100.000 03 -4.89×10-9 2.58×10-11 2.29×10-9
    下载: 导出CSV

    表  3  两种模型模拟定轨结果与给定初轨的差值

    Table  3.   Initial Orbit Differences Between POD Simulation Results and Given Results Using Two Models, Respectively

    模型 ΔPx/m ΔPy/m ΔPz/m ΔVx/m·s-1 ΔVy/m·s-1 ΔVz/m·s-1
    传统模型 0.097 809 -0.104 888 0.201 331 2.670×10-5 -6.317×10-7 -1.287×10-5
    星际模型 -1.764×10-5 1.934×10-5 -3.640×10-5 -4.890×10-5 2.581×10-11 2.294×10-9
    下载: 导出CSV

    表  4  两种模型模拟定轨的最后一次迭代的观测值残差统计信息/mm·s-1

    Table  4.   Statistical Results of POD Measurement Residuals During the Last Iteration Using Two Models, Respectively/mm·s-1

    模型 最后一次迭代观测值残差均值 最后一次迭代观测值残差RMS
    传统模型 -0.004 1 0.003 1
    星际模型 4.589 5×10-8 4.768 8×10-7
    下载: 导出CSV
  • [1] 叶茂, 李斐, 鄢建国, 等. GRAIL月球重力场模型定轨性能分析[J].武汉大学学报·信息科学版, 2016, 41(1):93-99 http://ch.whu.edu.cn/CN/abstract/abstract3439.shtml

    Ye Mao, Li Fei, Yan Jianguo, et al. The Analysis of GRAIL Gravity Model's Orbit Determination Ability[J]. Geomatics and Information Science of Wuhan University, 2016, 41(1):93-99 http://ch.whu.edu.cn/CN/abstract/abstract3439.shtml
    [2] Konopliv A S, Park R S, Folkner W M. An Improved JPL Mars Gravity Field and Orientation from Mars Orbiter and Lander Tracking Data[J]. Icarus, 2016, 274:253-260 doi:  10.1016/j.icarus.2016.02.052
    [3] Zannoni M, Tortora P. Numerical Error in Interplanetary Orbit Determination Software[J]. Journal of Guidance, Control, and Dynamics, 2013, 36(4):1008-1018 doi:  10.2514/1.59294
    [4] Asmar S W, Armstrong J W, Iess L, et al. Spacecraft Doppler Tracking:Noise Budget and Accuracy Achievable in Precision Radio Science Observations[J]. Radio Science, 2005, 40(2):1-9 http://cn.bing.com/academic/profile?id=258f5f118e171c5378b34c0ac6a30a86&encoded=0&v=paper_preview&mkt=zh-cn
    [5] Armstrong J W, Iess L, Tortora P, et al. Stochastic Gravitational Wave Background:Upper Limits in the 10-6 to 10-3 Hz Band[J]. The Astrophysical Journal, 2003, 599(2):806 doi:  10.1086/apj.2003.599.issue-2
    [6] Bertotti B, Iess L, Tortora P. A Yest of General Re-lativity Using Radio Links with the Cassini Spacecraft[J]. Nature, 2003, 425(6956):374-376 doi:  10.1038/nature01997
    [7] Pätzold M, Häusler B, Tyler G L, et al. Mars Express 10 Years at Mars:Observations by the Mars Express Radio Science Experiment (MaRS)[J]. Planetary and Space Science, 2016, 127:44-90 doi:  10.1016/j.pss.2016.02.013
    [8] Anderson J D, Lau E L, Schubert G, et al. Gravity Inversion Considerations for Radio Doppler Data from the JUNO Jupiter Polar Orbiter[C]. Bulletin of the American Astronomical Society, San Diego, America, 2004
    [9] Moyer T D. Formulation for Observed and Computed Values of Deep Space Network Data Types for Na-vigation[M]. Hoboken:John-Wiley & Sons, 2003
    [10] 李济生.人造卫星精密轨道确定[M].北京:解放军出版社, 1995

    Li Jisheng. Satellite Precision Orbit Determination[M]. Beijing:Chinese People's Liberation Army Publishing House, 1995
    [11] 叶茂.月球探测器精密定轨软件研制与四程中继跟踪测量模式研究[D].武汉: 武汉大学, 2016

    Ye Mao.Development of Lunar Spacecraft Precision Orbit Determination Software System and Research on a Four-Way Relay Tracking Measurement Model[D]. Wuhan: Wuhan University, 2016
    [12] Tommei G, Milani A, Vokrouhlický D. Light-time Computations for the BepiColombo Radio Science Experiment[J]. Celestial Mechanics and Dynamical Astronomy, 2010, 107(1/2):285-298 https://core.ac.uk/display/50025659
    [13] Widrow B. Statistical Analysis of Amplitude-quantized Sampled-data Systems[J]. Transactions of the American Institute of Electrical Engineers, Part Ⅱ:Applications and Industry, 1961, 79(6):555-568 doi:  10.1109/TAI.1961.6371702
    [14] 曹建峰, 黄勇, 刘磊, 等.深空探测器三程多普勒建模与算法实现[J].宇航学报, 2017, 38(3):304-309 doi:  10.3873/j.issn.1000-1328.2017.03.011

    Cao Jianfeng, Huang Yong, Liu Lei, et al. Modeling and Algorithm Realization of Three-Way Doppler for Deep Space Exploration[J]. Journal of Astronautics, 2017, 38(3):304-309 doi:  10.3873/j.issn.1000-1328.2017.03.011
    [15] 叶茂, 李斐, 鄢建国, 等.国内外深空探测器精密定轨软件研究综述及WUDOGS简介[J].飞行器测控学报, 2017, 36(1):1-12 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb201701008

    Ye Mao, Li Fei, Yan Jianguo, et al. Review of Deep-space Orbit Determination Software and Introduction to WUDOGS[J]. Journal of Spacecraft TT&C Technology, 2017, 36(1):1-12 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb201701008
    [16] 黄勇, 胡小工, 曹建峰, 等.上海天文台火星卫星定轨软件系统[J].飞行器测控学报, 2009, 28(6):83-89 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb200906020

    Huang Yong, Hu Xiaogong, Cao Jianfeng, et al. The Mars Satellite Orbit Determination Software at Shanghai Astronomical Observatory[J]. Journal of Spacecraft TT&C Technology, 2009, 28(6):83-89 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb200906020
    [17] 胡松杰, 唐歌实.北京中心深空探测器精密定轨与分析软件系统[J].飞行器测控学报, 2010, 29(5):69-74 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb201005016

    Hu Songjie, Tang Geshi. BACC Orbit Determination and Analysis Software for Deep-Space Explorers[J]. Journal of Spacecraft TT&C Technology, 2010, 29(5):69-74 http://d.old.wanfangdata.com.cn/Periodical/fxqckxb201005016
    [18] 叶茂, 李斐, 鄢建国, 等.深空探测器精密定轨与重力场解算软件系统WUDOGS及其应用分析[J].测绘学报, 2017, 46(3):288-296 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201703004.htm

    Ye Mao, Li Fei, Yan Jianguo, et al. Wuhan University Deep-Space Orbit Determination and Gravity Recovery System WUDOGS and Its Application Analysis[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(3):288-296 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201703004.htm
    [19] 杨轩, 鄢建国, 叶茂, 等.火星探测器精密定轨软件研制及实测数据处理[J].武汉大学学报·信息科学版, DOI: 10.13203/j.whugis20160494

    Yang Xuan, Yan Jianguo, Ye Mao, et al. Development of Precise Orbit Determination Software for Mars Probe and Data Processing for MEX[J]. Geomatics and Information Science of Wuhan University, DOI: 10.13203/j.whugis20160494
    [20] Budnik F, Morley T A, Mackenzie R A. ESOC's System for Interplanetary Orbit Determination: Implementation and Operational Experience[C]. 18th International Symposium on Space Flight Dynamics, Munich, Germany, 2004
    [21] ANSI/IEEE Std 754-1985. IEEE Standard for Binary Floating-Point Arithmetic[S]. IEEE, New York, 1985
    [22] Standish E M. The JPL Planetary Ephemerides[J].Celestial Mechanics and Dynamical Astronomy, 1982, 26(2):181-186 http://d.old.wanfangdata.com.cn/OAPaper/oai_arXiv.org_1104.0042
    [23] Acton C H. Ancillary Data Services of NASA's Navigation and Ancillary Information Facility[J]. Planetary and Space Science, 1996, 44(1):65-70 doi:  10.1016/0032-0633(95)00107-7
    [24] Pavlis D E, Wimert J, McCarthy J. GEODYN Ⅱ System Description, NASA-CR-139149[R]. Greenbelt: SGT Inc., 2013
  • [1] 谭涵, 吴家齐.  星间单差模糊度固定的低轨卫星精密定轨精度分析 . 武汉大学学报 ● 信息科学版, 2022, 47(9): 1460-1469. doi: 10.13203/j.whugis20200385
    [2] 赵之若, 王少宇, 王心宇, 钟燕飞.  一种改进的火星车多光谱影像深度新颖目标探测方法 . 武汉大学学报 ● 信息科学版, 2022, 47(8): 1328-1335. doi: 10.13203/j.whugis20220119
    [3] 杨轩, 鄢建国, 叶茂, 李斐, 郝卫峰, 金炜桐.  火星探测器精密定轨软件研制及实测数据处理 . 武汉大学学报 ● 信息科学版, 2019, 44(3): 385-391. doi: 10.13203/j.whugis20160494
    [4] 刘山洪, 鄢建国, 杨轩, 叶茂, 金炜桐, 李斐.  水星探测器精密定轨软件研制及应用 . 武汉大学学报 ● 信息科学版, 2019, 44(4): 510-517. doi: 10.13203/j.whugis20170211
    [5] 王琰, 郭睿, 张传定, 胡小工, 唐成盼, 李冉, 潘军洋.  Bernese ECOM光压模型在BDS卫星精密定轨中的应用 . 武汉大学学报 ● 信息科学版, 2018, 43(2): 194-200. doi: 10.13203/j.whugis20160024
    [6] 宋敏, 袁运斌.  火星探测巡航段自主导航方法研究 . 武汉大学学报 ● 信息科学版, 2016, 41(7): 952-957. doi: 10.13203/j.whugis20140278
    [7] 叶茂, 李斐, 鄢建国, 郝卫峰, 邵先远.  GRAIL月球重力场模型定轨性能分析 . 武汉大学学报 ● 信息科学版, 2016, 41(1): 93-99. doi: 10.13203/j.whugis20150008
    [8] 马洋, 欧吉坤, 袁运斌, 霍星亮, 丁文武.  导航卫星天线相位中心变化估计及对LEO精密定轨影响 . 武汉大学学报 ● 信息科学版, 2015, 40(7): 894-900. doi: 10.13203/j.whugis20130626
    [9] 郝卫峰, 叶茂, 李斐, 张胜凯, 朱婷婷.  南极长城站建立深空测控站的可行性 . 武汉大学学报 ● 信息科学版, 2015, 40(10): 1360-1365. doi: 10.13203/j.whugis20130720
    [10] 张小红, 李盼, 左翔.  固定模糊度的精密单点定位几何定轨方法及结果分析 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1009-1013.
    [11] 李文文, 李敏, 施闯, 赵齐乐.  利用DORIS相位数据的Jason-2卫星精密定轨 . 武汉大学学报 ● 信息科学版, 2013, 38(10): 1207-1211.
    [12] 陈明, 唐歌实, 曹建峰, 张宇.  嫦娥一号绕月探测卫星精密定轨实现 . 武汉大学学报 ● 信息科学版, 2011, 36(2): 212-217.
    [13] 张守建, 李建成, 邹贤才, 金涛勇.  GRACE卫星非差运动学精密定轨分析 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 679-682.
    [14] 周晓青, 胡志刚, 张新远.  低轨卫星星载GNSS精密定轨的精度检核方法 . 武汉大学学报 ● 信息科学版, 2010, 35(11): 1342-1345.
    [15] 王正涛, 靳祥升, 党亚民, 姜卫平.  低轨卫星精密定轨的初轨向量与力模型参数数值积分误差分析 . 武汉大学学报 ● 信息科学版, 2009, 34(6): 728-731.
    [16] 施闯, 李敏, 楼益栋, 邹蓉.  利用区域基准站进行导航卫星近实时精密定轨研究 . 武汉大学学报 ● 信息科学版, 2008, 33(7): 697-700.
    [17] 李敏, 赵齐乐, 葛茂荣.  GIOVE-A卫星精密定轨仿真研究 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 818-820.
    [18] 赵齐乐, 施闯, 柳响林, 葛茂荣.  重力卫星的星载GPS精密定轨 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 810-814.
    [19] 赵齐乐, 刘经南, 葛茂荣, 施闯.  CHAMP卫星cm级精密定轨 . 武汉大学学报 ● 信息科学版, 2006, 31(10): 879-882.
    [20] 刘经南, 赵齐乐, 张小红.  CHAMP卫星的纯几何定轨及动力平滑中的动力模型补偿研究 . 武汉大学学报 ● 信息科学版, 2004, 29(1): 1-6.
  • 加载中
图(3) / 表(4)
计量
  • 文章访问数:  947
  • HTML全文浏览量:  95
  • PDF下载量:  199
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-04-28
  • 刊出日期:  2018-10-05

行星际深空探测中双程测速的高精度计算方法

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

    国家自然科学基金 41374024

    国家自然科学基金 41604004

    国家自然科学基金 41174019

    湖北省自然科学基金重点项目创新群体类项目 2015CFA011

    宇航动力学国家重点实验室开放基金 2016ADL-DW0103

    中国博士后科学基金 2016M602360

    作者简介:

    金炜桐, 硕士生, 主要从事行星探测器精密定轨方面的研究。jwt_whdx@whu.edu.cn

    通讯作者: 叶茂, 博士。mye@whu.edu.cn
  • 中图分类号: P228

摘要: 对于行星际深空探测(距地球1亿km以上)任务,由于受到计算机字长的限制,传统双程测速模型的计算精度无法满足高精度定轨的需要,其最大误差源于多普勒频移周计数终点和始点上行几何距离之间和下行几何距离之间差分值的计算过程。对此建立行星际双程测速模型,高精度地计算了两个差分值,推导模型的计算公式并给出详细步骤,同时给出计算过程中需要的切比雪夫差分多项式递推公式的形式。将该模型在深空探测器精密定轨与重力场解算软件系统(Wuhan University deep-space orbit determination and gravity recovery system,WUDOGS)中进行了实现,并以欧空局火星快车号(Mars express,MEX)探测任务为背景,利用该软件进行仿真测试,从计算精度和定轨结果两个方面验证该模型的优越性。结果表明,该模型将双程测速的计算值在计算机中表达的精度提高2个数量级,同时避免了定轨过程中引入额外的数值误差,可以为后续高精度的行星际深空探测任务的定轨提供参考。

English Abstract

金炜桐, 李斐, 杨轩, 鄢建国, 叶茂, 郝卫峰, 曲春凯. 行星际深空探测中双程测速的高精度计算方法[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1483-1489. doi: 10.13203/j.whugis20170020
引用本文: 金炜桐, 李斐, 杨轩, 鄢建国, 叶茂, 郝卫峰, 曲春凯. 行星际深空探测中双程测速的高精度计算方法[J]. 武汉大学学报 ● 信息科学版, 2018, 43(10): 1483-1489. doi: 10.13203/j.whugis20170020
JIN Weitong, LI Fei, YANG Xuan, YAN Jianguo, YE Mao, HAO Weifeng, QU Chunkai. High Precision Computational Method of Two-Way Range-Rate in Long-Distance Deep Space Exploration[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1483-1489. doi: 10.13203/j.whugis20170020
Citation: JIN Weitong, LI Fei, YANG Xuan, YAN Jianguo, YE Mao, HAO Weifeng, QU Chunkai. High Precision Computational Method of Two-Way Range-Rate in Long-Distance Deep Space Exploration[J]. Geomatics and Information Science of Wuhan University, 2018, 43(10): 1483-1489. doi: 10.13203/j.whugis20170020
  • 双程多普勒测速(以下简称双程测速)是行星探测器精密定轨中最为常用的观测模式之一,也是无线电科学实验(radio science experiment)中最主要的数据源[1-2]

    随着深空网结构的优化、频段的升级和服务能力的提高,加之先进的校正技术[3]和对来自于不同误差源噪声的详细建模[4],双程测速的测量精度越来越高。国外的部分探测任务诸如“卡西尼-惠更斯”、“朱诺”和欧空局火星快车号(Mars express, MEX)等在其特定的任务时段,双程测速精度均可以达到10-2 ~10-3 mm/s[5-8]

    精密定轨的精度同时也与模型理论值的计算精度有关。传统的双程测速的理论值由多普勒频移周计数始点的双程平均距离到周计数终点的双程平均距离之平均变化率计算[9-11]。在行星际的探测任务中,由于探测器距离地面测站遥远(距地球1亿km以上), 且计算机浮点数模型的精度有限,距离量和必要的距离差分量在计算机中只能精确表达至毫米量级,从而导致测速值的计算精度无法满足高精度定轨的需要。

    在国外的研究中,文献[12]使用7节点的高斯求积公式数值求解距离差分量。对于典型的30 s多普勒计数间隔,该方法最终可以将双程测速值的计算精度提高至2×10-3 mm/s。但是该文献中并没有给出其详细的解法与计算公式。文献[3]基于Widrow提出的量化误差模型[13]对数值噪声进行建模并给出其统计特性,但并没有定量地讨论减少计算误差的策略。在国内的研究中,曹建峰等提出了一种间接使用中心天体、航天器的速度与加速度进行多普勒建模的方法[14],但相关文献仍旧较少。

    与此同时,国际上很多优秀定轨软件对国内用户限制[15-17],无法与其进行更加详细的交叉对比验证。一种通用的方法是将定轨软件中的浮点数均改为四精度存储,但是由于定轨软件架构复杂,更改浮点变量的精度操作繁琐,且会极大增加内存占用和时间复杂度[3]

    本文建立行星际双程测速模型(以下简称星际模型),高精度地计算双程测速的理论值,减少传统双程测速模型(以下简称传统模型)在计算过程中丢失的有效位,并将该模型在武汉大学自主研发的深空探测器精密定轨与重力场解算软件系统(Wuhan University deep-space orbit determination and gravity recovery system, WUDOGS)中进行了实现[18-19],并以MEX探测任务为背景,利用该软件进行仿真测试,从计算精度和模拟定轨结果两个方面验证该模型的优越性。这一工作可以应用到后续高精度的火星、小行星等探测任务中。

    • Te3Ts3分别表示多普勒周计数的结束时刻和起始时刻,若多普勒测速的时标位于Te3,则相应的载波信号被探测器在Te2时刻转发,地面测站在Te1时刻发射;同样,Ts3时刻对应的载波信号被探测器在Ts2时刻转发,地面测站在Ts1时刻发射。忽略探测器应答机的转发时延,则双程测速的理论计算值ρ(Te3)可以表示为:

      $$ \rho \left( {{T_{e3}}} \right) = {\rho ^g}\left( {{T_{e3}}} \right) + {\rho ^r}\left( {{T_{e3}}} \right) + {\rho ^a}\left( {{T_{e3}}} \right) + {\rho ^t}\left( {{T_{e3}}} \right) $$ (1)

      式中, ρg(Te3)为几何双程测速值; ρr(Te3)为广义相对论时延等效距离变率改正; ρa(Te3)为大气时延等效距离变率改正; ρt(Te3)为时间尺度不同造成的等效距离变率改正。且各项满足:

      $$ \left\{ \begin{array}{l} {\rho ^g}\left( {{T_{e3}}} \right) = \frac{1}{{2{t_c}}}\left[ {\left( {\rho _4^g - \rho _2^g} \right) + \left( {\rho _3^g - \rho _1^g} \right)} \right]\\ {\rho ^r}\left( {{T_{e3}}} \right) = \frac{1}{{2{t_c}}}\left[ {\left( {\rho _4^r - \rho _2^r} \right) + \left( {\rho _3^r - \rho _1^r} \right)} \right]\\ {\rho ^a}\left( {{T_{e3}}} \right) = \frac{1}{{2{t_c}}}\left[ {\left( {\rho _4^a - \rho _2^a} \right) + \left( {\rho _3^a - \rho _1^a} \right)} \right]\\ {\rho ^t}\left( {{T_{e3}}} \right) = \frac{1}{{2{t_c}}}\left( {c\Delta {T_e} - c\Delta {T_s}} \right) \end{array} \right. $$ (2)

      式中,ρ2gρ4g为下行几何距离;ρ1gρ3g为上行几何距离;c为光速;ρ4rρ2rρ3rρ1r为广义相对论时延的等效距离改正;ρ4aρ2aρ3aρ1a为大气介质时延的等效距离改正;ΔTe和ΔTs代表收发信号时刻(barycentric dynamical time, TDB)和世界标准时间(coordinated universal time, UTC)两种时间尺度之差; tc为多普勒积分间隔。

    • 模型的计算精度(model accuracy)至少需要比观测量的随机噪声水平小一个量级[20]。以MEX任务为例,10 s采样率的双程测速理论值的模型计算精度至少需要达到0.005 mm/s的量级才能够满足高精度定轨的需要。而根据IEEE754浮点数标准,最常用的双精度浮点数在计算机内存中分配的二进制数位为64位,1位存储符号位,11位存储指数位,52位存储有效数位,相当于十进制的[lg(252+1)]=15位有效数字。有效数位不能表达的部分按照最近舍入法则进行舍入或者截断[21]

      根据式(1)和式(2),若想计算ρ(Te3),首先需要计算7个差分值:ρ4gρ2gρ3gρ1gρ4rρ2rρ3rρ1rρ4aρ2aρ3aρ1acΔTecΔTs,其中误差最大的是前两个差分值的计算。以火星任务为例,由于地面测站至探测器的距离至少达到1011 m的量级[19]ρ4gρ2gρ3gρ1g的整数部分至少占用12个十进制有效位,小数部分至多占用3个十进制有效位,计算机至多能有效表达至毫米量级,将两数直接相减,则ρ4gρ2gρ3gρ1g两个差分值也至多能有效表达至毫米量级;若多普勒积分间隔为10 s,则ρg(Te3)至多能有效表达至0.1 mm/s量级。这样即使式(1)中的后3项可以在计算机中表达的十分精确,计算值ρ(Te3)也只能有效表达至0.1 mm/s量级,简单的采用式(1)和(2),显然无法满足行星探测器高精度定轨的需求。

    • 星际双程测速模型的实质在于规避两个量级在1011及以上的大数相减的情况,从而高精度地计算差分值ρ4gρ2gρ3gρ1g

    • ρ4gρ2g为例,给出如下高精度的计算方法。

      按照式(3)计算ρ2gρ4g[14]

      $$ \left\{ \begin{array}{l} \rho _2^g = \left| {\left[ {{\mathit{\boldsymbol{X}}_M}\left( {{T_{s2}}} \right) + {\mathit{\boldsymbol{X}}_S}\left( {{T_{s2}}} \right)} \right] - } \right.\\ \;\;\;\;\;\;\left. {\left[ {{\mathit{\boldsymbol{X}}_B}\left( {{T_{s3}}} \right) + {\mathit{\boldsymbol{X}}_E}\left( {{T_{s3}}} \right) + {\mathit{\boldsymbol{X}}_A}\left( {{T_{s3}}} \right)} \right]} \right|\\ \rho _4^g = \left| {\left[ {{\mathit{\boldsymbol{X}}_M}\left( {{T_{e2}}} \right) + {\mathit{\boldsymbol{X}}_S}\left( {{T_{e2}}} \right)} \right] - } \right.\\ \;\;\;\;\;\;\left. {\left[ {{\mathit{\boldsymbol{X}}_B}\left( {{T_{e3}}} \right) + {\mathit{\boldsymbol{X}}_E}\left( {{T_{e3}}} \right) + {\mathit{\boldsymbol{X}}_A}\left( {{T_{e3}}} \right)} \right]} \right| \end{array} \right. $$ (3)

      式中,XM(Te2)和XM(Ts2)为Te2Ts2时刻中心天体质心相对于太阳系质心的位置;XS(Te2)和XS(Ts2)为Te2Ts2时刻探测器相对于中心天体质心的位置;XB(Te3)和XB(Ts3)为Te3Ts3时刻地月质心相对于太阳系质心的位置;XE(Te3)和XE(Ts3)为Te3Ts3时刻地球质心相对地月质心的位置;XA(Te3)和XA(Ts3)为Te3Ts3时刻地面测站的地心天球坐标。

      将向量用坐标分量形式展开,引入以下表示:ΔXSXS, ΔYS, ΔZS)、ΔXMXM, ΔYM, ΔZM)、ΔXB2XB2, ΔYB2, ΔZB2)、ΔXE2XE2、ΔYE2, ΔZE2)和ΔXA2XA2、ΔYA2, ΔZA2)。对于ΔXS和ΔXM的坐标分量,分别由对应的Te2Ts2时刻的坐标分量相减得到;对于ΔXB2、ΔXE2、ΔXA2则分别由对应的Te3Ts3时刻的坐标分量相减得到。进一步令:

      $$ \left\{ \begin{array}{l} \Delta {X_2} = \Delta {X_M} + \Delta {X_S} - \Delta {X_{B2}} - \Delta {X_{E2}} - \Delta {X_{A2}}\\ \Delta {Y_2} = \Delta {Y_M} + \Delta {Y_S} - \Delta {Y_{B2}} - \Delta {Y_{E2}} - \Delta {Y_{A2}}\\ \Delta {Z_2} = \Delta {Z_M} + \Delta {Z_S} - \Delta {Z_{B2}} - \Delta {Z_{E2}} - \Delta {Z_{A2}}\\ {R_{x2}} = {X_M}\left( {{T_{s2}}} \right) + {X_S}\left( {{T_{s2}}} \right) - {X_B}\left( {{T_{s3}}} \right) - \\ \;\;\;\;\;\;\;\;{X_E}\left( {{T_{s3}}} \right) - {X_A}\left( {{T_{s3}}} \right)\\ {R_{y2}} = {Y_M}\left( {{T_{s2}}} \right) + {Y_S}\left( {{T_{s2}}} \right) - {Y_B}\left( {{T_{s3}}} \right) - \\ \;\;\;\;\;\;\;\;{Y_E}\left( {{T_{s3}}} \right) - {Y_A}\left( {{T_{s3}}} \right)\\ {R_{z2}} = {Z_M}\left( {{T_{s2}}} \right) + {Z_S}\left( {{T_{s2}}} \right) - {Z_B}\left( {{T_{s3}}} \right) - \\ \;\;\;\;\;\;\;\;{Z_E}\left( {{T_{s3}}} \right) - {Z_A}\left( {{T_{s3}}} \right) \end{array} \right. $$ (4)

      可以得到:

      $$ \rho _4^g = f\left( {\Delta {X_2},\Delta {Y_2},\Delta {Z_2}} \right) = \sqrt {{{\left( {\rho _2^g} \right)}^2} + 2\left( {{R_{x2}}\Delta {X_2} + {R_{y2}}\Delta {Y_2} + {R_{z2}}\Delta {Z_2}} \right) + \Delta X_2^2 + \Delta Y_2^2 + \Delta Z_2^2} $$ (5)

      由于距离较远,两时刻距离之间的差分值相对于距离本身为小量,故在Ts3时刻处,将式(5)沿Ts3时刻对应的下行支路方向进行泰勒展开。

      令:

      $$ \left\{ \begin{array}{l} {\mathit{\boldsymbol{D}}_2} = \Delta {X_2}{\bf{i}} + \Delta {Y_2}{\bf{j}} + \Delta {Z_2}{\bf{k}}\\ {a_2} = \frac{{\Delta {X_2}}}{{\left| {{\mathit{\boldsymbol{D}}_2}} \right|}},{b_2} = \frac{{\Delta {Y_2}}}{{\left| {{\mathit{\boldsymbol{D}}_2}} \right|}},{c_2} = \frac{{\Delta {Z_2}}}{{\left| {{\mathit{\boldsymbol{D}}_2}} \right|}} \end{array} \right. $$ (6)

      Ts3时刻对应的下行支路径向方向可以用向量D2的方向向量α2(a2, b2, c2)表示。

      f1f2f3分别代表Ts3时刻f沿α2方向的一阶、二阶和三阶方向导数值。则式(5)可以写成如下形式:

      $$ \begin{array}{*{20}{c}} {\rho _4^g - \rho _2^g = {f_1} \cdot \left| {{\mathit{\boldsymbol{D}}_2}} \right| + \frac{1}{{2!}}{f_2} \cdot {{\left| {{\mathit{\boldsymbol{D}}_2}} \right|}^2} + }\\ {\frac{1}{{3!}}{f_3} \cdot {{\left| {{\mathit{\boldsymbol{D}}_2}} \right|}^3}} \end{array} $$ (7)

      考虑到实际求解精度,式(7)中忽略了比三阶高的无穷小O(|D2|3)项。

      式(7)中3个方向导数的计算公式为:

      $$ \left\{ \begin{array}{l} {f_1} = \frac{{{w_2}}}{{\rho _2^g}}\\ {f_2} = \frac{{{{\left( {\rho _2^g} \right)}^2} - w_2^2}}{{{{\left( {\rho _2^g} \right)}^3}}}\\ {f_3} = \frac{{ - 3{{\left( {\rho _2^g} \right)}^2}{w_2} + 3w_2^3}}{{{{\left( {\rho _2^g} \right)}^5}}} \end{array} \right. $$ (8)

      式中,w2=a2Rx2+b2Ry2+c2Rz2,代入式(7),即可算出ρ4gρ2g的值。

      对于测距上行支路,将式(3)至式(8)相应的下行量替换为上行量即可计算出ρ3gρ1g的值。

    • 在计算几何距离差分值的过程中涉及到不同时刻天体星历差分值的计算。对于火星、地球和地月质心本身而言,其距离太阳系质心遥远,若直接将两时刻对应的位置相减,同样会引起上述提到的数值误差。因此,本文采用切比雪夫多项式差分的方式计算ΔXM、ΔXB和ΔXE

      根据历表插值的方法[22],两个时刻同一天体位置之差可以表示为:

      $$ \Delta {R_i} = \sum\limits_{k = 0}^{N - 1} {\left[ {{C_k}\left( {{\tau _2}} \right){P_{i,k}}\left( {{\tau _2}} \right) - {C_k}\left( {{\tau _1}} \right){P_{i,k}}\left( {{\tau _1}} \right)} \right]} $$ (9)

      式中, N代表切比雪夫多项式的阶数; τ为按照对应时刻所在的时间区间标准化至(-1, 1)范围内的时刻; Ck(τ)为τ对应的切比雪夫多项式系数; Pi, k(τ)为τ对应的各坐标分量的切比雪夫多项式。

      令$\left\{ \begin{array}{l} \Delta {C_k} = {C_k}\left( {{\tau _2}} \right) - {C_k}\left( {{\tau _1}} \right)\\ \Delta {P_{i,k}} = {P_{i,k}}\left( {{\tau _2}} \right) - {P_{i,k}}\left( {{\tau _1}} \right) \end{array} \right.$, 则有:

      $$ \Delta {R_i} = \sum\limits_{k = 0}^{N - 1} {\left[ {{C_k}\left( {{\tau _1}} \right)\Delta {P_{i,k}} + \Delta {C_k}{P_{i,k}}\left( {{\tau _2}} \right)} \right]} $$ (10)

      式中,ΔPi, k按照如下递推公式计算:

      $$ \left\{ \begin{array}{l} \Delta {P_{i,0}} = 0;\Delta {P_{i,1}} = {\tau _2} - {\tau _1}\\ \Delta {P_{i,k}} = 2{\tau _2}\Delta {P_{i,k - 1}} + 2\left( {{\tau _2} - {\tau _1}} \right)\Delta {P_{i,k - 1}}\left( {{\tau _1}} \right) - \\ \;\;\;\;\;\;\;\;\;\Delta {P_{i,k - 2}} \end{array} \right. $$ (11)

      式中, k=2, 3…N-1。

      t2t1两个时刻位于同一区间内,ΔCkPi, k(τ2)=0。式(10)可以写为:

      $$ \Delta {R_i} = \sum\limits_{k = 0}^{N - 1} {{C_k}\left( {{\tau _1}} \right)\Delta {P_{i,k}}} $$ (12)

      至此,从式(3)至式(12),严格避免了两个数量级在1011及以上的大数字相减,ρ4gρ2gρ3gρ1g的值可以通过该高精度算法计算得出。

    • 不考虑式(1)中后3项具体计算细节,给出星际模型计算ρg(Te3)的详细步骤,如图 1所示。

      图  1  星际模型详细计算步骤

      Figure 1.  Detailed Calculation Steps of Interplanetary Model

    • 双程测速模型计算值的具体策略描述如下:

      1) 利用SPICE[23]读取某时刻MEX探测器的位置和速度作为初始轨道根数;

      2) 给定力模型、积分弧段和积分步长,积分出MEX探测器的轨道;

      3) 将积分出的轨道分别代入传统和星际双程测速模型,生成两模型的计算值。

      为了定量分析两种模型的计算精度,需要一个“准确”的理论值作为参考。其中,美国戈达德宇航中心的GEODYN-Ⅱ软件系统广泛地应用于行星探测器的精密定轨、测地学参数的估计、卫星轨道预报、仪器校准等深空探测领域,经受过数次行星际实测任务的考验[15, 18-19, 24],可以认为GEODYN-Ⅱ的双程测速模型可靠准确,基本参数设置参照文献[19],将GEODYN-Ⅱ生成的计算值作为“准确”的理论计算值。

      将两模型的计算值分别与GEODYN-Ⅱ生成的计算值作差,得到计算值的残差。图 2分别给出了两模型的计算值残差图,表 1给出了残差的统计信息,RMS表示均方根误差。由于测试中均使用同一组轨道且策略严格保持一致,所以残差可以直接显示出两模型对于计算机浮点误差的敏感程度。对于传统模型,残差均值为-0.007 005 mm/s,超过计算值精度需要的0.005 mm/s,计算机的浮点误差最大达到0.111 22 mm/s,超过了MEX 10 s采样率的测量精度,残差的RMS为0.040 993 mm/s,这表明残差围绕均值在十分接近于测量精度的范围内分布,传统测速模型对浮点误差具有较强的敏感性。这也同样对应于浮点数运算规则的特性:在两个相近的大数相减时,有效位数将会大量缺失,而最近舍入法则会让这种数值噪声呈现一定的随机性[3]。对于星际模型,计算机浮点误差的量级最大为5×10-3 mm/s,均值为-0.003 43 mm/s,而RMS比均值还要小一个量级。在现有的观测精度下,该量级的数值噪声完全可以忽略不计,这是因为星际模型规避了对浮点误差最具敏感性的大数相减运算。距离差分值通过在多普勒计数起始点沿支路方向泰勒展开,从而转化为“小项”相加得到。星际模型能够将计算值的精度提高2倍,在一定程度上保证了计算值的“纯净度”。

      图  2  两种模型计算值残差图

      Figure 2.  Residual Plot of Calculation Values Using Two Models

      表 1  计算值残差统计信息

      Table 1.  Statistical Information of Calculation Residuals

      模型 计算值残差/mm·s-1
      最大值 均值 RMS
      传统双程测速 0.111 22 -0.007 005 0.040 993
      星际双程测速 5.000 0×10-3 -3.429 86×10-3 9.419 85×10-4

      为定量分析数值误差对定轨的影响,使用GEODYN-Ⅱ生成一组3 d弧段的模型计算值作为输入到WUDOGS中的真实观测值,在初始轨道根数的位置分量上分别加上-100 m、+100 m、+100 m的偏差,分别采用两种模型进行模拟定轨,得到最后一次迭代的观测值残差图,如图 3所示。表 2给出了两种模型模拟定轨的初轨改正数,表 3表 4分别给出了两种模型的模拟定轨结果与给定初始轨道根数的差值以及最后一次迭代的观测值残差统计信息。

      图  3  两种模型模拟定轨最后一次迭代残差图

      Figure 3.  Post-fit Residuals of Orbit Determination Using Two Models

      表 2  两种模型模拟定轨的初轨改正数

      Table 2.  Initial Orbit Corrections of POD Using Two Models, Respectively

      模型 ΔPx/m ΔPy/m ΔPz/m ΔVx/m·s-1 ΔVy/m·s-1 ΔVz/m·s-1
      传统模型 100.097 80 -100.104 90 -99.798 70 2.67×10-5 -6.32×10-7 -1.29×10-5
      星际模型 99.999 98 -99.999 98 -100.000 03 -4.89×10-9 2.58×10-11 2.29×10-9

      表 3  两种模型模拟定轨结果与给定初轨的差值

      Table 3.  Initial Orbit Differences Between POD Simulation Results and Given Results Using Two Models, Respectively

      模型 ΔPx/m ΔPy/m ΔPz/m ΔVx/m·s-1 ΔVy/m·s-1 ΔVz/m·s-1
      传统模型 0.097 809 -0.104 888 0.201 331 2.670×10-5 -6.317×10-7 -1.287×10-5
      星际模型 -1.764×10-5 1.934×10-5 -3.640×10-5 -4.890×10-5 2.581×10-11 2.294×10-9

      表 4  两种模型模拟定轨的最后一次迭代的观测值残差统计信息/mm·s-1

      Table 4.  Statistical Results of POD Measurement Residuals During the Last Iteration Using Two Models, Respectively/mm·s-1

      模型 最后一次迭代观测值残差均值 最后一次迭代观测值残差RMS
      传统模型 -0.004 1 0.003 1
      星际模型 4.589 5×10-8 4.768 8×10-7

      由于模拟定轨中没有引入其他形式的误差,理论上,若计算不存在浮点误差,最后一次迭代的观测值残差应接近于0,且初轨位置分量的改正数应接近于+100 m、-100 m和-100 m,但是数值噪声的存在会造成与理论结果的不一致性。据表 3表 4显示,传统模型的定轨结果与给定初始轨道位置的差值在分米量级,最后一次迭代收敛后残差的均值和RMS分别为-0.004 1 mm/s和0.003 1 mm/s,星际模型的定轨结果与给定初始轨道位置的差值在10-5 m量级,最后一次迭代收敛后残差的均值和RMS分别为4.589 5×10-8 mm/s和4.768 8×10-7 mm/s,且从图 3中可以看出,残差量级最大不超过2×10-6 mm/s,可以认为星际模型中的浮点误差成分极其微弱,在定轨过程中不会引入过大的数值噪声,即星际模型能够更好地提取观测值中的有用信息,使残差成分在最小二乘的意义下更“纯净”。

    • 针对距离相对遥远的行星际精密定轨任务,本文建立星际双程测速模型,高精度地计算了两个几何距离差分值,得到如下结论:

      1) 星际双程测速模型能避免两个相近的大数相减,从而高精度地计算多普勒周计数终点和始点上行与下行几何距离之间的差分值,避免传统双程测速模型计算机浮点数误差过大的问题。

      2) 各大天体的星历在短时间间隔内的差分值需由切比雪夫差分多项式及其相应系数的线性组合计算,同样可以避免计算机浮点数误差过大的问题,是星际双程测速模型的一个必要步骤。

      3) 由于计算机浮点数模型本身无法精确表达大部分实数,星际双程测速模型只能尽可能地减弱计算机浮点误差的影响,并不能完全消除。根据本文的仿真结果,星际双程测速模型将模型计算值在计算机中表达的精度提高2个量级,避免定轨过程中引入过大的数值误差,可以为后续行星际深空探测任务的高精度定轨提供参考。

参考文献 (24)

目录

    /

    返回文章
    返回