基于优化残差组合的北斗卫星短期钟差预报研究

周仕琦，蔡成林

(湘潭大学 自动化与电子信息学院,湖南 湘潭 411105)

0 引言

随着我国北斗卫星导航系统(BDS)率先在国际上提供精密单点定位(PPP)星上服务，如何让BDS 提供更加精准的服务变成了当下的热门研究.而可靠的卫星钟差短期预报结果将极大地提升PPP 的精度，能够推动实时精密单点定位(RT-PPP)技术的研究[1].

在卫星钟差预报方面，国内外学者已经做了大量工作，对卫星钟差预报提出了若干钟差预报模型.其中，灰度预测模型(GM(1,1))、二次多项式预测(QP)模型和自回归滑动平均模型(ARMA)是比较常用的卫星钟差预报模型[2-3].(GM(1,1))模型和QP 模型建模简单，ARMA 模型能很好地反映数据间的相关性，因此这几种常用模型已经被广泛应用于测绘导航之中.但上述方法仅仅是对卫星的原始钟差进行的研究，王宇谱等[4]研究表明使用钟差一次差分序列可以有效降低预报误差，提升预报精度.然而，经过一次差分后的序列呈非线性，多项式模型不适用于非线性数据的建模；(GM(1,1))模型受发展系数和灰作用量影响较大，对非线性序列预报效果也不理想；ARMA模型虽然在钟差预报中显示出优良的特性，但在拟合模型参数时，模型的识别和定阶的不准可能会影响ARMA 预报精度[5].针对这些问题，王润等[6]利用Elman 神经网络对卫星钟差一次差分的序列进行建模预报，取得了一定的效果.然而，传统的神经网络由于使用梯度算法，容易导致局部最小而影响预报精度.极限学习机(ELM)与传统的神经网络模型相比，泛化能力更强，并且能够解决传统基于梯度算法的神经网络模型易陷入局部最小、需要合适的学习率等问题[7].

在这些研究的基础上，本文提出一种利用一次差分数据优化残差组合预测模型，并与使用钟差一次差分的传统模型的预测效果进行比较，探究组合方法的性能.

1 单一预测模型

1.1 ARMA 模型

ARMA 模型是一种被广泛应用于钟差预报领域的线性平稳时间序列模型，模型公式为

式中：e(t)为 被预测变量；e(t-i)为与被预测变量相关的滞后阶；c为常数项；εt为白噪声(WN)序列，服从均值为0 的正态分布；αi为自相关系数(i=1，2，3，···，q)；p，q分别为自回归和移动平均中的滞后阶.

p，q参数的选取对ARMA 模型的拟合效果至关重要，若参数选取不当，则会导致模型精度降低.当训练钟差样本N固定时，使赤池信息量准则(AIC)达到最小值的(p，q)就是最佳组合参数，AIC 的定义为

式中：N为卫星钟差样本个数；L为似然函数.

1.2 ElM 模型

ELM 模型是一种基于单隐层的前向神经网络模型.通过随机搜索或人为给定的方式，确定输入层的权值和隐含层的偏置，再利用广义逆矩阵理论计算输出层权值，并根据此权值计算出模型的预测值[8].其拓扑结构如图1 所示.

图1 ELM 模型结构

图1 中，xj∈Rn,j=1,2,···,N为N个样本组成的样本空间；ai=[a1,a2,···,aN]为 输入权值；bi为隐含层节点的偏置；βi为隐含层到输出层的权重；Oj为模型输出.由图1 可知，一个具有L个隐含节点的ELM 模型可表示为

式中，g(x)为激活函数.当网络模型能够零误差逼近输出样本矩阵Tm=[t1,t2,···,tm]∈Rm时，即

此时模型就能表示为

公式(5)可以简化为

式中，H为隐藏层输出矩阵.这时式(6)可看作寻找线性系统的最小二乘解，由广义逆理论可得其解为

式中，H+为H的广义逆矩阵.

2 一次差分优化残差组合模型

2.1 钟差的一次差处理及预处理

