BDS-3卫星钟差特性分析

毛亚, 王潜心, 胡超, 杨鸿毅, 张铭彬

毛亚, 王潜心, 胡超, 杨鸿毅, 张铭彬. BDS-3卫星钟差特性分析[J]. 武汉大学学报 ( 信息科学版), 2020, 45(1): 53-61. DOI: 10.13203/j.whugis20180224
引用本文: 毛亚, 王潜心, 胡超, 杨鸿毅, 张铭彬. BDS-3卫星钟差特性分析[J]. 武汉大学学报 ( 信息科学版), 2020, 45(1): 53-61. DOI: 10.13203/j.whugis20180224
MAO Ya, WANG Qianxin, HU Chao, YANG Hongyi, ZHANG Mingbin. Analysis of the Characterization for BDS-3 Satellite Clock Error[J]. Geomatics and Information Science of Wuhan University, 2020, 45(1): 53-61. DOI: 10.13203/j.whugis20180224
Citation: MAO Ya, WANG Qianxin, HU Chao, YANG Hongyi, ZHANG Mingbin. Analysis of the Characterization for BDS-3 Satellite Clock Error[J]. Geomatics and Information Science of Wuhan University, 2020, 45(1): 53-61. DOI: 10.13203/j.whugis20180224

BDS-3卫星钟差特性分析

基金项目: 

国家自然科学基金 41730109

江苏省自然科学基金 BK20181361

江苏省研究生科研与实践创新计划 SJCX18_0668

详细信息
    作者简介:

    毛亚, 硕士, 主要从事精密卫星钟差估计与预报工作。maoya0428@foxmail.com

    通讯作者:

    王潜心, 博士, 教授。wqx@cumt.edu.cn

  • 中图分类号: P228

Analysis of the Characterization for BDS-3 Satellite Clock Error

Funds: 

The National Natural Science Foundation of China 41730109

the Natural Science Foundation of Jiangsu Province BK20181361

Postgraduate Research and Practice Innovation Program of Jiangsu Province SJCX18_0668

