基于奇异谱分析的BDS卫星钟差周期项提取.pdf

第 31 卷第 9 期 中国惯性技术学报 Vol.31 No.9 2023年09月 Journal of Chinese Inertial Technology Sep.2023 收稿日期：收稿日期：2023-04-19；修回日期：修回日期：2023-07-30 基金项目：基金项目：国家自然科学基金（11503031）；陕西省自然科学基础研究计划（2023-JC-YB-057,2022-JM031）作者简介：作者简介：雷雨（1983），男，讲师，博士，从事时间保持方法研究。文章编号：文章编号：1005-6734(2023)09-0909-09 doi.10.13695/ki.12-1222/o3.2023.09.008 基于奇异谱分析的基于奇异谱分析的 BDS 卫星钟差周期项提取卫星钟差周期项提取 雷 雨1，赵丹宁2（1.西安邮电大学 计算机学院，西安 710121；2.宝鸡文理学院 电子电气工程学院，宝鸡 721016）摘要：摘要：在建立星座自主时间基准时，必须扣除卫星钟差的周期波动，以避免将其引入系统时间。为准确地扣除周期波动，提出一种基于奇异谱分析（SSA）的 BDS 卫星钟差周期项提取方法。首先，经SSA 分解获得钟差信号的多个重构成分；然后，引入重标极差分析方法计算各重构成分的 Hurst 指数，根据 Hurst 指数辨识钟差的低频信号主导分量和高频信号随机分量，从而实现周期项的重建提取。对不同类型的 BDS 卫星钟进行了分析和研究，结果表明：所提方法能够准确提取钟差周期项，所提取周期项的频谱特征比钟差多项式拟合残差的频谱更为清晰；利用多项式拟合方法扣除周期项后，万秒稳和日稳平均分别提升 4.1%和 2.4%，而利用所提方法扣除周期项后，万秒稳和日稳平均分别提升 20.9%和 34.1%，提高了卫星钟频率稳定度，为星座自主时间基准的建立提供了基础。关 键 词：关 键 词：北斗卫星导航系统；卫星钟；周期项；奇异谱分析；频率稳定度 中图分类号：中图分类号：P228.1 文献标志码：文献标志码：A Periodic term extraction of BDS satellite clock offset based on singular spectrum analysis LEI Yu1,ZHAO Danning2(1.School of Computer Science&Technology,Xian University of Posts&Telecommunications,Xian 710121,China;2.School of Electrical&Electronic Engineering,Baoji University of Arts and Sciences,Baoji 721016,China)Abstract:In establishing the autonomous time reference of the navigation constellation,the periodic fluctuations of satellite clock offset must be deducted to avoid introducing them into the system time.In order to accurately remove the periodic fluctuations,a decomposition,identification and rebuilt method based on singular spectrum analysis(SSA)is proposed to extract the periodic terms of BeiDou satellite clock offset.Firstly,multiple reconstructed components of the clock offset signal are obtained by SSA decomposition.Then,rescaled range analysis is introduced to calculate the Hurst index of each reconstructed component.The dominant low-frequency components and random high-frequency components of clock offset subsequently are identified according to the Hurst index,so as to realize the reconstruction and extraction of the period term.The analysis and study of different types of BDS satellite clocks show that the proposed method can accurately extract the period terms of satellite clock offset,and the spectral features of the extracted periodic terms are clearer than that of the polynomial fitting residuals.While the mean improvements are 4.1%and 2.4%respectively for the frequency stability at the ten-thousand and daily sampling intervals via the polynomial fitting method,through the proposed method the improvements can reach up to 20.9%and 34.1%,respectively,which improves the frequency stability of the satellite clock and provides the reference for the establishment of an autonomous time reference for the navigation constellation.Key words:Beidou navigation satellite system;satellite clock;periodic terms;singular spectrum analysis;frequency stability 910 中国惯性技术学报 第 31 卷 卫星导航系统的定位、导航和授时服务依赖于准确的时间基准，目前北斗卫星导航系统（Beidou Navigation Satellite System,BDS）时间基准的建立和保持由地面监测站原子钟通过综合原子时算法计算产生，这种模式过度依赖地面站，一旦地面监测站出现异常，系统时间基准就会不连续。为保证系统时间基准的连续性和稳定性，必须降低对地面监测站的依赖。随着 BDS 卫星数量与星载钟性能的不断提高，利用星载钟建立和保持星座自主时间基准，对提高系统的生存能力具有重要的现实意义1,2。受星载钟自身因素和太空环境的影响，卫星钟差序列一般含有趋势项、周期项和随机项3，其中，星载钟受太空环境影响的主要表现为钟差的周期波动，导致星载钟在轨性能逊于地面钟。卫星钟差的周期波动不仅和其轨道周期相关，而且和日食周期、测量性能等有关4,5。在未扣除卫星钟差的周期波动之前，不应将星载钟纳入星座自主时间基准的建立和保持，以避免将周期波动引入系统时间。现有的卫星钟差特征分量提取方法有多项式拟合、傅里叶变换、小波变换和经验模态分解等，Huang 等构建多项式模型拟合提取钟差周期项，提升了 BDS 卫星钟差预报精度3；李骁逸等基于傅里叶变换的频谱分析方法校正 BDS 卫星钟的周期波动，校正后 BDS 各类星载钟的万秒稳均得到提高6；雷雨等利用小波变换将钟差序列分解成具有不同频率特征的信号分量，并根据各分量的特征构建预报模型，改进了钟差预报精度7；梁益丰等提出以完备集合经验模态分解为基础的钟差信噪分离方法，并将其应用于 BDS 卫星钟差周期项识别，扣除周期项后卫星钟万秒稳获得提升8。上述方法虽然在钟差特征分量分离和提取方面取得了较好的效果，但均存在一定的局限性。传统的多项式拟合方法无法提取钟差序列中的时变信号，提取的钟差趋势项、周期项和随机项不够准确；傅立叶变换将信号从时间域转化到频率域，以分离信号中的周期分量，但该方法不适用于非线性、非平稳信号；经验模态分解能将非线性、非平稳信号分解为一系列固有模态分量，并根据不同模态分量的特征识别周期成分，但这类方法存在端部效应和模态混叠的问题。奇异谱分析（Singular Spectrum Analysis,SSA）是20 世纪 90 年代兴起的一种研究非线性、非平稳信号的有效方法。根据所观测到的时间序列构造出轨迹矩阵，并对轨迹矩阵进行分解和重构，从而提取出代表原时间序列不同成分的信号，如长期趋势信号、周期信号和噪声信号，从而对时间序列的波动特征进行分析9,10。SSA 方法提取信号不需要先验信息，不受正弦波假定的约束，能够较好地从含噪声的时间序列中提取时变信号，目前已在 GPS 坐标时间序列分析、多路径效应消除、钟差分析和预报11,12等方面取得了成功应用。其中，肖胜红等11提出 SSA 和傅立叶带通滤波器相结合的卫星钟差周期项提取方法，有效抑制了重建周期项时频率混叠和边界效应；Xue 等12利用SSA 分离卫星钟差的趋势项和随机项，提高了钟差短期预报精度。上述研究表明，SSA 可以较好地从复杂卫星钟差序列中提取趋势项和周期项信息。本文在分析卫星钟差特性的基础上，提出一种基于 SSA 的“分解辨识重建”钟差周期项提取方法。首先应用 SSA 对钟差序列进行分解，获得时间序列的多个重构成分，然后利用重标极差分析（Rescaled Range Analysis,R/S）方法计算各重构成分的 Hurst 指数，最后根据 Hurst 指数对各重构成分的波动特征进行辨识，进而完成钟差趋势项、周期项和随机项的分离和提取。利用 BDS 卫星钟差数据分析和检验了本文方法的有效性与实用性。1 星载原子钟物理特性 1 星载原子钟物理特性 1.1 星载原子钟物理模型 受星载原子钟自身因素及空间外界环境的影响，星载原子钟信号不仅存在趋势性变化，而且表现出复杂的周期性波动。星载原子钟模型通常可表示为13：20001122301()sin(2)2()()d()tx txy tdtAf tW tW sst (1)其中，()x t为t时刻原子钟相位；0 x为原子钟初始相位；0y为原子钟钟速；d为原子钟频漂；A、0f和分别为原子钟相位周期分量的振幅、频率和初始相位；11()W t和220()dtW ss表示两种起主导作用的原子钟 噪声，1()W t和2()W s分别为两个独立的维纳过程，1和2分别为两个维纳过程的扩散系数，用来表示两种原子钟噪声强度；3()t 为原子钟测量噪声，3表示测量噪声强度。与趋势分量相比，周期分量和随机分量的数量级较小，但对频率稳定度的影响会逐渐累积。1.2 各分量对频率稳定度的影响 阿伦方差与原子钟频漂、相位白噪声、维纳过程的扩散系数、周期分量之间的关系可表示为13：2422222230122203sin()11()23()yfdAf (2)其中，2()y表示平滑时间为时的阿伦方差；21/、22 和2233/分别表示频率白噪声、频率随机游走噪声和相位白噪声对频率稳定度的影响。根据式(2)可知，当星载原子钟在轨运行时钟差周期特征显著时，第9期 雷雨等：基于奇异谱分析的BDS卫星钟差周期项提取 911 频率稳定度也将出现周期性波动，表现为时域频率稳定度的异常突变，从而影响原子钟噪声系数的拟合估计。因此，在建立星座自主时间基准时，必须消除卫星钟差的的周期性效应，以免将其引入系统时间。2 基于2 基于 SSA 的卫星钟差周期项分析方法的卫星钟差周期项分析方法 基于SSA的卫星钟差周期项分析方法主要包括钟差数据预处理、钟差分解与重构以及钟差重构成分辨识等过程，图1给出卫星钟差周期项分析流程。图 1 卫星钟差周期项分析流程 Fig.1 Analysis process of satellite clock offset periodic terms 已知卫星钟差序列为12,Nx xx，钟差周期项分析的具体步骤为：1）钟差数据预处理 运用中位数法检测和剔除钟差序列中的粗差，若钟差的一次差分数据ix(1,21iN)满足式(3)即认为ix为粗差14：MAD=medianMAD=median/0.6745iiixmnmxxm (3)其中，m表示钟差一次差分序列的中位数；n表示检测阈值，为保证粗差检测的准确性和可靠性，本文取=6n。对于粗差点，利用线性内插法对其进行插补。2）钟差序列分解 利用SSA对数据预处理后的钟差序列进行分解，主要计算过程为9,10：奇异值分解。根据钟差序列构建时滞矩阵121,NMXY YY，矩阵X的每行、每列为钟差序列的子序列，即T121,ii Mx xxY，11i NM，其中M表示时滞窗口长度。对时滞矩阵的协方差矩阵TWXX进行奇异值分解，获得特征值120M及矩阵W对应的特征向量是12M、UUU，其中jU称为时间正交函数。令T(1,2)jjjjMVX U，则时滞矩阵的奇异值分解可写为：12+IXXX+X (4)其中，max|0jIj；初等矩阵TjjjjXU V。对角平均。计算时滞矩阵X在jU上的投影：T,i jijaX U (5)其中，,i ja称为时间主成分。根据时间正交函数与时间主成分进行对角平均，得到钟差序列的重构成分（Reconstructed Component,RC）：*,1*,*11*,11 11 1 1ii k jj kkMi ji k jj kkNKi k jj kk i KiMizMiKMKi NNi aUaUaU (6)其中，*min(,1)MM NM，*max(,1)KM NM。分组。卫星钟差序列中通常包含多种周期性变化、趋势性变化和随机性变化，利用w-correlation分析RC成分之间的相关性，将信号特征相似的RC分组。将不同RC分别用()iZ、()jZ表示，则重构分量之间的w-correlation可表示为15：()()()(),()()()()(),1,ijiji jwwiijwNijkk ik jkw z zZZZZZZZZZ 1,i j M(7)其 中，kw为 权 重 系 数，其 定 义 为min(,)kwk M Nk；,i j的绝对值越接近于1，说明()iZ、()jZ两者之间的相关性越强。根据经验，,0.6i j即认为两者之间存在相关性，因此，将RC成分之间相关性大于0.6的两者视为是同一组信号。3）RC成分辨识 一个时间序列在一段时间内的波动如何随时间跨度大小而变化往往可以揭示该时间序列的特性，对于卫星钟差的趋势项和周期项，其时间序列当前或过去的取值以远超随机扰动所能达到的程度影响该序列在未来的取值，统计学上称为时间序列存在长期记忆性（Long-term Memory）和长时间相关效应，而随机项的未来值和当前或过去的取值相关性不强，长期记忆性较弱，属于均值回复过程。由英国水利学家赫尔斯特提出Hurst指数是用来衡量时间序列是否有长期记忆性的一个指标16,17，体现了时间序列的自相关性，尤其能够反映时间序列中隐藏的长期趋势。本文采用Hurst指数H来辨识RC成分序列特征，以定性分析钟差序列中的信号主导分量和噪声主导分量。若RC成分的Hurst指数0.51H，则RC成分属于趋势分量或周期信号主导分量，且H越大，规律性越强；若00.5H，则判定RC成分为随机噪声主导分量，912 中国惯性技术学报 第 31 卷 且H越小，随机性越强。基于R/S分析法计算Hurst指数的原理为16,17：对于时间序列12,Nx xx，将其分为C个长度为L的等长子区间，/LN C，则累积离差为：,1()kc kc icidxe，1,2kL (8)其中，ce为第c个子区间序列的平均值，1,2cC。令第c个子区间序列的标准差为cS，则极差cR为：,max()min()cc kc kRdd (9)定义重标极差L为：11CHcLccROLCS (10)其中，H为Hurst指数。对L和L进行双对数线性回归拟合，回归方程的截距就是式(11)中的常数O，而斜率就是H。3 算例分析 3 算例分析 采 用 德 国 地 球 科 学 研 究 中 心（Deutsches Geoforschungs Zentrum,GFZ）提供的BDS精密卫星钟差数据（ftp:/ftp.gfz-potsdam.de/GNSS/products/mgex），选取钟差序列连续性较好的C05、C09、C11、C43星载钟，包括GEO（C05）、IGSO（C09）与MEO（C11与C43）轨道类型，铷原子钟（C05、C09、C11）与氢原子钟（C43）类型，选择2022年7月29日至8月11日共14天4032个数据点，数据采样间隔为5 min。利用SSA对时间序列进行分析，重点是确定窗口长度M，一般1/2MN且为周期的最小公倍数。BDS卫星钟差序列中通常存在6 h、12 h和24 h的周期18-20，因此，选择窗口长度M为卫星钟差周期的最小公倍数，即288M。3.1 钟差序列分解和重建试验 由于篇幅限制，以C05卫星钟差为例展示钟差序列分解和重建的过程及效果。卫星钟差数据经SSA分解获得288个RC序列，频率按照RC阶次依次从低到高排列。原始钟差序列与所有RC之和的残差均方根为1.1710-19 s，表明原始钟差序列与重建序列之差接近0，证明SSA用于钟差序列分解和重建是可行的。C05卫星钟差前20阶RC成分之间的w-correlation分析结果如图2所示。根据SSA分组原理，当RC之间的相关系数大于一定的阈值，则判断属于同一周期信号。若阈值选择过大，会导致一些周期信号被过滤掉，反之则不能较好地排除一些周期性不明显的RC成分。经实验，本文选择相关性阈值为0.6。从图2可以看出，当阶次大于14时，各RC之间就不能很好地相互分离，说明非周期变化占较大部分。图 2 前 20 阶 RC 成分的 w-correlation Fig.2 W-correlation of the first 20 RCs 根据w-correlation分析结果，将,0.6i j的RC成分合并为同一周期信号，C05卫星钟差的分析结果如图3所示。图 3 C05 卫星钟差序列及其 RC1RC22序列 Fig.3 Time-series of the C05 satellite clock offset and the first twenty-two RCs 第9期 雷雨等：基于奇异谱分析的BDS卫星钟差周期项提取 913 从图3可以看出，第1阶RC1序列与原始钟差序列变化趋势非常相似，代表卫星钟差趋势项；RC2+RC3、RC4+RC5+RC6、RC7+RC8+RC9+RC10、RC11+RC12+RC13合成信号表现为显著的周期性变化，代表卫星钟差周期项；第1422阶RC14RC22序列表现为不规则的拟周期性变化，RC22以后的高阶RC序列则表现为随机性变化，限于篇幅未在图3中给出。为分析RC序列的波动特性，利用R/S分析法计算C05卫星钟差的各阶RC序列的Hurst指数，结果如图4所示。图 4 C05 卫星钟差各 RC 序列的 Hurst 指数 Fig.4 Hurst exponent of each RC sequence for C05 satellite clock offset 从图4可以看到，第1阶RC1序列的Hurst指数为0.99，与原始钟差的Hurst指数0.98几乎相等，说明RC1序列几乎不含波动分量，可直接作为趋势项；第213阶RC2RC13序列的Hurst指数在0.51至0.94之间，且RC4和RC5序列的Hurst指数大于RC3，这是由于RC4和RC5的周期性特征比RC3更显著，RC2RC13属于钟差低频信号主导分量；其余RC序列的Hurst指数小于0.5，属于钟差高频不规则信号分量。图4的Hurst指数计算结果与图2的RC序列变化特征具有很好的对应性，说明本文钟差分量辨识方法具有合理性。为进一步说明本文方法的合理性，利用快速傅立叶变换（Fast Fourier Transform,FFT）对C05卫星钟差的RC成分进行频谱分析，分析合并信号的周期性特征，结果如图5所示。从图5可以发现，RC2+RC3、RC4+RC5+RC6、RC7+RC8+RC9+RC10、RC11+RC12+RC13合 成 信 号 的 频 谱 具 有 周 期 性 特 征，而RC14+RC15+RC288合成信号的频谱是非周期性，说明RC14以后的高阶RC序列属于高频不规则分量，进一步验证了本文钟差分量辨识方法的有效性。图 5 C05 卫星钟差 RC 序列的频谱 Fig.5 Frequency spectrum of the RC sequence of C05 satellite clock offset 将C05卫星钟差的RC2RC13序列相加获得钟差低频周期分量，将第14阶RC14以后剩余的RC序列相加获得钟差高频随机分量，其他三颗卫星钟差也作类似处理，结果如图6所示。从图6可以看出，C05、C09、C11和C43四颗卫星钟差的变化趋势均为递增，表明SSA方法能准确地分解和重建卫星钟差序列中的趋势项、周期项和随机项，即使近似平稳的微小波动也能检测；卫星钟差呈复杂的多周期性变化，这是914 中国惯性技术学报 第 31 卷 由于卫星钟周期性变化不仅与卫星轨道周期相关，还与日食周期、测量性能等多种因素有关，卫星钟差周期项的数值普遍大于随机项，其中，C05和C09卫星钟差周期项的数值在ns量级，C11和C43卫星钟差周期项的数值在亚ns量级，说明GEO和IGSO卫星钟受周期项的影响更大；随机项呈不规则性变化，其中，C43卫星钟随机项的数值小于其他三颗卫星，反映了BDS-3星载氢钟的优良性能。图 6 四颗卫星钟差趋势项、周期项和随机项的重建结果 Fig.6 Reconstruction results of the trend term,period term and random term of the four satellite clock offset 3.2 周期分量特性分析 目前常用多项式拟合残差频谱分析方法对BDS卫星钟的周期特性进行分析，研究表明卫星钟的主周期与轨道周期耦合18-20（通常近似为其卫星轨道周期的1倍或1/2倍）。为更好地分析BDS卫星钟的周期特性，利用SSA+FFT方法对BDS卫星钟差序列中存在的周期分量进行检测，四颗卫星钟差SSA重建周期项的频谱分析结果如图7所示。为进行对比分析，二次多项式拟合残差频谱分析结果也在图7中给出。通过对比分析卫星钟差重建周期项和拟合残差的频谱分析结果，可以得到如下结论：1）多项式拟合方法无法完整地提取卫星钟差趋势项，导致频谱图中靠近右侧纵轴处出现残余的趋势信号。拟合残差中包含的剩余随机分量导致靠近左侧纵轴处出现大量不规则波动，而钟差SSA重建方法提取的周期信号低频处没有出现不规则峰值，高频部分没有出现异常波动，周期分量的频谱特征明显，证明SSA重建方法提取的周期信号不包含趋势分量和随机分量，周期信号提取能力优于多项式拟合方法。2）C43卫星钟差拟合残差和重建周期项频谱分析检测的周期存在区别，钟差拟合残差频谱分析检测的周期按幅值排序约为12 h、17 h、21 h、8 h、6 h、4.3 h和3.4 h等，重建周期项频谱分析检测的周期排序约为12 h、17 h、8 h和21 h。卫星钟差的周期性变化主要由卫星运行所致，拟合残差的高频周期约为主周期12 h的公约数，这些高频周期分量可能是在随机分量的干扰下耦合产生的，SSA信号重建方法能有效地分离和提取钟差周期分量和随机分量，避免了这种耦合现象。3）卫星钟周期与其卫星轨道周期密切相关，C05卫星钟周期为24 h和12 h，约为GEO卫星轨道周期24 h的1倍和1/2倍；C09卫星钟周期为24 h、16 h和12 h，约为IGSO卫星轨道周期24h的1倍、2/3倍和1/2倍；C11和C43卫星钟周期为24 h、17 h、12 h和8 h，近似为MEO卫星轨道周期12 h的2倍、4/3倍、1倍和2/3倍。其中，C11卫星钟差拟合残差频谱分析未检测到8 h周期信号，这是由于钟差拟合残差中的高频不规则信号干扰了频谱分析结果，周期分量第9期 雷雨等：基于奇异谱分析的BDS卫星钟差周期项提取 915 被干扰信号淹没，重建周期项频谱中部分周期信号的波峰附近也存在多余信号，这可能由多周期频谱旁瓣引起，但相比钟差拟合残差频谱有明显改善；MEO卫星钟主周期分量的幅值比其他轨道类型的卫星钟低一个数量级，说明MEO卫星钟受周期项的影响较小。4）本文方法与多项式拟合方法对于C05和C09卫星钟差的周期项分析结果相同，但对于C11和C43卫星钟差的分析结果却存在差异。具体而言，C11和C43卫星钟差的多项式拟合残差频谱不仅受高频信号干扰，而且存在残余趋势信号，这是由于多项式拟合方法仅能提取钟差序列中的“平均”趋势性效应，无法反映钟差信号的时变特性，对于变化复杂的C11和C43卫星钟差（图6）多项式拟合残余的趋势分量和高频分量势必影响钟差的频谱分析结果。图 7 四颗卫星钟差的频谱 Fig.7 Frequency spectrum of the four satellite clock offset 3.3 频率稳定度分析 为进一步检验本文方法分离和提取卫星钟差周期项的准确性，分别计算四颗卫星钟差扣除周期项前后的频率稳定度，定量分析卫星钟周期项对频率稳定度的影响。为消除铷钟和氢钟频漂对频率稳定度分析的影响，同时提高方差估计值的置信度，选择能较好消除频漂影响的重叠哈达玛方差计算卫星钟的频率稳定度14，同时对比本文方法和多项式拟合方法的周期性波动校正效果，结果如图8所示。从图8可以看出，由四颗卫星的原始钟差计算的重叠哈达玛方差曲线均存在不同程度的隆起“鼓包”现象，其中C05和C09卫星钟差的重叠哈达玛方差曲线隆起现象更为明显。根据式(2)易知，这种隆起异常程度与周期项大小密切相关，即周期项越大，隆起越明显。对比图7中的四颗卫星钟差的频谱不难发现，C05和C09卫星钟差的周期项幅值大于C11和C43的周期项幅值；重叠哈达玛方差曲线“鼓包”与重建周期项的重叠哈达玛方差曲线的凸起比较吻合，表明SSA重建钟差周期项能表征周期项对频率稳定度的影响；在扣除钟差周期项后，各卫星钟在不同取样时间内频率稳定度均有一定程度的提高。从卫星轨道类型而言，C05 GEO和C09 IGSO卫星钟频率稳定度受周期项的影响最为明显，峰值约为7.510-14，C11 MEO和C43 MEO卫星钟频率稳定度受周期项的影响较小，峰值分别约为610-14和3.510-14；从卫星钟类型而言，BDS-3星载氢钟频率稳定度受周期项的影响最小，这与卫星钟周期项幅值大小密切相关（图7）。多项式拟合校正方法对C05和C09卫星钟差的频率稳定度提高效果较好，但对C11和C43卫星钟差的提高效果有限，结合图7可以发现，这是由于C05和C09卫星钟差的多项式拟合残916 中国惯性技术学报 第 31 卷 差的频谱受残余趋势信号的影响较小，而C11和C41卫星钟差的多项式拟合残差中则包含较强的残余趋势信号；相对于多项式拟合校正方法，无论何种轨道类型和何种星载钟类型，本文方法扣除周期项后的重叠哈达玛方差曲线均有效地削去了原始钟差稳定度曲线的隆起“鼓包”，频率稳定度均得到提高。图 8 四颗卫星钟差的重叠哈达玛方差 Fig.8 Overlapping Hadamard variance of the four satellite clock offset 四颗卫星钟差扣除周期项前后的频率稳定度数值结果如表1所示。由表1可以看到，利用多项式拟合方法扣除周期项后，C05、C09、C11和C43卫星钟的万秒稳分别提高6.5%、6.4%、0.7%和2.8%，C11和C43卫星钟日稳分别提高7.7%和1.9%，扣除周期项后C05和C09卫星钟的日稳没有得到改善，平均而言，万秒稳和日稳平均分别提升4.1%和2.4%；而利用本文方法扣除周期项后，四颗卫星钟的万秒稳和日稳都得到一定的改善，且相对于多项式拟合方法对卫星钟的万秒稳和日稳的提高效果更为明显，C05、C09、C11和C43卫星钟的万秒稳分别提高21.0%、23.1%、17.7%和21.6%，日稳分别提高48.8%、54.7%、20.0%和13.1%，万秒稳和日稳平均分别提升20.9%和34.1%，其中，C09 IGSO卫星钟的提高幅度最大，万秒稳和日稳分别由4.9810-14和2.8710-14提高至3.8310-14和1.3010-14，MEO卫星钟的日稳提高幅度较小，这与MEO卫星钟的主周期（12 h）有关，但由图8可知，MEO卫星钟频率稳定度提高最明显的取样时间在312 h之间，对于以MEO卫星为主的BDS-3系统，这种稳定度提高非常有利于星座时间基准的建立和保持。表 1 四颗卫星钟的万秒稳和日稳 Tab.1 Ten thousand second stability and daily stability of four satellite clocks 卫星序号和类型 万秒稳(110-14)日稳(110-14)本文方法改善率%原始钟差 多项式方法本文方法原始钟差多项式方法本文方法 万秒稳 日稳 GEO C05（BDS-2 铷钟）6.14 5.74 4.85 2.52 2.52 1.29 21.0 48.8 IGSO C09（BDS-2 铷钟）4.98 4.66 3.83 2.87 2.87 1.30 23.1 54.7 MEO C11（BDS-2 铷钟）7.25 7.20 5.97 1.95 1.80 1.56 17.7 20.0 MEO C43（BDS-3 氢钟）4.95 4.81 3.88 2.13 2.09 1.85 21.6 13.1 4 结 论 4 结 论 本文分析了卫星钟周期项对频率稳定度的影响，阐述了周期项分离和提取的必要性，提出了融合SSA分解和R/S分析方法的卫星钟周期项分离和提取方法。利用不同类型的BDS卫星钟差数据进行分析和研究，结果表明，所提方法能有效地分离卫星钟差序列钟的趋势项、周期项和随机项，所提取的周期项频谱图比钟差拟合残差频谱图更加清晰，周期性特征更为第9期 雷雨等：基于奇异谱分析的BDS卫星钟差周期项提取 917 明显。通过对比分析扣除周期项前后卫星钟频率稳定度的差异发现，原始钟差的频率稳定度曲线存在隆起“鼓包”现象，且与周期项稳定度曲线的异常隆起吻合度很高，印证了所提方法的有效性；利用所提方法扣除周期项后，卫星钟的频率稳定度提高显著，万秒稳和日稳平均分别提升20.9%和34.1%，这非常有利于高稳定度星座时间基准的建立和保持。本文方法不仅可以应用于卫星钟周期项提取，而且在原子钟信号降噪、钟差预报和时间尺度算法等方面具有潜在的应用价值。参考文献（参考文献（References）：）：1 Chen G,Xing N,Tang C,et al.Clock ensemble algorithm test in the establishment of space-based time referenceJ.Remote Sensing,2023,15(5):1227.2 程梦飞，王宇谱，薛申辉，等.北斗二号卫星导航系统运行末期卫星钟性能评估J.中国惯性技术学报，2020,28(3):372-379.Cheng M,Wang Y,Xue S,et al.The performance evaluation of satellite clocks for BD-2 navigation satellite system in the final phaseJ.Journal of Chinese Inertial Technology,2020,28(3):372-379.3 Huang G,Cui B,Zhang Q,et al.An improved predicted model for BDS ultra-rapid satellite clock offsetsJ.Remote sensing,2018,10(1):60.4 Cao Y,Huang G,Xie S,et al.An evaluation method of GPS satellite clock in-orbit with periodic terms deductedJ.Measurement,2023,214(15):112765.5 Wang D,Li M.,Xue H,et al.Analysis of the J2 relativistic effect on the performance of on-board atomic clocksJ.GPS Solutions,2023,27(3):114.6 李骁逸，龚航，彭竞，等.北斗系统星载原子钟周期性波动的频谱分析校正方法J.国防科技大学学报，2021,43(5):86-92.Li X,Gong H,Peng J,et al.Correction method for periodic fluctuation of the BDS satellite atomic clocks based on spectrum analysisJ.Journal of National University of Defense Technology,2021,43(5):86-92.7 雷雨，赵丹宁，李变，等.基于小波变换和最小二乘支持向量机的卫星钟差预报J.武汉大学学报：信息科学版，2014,39(7)：815-819.Lei Y,Zhao D,Li B,et al.Prediction of satellite clock bias based on wavelet transform and least squares support vector machinesJ.Geomatics and Information Science of Wuhan University,2014,39(7):815-819.8 梁益丰，许江宁，吴苗，等.北斗卫星钟差的CEEMDAN分解与周期项提取方法J.中国惯性技术学报，2022,30(4):476-484.Liang Y,Xu J,Wu M,et al.CEEMDAN decomposition and periodic term extraction method of the BDS clock biasJ.Journal of Chinese Inertial Technology,2022,30(4):476-484.9 Vautard R,Yiou P,Ghil M.Singular-spectrum analysis:A toolkit for short,noisy chaotic signalsJ.Physica D:Nonlinear Phenomena,1992,58(1-4):95-126.10 Hassani H,Zhigljavsky A.Singular spectrum analysis:methodology and applicatio