对于同一颗卫星的钟差数据而言，相邻历元在无钟跳和粗差的情况下变化量较小，从整体上说呈线性趋势.然而，由于ELM 神经网络中非线性激活函数的影响，导致其对未作差分处理的钟差序列预报能力可能不及传统模型[9].因此本文中使用钟差的一次差分序列作为样本，而要使用一次差分序列，对钟差的预处理是十分有必要的.文献[10]提出了利用中位数探测(MAD)方法处理北斗卫星钟差一次差分后的差分序列，并使用GM(1,1)与ARMA 模型组合进行预测的方法，取得了一定的效果.然而，由于北斗卫星钟差数据存在数据质量较低、粗差较多尤其是存在较多微小的跳变等问题，并且同时由于MAD 法对粗差大小不敏感，且存在参数设置问题，会导致MAD 法的粗差探测效果不佳，容易漏过许多较小的粗差和跳变.文献[11]表明在GPS 原始时间序列的粗差剔除中，四分位法(IQR)对粗差的探测效果要优于MAD法.故本文将利用四分位法和拉格朗日插值法对北斗卫星钟差一次差分数据进行预处理，然后使用经过预处理后的北斗卫星钟差一次差分数据.

一次差分原理如下：设X={x(1),x(2),···,x(n-1)}为一组n维的卫星钟差序列，将卫星钟差数据进行历元间差分，可以得到一组未作预处理的一次差分序列：

对一次差分序列进行预处理后，可以得到经过预处理后的一次差分序列

2.2 一次差分优化残差组合模型

平稳时间序列一般可以通过下式表示[12]：

式中：Lt为趋势项；Rt为残差随机项.经过一次差分的卫星钟差序列也存在趋势项和残差随机项，将ARMA模型对卫星钟差序列拟合预测的结果作为趋势项进行提取，并与ELM 模型优良的构建非线性映射能力相结合，即组成一次差分优化残差组合模型，理论上可以提升模型的预测效果.一次差分优化残差组合模型原理如图2 所示.

具体流程为：

1)对原始一次差分序列做预处理，采用IQR 法进行粗差检测并剔除粗差，得到预处理后的序列{Δx′(1),Δx′(2),···,Δx′(n-1)}.

2)利用ARMA 模型对序列进行建模，得到建模序列 {Δ(1),Δ(2),···,Δ(n-1)}，并与预处理后的序列作差，得到残差序列：

5)利用一次差分还原公式将组合模型预报结果进行差分还原，得到一次差分优化残差组合预测结果，即最终预测结果，差分还原的公式为

式中：X(k)为 第k个历元的预报值；X(n-1)为原始钟差序列的最后一个钟差值；ΔX′′(k)为一次差分优化残差组合预测值.

3 实验结果与分析

为了探究一次差分优化残差组合模型在北斗卫星钟差短期预报的效果，本文使用了2022-03-16 来源于武汉大学IGS 数据中心提供的共24 h 事后精密钟差数据，采样间隔为30 s.并选取同在中轨道地球卫星(MEO)轨道的三颗搭载不同类型原子钟的北斗三号(BDS-3)卫星C20(RAFS)、C23(RAFS)和C28(PHM)前12 h 作为建模数据进行建模，三颗卫星原始相位如图3 所示.

图3 三颗卫星原始相位

3.1 数据预处理和模型构建

由于卫星钟差的原始数据通常因为卫星钟漂或者随机噪声，不可避免地会产生一定的系统误差，当系统误差过大时，数据就会出现突兀点和缺失点，不适合直接当作建模数据进行建模分析，因此，对数据的预处理是十分重要的.对一次差分后的频率数据进行粗差探测剔除可减少误差影响，从而提高预报精度.根据之前的讨论，本文采用IQR 进行粗差检测，并与中位数探测法进行比较，IQR 具体公式如下：

式中：C为判定系数；x为钟差的频率数据；M为中位数；Q1和Q2分别为第25百分数和第75百分数.当|C|≤2 时，表示数据是十分优质的；当|C|≥3时表示为异常结果，即存在问题数据，本文将|C|＜3作为正常数据的标准，以C20 号卫星的数据为例，实验结果如图3～5 所示.

由图3(a)和图4 可知，原始相位中存在着肉眼观测不到的粗差，将相位数据转换成频率数据后，粗差将变得十分明显，这时采用合适的粗差探测手段可以实现粗差的探测与剔除.

图4 C20 频率数据

由图5 可知，MAD 法会漏过一些微小的跳变，导致探测效果不佳，从而影响模型预测精度.与MAD法相比，IQR 探测效果更好，对微小跳变更为敏感，更适合对质量数据欠佳的北斗卫星钟差数据进行粗差检测.利用IQR 对数据进行粗差检测，并确定粗差发生的位置后.可以先将粗差处的数据置0，再利用拉格朗日插值法对数据进行插值处理，完成对粗差的剔除.

图5 粗差探测方法比较

