留言板

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

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

一种卫星钟差异常精确探测的两步法策略

吕宏春 卢晓春 武建锋

吕宏春, 卢晓春, 武建锋. 一种卫星钟差异常精确探测的两步法策略[J]. 武汉大学学报 ● 信息科学版, 2019, 44(2): 193-199. doi: 10.13203/j.whugis20170076
引用本文: 吕宏春, 卢晓春, 武建锋. 一种卫星钟差异常精确探测的两步法策略[J]. 武汉大学学报 ● 信息科学版, 2019, 44(2): 193-199. doi: 10.13203/j.whugis20170076
LÜ Hongchun, LU Xiaochun, WU Jianfeng. A Two-step Method for Precise Detection of Satellite Clock's Abnormal Data[J]. Geomatics and Information Science of Wuhan University, 2019, 44(2): 193-199. doi: 10.13203/j.whugis20170076
Citation: LÜ Hongchun, LU Xiaochun, WU Jianfeng. A Two-step Method for Precise Detection of Satellite Clock's Abnormal Data[J]. Geomatics and Information Science of Wuhan University, 2019, 44(2): 193-199. doi: 10.13203/j.whugis20170076

一种卫星钟差异常精确探测的两步法策略

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

中国科学院西部之光项目 2013YB06

详细信息
    作者简介:

    吕宏春, 博士, 主要从事卫星导航理论与方法、时间同步技术研究。lvhongchun123@163.com

  • 中图分类号: P228.4

A Two-step Method for Precise Detection of Satellite Clock's Abnormal Data

Funds: 

West Light Foundation of the Chinese Academy of Sciences 2013YB06

