中文核心期刊

中国科技核心期刊

中国科学引文数据库(CSCD)来源期刊

中国高校优秀科技期刊

中国宇航学会深空探测技术专业委员会会刊

高级检索

留言板

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

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

“嫦娥四号”中继星再生伪码测距数据定轨精度分析

牛东文 段建锋 欧阳琦 张宇 陈略 王美

牛东文, 段建锋, 欧阳琦, 张宇, 陈略, 王美. “嫦娥四号”中继星再生伪码测距数据定轨精度分析[J]. 深空探测学报(中英文), 2022, 9(1): 21-28. doi: 10.15982/j.issn.2096-9287.2022.20191213002
引用本文: 牛东文, 段建锋, 欧阳琦, 张宇, 陈略, 王美. “嫦娥四号”中继星再生伪码测距数据定轨精度分析[J]. 深空探测学报(中英文), 2022, 9(1): 21-28. doi: 10.15982/j.issn.2096-9287.2022.20191213002
NIU Dongwen, DUAN Jianfeng, OUYANG Qi, ZHANG Yu, CHEN Lue, WANG Mei. Regenerative Pseudo-Random Code Ranging Orbit Determination Accuracy Analysis for Chang’E-4 Relay Satellite[J]. Journal of Deep Space Exploration, 2022, 9(1): 21-28. doi: 10.15982/j.issn.2096-9287.2022.20191213002
Citation: NIU Dongwen, DUAN Jianfeng, OUYANG Qi, ZHANG Yu, CHEN Lue, WANG Mei. Regenerative Pseudo-Random Code Ranging Orbit Determination Accuracy Analysis for Chang’E-4 Relay Satellite[J]. Journal of Deep Space Exploration, 2022, 9(1): 21-28. doi: 10.15982/j.issn.2096-9287.2022.20191213002

“嫦娥四号”中继星再生伪码测距数据定轨精度分析

doi: 10.15982/j.issn.2096-9287.2022.20191213002
基金项目: 国家自然科学基金资助项目(11973015)
详细信息
    作者简介:

    牛东文(1986– ),男,工程师,主要研究方向:航天测控、航天器测定轨技术。通讯地址:北京市海淀区北清路26号院(100094)电话:(010)66360292E-mail:ndw05044050@163.com

    张宇(1979– ),男,高级工程师,主要研究方向:航天器测定轨技术。通讯地址:北京5130信箱(100094)电话:(010)66361456E-mail:zackyzy@163.com

  • ● China’s first regenerative pseudo-random code ranging experiment beyond the distance between the Earth and the Moon is introduced. ● Based on the actual track measurement data, the orbit determination accuracy of the regenerative pseudo-random code ranging is analyzed and evaluated. ● Based on the actual orbit measurement data, through the precise orbit determination analysis, it is quantitatively verified and further illustrated that the measurement accuracy and orbit determination accuracy of the regenerative pseudo-random code ranging are better than those of the sidetones ranging.
  • 中图分类号: V412.41