ARMA 模型参数组合的选取和ELM 神经网络隐含层神经元网络个数的确定是进行卫星钟差预报实验中极其重要的一环，本文实验中ARMA 模型参数组合如表1 所示，ELM 神经网络输入层设置为6，输出层设置为1，隐含层一般为[2,10]范围内的常数，具体根据实验确定，本文结构为6-2-1，激活函数为Sigmoid 函数.

表1 ARMA 模型参数组合

本文将采用均方根误差(RMSE)来评价预报结果的准确度，计算公式为[13]

式中：为IGS 提供的卫星钟差；xi为模型预测值.

3.2 钟差预报

为了比较模型的预报精度，本文利用GM(1,1)模型、ARMA 模型和文中的组合预报模型分别对C20、C23、C28 三颗卫星前12 h 的卫星钟差数据进行建模，预测后6 h、12 h 的卫星钟差，进行探究北斗卫星钟差短期预报中一次差分优化残差组合模型的性能的短期实验.GM(1,1)模型由于在卫星钟差一次差分前就进行了自加和自减，在进行一次差分处理后会导致精度不确定，因此本文对比试验GM(1,1)模型使用未经一次差分处理的原始钟差序列，其他模型使用经过一次差分处理后的钟差序列.实验结果如图6～7和表2 所示.

图6 各模型预报12 h 结果(局部)

图6(a)～(c)各模型预报各卫星12 h 钟差的局部结果对比，图7(a)～(c)分别为不同模型预报各卫星12 h钟差的残差的绝对值对比.表2为不同模型预报6 h和12 h 的预报结果统计，结合图表分析可知：

图7 各卫星不同模型预报残差(12 h)

表2 不同模型短期预报结果RMSE 均值统计

1)当卫星钟差原始序列呈较强的线性时，ARMA模型、ELM 模型、优化残差组合模型在使用其一次差分序列后的预测精度都高于GM(1,1)模型，在预报6 h和12 h 卫星钟差序列中都达到了亚纳米级；当卫星钟差原始序列呈较强的非线性时，使用一次差分序列的ELM 模型预测精度就会大幅下降，预测精度甚至不如GM(1,1)模型.而使用一次差分序列的ARMA模型和优化残差组合模型的预测能力明显优于其它两种模型，表现出了一定的优势.

2)在优化残差组合模型中，主要是ARMA 模型的影响，ELM 的残差预测项可看作对ARMA 模型的补偿与修正，而这种补偿效应的效果跟预处理后的一次差分序列中的趋势性有关.当卫星钟差呈线性时，其一次差分序列往往表现出比较好的趋势性，而ARMA 模型对趋势性较强的序列往往有较高的预测精度，此时ELM 模型残差预测结果的补偿具有一定的不确定性，以C20和C23 号卫星为例，组合模型在C20 上预测精度不及ARMA 模型，而在C23 号卫星上预测精度则高于ARMA 模型.当原始钟差序列呈非线性时，一次差分序列表现出较强的随机性，这时ELM 模型的残差预测结果就能在一定程度上补偿ARMA 模型，提升ARMA 模型的预测精度，以C28 号卫星为例，组合模型预报精度比ARMA 模型的预报精度在C28 上提升了38.2%，说明组合模型具有一定的效果.

3)由残差对比结果可知，使用原始卫星钟差数据的GM(1,1)模型预报精度随预报时长增加预报精度逐渐下降，稳定性最差；使用一次差分数据的ELM模型面对非线性数据时稳定性不佳；ARMA 模型和组合模型稳定性明显优于其他两种模型，且组合模型在数据呈非线性时预报精度更高，适合被用于北斗卫星钟差短期预报工作中.

4 总结

本研究将对ARMA 模型良好的趋势提取能力与构建非线性映射能力较强的ELM 神经网络模型相结合，并应用于BDS-3 钟差的预报之中，取得了一定的效果.本文首先利用IQR 代替传统的MAD 法对北斗卫星钟差一次差分序列进行粗差探测，结果表明，相比于MAD 法，IQR 更适用于对北斗卫星的卫星钟差一次差分序列的粗差探测.然后利用ARMA 模型将经过预处理后的卫星钟差一次差分序列分成趋势项和残差随机项，接着利用ELM 神经网络模型对残差部分进行建模预测，最后将ARMA模型的预测结果和ELM 神经网络的残差预测结果求和后进行差分还原.结果表明，当卫星钟差序列呈非线性时，ELM 模型的残差预测项对ARMA 模型具有一定的补偿效果，可以提升模型的预报精度.在之后的研究中，将采用优化算法对ELM 模型进行优化，并探究更多种不同模型间组合的效果.