More Information
    Author Bio:

    MAO Ya, master, specializes in precise satellite clock offset estimation and prediction. E-mail:maoya0428@foxmail.com

    Corresponding author:

    WANG Qianxin, PhD, professor. E-mail:wqx@cumt.edu.cn

  • 摘要: 北斗全球导航卫星系统(简称北斗三号系统,BDS-3)载有高精度的国产新型铷钟和被动型氢原子钟,相比于北斗二号系统(BDS-2),其性能得到了极大的提升。为了分析BDS-3卫星钟的性能,首先采用中国矿业大学北斗数据处理与分析中心定轨软件估计了BDS-3的卫星钟差,其重叠弧段标准差精度优于1 ns。然后针对常用卫星钟差异常值探测方法的缺点,改进了Baarda粗差探测法,使其能够有效地处理北斗卫星钟差数据中存在的粗差、钟跳等异常值。最后采用处理后的钟差数据,分析了北斗卫星钟差的周期特性和频率稳定性。结果表明,加入两个周期项后,拟合残差序列更加稳定,拟合精度得到了很大的提升;BDS-3卫星的频率稳定性和BDS-2相比,提升了20.85%~75.09%;北斗中圆地球轨道(medium earth orbit,MEO)卫星的频率稳定性要优于倾斜地球同步轨道(inclined geosynchronous orbit,IGSO)卫星;氢原子钟频率稳定性高于铷原子钟40.16%~65.60%,且氢原子钟的长期稳定性要明显优于铷原子钟。
    Abstract: BeiDou global navigation satellite system (BDS-3) satellites are equipped with high-precision domestic new rubidium clocks and passive hydrogen atomic clocks. Compared with the satellites of BeiDou-2 system (BDS-2), the performance of BDS-3 has great promotion. In order to analyze the performance of the BDS-3 satellite clock, firstly the satellite clock errors of BDS-3 were estimated using the orbit determination software from BeiDou Data Processing and Analysis Center, China University of Mining and Technology, and the standard deviation accuracy of the overlapping arc segments is better than 1 ns. Then, based on the shortcomings of the commonly used satellite clock error outlier detection method, the Baarda gross error detection method was improved so that it can effectively handle the outliers with gross errors and clock slips in the clock errors data. Finally, the processed clock errors data were used to analyze the period characteristics and frequency stability. Analysis results show that the fitting residuals are more stable and the fitting accuracy has been greatly improved after adding two periodic terms. Compared with BDS-2 satellites, the frequency stability of BDS-3 satellites has increased by 20.85%-75.09%. The frequency stability of medium earth orbit (MEO) satellites of BeiDou navigation satellite system is better than that of inclined geosynchronous orbit (IGSO) satellites. The frequency stability of the hydrogen atomic clock is higher than that of the europium atomic clock up to 40.16%-65.60%, and the long-term stability of the hydrogen atomic clock is significantly better than that of the europium atomic clock.
  • 全球导航卫星系统(global navigation satellite system, GNSS)旨在为全球用户提供全天候、全时段、高精度的导航、定位与授时服务[1]。北斗卫星作为北斗卫星导航系统(BeiDou navigation satellite system,BDS)空间端的主要组成部分,为整个系统提供重要空间基准,而星载原子钟身为北斗卫星的重要载荷之一,维持着导航系统星上时间基准,决定了北斗卫星在轨寿命,其性能不仅会影响北斗卫星自主导航能力,还会对其服务能力产生较大影响[2-3]

    GNSS星载原子钟在轨运行期间,会受到外部环境变化和设备老化等因素影响,导致其输出的钟差序列中存在粗差、钟跳等异常情况,而干净的卫星钟差数据是后续卫星钟性能研究与预报的前提[4-6],采用预处理后的钟差数据分析BDS卫星钟差特性可以取得较为可信的结果。常用的钟差数据质量控制方法主要有中位数法(median absolute deviation, MAD)和Baarda粗差探测法等。其中中位数法对粗差大小不敏感,简单有效,但很难探测出小的误差;Baarda法则是通过对验后残差进行处理,能够很好地识别隐藏在钟差数据中的小误差,但在进行多项式拟合钟差数据时,容易受到钟跳等异常值的影响无法准确获得拟合系数。针对常用钟差异常值探测方法的缺点,本文采用改进后的Baarda法对频率数据进行处理,能够准确探测出频率数据中的异常点,并采用多项式拟合值进行补充处理。

    2015年, 我国发射了第一颗北斗全球卫星导航系统(简称北斗三号系统,BDS-3)试验卫星。截止到2018年4月,BDS-3卫星共在轨5颗试验卫星(C31~C35)和8颗组网卫星(C19、C20、C21、C22、C27、C28、C29、C30)。文献[7]指出我国北斗星载原子钟天稳定度要比GPS和GLONASS系统的星载原子钟低一个数量级,因此提高我国原子钟的物理性能是提高BDS定位精度一个亟待解决的问题。文献[8]利用北斗精密钟差评估了北斗频率稳定度,发现北斗原子钟的千秒稳定性、万秒稳定性和日稳定性分别能达到(3~4)×10-13、(1~2)×10-13、(1~2)×10-14量级。文献[9-11]分析了BDS-3试验卫星的频率稳定度,结果表明BDS-3试验卫星频率稳定度比北斗二号系统(BDS-2)高一个量级。为分析BDS-3卫星钟差的性能,本文首先设计了3 d弧段的BDS-3精密卫星钟差解算实验;然后,根据北斗卫星钟差的噪声特性改进了Baarda法,为计算BDS-3卫星周期特性和频率稳定性指标提供了干净的钟差数据;最后,详细探讨了BDS-2/BDS-3之间、不同BDS卫星类型之间和不同星载原子钟类型之间的频率稳定性的差异。

    BDS-2卫星精密钟差可以从国际GNSS服务组织(International GNSS Service, IGS)和国际GNSS监测评估系统(International GNSS Monitoring & Assessment System, iGMAS)等途径获得,但是现有的途径均不提供BDS-3精密钟差产品。

    实验采用了iGMAS测站估计的3月1日—31日3 d弧段的BDS-2/BDS-3卫星钟差, 各卫星参数如表 1所示。

    表  1  北斗卫星参数
    Table  1.  Parameters of BDS satellites
    卫星 卫星号 卫星类型 原子钟
    BDS-2 C06 IGSO
    C07 IGSO
    C08 MEO
    C09 MEO
    C10 MEO
    C11 MEO
    C12 MEO
    C13 IGSO
    C14 MEO
    BDS-3 C19 MEO
    C20 MEO
    C27 MEO
    C28 MEO
    C31 IGSO
    C32 IGSO
    C33 MEO
    C34 MEO
    下载: 导出CSV 
    | 显示表格

    目前,已经建设好的iGMAS测站共有24个,其中能够接收到BDS-3信号的测站有17个,选择hmns和bjf1两个iGMAS测站分析跟踪BDS-3信号,结果如图 1所示。由图 1可以看出,iGMAS测站跟踪倾斜地球同步轨道(inclined geosynchronous orbit,IGSO)卫星的弧段长度要略长于跟踪中圆地球轨道(medium earth orbit, MEO)卫星的弧段长度,这也是估计得到的IGSO卫星钟差精度优于MEO卫星钟差精度的主要原因。图 2统计了iGMAS站接收BDS-3卫星信号的数据有效率,基本在50%左右。

    图  1  BDS-3卫星跟踪情况
    Figure  1.  BDS-3 Satellite Tracking Conditions
    图  2  数据有效率
    Figure  2.  Histogram of the Rate of Effective Data

    图 3表 2统计了BDS-3卫星钟差重叠弧段的均方根误差(root mean square, RMS)和标准差(standard deviation,STD)。从图 3表 2中可以发现,STD和RMS的精度相近,这说明连续两次定轨过程中不存在明显的系统性误差,而且大部分钟差精度优于1 ns。其中IGSO卫星的钟差精度优于MEO卫星,主要原因可能是iGMAS测站分布不均匀,MEO卫星的观测弧段较少。

    图  3  钟差解算精度
    Figure  3.  Accuracy of Clock Offset Estimation
    表  2  钟差精度统计/ns
    Table  2.  Accuracy Statistics of Clock Error/ns
    精度 C19 C20 C27 C28 C31 C32 C33 C34
    标准差 0.74 0.82 0.63 0.54 0.58 0.88 0.74 0.82
    均方差 0.61 0.80 0.89 0.80 0.74 0.63 0.97
    下载: 导出CSV 
    | 显示表格

    Baarda法能够有效地探测出钟差序列中的小粗差并将其剔除[12]。为了避免粗差、钟跳等异常值对二次多项式模型参数求解的影响,并且考虑到频率数据较钟差相位数据更能突出异常点,本文采用抗差最小二乘估计模型参数对频率数据进行建模,对探测出的异常点采用多项式拟合值进行填充,并将频率数据恢复成钟差相位数据。具体模型如下:

    $$ x(t) = {a_0} + {a_1}t + {a_2}{t^2} + {e_t} $$ (1)
    $$ {x^\prime }(t) = {a_1} + 2{a_2}t + {w_t} $$ (2)
    $$ y(t)=x^{\prime}(t)=\frac{x(t+\tau)-x(t)}{\tau} $$ (3)

    式中,x(t)为钟差相位数据;y(t)为频率数据;a0a1a2分别表示相位、钟速和钟漂;τ为采样间隔;etwt分别表示钟差和频率时间序列的拟合残差。

    为了确定频率数据的权,本文选用中国矿业大学分析中心解算的C19和C34两颗卫星在2018年年积日第60~90天共31 d的BDS-2/BDS-3精密钟差数据,来分析BDS星载原子钟的噪声水平,采样间隔为300 s。采用直方图和分位(quantile-quantile, Q-Q)图对多项式拟合残差进行分析,检验拟合残差是否符合标准正态分布,结果分别如图 4图 5所示。图 4显示了拟合残差的均值和标准差,从中可以看出拟合残差序列基本符合正态分布;图 5表明拟合残差序列在95%的置信区间内可以很好的拟合。从以上分析中可以得出:钟差残差序列基本符合正态分布。

    图  4  拟合残差直方图和正态分布图
    Figure  4.  Histograms and Normal Distribution of the Fitting Residuals
    图  5  拟合残差Q-Q图
    Figure  5.  Q-Q Distribution of the Fitting Residuals

    为了直观地展示本文改进的钟差数据异常值探测方法,图 6给出了本文进行数据预处理的详细流程图, 首先根据式(3),将读取的钟差数据转变为频率数据后,用线性模型对频率数据进行拟合得到钟速、钟漂和拟合残差,然后根据残差的标准差,归一化拟合残差得:

    $$ \tilde{v}_{i}=\frac{v_{i}}{\sigma_{v_{i}}} $$ (4)
    图  6  数据预处理流程图
    Figure  6.  Flowchart of the Data Pre-processing Method

    式中,$ {{\sigma _{{v_i}}}}$为残差vi对应的标准差。若归一化残差大于阈值k(一般取3),则剔除最大残差所对应的频率数据,并重新构建法方程,估计钟差参数得到残差, 直至没有残差大于阈值为止。

    根据改进的异常值探测方法,对C34号卫星钟差序列中存在的异常值进行探测,探测结果如图 7所示。根据图 7(a)中的频率数据可以看出,在钟差序列中存在些许较小的粗差,而图 7(b)中的频率数据没有发现显著的跳变,这表明本文所改进的Baarda粗差探测法能够有效探测出钟差序列中的小粗差。尤其是针对频漂特性比较明显的钟差序列,采用传统的Baarda粗差探测法和中位数法均不能得到很好的处理,而采用本文所改进的Baarda粗差探测法效果十分明显。另外,由于钟跳的出现会对钟差性能指标的计算产生较大的影响,有必要对钟差序列中存在的钟跳进行处理,将钟跳当作频率数据中的粗差进行处理,能够较好地处理钟差序列中存在的钟跳情况。

    图  7  C34钟差数据异常值探测结果
    Figure  7.  Results of Outlier Detection for C34

    卫星轨道和钟差高度相关,尤其是轨道径向的误差会被卫星钟差所吸收,因此在钟差数据中会有一定的周期特性[13]。本实验采用§2中处理好的卫星钟差数据,在去除二次趋势项之后,得到相应的钟差残差数据,再用快速傅里叶变换的方法获得卫星钟差的主要频率,以此来分析北斗卫星钟差的周期特性[14-15]。在采用傅里叶变换将时间域信号转换到频率域信号的过程中,原始的复杂信号可以分解为多个简单的正弦和余弦信号的叠加,在这些简单的正弦和余弦信号中, 分别有其各自的频率, 并且都可以通过功率谱体现出来,功率谱越大, 则表明对应频率的信号在原始信号中的作用越大, 否则越小[16-19]。因此,可以根据频谱分析结果中振幅的大小,选择钟差数据中的显著周期项。BDS卫星钟差主频率分析结果如图 8图 9所示,图中横轴给出的是一天中出现某个周期的次数。

    图  8  BDS-2卫星频谱分析
    Figure  8.  Frequency Analysis of BDS-2 Satellites
    图  9  BDS-3卫星频谱分析
    Figure  9.  Frequency Analysis of BDS-3 Satellites

    图 8图 9可以看出:(1)BDS-2和BDS-3的IGSO卫星以及BDS-2的C14卫星有明显的12 h和24 h周期项,另外除了C08,其他卫星还存在明显的8 h或6 h周期项,并且随着周期的缩短,振幅有逐渐降低的趋势,与文献[16]所得的结论一致。在BDS-2和BDS-3的MEO卫星中,C11、C28主要周期项不明显。(2)造成周期性不一致的原因可能是卫星钟差对外界环境变化比较敏感,导致估计的卫星钟差除了跟轨道周期耦合外,还会受到昼夜变化的影响。卫星钟差本身也存在一定的误差,导致拟合残差中也存在误差。另外能够接收到BDS-3信号的测站少、分布不均匀,导致钟差数据存在部分历元的缺失,也可能会造成周期性不一致。(3)不同的卫星主要周期项不同,钟差数据中的周期项将会对预报精度产生较大的影响。在对钟差数据进行建模预报时,必须考虑到周期特性的影响,根据文献[17]所得的研究结果,每颗卫星选择两个主要周期项加到式(1)中,构建附加周期项的钟差模型如下:

    $$ \begin{array}{*{20}{c}} {x(t) = {a_0} + {a_1}t + {a_2}{t^2} + \sum\limits_{i = 1}^P {\left( {{A_i}\sin \left( {\frac{{2\pi }}{T}t} \right) + } \right.} }\\ {\left. {{B_i}\cos \left( {\frac{{2\pi }}{T}t} \right)} \right) + {\mathit{\Delta} _i}} \end{array} $$ (5)

    式中,P为周期项的个数; AiBiT分别为正、余弦函数的振幅和周期;${{\mathit{\Delta} _i}} $表示附加周期项拟合残差。

    图 10为C14、C19、C34 3颗卫星不加周期项和加入两个主要周期项之后的拟合残差序列图,可以看出,加入两个周期项后,拟合残差序列更加稳定,拟合精度得到了很大的提升,这也证明了不同卫星钟差序列中存在几何周期项[17]

    图  10  C14、C19、C34卫星钟差拟合残差
    Figure  10.  Fitting Residuals of C14, C19, C34 Clock Offsets

    本实验仍采用§2处理后的数据,将钟差序列中的粗差、钟跳、天跳变以及频率跳变进行细密的处理后,计算卫星钟Allan方差和Hadamard方差(下文简称Had),结果如图 11图 12所示。表 3中给出了BDS-2的IGSO、MEO卫星和BDS-3卫星的千秒稳定性、万秒稳定性和日稳定性。可以看出,采用Allan方差和Had方差计算的频率稳定性指标基本一致。表 4中还统计了两种算法计算的频率稳定性指标之间的相似度,结果发现,相似度达到95%以上,且平滑时间在600~2 000 s左右时,Allan方差和Had方差序列均接近$ \tau^{-\frac{1}{2}}$,主要受到频率白噪声的影响[13, 16, 18-19]。超过2 000 s时,序列表现出复杂的、非幂率变化,并在10 000 s附近出现一个“凸起”,这是由于周期特性引起的,可见钟差序列中的周期特性极大地影响了钟差的频率稳定性。

    图  11  Allan方差计算结果
    Figure  11.  Calculated Results of Allan Variance
    图  12  Hadamard方差计算结果
    Figure  12.  Calculated Results of Hadamard Variance
    表  3  北斗卫星钟差频率稳定性计算结果
    Table  3.  Calculated Results of BeiDou Satellite Clock Offset Frequency Stability
    卫星 卫星号 千秒稳定性/10-13 万秒稳定性/10-14 日稳定性/10-15
    Allan Had Allan Had Allan Had
    BDS-2 C06 6.73 6.68 1.75 15.6 18.9 22.2
    C07 3.09 3.07 8.35 7.97 5.81 9.48
    C08 2.53 2.51 7.51 7.62 1.14 1.01
    C09 1.97 1.97 5.53 4.94 6.46 6.81
    C10 3.07 3.07 7.51 8.09 5.28 7.12
    C11 2.83 2.58 4.79 4.94 3.76 5.55
    C12 1.08 1.10 4.73 4.61 0.88 0.92
    C13 2.51 2.51 6.33 6.31 7.88 7.64
    C14 0.80 0.80 2.32 2.17 4.54 6.56
    BDS-3 C19 0.78 0.78 3.59 3.53 3.27 3.25
    C20 0.51 0.51 2.83 2.13 1.85 1.00
    C27 0.41 0.41 2.46 2.23 4.32 5.88
    C28 0.74 0.74 8.65 8.46 4.58 2.43
    C31 1.50 1.50 6.35 6.66 4.26 4.79
    C32 0.43 0.43 2.41 2.57 1.37 1.31
    C33 0.58 0.57 3.51 3.64 2.22 3.66
    C34 0.51 0.51 4.50 4.14 4.61 5.65
    下载: 导出CSV 
    | 显示表格
    表  4  两种算法频率稳定性统计结果
    Table  4.  Statistical Results of Frequency Stability for Two Algorithms
    稳定性类型 Allan Had 相似度/%
    千秒稳定性/10-13 1.77 1.75 99
    万秒稳定性/10-14 4.89 5.62 98
    日稳定性/10-15 4.77 5.60 96
    下载: 导出CSV 
    | 显示表格

    表 5给出了BDS-3和BDS-2卫星的频率稳定性指标的对比结果,统计结果表明,BDS-3稳定性整体优于BDS-2卫星,BDS-3的频率稳定性相对于BDS-2频率稳定性最高提升了75.09%。BDS-3卫星性能的提升主要得益于载荷有新型国产铷钟和被动型氢原子钟。

    表  5  BDS-2/BDS-3频率稳定性统计结果
    Table  5.  Statistical Results of Frequency Stability for BDS-2/BDS-3
    稳定性类型 算法 卫星系统 提升幅度/%
    BDS-2 BDS-3
    千秒稳定性/10-13 Allan 2.73 0.68 75.09
    Had 2.70 0.68 74.81
    万秒稳定性/10-14 Allan 5.42 4.29 20.85
    Had 6.92 4.17 39.74
    日稳定性/10-15 Allan 6.07 3.31 45.47
    Had 7.48 3.50 53.21
    下载: 导出CSV 
    | 显示表格

    表 6给出了BDS系统IGSO和MEO卫星频率稳定性的统计结果。从表 6中可以看出,MEO卫星的频率稳定性整体优于IGSO卫星,而BDS-3 IGSO卫星的日稳定性优于MEO卫星的日稳定性,主要由于BDS-3 C32卫星载荷有氢原子钟,侧面反映出氢原子钟的长期稳定性要优于铷原子钟。

    表  6  不同卫星类型频率稳定性统计结果
    Table  6.  Statistical Results of Frequency Stability For Different Satellites Types
    稳定性类型 算法 BDS-2 BDS-3
    IGSO MEO 提升幅度/% IGSO MEO 提升幅度/%
    千秒稳定性/10-13 Allan 3.32 1.57 52.66 0.97 0.59 39.03
    Had 3.30 1.49 54.77 0.97 0.59 39.21
    万秒稳定性/10-14 Allan 6.16 3.95 35.97 4.38 4.26 2.82
    Had 8.42 3.91 53.61 4.62 4.02 12.86
    日稳定性/10-15 Allan 7.58 3.06 59.62 2.82 3.48 -23.45
    Had 9.04 4.34 51.97 3.05 3.65 -19.51
    下载: 导出CSV 
    | 显示表格

    表 7给出了BDS系统搭载的不同类型原子钟的频率稳定性的统计结果(实验期间仅C32卫星启用了氢原子钟),结果表明,氢原子钟千秒稳定性、万秒稳定性和日稳定性分别高于铷原子钟40.16%、47.10%、65.60%,并且氢原子钟的长期稳定性要明显优于铷原子钟。

    表  7  不同原子钟类型频率稳定性统计结果
    Table  7.  Statistical Results of Frequency Stability for Different Atomic Clocks Types
    稳定性类型 算法 铷钟 氢钟 提升幅度/%
    千秒稳定性/10-13 Allan 0.72 0.43 40.16
    Had 0.72 0.43 40.04
    万秒稳定性/10-14 Allan 4.56 2.41 47.10
    Had 4.40 2.57 41.57
    日稳定性/10-15 Allan 3.59 1.37 61.81
    Had 3.81 1.31 65.60
    下载: 导出CSV 
    | 显示表格

    本文详细分析了BDS-2 IGSO、MEO卫星和BDS-3卫星的周期特性和频率稳定性,采用不同算法讨论了BDS系统不同卫星类型之间、不同星载原子钟类型之间的频率稳定性的差异。分析实验结果得到以下结论:(1)采用Allan方差和Had方差估计的BDS卫星钟差频率稳定性基本一致,相似度达到95%以上。(2)BDS-3卫星载荷了新一代的国产铷钟和被动型氢原子钟,其性能相对于BDS-2系统来说具有较大的提升,BDS-3卫星的频率稳定性相对于BDS-2卫星提升了20.85%~75.09%。(3)BDS卫星中MEO卫星的频率稳定性要优于IGSO卫星。(4)氢原子钟千秒稳定性、万秒稳定性和日稳定性分别高于铷原子钟40.16%、47.10%、65.60%,且氢原子钟的长期稳定性要明显优于铷原子钟。

    致谢: 感谢江苏省双创团队项目支持和国际GNSS监测评估系统(iGMAS)提供的数据产品。
  • 图  1   BDS-3卫星跟踪情况

    Figure  1.   BDS-3 Satellite Tracking Conditions

    图  2   数据有效率

    Figure  2.   Histogram of the Rate of Effective Data

    图  3   钟差解算精度

    Figure  3.   Accuracy of Clock Offset Estimation

    图  4   拟合残差直方图和正态分布图

    Figure  4.   Histograms and Normal Distribution of the Fitting Residuals

    图  5   拟合残差Q-Q图

    Figure  5.   Q-Q Distribution of the Fitting Residuals

    图  6   数据预处理流程图

    Figure  6.   Flowchart of the Data Pre-processing Method

    图  7   C34钟差数据异常值探测结果

    Figure  7.   Results of Outlier Detection for C34

    图  8   BDS-2卫星频谱分析

    Figure  8.   Frequency Analysis of BDS-2 Satellites

    图  9   BDS-3卫星频谱分析

    Figure  9.   Frequency Analysis of BDS-3 Satellites

    图  10   C14、C19、C34卫星钟差拟合残差

    Figure  10.   Fitting Residuals of C14, C19, C34 Clock Offsets

    图  11   Allan方差计算结果

    Figure  11.   Calculated Results of Allan Variance

    图  12   Hadamard方差计算结果

    Figure  12.   Calculated Results of Hadamard Variance

    表  1   北斗卫星参数

    Table  1   Parameters of BDS satellites

    卫星 卫星号 卫星类型 原子钟
    BDS-2 C06 IGSO
    C07 IGSO
    C08 MEO
    C09 MEO
    C10 MEO
    C11 MEO
    C12 MEO
    C13 IGSO
    C14 MEO
    BDS-3 C19 MEO
    C20 MEO
    C27 MEO
    C28 MEO
    C31 IGSO
    C32 IGSO
    C33 MEO
    C34 MEO
    下载: 导出CSV

    表  2   钟差精度统计/ns

    Table  2   Accuracy Statistics of Clock Error/ns

    精度 C19 C20 C27 C28 C31 C32 C33 C34
    标准差 0.74 0.82 0.63 0.54 0.58 0.88 0.74 0.82
    均方差 0.61 0.80 0.89 0.80 0.74 0.63 0.97
    下载: 导出CSV

    表  3   北斗卫星钟差频率稳定性计算结果

    Table  3   Calculated Results of BeiDou Satellite Clock Offset Frequency Stability

    卫星 卫星号 千秒稳定性/10-13 万秒稳定性/10-14 日稳定性/10-15
    Allan Had Allan Had Allan Had
    BDS-2 C06 6.73 6.68 1.75 15.6 18.9 22.2
    C07 3.09 3.07 8.35 7.97 5.81 9.48
    C08 2.53 2.51 7.51 7.62 1.14 1.01
    C09 1.97 1.97 5.53 4.94 6.46 6.81
    C10 3.07 3.07 7.51 8.09 5.28 7.12
    C11 2.83 2.58 4.79 4.94 3.76 5.55
    C12 1.08 1.10 4.73 4.61 0.88 0.92
    C13 2.51 2.51 6.33 6.31 7.88 7.64
    C14 0.80 0.80 2.32 2.17 4.54 6.56
    BDS-3 C19 0.78 0.78 3.59 3.53 3.27 3.25
    C20 0.51 0.51 2.83 2.13 1.85 1.00
    C27 0.41 0.41 2.46 2.23 4.32 5.88
    C28 0.74 0.74 8.65 8.46 4.58 2.43
    C31 1.50 1.50 6.35 6.66 4.26 4.79
    C32 0.43 0.43 2.41 2.57 1.37 1.31
    C33 0.58 0.57 3.51 3.64 2.22 3.66
    C34 0.51 0.51 4.50 4.14 4.61 5.65
    下载: 导出CSV

    表  4   两种算法频率稳定性统计结果

    Table  4   Statistical Results of Frequency Stability for Two Algorithms

    稳定性类型 Allan Had 相似度/%
    千秒稳定性/10-13 1.77 1.75 99
    万秒稳定性/10-14 4.89 5.62 98
    日稳定性/10-15 4.77 5.60 96
    下载: 导出CSV

    表  5   BDS-2/BDS-3频率稳定性统计结果

    Table  5   Statistical Results of Frequency Stability for BDS-2/BDS-3

    稳定性类型 算法 卫星系统 提升幅度/%
    BDS-2 BDS-3
    千秒稳定性/10-13 Allan 2.73 0.68 75.09
    Had 2.70 0.68 74.81
    万秒稳定性/10-14 Allan 5.42 4.29 20.85
    Had 6.92 4.17 39.74
    日稳定性/10-15 Allan 6.07 3.31 45.47
    Had 7.48 3.50 53.21
    下载: 导出CSV

    表  6   不同卫星类型频率稳定性统计结果

    Table  6   Statistical Results of Frequency Stability For Different Satellites Types

    稳定性类型 算法 BDS-2 BDS-3
    IGSO MEO 提升幅度/% IGSO MEO 提升幅度/%
    千秒稳定性/10-13 Allan 3.32 1.57 52.66 0.97 0.59 39.03
    Had 3.30 1.49 54.77 0.97 0.59 39.21
    万秒稳定性/10-14 Allan 6.16 3.95 35.97 4.38 4.26 2.82
    Had 8.42 3.91 53.61 4.62 4.02 12.86
    日稳定性/10-15 Allan 7.58 3.06 59.62 2.82 3.48 -23.45
    Had 9.04 4.34 51.97 3.05 3.65 -19.51
    下载: 导出CSV

    表  7   不同原子钟类型频率稳定性统计结果

    Table  7   Statistical Results of Frequency Stability for Different Atomic Clocks Types

    稳定性类型 算法 铷钟 氢钟 提升幅度/%
    千秒稳定性/10-13 Allan 0.72 0.43 40.16
    Had 0.72 0.43 40.04
    万秒稳定性/10-14 Allan 4.56 2.41 47.10
    Had 4.40 2.57 41.57
    日稳定性/10-15 Allan 3.59 1.37 61.81
    Had 3.81 1.31 65.60
    下载: 导出CSV
  • [1]

    Yun Q, Lou Y, Dai X, et al. Benefits of Satellite Clock Modeling in BDS and Galileo Orbit Determination[J]. Advances in Space Research, 2017, 60(12):2 550-2 560 doi: 10.1016/j.asr.2017.03.040

    [2]

    Wang D, Guo R, Xiao S, et al. Atomic Clock Performance and Combined Clock Error Prediction for the New Generation of BeiDou Navigation Satellites[J]. Advances in Space Research, 2018, 63(9): 2 889-2 898

    [3]

    Huang G, Zhang Q, Li H, et al. Quality Variation of GPS Satellite Clocks on Orbit Using IGS Clock Products[J]. Advances in Space Research, 2013, 51(6):978-987 doi: 10.1016/j.asr.2012.09.041

    [4] 黄观文. GNSS星载原子钟质量评价及精密钟差算法研究[D].长安: 长安大学, 2012 http://cdmd.cnki.com.cn/Article/CDMD-10710-1013017231.htm

    Huang Guanwen.GNSS Satellite Atomic Clock Quality Evaluation and Precision Clock Algorithm[D]. Chang'an: Chang'an University, 2012 http://cdmd.cnki.com.cn/Article/CDMD-10710-1013017231.htm

    [5] 田婕, 黄观文, 王利, 等.一种北斗二代卫星钟差的数据质量控制方法[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 Geodynamics, 2017, 37(2):215-220 http://d.old.wanfangdata.com.cn/Periodical/dkxbydz201702022

    [6]

    Zhai Z C. The New Creations and Developments of On-board Atomic Clocks for the 3rd Generation of GNSS[J]. Progress in Astronomy, 2008, 26(4):301-310 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=twxjz200804001

    [7] 王宁, 王宇谱, 李林阳, 等. BDS星载原子钟频率稳定性分析[J].武汉大学学报·信息科学版, 2017, 42(9):1 256-1 263 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201709011

    Wang Ning, Wang Yupu, Li Linyang, et al. Stability Analysis of the Space-Borne Atomic Clock Frequency for BDS[J]. Geomatics and Information Science of Wuhan University, 2017, 42(9):1 256-1 263 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201709011

    [8] 胡志刚.北斗卫星导航系统性能评估理论与试验验证[D].武汉: 武汉大学, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10486-1014136513.htm

    BeiDou Navigation Satellite System Performance Assessment Theory and Experimental Verification[D]. Wuhan: Wuhan University, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10486-1014136513.htm

    [9]

    Steigenberger P, Hugentobler U, Hauschild A, et al. Orbit and Clock Analysis of Compass GEO and IGSO Satellites[J]. Journal of Geodesy, 2013, 87(6):515-525 doi: 10.1007/s00190-013-0625-4

    [10]

    Lou Y, Liu Y, Shi C, et al. Precise Orbit Determination of BeiDou Constellation: Method Comparison[J]. GPS Solutions, 2016, 20(2):259-268 doi: 10.1007/s10291-014-0436-y

    [11]

    Chen J P, Hu X G, Tang C P, et al. Orbit Determination and Time Synchronization for New-Generation BeiDou Satellites:Preliminary Results[J]. Scientia Sinica, 2016, 46(11):85-95

    [12]

    Zhao Q, Wang C, Guo J, et al. Precise Orbit and Clock Determination for BeiDou-3 Experimental Satellites with Yaw Attitude Analysis[J]. GPS Solutions, 2018, 22(1):4-9 doi: 10.1007/s10291-017-0673-y

    [13] 李跃华, 蔺玉亭, 张达, 等.钟差数据处理方法研究[C].中国卫星导航学术年会, 上海, 2011

    Li Yuehua, Lin Yuting, Zhang Da, et al. Method Research on Clock Deviation Process[C]. China Satellite Navigation Conference, Shanghai, 2011

    [14] 黄观文, 余航, 郭海荣, 等.北斗在轨卫星钟中长期钟差特性分析[J].武汉大学学报·信息科学版, 2017, 42(7):982-988 http://ch.whu.edu.cn/CN/abstract/abstract5784.shtml

    Huang Guanwen, Yu Hang, Guo Hairong, et al. Analysis of the Mid-long Term Characterization for BDS On-orbit Satellite Clocks[J]. Geomatics and Information Science of Wuhan University, 2017, 42(7):982-988 http://ch.whu.edu.cn/CN/abstract/abstract5784.shtml

    [15] 毛亚, 王潜心, 胡超, 等.北斗三号试验卫星的钟差评估及预报[J].天文学报, 2018, 59(1):56-69 http://d.old.wanfangdata.com.cn/Periodical/twxb201801006

    Mao Ya, Wang Qianxin, Hu Chao, et al. Performance Analysis and Prediction of Clock Offsets for the BDS-3 Test Satellites[J]. Acta Astronomica Sinica, 2018, 59(1):56-69 http://d.old.wanfangdata.com.cn/Periodical/twxb201801006

    [16]

    Wang B, Lou Y, Liu J, et al. Analysis of BDS Satellite Clocks in Orbit[J]. GPS Solutions, 2015, 20(4):112-119 http://d.old.wanfangdata.com.cn/Periodical/whchkjdxxb201707018

    [17]

    Huang G, Cui B, Zhang Q, et al. An Improved Predicted Model for BDS Ultra-rapid Satellite Clock Offsets[J]. Remote Sensing, 2018, 10(1):60-67 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=remotesensing-10-00060

    [18]

    Xie X, Geng T, Zhao Q, et al. Performance of BDS-3: Measurement Quality Analysis, Precise Orbit and Clock Determination[J]. Sensors, 2017, 17(6):1 233-1 240 doi: 10.3390/s17061233

    [19]

    Montenbruck O, Steigenberger P, Prange L, et al. The Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS) Achievements, Prospects and Challenges[J]. Advances in Space Research, 2017, 59(7):1 671-1 697 doi: 10.1016/j.asr.2017.01.011

  • 期刊类型引用(31)

    1. 王中元,周圣淇,胡超,王瑞光,姚雪燕,来闯. 顾及钟差短时相关性的导航卫星实时钟差估计. 中国惯性技术学报. 2025(01): 27-35 . 百度学术
    2. 周圣淇,王中元,胡超,王瑞光,姚雪燕. 一种顾及参数稀疏的GNSS瞬时钟差模型. 大地测量与地球动力学. 2024(03): 263-266+320 . 百度学术
    3. 潘雄,黄伟凯,赵万卓,张思莹,张龙杰,金丽宏. 基于BiLSTM模型的BDS-3短期钟差预报精度研究. 测绘学报. 2024(01): 65-78 . 百度学术
    4. 曹新运,沈飞,李建成,张守建. BDS-3/GNSS非组合精密单点定位. 武汉大学学报(信息科学版). 2023(01): 92-100 . 百度学术
    5. 丁慧敏,王文贯,李浩军. BDS-3卫星钟差模型化特性分析. 大地测量与地球动力学. 2023(02): 153-157 . 百度学术
    6. 阚昊宇,胡志刚,吕逸飞,谢新,周仁宇,赵齐乐. 利用不同时间同步体制钟差评估北斗三号星载原子钟性能. 武汉大学学报(信息科学版). 2023(04): 604-610 . 百度学术
    7. 王西龙,许小龙,赵齐乐. 北斗三号系统信号质量分析及轨道精度验证. 武汉大学学报(信息科学版). 2023(04): 611-619 . 百度学术
    8. 马成,凡艳伟,张澄. 北斗三号双频长基线单历元解算性能评估. 地理空间信息. 2023(05): 91-95 . 百度学术
    9. 刘明亮,安家春,王泽民,张保军,宋翔宇. BDS-3多频伪距定位性能分析. 武汉大学学报(信息科学版). 2023(06): 902-910 . 百度学术
    10. 赖金富,史俊波,胡翾,田坤,欧阳晨皓. 北斗三号卫星导航系统在非洲地区控制测量的应用研究. 测绘通报. 2023(09): 25-29+45 . 百度学术
    11. 李坤,王潜心,闵扬海,龚佑兴,苗伟,程彤. 附有周期项的二次多项式LASSO钟差预报模型. 河南理工大学学报(自然科学版). 2022(03): 74-80 . 百度学术
    12. 王德盛,崔太岷,胡燕,杨玉锋. 多核相关向量机的BDS卫星钟差预报算法. 测绘科学. 2022(01): 40-48 . 百度学术
    13. 张熙,刘长建,章繁,吴庆,胡小华. 四大GNSS广播星历精度评估与对比分析. 武汉大学学报(信息科学版). 2022(02): 208-218 . 百度学术
    14. 马会林,吴文坛,吴秀丽. 亚太区域北斗三号对北斗二号精密单点定位的增强效果分析. 测绘通报. 2022(01): 56-60 . 百度学术
    15. 简濠骏,王逸石,张元泰,魏英韬,周裕欣. 北斗卫星广播星历精度评估与单点定位优化模型. 大地测量与地球动力学. 2022(03): 291-297 . 百度学术
    16. 朱绍涛,孙建霞. BDS-3完整星座多频标准单点定位精度分析. 地理空间信息. 2022(03): 84-87+96 . 百度学术
    17. 王棣星,李江伟,陶清瑞. 部分相对论效应对北斗原子钟性能影响分析. 测绘科学. 2022(03): 29-36+64 . 百度学术
    18. 刘路,郭金运,周茂盛,鄢建国,纪兵,赵春梅. GNSS广播星历轨道和钟差精度分析. 武汉大学学报(信息科学版). 2022(07): 1122-1132 . 百度学术
    19. 樊礼谦,焦文海,蔡洪亮,周巍,徐颖,周舒涵. 北斗三号卫星钟长期稳定性分析. 导航定位学报. 2022(04): 11-19 . 百度学术
    20. 姚文豪,李方超,高文亮,王明君. BDS-3在轨卫星广播星历精度评估. 测绘科学. 2022(07): 1-8+26 . 百度学术
    21. 伏军胜,贾小林,刘家龙,许瑾,贺延伟,张奋. BDS-3卫星与其他GNSS系统卫星原子钟性能分析. 真空与低温. 2022(05): 615-622 . 百度学术
    22. 李送强,赵兴旺,胡豪杰,刘超. BDS-3广播星历轨道、钟差精度分析. 合肥工业大学学报(自然科学版). 2021(03): 407-412 . 百度学术
    23. 刘生锋,杨文龙,谷涛. BDS-3双频组合精密单点定位精度分析. 测绘通报. 2021(03): 7-11+17 . 百度学术
    24. 王汉民. 北斗二号对北斗三号伪距单点定位精度影响分析. 全球定位系统. 2021(02): 44-48 . 百度学术
    25. 押少帅,赵兴旺,胡豪杰,刘超,陈健. 北斗三号卫星钟差短期预报与稳定性分析. 全球定位系统. 2021(03): 39-46 . 百度学术
    26. 苗伟,王潜心,程彤,李萌萌,李坤,余志浩. 基于多项式附加周期项模型的iGMAS钟差产品综合方法. 浙江理工大学学报(自然科学版). 2021(05): 649-656 . 百度学术
    27. 李骁逸,龚航,彭竞,于美婷,孙广富. 北斗系统星载原子钟周期性波动的频谱分析校正方法. 国防科技大学学报. 2021(05): 86-92 . 百度学术
    28. 胡燕,王德盛,杨玉锋. 基于EM算法优化相关向量机的BDS-3超快速钟差预报. 大地测量与地球动力学. 2021(12): 1230-1234 . 百度学术
    29. 张澄. BDS-3/GPS/Galileo/QASS兼容频率伪距单点定位精度分析. 地理信息世界. 2020(04): 135-138 . 百度学术
    30. 彭劲松. GEO卫星对BDS-3伪距单点定位性能定量提升分析. 全球定位系统. 2020(05): 62-66 . 百度学术
    31. 程彤,王潜心,胡超,李萌萌,苗伟,龚佑兴. BDS-2和BDS-3卫星原子钟特性分析. 测绘科学. 2020(12): 69-76 . 百度学术

    其他类型引用(18)

图(12)  /  表(7)
计量
  • 文章访问数:  1747
  • HTML全文浏览量:  231
  • PDF下载量:  334
  • 被引次数: 49
出版历程
  • 收稿日期:  2018-10-20
  • 发布日期:  2020-01-04

目录

/

返回文章
返回