Regenerative Pseudo-Random Code Ranging Orbit Determination Accuracy Analysis for Chang’E-4 Relay Satellite

  • 摘要: “嫦娥四号”中继星目前正稳定运行在地月L2点使命轨道,期间进行了再生伪码测距试验,试验期间交叉进行了再生伪码测距和侧音测距。利用“嫦娥四号”中继星在地月L2点使命轨道的再生伪码测距数据和侧音测距数据,进行了定轨精度评估。结果显示再生伪码测距数据的均方根误差(Root Mean Square Error,RMS)优于侧音测距数据RMS一个量级,同时再生伪码测距数据定轨和预报精度优于侧音测距数据一倍。这对再生伪码测距在中国后续深空探测任务中的应用具有参考意义。
    Highlights
    ● China’s first regenerative pseudo-random code ranging experiment beyond the distance between the Earth and the Moon is introduced. ● Based on the actual track measurement data, the orbit determination accuracy of the regenerative pseudo-random code ranging is analyzed and evaluated. ● Based on the actual orbit measurement data, through the precise orbit determination analysis, it is quantitatively verified and further illustrated that the measurement accuracy and orbit determination accuracy of the regenerative pseudo-random code ranging are better than those of the sidetones ranging.
  • 图  1  再生伪码测距系统结构图

    Fig.  1  Chart of regenerative pseudo-random code ranging

    图  2  再生测距陶思沃斯码

    Fig.  2  Regenerative ranging Tausworthe code

    图  3  下行链路再生伪码测距与侧音测距信号功率比

    Fig.  3  Power ratio of downlink regenerative pseudo-random code ranging to sidetone ranging signal

    表  1  “鹊桥”中继星再生伪码测距试验时段(北京时间)

    Table  1  Regenerative pseudo-random code ranging experiment schedule(BJT)

    日期测站观测时长(hh:mm)
    2019.06.14纳米比亚站2:49
    阿根廷深空站2:59
    佳木斯深空站3:03
    喀什深空站3:12
    2019.06.15喀什深空站1:57
    2019.06.16纳米比亚站2:59
    阿根廷深空站2:54
    2019.06.17佳木斯深空站2:46
    2019.06.18喀什深空站2:58
    纳米比亚站3:23
    2019.06.19纳米比亚站2:23
    阿根廷深空站2:59
    佳木斯深空站1:49
    2019.06.20佳木斯深空站1:10
    喀什深空站2:50
    下载: 导出CSV

    表  2  定轨分析策略

    Table  2  Strategy for the orbit determination

    动力学模型参数
    地球质点J2000地心惯性系
    地球非球形JGM-38×8阶次
    月球质点第三体引力JPL-DE421历表
    月球非球形Grail6608×8阶次
    太阳质点第三体引力JPL-DE421历表
    其它大行星第三体引力JPL-DE421历表
    太阳光压球形模型
    固定面质比
    锥形阴影,Cr=1.7
    光压等效面积8.3 m2
    下载: 导出CSV

    表  3  不同策略定轨数据质量评估

    Table  3  Orbit determination data quality evaluation for different strategies

    策略测控站数据类型RMS/m系统差/m
    1. 侧音测距+测速+VLBI纳米比亚站测距和1.78018.750
    佳木斯深空站测距0.612
    阿根廷深空站测距和1.56019.797
    喀什深空站测距和0.87628.367
    2. 侧音测距+测速纳米比亚站测距和1.830–22.718
    佳木斯深空站测距2.360
    阿根廷深空站测距和3.70025.121
    喀什深空站测距和0.90736.700
    3. 伪码测距+测速+VLBI纳米比亚站测距和0.47219.692
    佳木斯深空站测距0.470
    阿根廷深空站测距和0.37518.230
    喀什深空站测距和0.368310.10
    4. 伪码测距+测速纳米比亚站测距和0.46119.262
    佳木斯深空站测距0.452
    阿根廷深空站测距和0.394–17.50
    喀什深空站测距和0.386310.10
    下载: 导出CSV

    表  4  不同策略定轨结果与精密轨道比较统计

    Table  4  Statistics on comparisons of orbit determination results of different strategies with the precision orbit

    策略误差项定轨误差预报误差
    位置/m速度/(m·s–1位置/m速度/(m·s–1
    MAXRMSMAXRMSMAXRMSMAXRMS
    1. 侧音测距+测速+VLBIR–10.3663.223–0.001 480.000 65974.01322.6240.002 8270.000 475
    T455.685262.395–0.001 370.000 083–680.08471.394–0.00 1540.000 059
    N–627.289184.754–0.002 430.001 427–370.460114.9600.002 6260.000 200
    Total637.548320.9290.003 060.001 574684.546137.2030.004 1560.000 519
    2. 侧音测距+测速R–26.71112.4260.000 8750.000 39460.67221.8360.002 0050.000 323
    T–271.719152.388–0.001 0790.000 087–494.13864.174–0.001 490.000 124
    N590.413371.482–0.002 340.000 690–583.08616.424–0.001 120.000 590
    Total643.014401.7160.002 4960.000 800746.17469.7490.002 6640.000 684
    3. 伪码测距+测速+VLBIR–14.0804.3330.000 2820.000 07844.41610.4010.000 5210.000 067
    T–61.84923.284–0.000 1470.000 048–70.3672.3010.000 0130.000 006
    N–307.383112.071–0.001 2030.000 648–226.67153.6580.001 3010.000 143
    Total312.021114.5470.001 2120.000 655235.09454.7050.001 4020.000 158
    4. 伪码测距+测速R–13.7784.4260.000 3480.000 10051.47812.1510.000 6820.000 093
    T–82.64435.100–0.000 2050.000 047–105.7216.444–0.000 1030.000 011
    N–335.438121.933–0.0013 110.000 709–246.35558.6000.001 4230.000 156
    Total341.095126.9620.0013 210.000 717260.00160.1930.001 5810.000 182
    下载: 导出CSV
  • [1] 段建锋,陈明,张宇,等. “嫦娥四号”中继星使命轨道段定轨计算与分析[J]. 深空探测学报(中英文),2019,6(3):247-253.

    DUAN J F,CHEN M,ZHANG Y,et al. Orbit determination and analysis of Chang’E-4 relay satellite on mission orbit[J]. Journal of Deep Space Exploration,2019,6(3):247-253.
    [2] 唐歌实. 深空测控无线电测量技术[M]. 北京: 国防工业出版社, 2012.
    [3] 巨兰,张碧雄,张晓明. 深空网伪码测距技术研究[J]. 飞行器测控学报,2005,24(6):9-13.

    JU L,ZHANG B X,ZHANG X M. Research on pseudo-noise ranging technologies for DSN network[J]. Journal of Spacecraft TT&C Technology,2005,24(6):9-13.
    [4] 王琦,吴斌. 航天测控系统中伪码测距精度分析[J]. 无线电工程,2009,39(1):39-44. doi:  10.3969/j.issn.1003-3106.2009.01.013

    WANG Q,WU B. Analysis on precision of pseudo noise code ranging in space TT&C systems[J]. Radio Engineering,2009,39(1):39-44. doi:  10.3969/j.issn.1003-3106.2009.01.013
    [5] 胡松杰,唐歌实. 北京中心深空探测器精密定轨预分析软件系统[J]. 飞行器测控学报,2010,29(5):69-74.

    HU S J,TANG G S. BACC orbit determination and analysis software for deep-space explorers[J]. Journal of Spacecraft TT&C\; Technology,2010,29(5):69-74.
    [6] RAPPAPORT T S. 无线通信原理与应用[M]. 蔡涛, 李旭, 杜振民, 译. 北京: 电子工业出版社, 1999.
    [7] BERNER J B, LAYLAND J M, KINMAN P W, et al. Regenerative pseudo-noise ranging for deep-space applications: TMO Progress Report 42-137[R]. [S. l.]: Communications Ground Systems Section, 1999.
    [8] CCSDS. Pseudo-Noise(PN) ranging system informational report: CCSDS 414.0-G-2[R]. Washington, DC : CCSDS, 2014.
    [9] 董光亮, 李国民, 王新永, 等. 中国深空网: 系统设计与关键技术(中)S/X/Ka三频段深空测控通信系统[M]. 北京: 清华大学出版社, 2016.
    [10] 金小军. 再生伪码测距技术及实现研究[D]. 杭州: 浙江大学, 2007.
    [11] 董光亮, 耿虎军, 李国民, 等. 中国深空网: 系统设计与关键技术(下)深空干涉测量系统[M]. 北京: 清华大学出版社, 2016.
    [12] 董光亮, 耿虎军, 李国民, 等. 中国深空网: 系统设计与关键技术(上)S/X双频段深空测控通信系统[M]. 北京: 清华大学出版社, 2016.
    [13] 刘磊,曹建峰,胡松杰,等. 地月平动点中继应用轨道维持[J]. 深空探测学报(中英文),2015(4):318-324.

    LIU L,CAO J F,HU S J. Maintenance of relay orbit about the Earth-Moon collinear libration points[J]. Journal of Deep Space Exploration,2015(4):318-324.
    [14] 汤锡生, 陈贻迎, 朱民才. 载人飞船轨道确定和返回控制[M]. 北京: 国防工业出版社, 2002.
    [15] 李济生. 人造卫星精密轨道确定[M]. 北京: 解放军出版社, 1995.
    [16] 陈明,张宇,曹建峰,等. 嫦娥二号探测器轨道确定与测轨技术[J]. 科学通报,2012,57(9):689-696. doi:  10.1360/972011-818

    CHEN M,ZHANG Y,CAO J F,et al. Orbit determination and tracking technology of CE-2 satellite[J]. Chinese Science Bulletin,2012,57(9):689-696. doi:  10.1360/972011-818
    [17] 曹建峰,黄勇,胡小工,等. 月球重力场对“嫦娥一号”近月轨道的影响[J]. 宇航学报,2010,31(4):998-1004. doi:  10.3873/j.issn.1000-1328.2010.04.010

    CAO J F,HUANG Y,HU X G,et al. The effect of lunar gravity field on the low orbit of lunar satellite CE-1[J]. Journal of Astronautics,2010,31(4):998-1004. doi:  10.3873/j.issn.1000-1328.2010.04.010
    [18] 陈明,唐歌实,曹建峰,等. 嫦娥一号绕月探测探测器精密定轨实现[J]. 武汉大学学报. 信息科学版,2011,36(2):212-217.

    CHEN M,TANG G S,CAO J F,et al. Precision orbit determination of CE-1 lunar satellite[J]. Geomatics and Information Science of Wuhan University,2011,36(2):212-217.
    [19] 牛东文, 段建锋, 王美, 等. 嫦娥四号探测器中国深空网干涉测量数据定轨精度分析[C]//第一届中国空间科学大会论文集. 厦门: 中国空间科学学会, 2019.

    NIU D W, DUAN J F, WANG M. A study of the orbit determination ability for China's Deep Space Network VLBI data for Chang 'E-4 prober[J]. The 1st China Space Science Assembly. Xiamen: Chinese Society of Space Research, 2019.
    [20] 段建锋,张宇,曹建峰,等. 中国月球探测任务轨道确定技术及发展综述[J]. 深空探测学报(中英文),2019,6(3):203-209.

    DUAN J F, ZHANG Y, CAO J F, et al. A summary of orbit determination techndogy for Chinese lunar exploration project[J]. Journal of Deep Space Exploration,2019,6(3):203-209.
  • [1] 王赤, 李磊, 张爱兵, 张珅毅, 侯东辉, 徐子贡, 谢良海, 王慧姿, 罗朋威, 郭静楠, 史全岐, 张小平.  月表太阳风和粒子辐射环境—“嫦娥四号”观测新结果 . 深空探测学报(中英文), 2022, 9(3): 239-249. doi: 10.15982/j.issn.2096-9287.2022.20220020
    [2] 刘思语, 黄勇, 李培佳, 杨鹏.  稀疏观测模式的“嫦娥四号”中继星轨道确定 . 深空探测学报(中英文), 2022, 9(1): 14-20. doi: 10.15982/j.issn.2096-9287.2022.20210071
    [3] 李康康, 谷德峰, 安子聪, 宋佳凝, 朱炬波.  空间引力波探测器轨道确定及仿真分析 . 深空探测学报(中英文), 2021, 8(2): 182-189. doi: 10.15982/j.issn.2096-9287.2021.20200081
    [4] 李海涛, 陈少伍, 李赞, 樊敏, 程承.  “嫦娥五号”任务再入返回段测控布站区域确定方法 . 深空探测学报(中英文), 2021, 8(3): 284-289. doi: 10.15982/j.issn.2096-9287.2021.20200053
    [5] 王明远, 王美, 平劲松, 韩松涛.  月球空间环境研究进展 . 深空探测学报(中英文), 2021, 8(5): 486-494. doi: 10.15982/j.issn.2096-9287.2021.20200013
    [6] 牛东文, 段建锋, 王美, 陈略, 路伟涛, 任天鹏.  中国深空网VLBI定轨精度及分时采集模式分析 . 深空探测学报(中英文), 2021, 8(5): 495-502. doi: 10.15982/j.issn.2096-9287.2021.20191120001
    [7] 高宇辉, 师明, 蔡敦波, 张弓.  一种通用型遥操作任务智能规划方法研究 . 深空探测学报(中英文), 2021, 8(2): 140-146. doi: 10.15982/j.issn.2096-9287.2021.20200071
    [8] 熊亮, 张磊.  月球南极中继通信轨道仿真分析 . 深空探测学报(中英文), 2020, 7(3): 264-270. doi: 10.15982/j.issn.2095-7777.2020.20191108006
    [9] 陈略, 谢剑锋, 韩松涛, 曹建峰, 平劲松.  “嫦娥4号”中继星开环测速方案设计与试验验证 . 深空探测学报(中英文), 2019, 6(3): 236-240. doi: 10.15982/j.issn.2095-7777.2019.03.006
    [10] 马传令, 刘勇, 梁伟光, 张尧.  “嫦娥4号”中继星应急轨道控制策略设计与分析 . 深空探测学报(中英文), 2019, 6(3): 269-276. doi: 10.15982/j.issn.2095-7777.2019.03.011
    [11] 段建锋, 张宇, 曹建峰, 陈略, 陈明, 谢剑锋.  中国月球探测任务轨道确定技术及发展综述 . 深空探测学报(中英文), 2019, 6(3): 203-209. doi: 10.15982/j.issn.2095-7777.2019.03.001
    [12] 欧阳琦, 牛东文.  “龙江2号”月球轨道微卫星定轨分析 . 深空探测学报(中英文), 2019, 6(3): 254-260. doi: 10.15982/j.issn.2095-7777.2019.03.009
    [13] 曹建峰, 陈略, 董卫华, 段建锋, 韩松涛, 张宇.  双站跟踪模式下“嫦娥4号”中继星定轨仿真分析 . 深空探测学报(中英文), 2019, 6(3): 241-246. doi: 10.15982/j.issn.2095-7777.2019.03.007
    [14] 段建锋, 李勰, 李翠兰, 王兆魁.  “嫦娥4号”中继星使命轨道段定轨计算与分析 . 深空探测学报(中英文), 2019, 6(3): 247-253. doi: 10.15982/j.issn.2095-7777.2019.03.008
    [15] 陈国辉, 华岳, 王波, 郑士昆, 肖勇, 马小飞, 张爱兵, 陈岚.  “嫦娥4号”中继星伞状可展开天线设计与验证 . 深空探测学报(中英文), 2018, 5(6): 524-530. doi: 10.15982/j.issn.2095-7777.2018.06.004
    [16] 徐进, 张爱兵, 贺一峰, 杨光文.  “嫦娥4号”中继星中继通信方案研究与分析 . 深空探测学报(中英文), 2018, 5(6): 506-514. doi: 10.15982/j.issn.2095-7777.2018.06.002
    [17] 张立华, 熊亮, 王鹏, 孙骥, 周文艳, 高珊, 刘适, 王晓磊, 关轶峰, 张爱兵, 徐进, 陈国辉.  “嫦娥4号”中继星任务分析与系统设计 . 深空探测学报(中英文), 2018, 5(6): 515-523. doi: 10.15982/j.issn.2095-7777.2018.06.003
    [18] 段建锋, 刘勇, 李勰, 陈明, 王兆魁.  “嫦娥4号”中继星任务轨道确定问题初探 . 深空探测学报(中英文), 2018, 5(6): 531-538. doi: 10.15982/j.issn.2095-7777.2018.06.005
    [19] 高珊, 周文艳, 梁伟光, 刘德成, 唐玉华, 杨维廉.  地月拉格朗日L2点中继星轨道分析与设计 . 深空探测学报(中英文), 2017, 4(2): 122-129. doi: 10.15982/j.issn.2095-7777.2017.02.004
    [20] 董光亮.  深空测控新技术研究进展 . 深空探测学报(中英文), 2014, 1(4): 243-249. doi: 10.15982/j.issn.2095-7777.2014.04.001
  • 加载中
图(3) / 表 (4)
计量
  • 文章访问数:  136
  • HTML全文浏览量:  44
  • PDF下载量:  24
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-13
  • 修回日期:  2021-02-02
  • 网络出版日期:  2022-03-24
  • 刊出日期:  2022-04-13

“嫦娥四号”中继星再生伪码测距数据定轨精度分析

doi: 10.15982/j.issn.2096-9287.2022.20191213002
    基金项目:  国家自然科学基金资助项目(11973015)
    作者简介:

    牛东文(1986– ),男,工程师,主要研究方向:航天测控、航天器测定轨技术。通讯地址:北京市海淀区北清路26号院(100094)电话:(010)66360292E-mail:ndw05044050@163.com

    张宇(1979– ),男,高级工程师,主要研究方向:航天器测定轨技术。通讯地址:北京5130信箱(100094)电话:(010)66361456E-mail:zackyzy@163.com

  • ● China’s first regenerative pseudo-random code ranging experiment beyond the distance between the Earth and the Moon is introduced. ● Based on the actual track measurement data, the orbit determination accuracy of the regenerative pseudo-random code ranging is analyzed and evaluated. ● Based on the actual orbit measurement data, through the precise orbit determination analysis, it is quantitatively verified and further illustrated that the measurement accuracy and orbit determination accuracy of the regenerative pseudo-random code ranging are better than those of the sidetones ranging.
  • 中图分类号: V412.41

摘要: “嫦娥四号”中继星目前正稳定运行在地月L2点使命轨道,期间进行了再生伪码测距试验,试验期间交叉进行了再生伪码测距和侧音测距。利用“嫦娥四号”中继星在地月L2点使命轨道的再生伪码测距数据和侧音测距数据,进行了定轨精度评估。结果显示再生伪码测距数据的均方根误差(Root Mean Square Error,RMS)优于侧音测距数据RMS一个量级,同时再生伪码测距数据定轨和预报精度优于侧音测距数据一倍。这对再生伪码测距在中国后续深空探测任务中的应用具有参考意义。

注释:
1)  ● China’s first regenerative pseudo-random code ranging experiment beyond the distance between the Earth and the Moon is introduced. ● Based on the actual track measurement data, the orbit determination accuracy of the regenerative pseudo-random code ranging is analyzed and evaluated. ● Based on the actual orbit measurement data, through the precise orbit determination analysis, it is quantitatively verified and further illustrated that the measurement accuracy and orbit determination accuracy of the regenerative pseudo-random code ranging are better than those of the sidetones ranging.

English Abstract

牛东文, 段建锋, 欧阳琦, 张宇, 陈略, 王美. “嫦娥四号”中继星再生伪码测距数据定轨精度分析[J]. 深空探测学报(中英文), 2022, 9(1): 21-28. doi: 10.15982/j.issn.2096-9287.2022.20191213002
引用本文: 牛东文, 段建锋, 欧阳琦, 张宇, 陈略, 王美. “嫦娥四号”中继星再生伪码测距数据定轨精度分析[J]. 深空探测学报(中英文), 2022, 9(1): 21-28. doi: 10.15982/j.issn.2096-9287.2022.20191213002
NIU Dongwen, DUAN Jianfeng, OUYANG Qi, ZHANG Yu, CHEN Lue, WANG Mei. Regenerative Pseudo-Random Code Ranging Orbit Determination Accuracy Analysis for Chang’E-4 Relay Satellite[J]. Journal of Deep Space Exploration, 2022, 9(1): 21-28. doi: 10.15982/j.issn.2096-9287.2022.20191213002
Citation: NIU Dongwen, DUAN Jianfeng, OUYANG Qi, ZHANG Yu, CHEN Lue, WANG Mei. Regenerative Pseudo-Random Code Ranging Orbit Determination Accuracy Analysis for Chang’E-4 Relay Satellite[J]. Journal of Deep Space Exploration, 2022, 9(1): 21-28. doi: 10.15982/j.issn.2096-9287.2022.20191213002
    • “嫦娥四号”月球探测器的中继卫星“鹊桥”,是人类历史上第一颗环绕地月L2点的专用中继通信卫星,于2018年5月21日利用“长征4号丙”火箭在西昌卫星发射中心成功发射,当前已成功运行在地月L2点的使命轨道。使命轨道为南簇Halo轨道,平均周期约为14 d,Z方向的轨道振幅约为1.3万km[1]。在地面应用系统和测控系统的支持下,“鹊桥”中继星在使命轨道一方面要定期维持其稳定运行在Halo轨道,另一方面还要为着陆器的月面工作以及“玉兔二号”巡视器的遥操作提供中继服务。根据任务规划,在“鹊桥”中继星进入地月L2点使命轨道段并正常运行后,将在着陆器、巡视器休眠期间择机开展再生伪码测距试验,验证地面测控设备与“鹊桥”中继星联合开展再生伪码测距的功能,评估再生伪码测距性能。

      在深空探测任务中,深空航天器的导航定位需要高精度测距。地面测控站接收到的信号具有时延巨大、功率微弱的显著特点,所以相比近地测距,深空测距的技术难度更大。测距方法通常有两种:第一种是侧音测距;第二种是伪噪声(Pseudo Noise,PN)码序列测距,简称伪码测距。这两种方法的测距原理都是从地面测控站发出带有特殊标记的信号,航天器接收测距信号将其转发回地面测量站,然后比较发射的测距信号与接收的测距信号间的时间差或相位差,从中得到测距信息。侧音测距的优势是设备相对简单、测距精度高、捕获时间短,主要缺点是在测量中需要多次解模糊,解模糊能力差,同时硬件实现频率小于0.5 Hz的次侧音有困难,而最大无模糊测量距离又是由最小次侧音决定的[2]。因此,侧音测距多用于近地航天器,而伪随机码测距则多用于深空航天器。在伪随机码测距中容易采用长的码周期,可以增大无模糊距离,多次解距离模糊的复杂问题也能够避免,且可以增强抗干扰性和保密性,同时调制载波后可以和其它信号共同占用一个射频带宽。

      随着深空探测任务的增多,近年来伪码测距在国际深空探测领域的发展与应用十分迅速,已经成为深空测距领域的一个研究热点。而再生伪码测距已经在美国航空航天局(National Aeronautics and Space Administration,NASA)的喷气推进实验室(Jet Propulsion Laboratory,JPL)设计的新航天器转发模块(Spacecraft Transponding Module,STM)中得到应用,此转发模块已经应用于深空与近地空间任务中。采用再生伪码测距后,地面测控站接收到的测距信号功率谱密度增加了30 dB。NASA已经开发出无线电通讯系统中的再生伪随机码测距跟踪系统,并应用于探测冥王星的“新地平线号”(New Horizons)任务中[2]。国内学者也对伪码测距的相关理论技术进行了研究,文献[3]介绍了NASA深空网采用的伪码测距基本原理,分析了再生伪码测距系统具体序列的产生、相关和距离解算的方法与途径,并计算分析了计算精度和捕获时间等系统关键技术指标。文献[4]结合工程经验数据,给出了伪码测距系统的具体实现结构,深入研究了相位噪声、群时延、电平抖动、时钟抖动、跟踪环路、热噪声、频率源稳定性和介质折射等引入的随机误差和系统误差,确定了伪码测距系统的主要误差源,为进一步提高伪码测距精度指出了方向。

      本文在介绍“嫦娥四号”“鹊桥”中继星再生伪码测距试验的基础上,基于北京航天飞行控制中心的轨道计算与分析软件平台(BACC Orbit Determination and Analysis System,BODAS)[5],利用再生伪码测距试验获得的测量数据,通过与精密轨道的对比,对再生伪码测距数据的定轨精度进行评价分析。

    • 伪随机码测距分为透明测距(Transparent Ranging)和再生测距(Regenerative Ranging)。所谓透明转发测距,是指航天器接收到地面发送的信号后从中提取出测距信号,经过滤波、下行调制处理,然后与遥测调制部分相加送到调制器下行发射。自由空间中距离发射机$ r $处天线的接收功率由Friis公式给出[6]

      $$ {P_\text r} = \frac{{{P_\text t}{G_\text t}{G_\text r}{\lambda ^2}}}{{{{\left( {4 \text{π} } \right)}^2}{r^2}\eta }} $$ (1)

      其中:${P_\text t}$为发射功率;${P_\text r}$为接收功率;${G_\text t}$为发射天线增益;${G_\text r}$为接收天线增益;$r$为发射和接收之间的距离;$\eta $为与传播无关的系统损耗因子;$ \lambda $为波长。由式(1)可知,对于透明转发测距,地面站将测距信号发出后,由航天器转发后地面站再接收到的信号信噪比正比于$ 1/{r^4} $,其中$ r $是需要测量的航天器到测站的距离。

      在再生测距系统中,伪随机测距码经过调相到上行载波中,然后从地面测站传输至航天器上。再生测距中用到的伪码信号是由测距时钟与一些伪随机码经逻辑组合而成。测距信号经航天器接收后,星上转发器将其进行解调,这样测距信号被获取。航天器重新生成与地面发射的测距信号完全一致的伪码测距信号,然后将此伪码信号通过调相方式调制到下行载波中。地面站接收到下行信号后解调出测距信号,将此测距信号与上行发射的伪码信号进行互相关处理,得到伪码信号的时延值,该时延值就对应了目标与地面站间的距离信息。由于航天器再生了测距信号,即测距信号由航天器再生后发给地面站,由式(1)可知,采用再生测距的方式地面站接收到信号的信噪比与$ 1/{r^2} $成正比。再生伪码测距系统结构如图1所示。

      图  1  再生伪码测距系统结构图

      Figure 1.  Chart of regenerative pseudo-random code ranging

      再生测距利用专门的电路对测距信号进行捕获和跟踪,待整个测距信号锁定以后,跟透明转发测距一样,进行下行调制度处理,再与遥测调制部分相加。采用再生测距,一方面由于再生了信号,所以不存在因为测距通道的前端滤波器滤掉测距信号高次谐波而引起信号质量下降的问题,另一方面再生测距码跟踪环的环路带宽非常小(一般为Hz量级),因此噪声大为减小(与透明转发测距的MHz量级相比)。根据上行信号强度不同,再生伪码测距能将返回到测距设备的测距信号的信噪比提高达30 dB,这部分增益可以以3种方式提高测量性能:一是下行测距信号调制度和测距时间不变,测距精度可以得到提高;二是减小测距信号的调制度,而遥测信号的功率增加,进而减少遥测误码率,有利于对航天器的状态和控制情况进行判断;三是测距时间可以减少[7]

    • 测距精度是测距的核心问题。空间数据系统咨询委员会(Consultative Committee for Space Data Systems,CCSDS)建议选择两种伪随机码,称为权重平衡陶思沃斯码(Weighted-Voting Balanced Tausworthe Code),$ \nu = 4 $(T4B)与$ \nu = 2 $(T2B)[2],即T4B权重平衡陶思沃斯码和T2B权重平衡陶思沃斯码。

      权重平衡陶思沃斯码是由6个周期的二进制序列通过组合算法组合而成,将权重因子$ \nu = 4 $$ \nu = 2 $赋予测距时钟序列$ {C_1} $。再生测距使用的伪随机码的测距序列是由测距时钟与数个PN分序列的逻辑组合而成的。测距时钟序列是周期为2的周期性二进制($ \pm 1 $)序列。PN码是周期性的二进制($ \pm 1 $)序列,在一个长度L的周期上,其周期性自相关函数有峰值+L,其余所有L–1个峰值外的值等于–1。PN码序列长度分别为2、7、11、15、19、23。T4B陶思沃斯码与T2B陶思沃斯码在结构上是一致的,其组成的序列结构如图2所示。

      图  2  再生测距陶思沃斯码

      Figure 2.  Regenerative ranging Tausworthe code

      对于T4B陶思沃斯码的最终组合序列

      $$ C = {\text{sign}}\left( {4{C_1} + {C_2} - {C_3} - {C_4} + {C_5} - {C_6}} \right) $$ (2)

      对于T2B陶思沃斯码的最终组合序列

      $$ C = {\text{sign}}\left( {2{C_1} + {C_2} - {C_3} - {C_4} + {C_5} - {C_6}} \right) $$ (3)
    • 地面站发射机将图2所示的伪随机码序列组成一个总长度为1 009 470位的唯一序列调制到上行载波中,采用的是线性调制方式。测距信号的码速率应该与上行载波的频率相关,表达式如下,其中$ k{\text{ = }}6 $l={1,2,4,8,16,32,64},或者l=2,$ k{\text{ = }}\left\{ {8,9,{\text{10}}} \right\} $

      $$ {F_{{\text{chip}}}} = 2{F_{{\text{clock}}}} = l\frac{{{f_{{\text{S-band}}}}}}{{128 \times {2^k}}} $$ (4)
      $$ {F_{{\text{chip}}}}{\text{ = }}2{F_{{\text{clock}}}} = l\left( {\frac{{221}}{{749}}} \right)\frac{{{f_{{\text{X-band}}}}}}{{128 \times {2^k}}} $$ (5)
      $$ {F_{{\text{chip}}}}{\text{ = }}2{F_{{\text{clock}}}} = l\left( {\frac{{221}}{{3\;599}}} \right)\frac{{{f_{{\text{Ka-band}}}}}}{{128\times 2^k}} $$ (6)

      其中:$ {F_{{\text{chip}}}} $为码速率,单位Mchip/s;$ {F_{{\text{clock}}}} $为测距时钟,单位MHz;$ {f_{{\text{S - band}}}} $$ {f_{{\text{X - band}}}} $$ {f_{{\text{Ka - band}}}} $分别为S波段、X波段、Ka波段的上行载波频率,单位MHz。双程伪码测距模糊度为

      $$ AMB = L \times c/{F_{{\text{chip}}}} $$ (7)

      其中:$ L $为伪码序列总长度;c为光速。

      地面站为了更好地协同工作,最少需要满足两种码速率的值:一种码速率的值约为2 Mchip/s($ l{\text{ = }}8 $$ k{\text{ = }}6 $),另一种码速率的值约为1 Mchip/s($ l{\text{ = }}4 $$ k{\text{ = }}6 $)。

    • 采用再生伪码测距时,下行链路测距信号的功率为[8]

      $$ {P_{{\text{PN}}}} = {P_{\text{T}}} \times {\cos ^2}{m_{\text {TM}}} \times {\sin ^2}{m_{{\text{PN}}}} $$ (8)

      而采用侧音测距时,下行链路测距信号的功率为[9]

      $$ {P_{{\text{ST}}}} = {P_{\text{T}}} \times {\cos ^2}{m_{{\text{TM}}}} \times J_1^2({m_{{\text{ST}}}}) \times {\text e^{ - m_{\text {N}} ^2}} $$ (9)

      其中:$ {P_{{\text{PN}}}} $是伪码测距信号功率;$ {P_{\text{T}}} $是下行信号总功率;$ {m_{{\text{TM}}}} $是遥测数据调制指数;$ {m_{{\text{PN}}}} $是伪码测距信号调制指数;$ {m_{{\text{ST}}}} $是下行测距音的有效调制指数;$ {m_{\text{N}}} $是下行测距信道中热噪声的有效调制指数;$ {J_1}({m_{{\text{ST}}}}) $是变量$ {m_{{\text{ST}}}} $的一阶贝塞尔函数。

      当调制指数$ {m_{{\text{ST}}}} < 1 $时,$ {m_{{\text{ST}}}} $的一阶贝塞尔函数的拟合多项式为

      $$ {J_1}({m_{{\text{ST}}}}) = - 0.089\;1 m_{{\text{ST}}}^{\text{2}} + 0.533\;8 {m_{{\text{ST}}}} - 0.002\;4 $$ (10)

      需要说明的是,在侧音测距中,调制指数$ {m_{{\text{ST}}}} $$ {m_{\text{N}}} $都不等于下行测距调制指数的设计值$ {m_{{\text{RD}}}} $,这是因为测距信道中的自动增益控制电路加强了如下约束

      $$ m_{\text{ST} }^2 + m_{\text{N}}^{\text{2}} = m_{{\text{RD}}}^{\text{2}} $$ (11)

      即下行测距信道的总功率(测距信号功率加上信道带宽内的噪声功率)是一个恒定值,有效的下行调制指数为

      $$ {m_{\text {ST}}} = {m_{{\text{RD}}}} \sqrt {\frac{\varGamma }{{1 + \varGamma }}} $$ (12)
      $$ {m_{\text{N}}} = \frac{{{m_{{\text{RD}}}}}}{{\sqrt {1 + \varGamma } }} $$ (13)

      其中:$\varGamma$为应答机测距信道的信噪比。

      而在再生伪码测距中,由于上行链路中测距通道的热噪声不带到下行测距通道中,其下行链路设计的伪码调制指数$ {m_{{\text{PN}}}} $就是实际下行链路的测距调制指数,同时$ {m_{{\text{PN}}}} $通常是按档位设计为固定值。

    • 2019年6月,“鹊桥”中继星进行了再生伪码测距试验,试验内容包括数字应答机功能验证测试、伪码测距功能测试、伪码测距特性测试和伪码测距精密定轨测试等,获取了数字应答机的捕获时间、捕获概率和测距结果等数据,测试结果正常。在伪码测距精密定轨测试中,采用了T4B陶思沃斯码,其码速率值约为2 Mchip/s($ l{\text{ = }}8 $$ k{\text{ = }}6 $),对再生伪码测距的距离捕获时间、概率和测距性能进行了测试。参与试验的测控站主要为佳木斯深空站、喀什深空站、阿根廷深空站、纳米比亚站。测试期间“鹊桥”中继星没有进行轨道维持或动量轮卸载,这样在试验后通过精密定轨来评估再生伪码测距性能时,避免了对轨道维持或动量轮卸载产生的速度增量的求解估计,即为通过精密定轨评估再生伪码测距性能提供了良好条件。再生伪码测距时段如表1所示,在其它部分时段进行了侧音测距,以便于进行定轨精度对比分析。

    • 为了便于比较下行链路伪码测距信号功率和侧音测距信号功率,在1.3节的式(8)和式(9)中,取伪码测距和侧音测距的下行信号总功率$ {P_{\text{T}}} $及遥测数据调制指数$ {m_{{\text{TM}}}} $相等,则有

      $$ \left( {\frac{{{P_{{\text{PN}}}}}}{{{P_{{\text{ST}}}}}}} \right)\text {dB} = 10 \times {\lg}\left( {\frac{{{{\sin }^2}{m_{{\text{PN}}}}}}{{J_1^2({m_{{\text{ST}}}}) \times {\text e^{ - m_{\text{N}}^2}}}}} \right) $$ (14)

      根据本次“鹊桥”中继星伪码测距试验的情况,在式(13)中取$ {m_{{\text{PN}}}} $的设计值为固定值0.5 rad,取侧音测距设计调制指数范围为$ {m_{{\text{RD}}}} = \left( {0:1} \right) $,同时取应答机测距信道的信噪比$\varGamma$为20 dB,侧音测距有效调制指数$ {m_{{\text{ST}}}} $和测距信道中热噪声的有效调制指数$ {m_{\text{N}}} $分别利用式(11)和式(12)进行计算,则${P_{{\text{PN}}}}/{P_{{\text{ST}}}}\left( \text {dB} \right)$理论计算结果如图3所示。从图3可看出,在侧音测距设计调制指数较小时,下行链路的伪码测距信号的功率要高出侧音测距信号功率很多,而随着侧音测距设计调制指数的增加,下行链路的伪码测距信号功率将与侧音测距信号功率趋于一致。

      表 1  “鹊桥”中继星再生伪码测距试验时段(北京时间)

      Table 1.  Regenerative pseudo-random code ranging experiment schedule(BJT)

      日期测站观测时长(hh:mm)
      2019.06.14纳米比亚站2:49
      阿根廷深空站2:59
      佳木斯深空站3:03
      喀什深空站3:12
      2019.06.15喀什深空站1:57
      2019.06.16纳米比亚站2:59
      阿根廷深空站2:54
      2019.06.17佳木斯深空站2:46
      2019.06.18喀什深空站2:58
      纳米比亚站3:23
      2019.06.19纳米比亚站2:23
      阿根廷深空站2:59
      佳木斯深空站1:49
      2019.06.20佳木斯深空站1:10
      喀什深空站2:50

      图  3  下行链路再生伪码测距与侧音测距信号功率比

      Figure 3.  Power ratio of downlink regenerative pseudo-random code ranging to sidetone ranging signal

      需要说明的是,在此次“鹊桥”中继星的伪码测距试验中,伪码测距的信噪比相对于侧音测距的信噪比提高了大约3~6 dB,并没有达到国外文献所说的30 dB。通过分析“鹊桥”中继星的试验条件和国外文献资料,“鹊桥”中继星的伪码测距是试验状态,是功能测量,当时采用的上行发送功率均较高,这样对于侧音测距来讲通过器上自动增益控制电路获得的测距调制指数较高,从而试验结果的3~6 dB是与图3的理论计算结果相符的。而国外文献描述的30 dB是在上行功率微弱的状态下取得的结果,这样在侧音测距中通过器上自动增益控制电路获得的测距调制指数也会较低,从图3的理论计算结果也可以看出这一结论。而再生伪码测距的下行测距信号功率在器上能正确解调出上行测距伪码的情况下,是与上行信号功率无关的,而与器上具体设计有关。所以“鹊桥”中继星的伪码测距试验中没有达到国外文献描述的比侧音测距信号的信噪比高出30 dB的情况。

    • 再生测距的误差源包括频率基准源、传播介质、热噪声等。频率基准源是测距系统的核心部件,基准频率不稳是测距误差的重要来源,高精度应用时一般选择超稳振荡器(Ultra-Stable Oscillator,USO)[10]。传播介质是测距中另一个重要误差源,主要包括对流层误差、电离层误差、太阳等离子误差等。对于再生伪码测距的这些误差,需要在进行精密定轨之前进行数据处理,以消除或者修正误差,提高后续精密定轨精度。

      再生伪码测距数据处理包括:

      1)数据合理性检验

      在数据合理性检验中,主要是对测距(测距和)、测速等外测数据和测站高程、温度、湿度、压强等气象数据设定一定的范围,超过该范围的数据即认为是不合理数据,予以剔除。

      2)数据平滑

      数据平滑中采用3阶最小二乘法多项式拟合,以减小随机误差。

      3)传播介质修正

      对于传播介质误差修正,利用测站配备的环境参数采集设备测量相关的气象参数,对电离层和对流层等传播介质误差项进行修正,采集设备包括水汽微波辐射计、GNSS双频接收机、气象仪等[11]。同时完成应答机零值修正和软件解模糊,双程模糊度计算公式即式(7),而测站测控设备的零值和轴系误差等修正已由测站设备级计算机完成[12]

    • “鹊桥”中继星在地月L2点长期运行期间,因其不稳定的轨道特性,需要进行经常性的轨道维持控制使其轨道得以保持[13]。根据第2.1节的再生伪码测距试验时段及“鹊桥”中继星轨道维持情况,精密定轨使用2019-06-13T19:00:00—2019-06-23T11:00:00共约10 d的侧音测距、伪码测距、测速和VLBI数据进行计算,轨道改进历元时刻为2019-06-13T19:00:00,其中2019-06-14、2019-06-16、2019-06-18、2019-06-19每天有约1 h的VLBI测量数据。分析采用的动力学模型如表2所示[14-19],月球非球形采用8×8的Grail660模型,积分中心为地心。再生伪码测距试验期间,“鹊桥”中继星没有过阴影,姿态对日稳定,在精密定轨的时段中采用固定光压系数并作为一个待估参数,光压系数初值(Cr)和光压等效面积如表2所示。

      表 2  定轨分析策略

      Table 2.  Strategy for the orbit determination

      动力学模型参数
      地球质点J2000地心惯性系
      地球非球形JGM-38×8阶次
      月球质点第三体引力JPL-DE421历表
      月球非球形Grail6608×8阶次
      太阳质点第三体引力JPL-DE421历表
      其它大行星第三体引力JPL-DE421历表
      太阳光压球形模型
      固定面质比
      锥形阴影,Cr=1.7
      光压等效面积8.3 m2
    • 本节采用表2的动力学模型,按照4种策略进行再生伪码测距与侧音测距定轨及预报精度分析,使用2019-06-13T19:00:00—2019-06-21T12:00:00共约8 d的测量数据,轨道改进历元时刻为2019-06-13T19:00:00,预报至2019-06-23T11:00:00。4种策略分别为:侧音测距、测速和VLBI数据联合定轨,侧音测距和测速联合定轨,伪码测距、测速和VLBI联合定轨,伪码测距和测速联合定轨。在定轨中均求解测距和的系统差。4种策略测距数据质量评估如表3所示,与精密轨道比较的RTN方向误差统计结果如表4所示。从表3可以看出,再生伪码测距数据的RMS要优于侧音测距数据RMS一个量级。从表4与精密轨道的比较结果可以看出,再生伪码测距数据定轨和预报精度要优于侧音测距数据一倍。通常来讲,定轨的误差来源通常有测量数据误差、动力学模型误差、积分方法误差、观测几何等几个方面,是综合了测量、模型、几何及数学方法等各方面的误差后的结果。而本小节是在3.2节的定轨基本策略基础上,在动力学模型和积分方法等其它方面均一致的条件下,分析得到了分别使用伪码测距数据和侧音测距数据与其它测量数据联合定轨和预报的误差结果。

      表 3  不同策略定轨数据质量评估

      Table 3.  Orbit determination data quality evaluation for different strategies

      策略测控站数据类型RMS/m系统差/m
      1. 侧音测距+测速+VLBI纳米比亚站测距和1.78018.750
      佳木斯深空站测距0.612
      阿根廷深空站测距和1.56019.797
      喀什深空站测距和0.87628.367
      2. 侧音测距+测速纳米比亚站测距和1.830–22.718
      佳木斯深空站测距2.360
      阿根廷深空站测距和3.70025.121
      喀什深空站测距和0.90736.700
      3. 伪码测距+测速+VLBI纳米比亚站测距和0.47219.692
      佳木斯深空站测距0.470
      阿根廷深空站测距和0.37518.230
      喀什深空站测距和0.368310.10
      4. 伪码测距+测速纳米比亚站测距和0.46119.262
      佳木斯深空站测距0.452
      阿根廷深空站测距和0.394–17.50
      喀什深空站测距和0.386310.10

      表 4  不同策略定轨结果与精密轨道比较统计

      Table 4.  Statistics on comparisons of orbit determination results of different strategies with the precision orbit

      策略误差项定轨误差预报误差
      位置/m速度/(m·s–1位置/m速度/(m·s–1
      MAXRMSMAXRMSMAXRMSMAXRMS
      1. 侧音测距+测速+VLBIR–10.3663.223–0.001 480.000 65974.01322.6240.002 8270.000 475
      T455.685262.395–0.001 370.000 083–680.08471.394–0.00 1540.000 059
      N–627.289184.754–0.002 430.001 427–370.460114.9600.002 6260.000 200
      Total637.548320.9290.003 060.001 574684.546137.2030.004 1560.000 519
      2. 侧音测距+测速R–26.71112.4260.000 8750.000 39460.67221.8360.002 0050.000 323
      T–271.719152.388–0.001 0790.000 087–494.13864.174–0.001 490.000 124
      N590.413371.482–0.002 340.000 690–583.08616.424–0.001 120.000 590
      Total643.014401.7160.002 4960.000 800746.17469.7490.002 6640.000 684
      3. 伪码测距+测速+VLBIR–14.0804.3330.000 2820.000 07844.41610.4010.000 5210.000 067
      T–61.84923.284–0.000 1470.000 048–70.3672.3010.000 0130.000 006
      N–307.383112.071–0.001 2030.000 648–226.67153.6580.001 3010.000 143
      Total312.021114.5470.001 2120.000 655235.09454.7050.001 4020.000 158
      4. 伪码测距+测速R–13.7784.4260.000 3480.000 10051.47812.1510.000 6820.000 093
      T–82.64435.100–0.000 2050.000 047–105.7216.444–0.000 1030.000 011
      N–335.438121.933–0.0013 110.000 709–246.35558.6000.001 4230.000 156
      Total341.095126.9620.0013 210.000 717260.00160.1930.001 5810.000 182

      文献[3]给出了伪码测距系统各种误差源的随机误差计算公式,并给出正常工作条件下码速率为1 Mchip/s时,伪码测距系统的RMS可以控制在1 m[3]。“鹊桥”中继星在进行伪码测试时采用约2 Mchip/s的码速率,按照文献[3]的计算公式伪码测距数据的RMS可以控制在0.5 m,表3的数据质量评估结果与这一结论相符。需要说明的是,定轨数据的RMS是数据的残差均方根差,其表征了测量数据的噪声水平,本文测量数据的RMS是在迭代定轨过程中,根据定轨结果计算的观测数据理论值与对应时刻实际观测值进行对比,得到的定轨数据的均方根差。

    • 在深空探测任务中,高精度测距为深空航天器的导航定位提供了测量元素[2]。利用 “鹊桥”中继星进行的再生伪码测距,是中国首次在地月距离以外进行的再生伪码测距试验[20],可以充分体现再生伪码测距在深空任务导航中的优势。本文主要对 “鹊桥”中继星再生伪码测距数据定轨精度进行了评估,得出再生伪码测距数据的RMS要优于侧音测距数据RMS一个量级,同时再生伪码测距数据定轨和预报精度要优于侧音测距数据一倍。从测量原理上易知,再生伪码测距的测量精度和定轨精度要优于侧音测距,但在中国的实际任务中再生伪码测距还没有真正的实际应用。本文利用 “鹊桥”中继星再生伪码测距试验的实际轨道测量数据,通过精密定轨分析从定量上验证和进一步说明了这一结论,对于再生伪码测距在中国后续深空探测任务中的实际应用可以提供参考。

参考文献 (20)

目录

    /

    返回文章
    返回