More Information
    Author Bio:

    LÜ Hongchun, PhD, specializes in the theories and methods of satellite navigation and time synchronization. E-mail: lvhongchun123@163.com

  • 摘要: 针对卫星钟差数据中的粗差和钟跳两类异常进行精确探测,提出使用历元单差算法和拟准检定法相结合的两步法策略,第1步,依据钟跳和粗差在历元单差序列中不同表现特点,采用历元单差法探测疑似钟跳点并剔除部分粗差;第2步,对按疑似钟跳分段的数据使用拟准检定法精确探测粗差。通过模拟数据算例阐述钟差异常探测过程并进行性能分析,通过实测数据进行算法有效性验证,结果表明,仅使用历元单差法能有效探测疑似钟跳点和部分粗差,其探测的疑似钟跳点包含了真实钟跳;使用两步法既能有效规避钟跳影响,又能精确探测小粗差,探测能力接近底噪的2倍中误差。
  • 图  1  原始卫星钟差数据

    Figure  1.  Primary Satellite Clock Bias Data

    图  2  采用方案4对实测数据的处理结果

    Figure  2.  Result of Method 4 of Practical Data

    图  3  钟跳点的局部特写

    Figure  3.  Features of Clock Jump Points

    表  1  单粗差

    Table  1.   Single Gross Error

    项目 理论值
    5 6 7 8 9 10 11
    增量 0 0 0.5 0 0 0 0
    观测值 5 6 7.5 8 9 10 11
    历元单差 - 1 1.5 0.5 1 1 1
    单差真误差 - 0 0.5 -0.5 0 0 0
    下载: 导出CSV

    表  2  一般连续粗差情况

    Table  2.   General Continuous Gross Errors

    项目 理论值
    5 6 7 8 9 10 11
    增量 0 0 0.5 -0.3 0 0 0
    观测值 5 6 7.5 7.7 9 10 11
    历元单差 - 1 1.5 0.2 1.3 1 1
    单差真误差 - 0 0.5 -0.8 0.3 0 0
    下载: 导出CSV

    表  3  相近连续粗差

    Table  3.   Continuous Gross Errors with Similar Value

    项目 理论值
    5 6 7 8 9 10 11
    增量 0 0 0.5 0.55 0.52 0 0
    观测值 5 6 7.5 8.55 9.52 10 11
    历元单差 - 1 1.5 1.05 0.97 0.48 1
    单差真误差 - 0 0.5 0.05 -0.03 -0.52 0
    下载: 导出CSV

    表  4  钟跳

    Table  4.   Clock Jump

    项目 理论值
    5 6 7 8 9 10 11
    增量 0 0 1 1 1 1 1
    观测值 5 6 8 9 10 11 12
    历元单差 - 1 2 1 1 1 1
    单差真误差 - 0 1 0 0 0 0
    下载: 导出CSV

    表  5  钟跳和粗差在相邻历元出现

    Table  5.   Clock Jump and Gross Errors Happened in Adjacent Epochs

    项目 理论值
    5 6 7 8 9 10 11
    增量 0 0 1 0.6 1 1 1
    观测值 5 6 8 8.6 10 11 12
    历元单差 - 1 2 0.6 1.4 1 1
    单差真误差 - 0 1 -0.4 0.4 0 0
    下载: 导出CSV

    表  6  模拟数据分析结果

    Table  6.   Analysis Result of Simulation Data

    粗差点位置 粗差模拟值/ns 方案1 方案2 方案3 方案4
    4 -52 000
    5 92 000
    9 2.5 ×
    13 2.34 ×
    16 -5 ×
    20 24 000
    21 -72 000
    49 0.77 × × ×
    51 -0.89 × ×
    60 10 × ×
    61 9 × ×
    62 9.5 × ×
    68 -68 000
    71 -40 000
    99 1.5 × ×
    100 1.6 × ×
    109 -0.66 × ×
    121 -156 000
    123 0.84 × × ×
    128 -88 000
    133 0.68 × ×
    134 0.54 × × × ×
    141 -184 000
    142 64 000
    145 160 000
    148 1.6 × ×
    150 -48 000
    154 -0.99 × ×
    157 -96 000
    181 148 000
    184 96 000
    185 96 040
      注:“√”表示粗差被探测;“×”表示粗差未被探测
    下载: 导出CSV
  • [1] 黄飞江, 卢晓春, 刘光灿, 等.星地动态双向时间同步与测距算法[J].宇航学报, 2014, 35(9):1050-1057 doi:  10.3873/j.issn.1000-1328.2014.09.010

    Huang Feijiang, Lu Xiaochun, Liu Guangcan, et al. An Algorithm for Dynamic Satellite-Ground Two-way Time Synchronization and Ranging[J]. Journal of Astronautics, 2014, 35(9):1050-1057 doi:  10.3873/j.issn.1000-1328.2014.09.010
    [2] 李瑞锋, 王元明.星地双向时间同步技术在COMPASS中的应用[J].武汉大学学报·信息科学版, 2013, 38(7):842-844 http://ch.whu.edu.cn/CN/abstract/abstract2685.shtml

    Li Ruifeng, Wang Yuanming. Application of Two-way Satellite-Ground Time Synchronization in COMPASS[J].Geomatics and Information Science of Wuhan University, 2013, 38(7):842-844 http://ch.whu.edu.cn/CN/abstract/abstract2685.shtml
    [3] 唐升, 刘娅, 李孝辉.星载原子钟自主完好性监测方法研究[J].宇航学报, 2013, 34(1):39-45 doi:  10.3873/j.issn.1000-1328.2013.01.006

    Tang Sheng, Liu Ya, Li Xiaohui. A Study on Onboard Satellite Atomic Clock Autonomous Integrity Monitoring[J]. Journal of Astronautics, 2013, 34(1):39-45 doi:  10.3873/j.issn.1000-1328.2013.01.006
    [4] 杨玲, 沈云中, 楼立志.基于中位参数初值的等价权抗差估计方法[J].测绘学报, 2011, 40(1):28-32 http://d.old.wanfangdata.com.cn/Periodical/chxb201101006

    Yang Ling, Shen Yunzhong, Lou Lizhi. Equivalent Weight Robust Estimation Method Based on Median Parameter Estimates[J].Acta Geodaetica et Cartographica Sinica, 2011, 40(1):28-32 http://d.old.wanfangdata.com.cn/Periodical/chxb201101006
    [5] 王潜心, 徐天河, 许国昌.粗差检测与抗差估计相结合的方法在动态相对定位中的应用[J].武汉大学学报·信息科学版, 2011, 36(4):476-480 http://ch.whu.edu.cn/CN/abstract/abstract508.shtml

    Wang Qianxin, Xu Tianhe, Xu Guochang. Application of Combining Method of Outlier Detection and Robust Estimation to GPS Kinematic Relative Positioning[J].Geomatics and Information Science of Wuhan University, 2011, 36(4):476-480 http://ch.whu.edu.cn/CN/abstract/abstract508.shtml
    [6] 吴静.利用中位数的GPS卫星钟跳探测方法[J].测绘科学, 2015, 40(6):36-41 http://d.old.wanfangdata.com.cn/Periodical/chkx201506008

    Wu Jing. Study on Detection of GPS Clock Jump Using Median Absolute Deviation[J].Science of Surveying and Mapping, 2015, 40(6):36-41 http://d.old.wanfangdata.com.cn/Periodical/chkx201506008
    [7] 田婕, 黄观文, 王利, 等.一种北斗二代卫星钟差的数据质量控制方法[J].大地测量与地球动力学, 2017, 37(2):215-220 http://d.old.wanfangdata.com.cn/Periodical/dkxbydz201702022

    Tian Jie, Huang Guanwen, Wang Li, et al. A Data Quality Control Method for the BeiDou-2 Satellite Clock Bias[J].Journal of Geodesy and Geodyna-mics, 2017, 37(2):215-220 http://d.old.wanfangdata.com.cn/Periodical/dkxbydz201702022
    [8] 方杨, 王广兴, 杜玉军, 等.三种粗差检测方法的比较及分析[J].测绘通报, 2009(9):3427-3433 http://d.old.wanfangdata.com.cn/Periodical/chtb200909002

    Fang Yang, Wang Guangxing, Du Yujun, et al. Comparison and Analysis of Three Methods of Gross Error Detection[J].Bulletin of Surveying and Mapping, 2009(9):3427-3433 http://d.old.wanfangdata.com.cn/Periodical/chtb200909002
    [9] 费业泰.误差理论与数据处理.[M].6版.北京:机械工业出版社, 2010

    Fei Yetai. Error Theory and Data Processing[M]. 6th ed. Beijing:China Machine Press, 2010
    [10] 欧吉坤.一种三步抗差方案的设计[J].测绘学报, 1996, 25(3):173-179 doi:  10.3321/j.issn:1001-1595.1996.03.003

    Ou Jikun. Design of a New Scheme of Robust Estimation by Three Steps[J]. Acta Geodaetica et Cartographica Sinica, 1996, 25(3):173-179 doi:  10.3321/j.issn:1001-1595.1996.03.003
    [11] 胡志刚.北斗卫星导航系统性能评估理论与试验验证[D].武汉: 武汉大学, 2013

    Hu Zhigang. BeiDou Navigation Satellite System Performance Assessment Theory and Experimental Verification[D]. Wuhan: Wuhan University, 2013
    [12] 欧吉坤.一种检测粗差的新方法——拟准检定法[J].科学通报, 1999, 44(16):1777-1781 doi:  10.3321/j.issn:0023-074X.1999.16.021

    Ou Jikun. A New Methool of Gross Error Detection-Quasi-accurate Detection Method[J]. Chinese Sciencc Bulletin, 1999, 44(16):1777-1781 doi:  10.3321/j.issn:0023-074X.1999.16.021
    [13] 欧吉坤.再论拟准检定法的原理、实施和应用[J].测绘工程, 2002, 11(4):3-6 doi:  10.3969/j.issn.1006-7949.2002.04.001

    Ou Jikun. Further on the Principle Implementation and Application of Quasi-accurate Detection Method[J].Engineering of Surveying and Mapping, 2002, 11(4):3-6 doi:  10.3969/j.issn.1006-7949.2002.04.001
    [14] 欧吉坤.拟准观测的选取和真误差估值的"分群"现象[J].测绘学报, 2000, 29(Sup):7-12 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=chxb2000z1011

    Ou Jikun. Selection of Quasi-accurate Observations and Hive off Phenomena About the Estimators of Real Errors[J].Acta Geodaetica et Cartographica Sinica, 2000, 29(Sup):7-12 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=chxb2000z1011
    [15] 柴艳菊.拟准检定法的理论、应用及程序设计[D].武汉: 中国科学院研究生院(测量与地球物理研究所), 2002

    Chai Yanju. Theory, Application and Program Design of Quasi-Accurate Detection[D]. Wuhan: Institute of Geodesy and Geophysics, University of Chinese Academy of Sciences, 2002
    [16] Kaplan, E D, Christopher J. Hegarty, Understanding GPS Principles and Applications0[M]. 2nd ed.Boston:Artech House Inc, 2005
  • [1] 赵祥鸿, 暴景阳, 欧阳永忠, 黄贤源, 黄辰虎, 陆秀平.  利用BP神经网络剔除多波束测深数据粗差 . 武汉大学学报 ● 信息科学版, 2019, 44(4): 518-524. doi: 10.13203/j.whugis20160336
    [2] 鲁铁定, 杨元喜, 周世健.  均值漂移模式几种粗差探测法的MDB比较 . 武汉大学学报 ● 信息科学版, 2019, 44(2): 185-192, 199. doi: 10.13203/j.whugis20140330
    [3] 邱春平, 秦志远, 熊新, 张红敏, 靳国旺, 李贺.  机载SAR影像定向中像点粗差的拟准检定方法 . 武汉大学学报 ● 信息科学版, 2019, 44(2): 283-288, 295. doi: 10.13203/j.whugis20150338
    [4] 舒红, 史文中.  浅谈测量平差到空间数据分析的可靠性理论延伸 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 1979-1985, 1993. doi: 10.13203/j.whugis20180339
    [5] 单杰.  光束法平差简史与概要 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 1797-1810. doi: 10.13203/j.whugis20180331
    [6] 赵俊, 郭飞霄.  利用L1范数和中位数选取拟准观测值 . 武汉大学学报 ● 信息科学版, 2018, 43(8): 1242-1248. doi: 10.13203/j.whugis20160156
    [7] 潘雄, 吕玉婷, 汪耀, 罗静, 徐景田.  基于半参数平差模型的粗差定位与定值研究 . 武汉大学学报 ● 信息科学版, 2016, 41(11): 1421-1426. doi: 10.13203/j.whugis20140493
    [8] 张凯, 赵建虎, 张红梅.  一种基于M估计的水下地形抗差匹配算法 . 武汉大学学报 ● 信息科学版, 2015, 40(4): 558-562. doi: 10.13203/j.whugis20130011
    [9] 薛树强, 杨元喜.  抗差高斯-雅克比组合平差法 . 武汉大学学报 ● 信息科学版, 2015, 40(7): 932-937. doi: 10.13203/j.whugis20130247
    [10] 王振杰, 聂志喜, 欧吉坤.  一种基于TurboEdit改进的GPS双频观测值周跳探测方法 . 武汉大学学报 ● 信息科学版, 2014, 39(9): 1017-1021. doi: 10.13203/j.whugis20130021
    [11] 王海涛, 欧吉坤, 袁运斌, 郭建锋.  估计观测值粗差三种方法的等价性讨论 . 武汉大学学报 ● 信息科学版, 2013, 38(2): 162-166.
    [12] 蒋庆仙, 王成宾, 马小辉, 陈晓璧.  利用中位数进行光纤陀螺信号抗差估计 . 武汉大学学报 ● 信息科学版, 2011, 36(6): 656-659.
    [13] 赵长胜, 陶本藻.  有色噪声作用下的抗差卡尔曼滤波 . 武汉大学学报 ● 信息科学版, 2007, 32(10): 880-882.
    [14] 韩保民, 欧吉坤, 柴艳菊.  用拟准检定法探测和修复GPS数据中的粗差和周跳 . 武汉大学学报 ● 信息科学版, 2002, 27(3): 246-250.
    [15] 王金岭, 陈永奇.  粗差数学期望平移模型的理论研究 . 武汉大学学报 ● 信息科学版, 1995, 20(2): 146-150.
    [16] 李明峰, 於宗俦, 李延兴, 于正林.  动态平差概括模型与假设检验 . 武汉大学学报 ● 信息科学版, 1995, 20(3): 234-239.
    [17] 李明峰, 於宗俦, 于正林.  动态平差模型下形变与粗差的可区分性研究 . 武汉大学学报 ● 信息科学版, 1995, 20(1): 51-56.
    [18] 郑肇葆, 谭春健, 张润根.  按残差绝对值和最小原理的光束法区域网平差 . 武汉大学学报 ● 信息科学版, 1990, 15(2): 34-40.
    [19] 金新祥.  观测值残差和未知参数一次范数最小平差 . 武汉大学学报 ● 信息科学版, 1989, 14(3): 59-67.
    [20] 张方仁, 王勇.  主成份检验和粗差定位 . 武汉大学学报 ● 信息科学版, 1986, 11(4): 44-56.
  • 加载中
