一种北斗卫星钟差短期预报组合方法

赵金霞，张小宁

( 1. 重庆三峡职业学院, 重庆 404155；2. 甘肃煤田地质局综合普查队, 甘肃 天水741000 )

0 引 言

2020年6月23日9时43分，我国在西昌卫星发射中心用长征三号乙运载火箭，成功发射第55颗北斗导航卫星，暨北斗三号(BDS-3)最后一颗全球组网卫星. 至此，我国BDS-3全球组网基本系统空间星座部署任务完成. 2020年7月31日上午，BDS-3全球卫星导航系统(GNSS)正式开通，北斗大规模应用已经开始. 北斗卫星实时精密单点定位(PPP)中需要提供高精度钟差值，但是国际GNSS服务(IGS)中心提供的精密钟差文件滞后14天，所以实时预报卫星钟差显得至关重要. 关于卫星钟差预报有很多方法，如二次多项式拟合模型、灰色理论预报模型、BP神经网络模型、小波变换模型等. 但通常单一模型很难准确地预测钟差文件，有些专家学者将各单一模型组合成新的模型，如文献[1]中将灰色理论模型与BP神经网络模型组合来实现卫星钟差的预报，通过灰色理论预测差值作为BP神经网络的输入值，将灰色预报的数值与BP神经网络预报的差值结合，作为最终预报值得到了较好的预报效果；文献[2]中提出将灰色模型与差分整合移动平均自回归(ARIMA)模型结合，利用灰色模型分离出钟差的趋势项，随后利用ARIMA模型对灰色模型预报残差建模分析，最后将二者预报的结果相加即可得到最终预报值；这些方法均得到了良好的预报效果.

本文探讨基于灰色理论模型与BP神经网络模型，通过科学定权策略确定二者的权重，实现最优权组合钟差预报，并结合北斗钟差数据计算与分析，进行应用验证，得出一些有益结论.

1 最优权组合模型钟差预报

1.1 基于最优权组合模型的建立

组合预测模型有很多种，常见的如等权组合模型、非线性组合模型、变权组合模型等[3-6]. 基于这些研究基础和实践认识，本文建立灰色模型与BP神经网络组成的北斗卫星钟差预报最优组合模型.

且wt(i)≥0 .

基于以上推导得出n钟单模型的最优权组合模型

式(2)中，eˆ(i) 为 最优权组合模型的第i个历元钟差预报值.

为获得该组合模型的最优预报结果，最优加权系数应满足下式到达最小：

由

可知

同时

因此得式(5)的最优权的解为

式(7)中，R为N维的列向量矩阵[7].

2 实验与分析

2.1 实验一

为分析最优权组合模型预报钟差的有效性与精度评定，采用武汉大学卫星导航定位技术研究中心通过多星联合解算的精密钟差数据，采样间隔为30 s，考虑到时间段内采样率过大，故从数据序列中采样间隔为5 min的钟差数据为例进行分析，选择2020年2月10日当天数据，年积日为41. 分别选择各类代表卫星，选取地球静止轨道(GEO)卫星C03，倾斜地球同步轨道(IGSO)卫星C13，中圆地球轨道(MEO)卫星C11，这3颗卫星均搭载铷原子钟，且均是服役的北斗二号(BDS-2)在轨卫星[8-9]. 分别将灰色模型(GM(1,1))、BP1神经网络模型、最优权组合预测模型分别运用在这三颗卫星钟差预报中. 其中每种模型预报时，均是采用前20个历元的卫星钟差数据建立模型求算相应的预报参数，预报后20历元的卫星钟差数据，与武汉大学北斗数据处理中心提供的精密钟差数据作比较，以此来分析每种预报模型的有效性.

为定量比较每种预报模型的优劣，本文采用均方根误差(RMSE)作为衡量指标，其计算公式为

式中：ai是由IGS/MGEX(Multi-GNSS Experiment)中心提供的钟差值；aˆi是各模型预报值[10].

1) 由表1可知，三种卫星的神经网络预报模型权系数均高于GM(1,1)模型，BP1神经网络模型是一种优良的单模型预报钟差方法；不同种类的北斗卫星具有不同的权重系数.

表1 各卫星组合模型的权重系数表

2) 由图1～4、表2可知，灰色理论模型预报的钟差模型呈线性，灰色理论模型预报前期预报精度较神经网络模型预报精度好，但随着预报时间的延长，预报精度降低. 预报后期精度均较BP神经网络模型差.BP神经网络模型预报精度前期不如灰色理论模型，但随着预报时间的延长，预报精度高于灰色理论模型. 并且预报趋势与钟差真值吻合较好. BP神经网络是一种优秀的非线性预报模型.

图1 C03卫星钟差预报结果

3) 最优权组合预报模型结合了灰理论模型与BP神经网络模型的优点，不论在预报前期或者预报后期均能得到满意的预报效果，该组合模型预报值与钟差真值吻合较好. 由表2可知，三种北斗卫星的最优权组合预报精度均优于单一模型.

图2 C13卫星钟差预报结果

图3 C11卫星钟差预报结果

图4 C11卫星钟差预报误差对比图

4) 表2中，MEO卫星预报精度最好，IGSO预报精度次之，GEO预报精度最差. 这三种卫星均搭载铷原子钟，预报精度的不同受到卫星轨道误差的影响.

表2 三种模型预报精度统计对比表

2.2 实验二

由于受到观测条件的限制，为验证最优权组合模型对BDS-3的钟差预报精度情况. 本文在撰写时BDS-3部分卫星并未投入使用，选取BDS-3 MEO C20、IGSO C31两颗具有代表性的卫星. 预报方案如同实验一，采用两种单一模型与最优权组合预测模型作比较，精度统计采用实验一的均方根误差.

1) 在获取数据时，部分BDS-3卫星并未启用. 选择MEO与IGSO各代表卫星. 由表3可知，在最优权组合模型中，BP神经网络的权重系数大于灰色理论模型，又一次证明了BP神经网络模型作为单一模型预测的准确性. 由图5～6可知，灰色模型，BP神经网络模型，最优权组合模型在BDS-3钟差预报中，预测效果如BDS-2卫星.

表3 权重系数

图5 PRN C20预报钟差精度

图6 PRN C31预报钟差精度

2) 由表2、表4中可以看出，不管是最优权组合模型还是单一模型，相同类别的BDS-3卫星钟差预测精度要高于BDS-2卫星钟差预测精度，表明了BDS-3卫星所搭载的原子钟更加稳定.

表4 三种模型预报精度统计表

结合实验一与实验二，最优权组合预测思想是综合BP神经网络与灰色理论预测模型，吸取各单一模型的预测长处与预测不确定性难点以达到取长补短的目的. 利用较好的单一预测模型对较差单一预测模型的修正，在充分顾及单一预测模型优势与不足的同时，达到对预测结果有效的融合. 最优权组合预测钟差方法是对卫星钟差预报的一个大胆尝试，这为组合预测思想应用提供实践探索.

3 结束语

基于组合预测思想，利用最优权组合方法将BP神经网络与灰理论两种模型科学合理结合. 并通过两组实验，包括BDS-2与BDS-3卫星钟差数据，充分证明了最优权组合预测思想的适用性与可行性；同时间接证明了BDS-3所携带的原子钟较BDS-2原子钟更加稳定. 组合预测的思想不仅局限于本文提及的两种单一模型，很多优秀的单一模型都可以通过组合预测的方法来提高钟差预测的精度，这也是本文下一步工作重点.