图(3) / 表(6)
计量
  • 文章访问数:  856
  • HTML全文浏览量:  100
  • PDF下载量:  299
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-09-29
  • 刊出日期:  2019-02-05

一种卫星钟差异常精确探测的两步法策略

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

    中国科学院西部之光项目 2013YB06

    作者简介:

    吕宏春, 博士, 主要从事卫星导航理论与方法、时间同步技术研究。lvhongchun123@163.com

  • 中图分类号: P228.4

摘要: 针对卫星钟差数据中的粗差和钟跳两类异常进行精确探测,提出使用历元单差算法和拟准检定法相结合的两步法策略,第1步,依据钟跳和粗差在历元单差序列中不同表现特点,采用历元单差法探测疑似钟跳点并剔除部分粗差;第2步,对按疑似钟跳分段的数据使用拟准检定法精确探测粗差。通过模拟数据算例阐述钟差异常探测过程并进行性能分析,通过实测数据进行算法有效性验证,结果表明,仅使用历元单差法能有效探测疑似钟跳点和部分粗差,其探测的疑似钟跳点包含了真实钟跳;使用两步法既能有效规避钟跳影响,又能精确探测小粗差,探测能力接近底噪的2倍中误差。

English Abstract

吕宏春, 卢晓春, 武建锋. 一种卫星钟差异常精确探测的两步法策略[J]. 武汉大学学报 ● 信息科学版, 2019, 44(2): 193-199. doi: 10.13203/j.whugis20170076
引用本文: 吕宏春, 卢晓春, 武建锋. 一种卫星钟差异常精确探测的两步法策略[J]. 武汉大学学报 ● 信息科学版, 2019, 44(2): 193-199. doi: 10.13203/j.whugis20170076
LÜ Hongchun, LU Xiaochun, WU Jianfeng. A Two-step Method for Precise Detection of Satellite Clock's Abnormal Data[J]. Geomatics and Information Science of Wuhan University, 2019, 44(2): 193-199. doi: 10.13203/j.whugis20170076
Citation: LÜ Hongchun, LU Xiaochun, WU Jianfeng. A Two-step Method for Precise Detection of Satellite Clock's Abnormal Data[J]. Geomatics and Information Science of Wuhan University, 2019, 44(2): 193-199. doi: 10.13203/j.whugis20170076
  • 星地双向时间同步[1-2]是目前精度最高的时间同步方法之一,可计算高精度的星地钟差,然而其直接计算得到的卫星钟差未经平滑,受到卫星钟、卫星收发设备、地面钟、地面收发设备、信号路径以及解算软件等众多因素影响,可能出现粗差或钟跳。针对卫星钟出现的异常可以通过监测的方式及时告警[3],减少卫星钟差异常值数量或对异常情况进行记录。但产生异常值的原因很多,仅在物理层面对卫星钟的监测难以杜绝卫星钟差异常的产生。在卫星数据处理领域,针对粗差未完全剔除的情况,常采用抗差估计以减弱异常值对数据处理结果的影响,以期得到一些不是很好但仍然可用的观测值[4-5],但抗差估计方法难以完全去除异常数据的影响。

    星地双向时间同步方法常常用于星地时间同步试验中,得到的钟差数据未经后处理,卫星钟差观测值可能出现缺失数据段,同时有效数据段内各类异常也可能同时出现。当钟跳和粗差两类异常同时段出现时,将对对方的探测产生不利影响,必须同时解决钟跳和粗差的探测问题;在卫星测试阶段可能出现粗差量级远大于钟差数据本身的异常,将使最小二乘残差算法的异常值处理方法失效,钟跳和粗差难以分离。目前同时解决钟跳和粗差探测的研究较少。文献[6]针对GPS卫星钟差进行钟跳判定和粗差剔除,采用对钟差按历元单差,以单差序列的10倍中位数参数为阈值,利用单个粗差使单差序列相应历元形成正负两个粗大值,而单个钟跳使单差序列相应历元形成一个粗大值的特点,同时探测钟跳和剔除粗差,但该方法未考虑实际情况下卫星钟差多个粗差连续出现以及钟跳邻近历元出现粗差的情况,这些情况下钟跳和粗差无法被有效区分,且该算法难以探测小粗差。文献[7]利用频率序列识别钟跳,采用阈值法进行异常数据段剔除以及结合绝对中位数法(median absolute deviation,MAD)和Baarda数据探测法共同进行粗差探测与剔除,对北斗系统多种钟差异常同时出现的情况进行有效控制,但该方法采用恒定的经验阈值(文中取1×10-5 s)判定钟跳,对小于经验阈值的异常钟跳无法判定,同时也将使包含这些钟跳的数据段内的粗差探测失效。

    本文对卫星钟差数据异常特点进行分析,顾及各类异常情况,据此探讨历元单差算法,使用中位数或中误差动态计算阈值,结合拟准检定法做进一步粗差探测,达到精准探测钟跳和粗差的目的,使用钟差模拟数据进行算法校验,利用卫星实测数据进行验证,分析了在复杂情况下剔除小粗差和探测钟跳的能力,另外还在算例中给出了探测成功后的进一步处理示例。

    • 粗差是相对于观测精度而言的,在经典的最小二乘平差中,通常认为大于3倍中误差的真误差为粗差[8]。常使用3倍中误差或残差的绝对中位数作为阈值。

      3倍中误差准则,即3σ准则,是最常用也是最简单的判别粗差的准则,其以测量次数充分大为前提,在实际测量中是一个近似的准则[9]。3倍中误差准则计算全样本的单位权中误差,当粗差较大时, 单位权中误差数值会相应较大,需多次应用3倍中误差准则才能有效消除粗差。

      残差的绝对中位数法在卫星数据处理中应用广泛,具有较强的抗差性,其判定准则为[10-11]

      $$ \left| {{V_j}-{\text{median}}\left( V \right)} \right| > n \cdot {\sigma _{{\text{MAD}}}} $$ (1)

      式中,V为残差序列;Vj为第j次测量的残差;σMAD为中位数误差;n为放大因子,一般取3~5,有时为降低误判概率,可取10;函数median·表示求取中位数。若式(1)成立,则认为该次观测含有粗差,应予剔除。残差的绝对中位数受粗差的影响很小,一般使用1次就可有效消除粗差。

    • 历元单差公式表示为:

      $$ {\mathit{\boldsymbol{D}}_j} = \left( {{\mathit{\boldsymbol{L}}_{j + 1}}-{\mathit{\boldsymbol{L}}_j}} \right)/\left( {{\mathit{\boldsymbol{T}}_{j + 1}}-{\mathit{\boldsymbol{T}}_j}} \right) $$ (2)

      式中,测量值序列记作L;相应的时间序列记作T;历元单差结果记作DLjTjDj分别表示各序列第j个元素。卫星钟差历元单差值Dj物理意义上可看成卫星钟速,在序列L没有异常值的情况下,D序列变化极小。

      依据§1.1方法设置阈值,若某个历元单差值超过阈值,则认为测量值相应位置发生了值的突变,可判定发生了粗差或者钟跳。

      典型的卫星钟差异常情况包括5种:单粗差情况、一般连续粗差情况、相近连续粗差值情况、钟跳情况以及钟跳相邻历元出现粗差的情况,分别如表 1-5所示。表中采用无量纲的数进行模拟分析,对历元单差真误差设置阈值0.1,真误差绝对值大于0.1的判定为异常。

      表 1  单粗差

      Table 1.  Single Gross Error

      项目 理论值
      5 6 7 8 9 10 11
      增量 0 0 0.5 0 0 0 0
      观测值 5 6 7.5 8 9 10 11
      历元单差 - 1 1.5 0.5 1 1 1
      单差真误差 - 0 0.5 -0.5 0 0 0

      表 2  一般连续粗差情况

      Table 2.  General Continuous Gross Errors

      项目 理论值
      5 6 7 8 9 10 11
      增量 0 0 0.5 -0.3 0 0 0
      观测值 5 6 7.5 7.7 9 10 11
      历元单差 - 1 1.5 0.2 1.3 1 1
      单差真误差 - 0 0.5 -0.8 0.3 0 0

      表 3  相近连续粗差

      Table 3.  Continuous Gross Errors with Similar Value

      项目 理论值
      5 6 7 8 9 10 11
      增量 0 0 0.5 0.55 0.52 0 0
      观测值 5 6 7.5 8.55 9.52 10 11
      历元单差 - 1 1.5 1.05 0.97 0.48 1
      单差真误差 - 0 0.5 0.05 -0.03 -0.52 0

      表 4  钟跳

      Table 4.  Clock Jump

      项目 理论值
      5 6 7 8 9 10 11
      增量 0 0 1 1 1 1 1
      观测值 5 6 8 9 10 11 12
      历元单差 - 1 2 1 1 1 1
      单差真误差 - 0 1 0 0 0 0

      表 5  钟跳和粗差在相邻历元出现

      Table 5.  Clock Jump and Gross Errors Happened in Adjacent Epochs

      项目 理论值
      5 6 7 8 9 10 11
      增量 0 0 1 0.6 1 1 1
      观测值 5 6 8 8.6 10 11 12
      历元单差 - 1 2 0.6 1.4 1 1
      单差真误差 - 0 1 -0.4 0.4 0 0

      表 1所示,若第N历元为单粗差情况,则第N-1和第N历元单差真误差判定为异常。

      表 2所示,若第N和第N+1历元出现连续粗差且粗差值相差较大,则第N-1至N+1这3个历元单差真误差判定为异常。

      表 3所示,第NN+2历元出现连续粗差且粗差值相近,则第N-1和N+2这两个历元单差真误差判定为异常。

      表 4所示,第N个历元开始出现钟跳,则第N-1个历元单差真误差判定为异常。

      表 5所示,第N个历元出现钟跳,并且第N+1个历元出现粗差,则第N-1至N+1这3个历元单差真误差判定为异常。

      探测粗差和钟跳应考虑这5种典型情况特点。利用历元单差法探测粗差和钟跳算法分为3步:

      1) 历元单差序列阈值计算。利用原始观测值计算历元单差序列,利用§1.1阈值设定方法,以历元单差序列为观测量求取阈值。

      2) 粗差剔除。若连续历元的单差值DjDj+1均超过阈值,则判定历元Tj+1对应的原始观测值含粗差,予以剔除。

      3) 疑似钟跳点判定。进行了第2)步粗差剔除之后,对剩余的数据计算历元单差序列,依第1)步所得的阈值Dth,若单差值Dj超过此阈值,则判定历元TjTj+1可能发生钟跳,记疑似钟跳历元为Tj+1

      上述历元单差法探测粗差和钟跳算法成功剔除了表 1表 2表 5形式的粗差,准确判定表 4钟跳位置,能够判定表 5钟跳,但也存在缺陷:表 3粗差异常被判为疑似钟跳,表 5钟跳初始的正常观测值被剔除。值得注意的是,在有大量正常观测样本的情况下,有个别正常值被误剔除并不影响系统观测特性;相近连续粗差的出现一般有其系统性原因,对其作为疑似钟跳标记位置是有益的。

      历元单差法能有效剔除粗差,其判定的疑似钟跳包含了真实钟跳情况。若在疑似钟跳观测值处对观测数据进行分段并做进一步处理,则可有效规避钟跳影响。

      一般依次考虑以下3个因素对数据进行分段:①由历元单差法判定的疑似钟跳点;②数据中断点;③根据前两个因素分段后,若某段数据量过大,计算机难以处理,可对该段进行数据分割。应注意并非所有的数据中断都需要进行数据分段,只有当中断达到一定的时长时才进行数据分段,这个时长通常由经验给出,本文后面将介绍的实测数据算例中,1 s采样间隔的数据当中断时长达100 s时即进行分段。

    • 对已经分段的钟差序列要进行更精准的粗差剔除,可分别对每段采用拟准检定法。

      1) 拟准检定法原理[12-13]

      拟准检定法认为实际观测值大部分正确,粗差是少数。则可以选定r个拟准观测,此处直接给出真误差估值如下:

      $$ \mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} = }}\left( \begin{gathered} {{\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}}_1} \hfill \\ {{\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}}_r} \hfill \\ \end{gathered} \right) =-{\left( {{\mathit{\boldsymbol{R}}^{\text{T}}}\mathit{\boldsymbol{R}} + {\mathit{\boldsymbol{P}}_Q}\mathit{\boldsymbol{G}}_Q^{\text{T}}{\mathit{\boldsymbol{G}}_Q}{\mathit{\boldsymbol{P}}_Q}} \right)^{-1}}{\mathit{\boldsymbol{R}}^{\text{T}}}\mathit{\boldsymbol{RL}} $$ (3)

      式中,${{\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}}_r}$是拟准观测的真误差估值;${{\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}}_1}$是非拟准观测的真误差估值;R表示平差因子阵;L是观测值;GQ=(0   ATr)取决于真误差估值对应的系数阵Ar;${\mathit{\boldsymbol{P}}_Q} = \left( {\begin{array}{*{20}{c}} 0&0 \\ 0&{{\mathit{\boldsymbol{P}}_r}} \end{array}} \right)$取决于真误差估值对应的权阵Pr

      拟准观测选取正确时,式(3)计算的真误差估值出现明显分群现象,即含粗差的真误差估值明显大于其他非正常的观测值的真误差估值。

      拟准观测的选取非常关键,在没有特别先验信息的情况下,可采用初选和复选两步方案。

      2) 拟准观测的初选、复选及分群准则[14-15]

      初选时,将观测值按一定规则分为4类:0类很可能含粗差;1类观测结构不好;3类含粗差可能性较小;其余为2类。实际操作中根据观测值计算指标对观测值分类,取第3类和标准化残差较小的2类观测值作为初选拟准观测值,参考文献[14-15]中对初选流程有详细阐述,此处不再赘述。

      通过初选的拟准观测,可以计算真误差估值$\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}$,进而计算相应指标进行复选观测。

      初次复选时, 计算指标为:

      $$ W_i^{\text{Ⅰ}} = {\hat \Delta _i}/{C_1}\left( {i = 1, 2 \cdots n} \right) $$ (4)

      式中,${C_1} = 1.483{\text{median}}\left( {\left| {\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}} \right|} \right)$。

      将$\left| {{{\hat \Delta }_i}} \right|$或|Wi|各自按大小排序,排序后对相邻数据求差,若$\left| {{{\hat \Delta }_i}} \right|$和|Wi|分群明显,则差值序列将有明显跳值, 且跳值点之下的$\left| {{{\hat \Delta }_i}} \right|$和|Wi|值都比较小,此处分界,分界线以下复选为拟准观测。

      初次复选后的分界线往往不够明显,需要迭代计算更为准确的复选指标,并且采用量化的分群准则,一般计算如下指标:

      $$ W_i^{{\text{Ⅱ}}} = {\hat \Delta _i}/{C_2}\left( {i = 1, 2 \cdots n} \right) $$ (5)

      式中,${C_2} = {\left[{\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}_r^{\text{T}}{{\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}}_r}/\left( {r-1} \right)} \right]^{1/2}}$,其中${{{\mathit{\boldsymbol{ \boldsymbol{\hat \varDelta} }}}_r}}$为上一次复选得到的拟准观测的真误差估值。可采用拟准观测的3倍中误差为阈值进行分群,判定粗差,即将指标Wi大于3倍中误差的观测判定为含有粗差。

      也可对真误差估值按大小排序,计算以下指标:

      $$ W_i^{{\text{Ⅲ}}} = \left( {{P_i}\left| {{{\hat \Delta }_i}} \right|-{P_{i + 1}}\left| {{{\hat \Delta }_{i + 1}}} \right|} \right)/{C_2} $$ (6)

      式中,P为观测权。可对按大小排序的真误差估值差分后,对差分序列使用2倍的中误差分群,即从指标Wi大于2倍中误差的位置进行分群。

      多次复选后,若分群仍不明显,说明没有粗差或选择的拟准观测不正确,若拟准观测的单位权中误差比较大,说明可能有粗差,更换拟准观测后返回拟准复选过程,否则认为没粗差。

    • 通过上述的处理过程,残余粗差的量值不大于拟准观测的3倍中误差,保证绝大多数的粗差被精确探测和剔除,此后可以对观测值的钟跳进行判定。判定方法有两种,第1种是对已经经过两步法剔除了粗差的数据再次计算历元单差序列,依据两步法中的第1步历元单差法计算的阈值对单差序列进行疑似钟跳的探测,探测结果即为真实钟跳点;第2种方法是分别利用疑似钟跳点之前和之后的数据拟合钟跳点钟差,并设定一个阈值,若拟合的两个值之差超过阈值,则判定此为钟跳。

    • 模拟数据能更好地模拟大粗差和各异常情况,更完整地体现卫星钟差异常探测过程,对卫星钟差实测数据的异常探测是算法的实践,以下分别进行算例分析。

    • 卫星钟差一般可表示成下列二次多项式[16]

      $$ {\Delta _{{\text{clk}}}} = {a_0} + {a_1}\left( {t-{t_{{\text{oc}}}}} \right) + {a_2}{\left( {t-{t_{{\text{oc}}}}} \right)^2} $$ (7)

      式中,toc表示参考时刻;t为计算时刻;a0a1a2为系数。

      给4个参数toca0a1a2分别依次取0、4×105 ns、0.86 ns/s、-6×10-6 ns/s2,每10 s得到一个钟差值,共计算200个数据。进行以下步骤, 产生模拟观测数据序列L

      1) 给所有历元加-0.5~0.5 ns的随机噪声,可将此随机噪声序列视为底噪,其中误差为0.283 ns;

      2) 给历元9、13、16、49、51、99、100、109、123、133、134、148、154分别加入2.5、2.34、-5、0.77、-0.89、1.5、1.6、-0.66、0.84、0.68、0.54、1.6、-0.99 ns的较小误差;

      3) 给历元4、5、20、21、68、71、121、128、141、142、145、150、157、181、184、185分别加入(-0.13, 0.23,0.06,-0.18,-0.17,-0.1,-0.39,-0.22,-0.46,0.16,0.4,-0.12,-0.24,0.37,0.24,0.240 1)×4×105 ns的较大量级误差;

      4) 给历元60、61、62分别加入10、9、9.5 ns的误差;

      5) 给第120及以后历元加-8×104 ns的偏差,模拟钟跳。

      使用4种方案处理模拟数据:

      方案1:使用历元单差法探测粗差和钟跳,使用绝对中位数法设置阈值,放大因子取10。

      方案2:与方案1方法相同,但放大因子取3。

      方案3:第1步,使用历元单差法探测粗差和钟跳,使用绝对中位数法设置阈值,放大因子取10;第2步,对各分段数据利用拟准检定法进行粗差探测,使用拟准观测的3倍中误差进行分群。

      方案4:第1步,使用历元单差法探测粗差和钟跳,使用绝对中位数法进行判定,放大因子取10;第2步,对各分段数据利用拟准检定法进行粗差探测,对真误差估值按大小排序并作差分,对差分序列使用拟准观测的2倍中误差分群。

      各方案处理结果如表 6所示。

      表 6  模拟数据分析结果

      Table 6.  Analysis Result of Simulation Data

      粗差点位置 粗差模拟值/ns 方案1 方案2 方案3 方案4
      4 -52 000
      5 92 000
      9 2.5 ×
      13 2.34 ×
      16 -5 ×
      20 24 000
      21 -72 000
      49 0.77 × × ×
      51 -0.89 × ×
      60 10 × ×
      61 9 × ×
      62 9.5 × ×
      68 -68 000
      71 -40 000
      99 1.5 × ×
      100 1.6 × ×
      109 -0.66 × ×
      121 -156 000
      123 0.84 × × ×
      128 -88 000
      133 0.68 × ×
      134 0.54 × × × ×
      141 -184 000
      142 64 000
      145 160 000
      148 1.6 × ×
      150 -48 000
      154 -0.99 × ×
      157 -96 000
      181 148 000
      184 96 000
      185 96 040
        注:“√”表示粗差被探测;“×”表示粗差未被探测

      方案1探测出第60、63和120历元的疑似钟跳点,以及使单差序列残差绝对值大于6.24 ns的钟差非连续粗差模拟值,本例中32个模拟粗差探测出16个。

      方案2探测出第60、63、101和120历元的疑似钟跳点,以及大部分使单差序列残差绝对值大于1.87 ns的钟差粗差模拟值,本例中32个模拟粗差探测出19个。

      方案3采用两步法,探测出疑似钟跳点后可探测小粗差,本例中32个模拟粗差探测出29个,探测出大于3倍底噪中误差的所有粗差模拟值和部分小于3倍底噪中误差但大于2倍底噪中误差的粗差模拟值。

      方案4采用两步法,探测出疑似钟跳点后可探测小粗差,本例中32个模拟粗差探测出31个,探测出粗差模拟值中大于2倍底噪中误差的所有粗差模拟值。

      粗差探测中难以避免地将部分正常观测值误判为粗差,本例误判统计:方案1为1个,方案2为1个,方案3为3个,方案4为4个。

      由算例可见,方案1和方案2中,历元单差法能有效探测真实钟跳点,并消除大部分较大粗差,同时选取较小的放大因子能提高粗差探测能力,但疑似钟跳点数量也随之增加;采用两步法策略有良好的粗差探测能力,方案4探测能力最好。

      粗差探测时不可避免地将某些正常观测值误判为粗差,本算例中方案4共判定粗差35个,包括31个模拟粗差和4个误判粗差。误判占总观测样本的2%,其减少了正常观测值个数,但在有大量正常观测样本的情况下,可认为少量的误判不影响正常观测特性。

      剔除粗差后可以对钟跳点进行判定,本例方案4判定出唯一钟跳点在第122历元处。这和本例模拟的钟跳点第120历元有差异,原因在于第121历元为真实粗差被探测并剔除,第120历元被误判为粗差而被剔除,剩余的不含粗差的观测数据中钟跳将出现在第122历元。

    • 本例针对某次星地双向时间同步钟差测试数据进行粗差剔除和钟跳探测,原始数据采样间隔1 s,数据时长约11 h,如图 1所示。

      图  1  原始卫星钟差数据

      Figure 1.  Primary Satellite Clock Bias Data

      粗差剔除和钟跳探测使用§2.1中的方案4,并对钟跳点进行判定,处理结果如图 2所示,图 2中钟跳点局部特写如图 3所示。

      图  2  采用方案4对实测数据的处理结果

      Figure 2.  Result of Method 4 of Practical Data

      图  3  钟跳点的局部特写

      Figure 3.  Features of Clock Jump Points

      图 2图 3可以看出无明显粗差存在,钟跳被良好探测和标记。卫星钟的正常调钟引起的钟跳往往持续时间较长且跳值较大,而从图 2图 3钟跳的分布特点和跳值可知这些钟跳并不是由正常的卫星调钟引起,对这些异常钟跳进行标记在实践中是有益的。

      至此对粗差的剔除和钟跳的判定已经完成,对探测出的钟跳值的后续处理应根据实际需求进行处理。

    • 本文依据历元单差法和拟准检定法的算法特点,分析了数据异常的典型形式,探讨使用历元单差法探测钟跳和粗差的方法和步骤,总结拟准检定法探测粗差的过程,提出将二者结合应用于卫星钟差数据异常探测的两步法策略,即先利用历元单差法探测疑似钟跳并剔除部分粗差,再对按疑似钟跳分段的数据使用拟准检定法。分别通过模拟数据和实测数据进行算例分析,得出以下结论:

      1) 历元单差算法简单,通过对历元单差序列使用残差的绝对中位数法检定,可有效探测粗差和疑似钟跳点,其探测的疑似钟跳点包含真实钟跳;

      2) 使用两步法策略,按历元单差法确定的疑似钟跳点对数据分段,每个段内使用拟准检定法进行粗差探测,结果表明,该方法可精确探测小粗差,探测能力接近底噪的2倍中误差。

      两步法策略步骤清晰,通过模拟数据和卫星钟差实测数据的分析,表明两步法策略能够精确剔除粗差和探测钟跳。

参考文献 (16)

目录

    /

    返回文章
    返回