基于MEABP神经网络的卫星钟差预报.docx

1、基于MEABP神经网络的卫星钟差预报摘要：卫星钟差是影响导航定位精度的重要因素之一，建立高精度的钟差预报模型对高精度定位有重要意义。针对常用模型卫星钟差在短期预报中随时间增加误差积累，以及传统BP神经网络不稳定，容易出现过拟合等问题，本文提出一种基于思维进化算法（MEA）优化的BP神经网络钟差预报模型和算法。首先对原始钟差数据进行一次差处理；然后利用思维进化算法对BP神经网络的初始权值和阈值进行优化，给出该模型进行钟差预报的具体步骤；选用IGS站提供的多天GPS精密钟差产品数据进行试验分析，使用GPS一天中前12 h数据建模，进行2、3、6和12 h的钟差预报。结果表明：利用MEA-BP模型得

2、到的上述4种时段的预报精度分别优于0.36、0.38、0.62和1.56 ns，预报误差曲线变化起伏较小，说明新模型的预报性能优于3种传统模型，新模型在钟差预报短期预报中的实用性及稳定性是较佳的。关键词：卫星钟差一次差思维进化算法BP神经网络钟差预报Prediction of the satellite clock bias based on MEA-BP neural networkAbstract: The satellite clock bias is one of the important factors that affect the accuracy of navigation

3、and positioning, so establishing a high-precision clock bias prediction model is of great significance to high-precision positioning. Aiming at the problem that satellite clock bias error accumulates by common models over time in short-term prediction, and the easy overfitting and instability of the

4、 traditional BP neural network, this paper proposed a model and algorithm of clock bias prediction based on BP neural network optimized by the mind evolutionary algorithm(MEA). First, original clock bias data made once difference to obtain the corresponding once difference sequences. Then, the initi

5、al weights and thresholds of the BP neural network were optimized by the mind evolutionary algorithm, the specific steps of using this model for the clock bias prediction were given. The multi-day GPS precision clock bias product data provided by the IGS station is used for experimental analysis. Th

6、e article used the GPS data for the first 12 h of the day for modeling were listed, and made short-term clock bias prediction within 2, 3, 6 and 12 h. The results showed that the above four periods of prediction precision obtained by using the MEA-BP model were better than 0.36, 0.38, 0.62 and 1.56

7、ns, respectively. The fluctuation of the prediction error curve was small, and the prediction performance of the new model was better than the three traditional models, which showed the new model is better in practicability and stability in the short-term prediction of clock bias.Key words:satellite

8、 clock biasonce differencemind evolutionary algorithm(MEA)BP neural networkclock bias prediction全球导航卫星系统(GNSS)是一种利用无线信号传播提供导航、定位及授时服务(PNT)无线电导航定位系统1-3。由于导航卫星的飞速发展及广泛使用，人们对精度的要求越来越严格，影响导航定位精度的因素较多，其中时间误差的影响不容忽视4。1 ns的时间误差相对应的距离误差为3 dm，严重影响当今导航系统的定位精度。为满足当今用户厘米级定位的需求，高精度的时间系统至关重要5。IGS发布的钟差产品为事后精密钟差文件，

9、精度较高，但获取时间较长，无法满足实时定位的需求。卫星钟差的预报有助于获取卫星自主导航所需要的先验信息及提高实时动态定位精度6-7，为提高实时定位精度，卫星钟差数据的精度和时效性亟待解决4。目前关于卫星钟差预报的模型较多，有多项式模型8、灰色模型9、多项式+周期项模型10、ARMA模型11、Kalman滤波模型12等，还有学者针对单一预报模型的局限性提出了线性组合模型13-14等。由于卫星钟易受外界环境的影响，且自身特征比较复杂，卫星钟差呈现非线性特性，线性预报模型往往难以精确地表示卫星钟差变化15。针对该问题，许多学者利用适用于非线性处理的神经网络在钟差预报方面进行了大量研究。文献16提出E

10、MD-SVM算法进行卫星钟差预测，解决了氢原子钟钟差预测易受非线性和不稳定性影响的问题，预测精度显著提高并不易发散；文献17利用径向基函数神经网络进行卫星钟差预报，得到了较高精度的钟差结果；文献18提出基于一次差的小波神经网络模型，在6、12和24 h的预测中比IGU-P钟差产品的精度分别提高了约13.53%、31.56%和49.46%，且预测稳定性更高；文献19利用ART网络确定ELM网络的隐含层节点数，在长达30 d的预报中，预报精度优于二次多项式模型和灰色模型；文献20通过对四阶谐波模型残差的分析得到了钟差数据的额外周期数，并实现了8阶谐波函数建模，预测得到的钟差精度高于IGU产品精度；

11、文献21在周期项模型的基础上添加循环项，并利用BP神经网络补偿非线性系统误差，提高了超快速钟差预报的精度；还有学者针对BDS数据建立超快速钟差预报模型，均得到了较好的钟差预报精度22-24。BP神经网络在1986年由文献25提出后使用广泛，具有良好的自适应性、抗差性和联想记忆功能，但在使过程中容易陷入局部最优，且收敛速度较慢，影响最终收敛精度。针对卫星钟差的特性和BP神经网络的局限性，利用思维进化算法(MEA)良好的全局搜索能力和较强的收敛性，本文引入思维进化算法来选择BP神经网络所需要的最优权值和阈值，将两者结合建立了基于思维进化算法优化BP神经网络(MEA-BP)的钟差预报模型，计算中对卫

12、星钟差数据作一次差处理，利用一次差后的数据进行建模，预报得到新的一次差值，再将一次差值还原得到预测钟差数据。大量算例表明，按照本文提出的网络结构得到了较高精度、稳定性的MEA-BP钟差预报算法。1 钟差一次差处理及常规预报模型1.1 钟差一次差处理同一卫星的钟差数据在前后相邻两个历元的变化量较小，从整体看呈现线性趋势。BP神经网络对于非线性数据有良好的适用性，对原始钟差数据序列不够敏感，对未作处理的数据进行建模预报时严重影响预报精度，故本文先对原始钟差数据作一次差处理，将相位数据转为频率数据，处理后的数据具有很好的非线性特征，适用于神经网络的建模与预报。设X=x(1),x(2), ,x(n)为

13、一组n维度的卫星钟差序列，其中x(i),i=1, 2, ,n为不同历元的钟差数据，对相邻历元的钟差数据作差，得到新的一次差数据序列为(1)式中，x(i)=x(i+1)x(i)。用神经网络对构成的数据序列X进行建模，对之后m历元的卫星钟差进行预报：对该n1个数据进行建模，预报后m个数据序列X=x(n), x(n+1), , x(m+n1)，利用预报的一次差值与建模最后一个钟差数据x(n)恢复待求历元的钟差，即(2)式中，nkn+m1。1.2 3种预报模型1.2.1 二次多项式模型(QP模型)8卫星钟差预报的二次多项式模型为(3)式中，i=1, 2, ,n；ti时刻的卫星钟差为xi；t0为卫星钟参

14、考时刻；a0、a1、a2为参考时刻的卫星相位、频偏和频漂；为模型随机误差。利用最小二乘原理即可求得参数估值，将其代入式(3)中可得预报历元的钟差。但多项式模型并不适合非平稳钟差数据，并且随着预报时间的增加，多项式模型的预报误差随之加大，预报精度会显著降低。1.2.2 二次多项式+周期项模型(SA-QP模型)10在二次多项式模型的基础上，附加周期项的钟差预报模型为(4)式中，Ak、fk、k分别对应周期项的振幅、频率和相位；p为主周期的个数；其余参数同式(3)。与二次多项式模型计算相似，利用最小二乘原理求得参数估值，代入式(4)中可求得预报历元的钟差值。本文针对卫星钟差数据的特性，将主周期个数设为

15、4，分别为12、6、4和3 h25。1.2.3 灰色模型(GM(1, 1)9灰色模型预测钟差的表达式为(5)式中，k=1, 2, ,n。参数估值可利用最小二乘原理求解，将得到的参数代入式(5)可得需要预测历元的卫星钟差值。但灰色模型易受指数参数的影响，使用最小二乘求解参数时有时会陷入局部最优，得到的卫星钟差精度较差。2 MEA-BP算法基本原理2.1 BP神经网络算法神经网络在非线性时间序列的拟合和预测方面效果显著，其中BP神经网络使用最为广泛。BP神经网络是一种误差反向传播的多层前馈神经网络，包括输入层、隐含层及输出层3部分，学习能力和记忆功能较强，属于自适应的非线性动力学系统25。当输出值

16、不是期望值时，该网络会根据误差反向传播调整权值和阈值，进而不断逼近期望输出值。其拓扑结构如图 1所示。图 1BP神经网络拓扑图Fig. 1Topologic diagram of BP neural network图选项BP神经网络的隐含层可为多层，为减少运算量、提高计算效率，本文将隐含层设为1层。在训练过程中，每层节点的输入层计算为(6)激活函数常采用sigmoid函数，即(7)该函数连续可导，可将数据从, +投影到区间0, 1中。则输入层计算后，对应的输出层计算为(8)式中，Outn,j为BP神经网络得到的输出值。BP神经网络通过误差后向传播进行权值更新，更新公式为(9)式中，t为迭代次数

17、；(t)为期望输出与实际输出之差；为神经网络的学习效率；y(t)为神经元的输出值。2.2 MEA算法思维进化算法(mind evolutionary algorithm，MEA)由文献26于1998年针对遗传算法的局限性提出的一种进化算法。MEA算法保留了遗传算法的“群体”、“个体”、“环境”和“进化”等思想，和遗传算法的“交叉”和“变异”不同的是，提出了“趋同”和“异化”等新的概念。利用MEA算法优化神经网络比使用遗传算法优化有更快的训练速度，大大缩短了神经网络的训练时间，实用性更广27-28。首先利用趋同操作对子群体中的个体进行优选，然后成熟的各子群体再通过异化操作进行全局竞争，大大提高了

18、优化的效率。MEA算法的结构如图 2所示。图 2MEA算法结构Fig. 2MEA algorithm图选项MEA算法是通过迭代进行优化的学习方法，其基本思路与步骤如下：(1) 在解空间内随机产生一定数量的个体，计算训练个体均方误差的倒数得到各个体的分数，从中选取得分相对较高的若干个体得到优胜个体和临时个体。(2) 以选出的个体为中心，在各个体周围产生新的个体，分别得到优胜子群体和临时子群体。(3) 对各个子群体分别进行趋同操作，即在各子群体范围内，个体进行竞争选出优胜者。当胜者不再变化时，表示该子群体已经成熟，将优胜者的得分作为该子群体的得分，并张贴在全局公告板上。(4) 各子群体完成趋同操作

19、后，在全局范围内对所有子群体进行异化操作，通过比较各子群体间的得分大小，完成优胜子群体和临时子群体的替换、废弃和子群体中个体释放的过程，最后得到全局范围内的最优个体及其得分。(5) 子群体中释放的个体在解空间中产生新的临时子群体，保持临时子群体总数不发生变化。判断是否满足精度要求，若不满足继续重复步骤(3)和步骤(4)操作，直至迭代结束或最优个体分数不发生变化，输出最优个体。2.3 MEA优化BP的钟差建模与预报算法BP神经网络是通过反复训练来调整网络中的各权值和阈值，使得输出值逼近期望值。但是初始权值和阈值的选择严重影响BP神经网络的收敛和精度，并且得到的最终结果容易陷入局部最优。针对BP神

20、经网络使用随机参数可能会导致算法结果异常的情况，本文使用MEA算法得到了神经网络中较好的初始化参数，可避免BP神经网络陷入局部最小值，大大提高卫星钟差预测的精度。设有一组卫星钟差数据为x1,x2, ,xn，使用该n个时刻的数据来进行建模，预报以后时刻的钟差值。通常输出层的个数要和输出数据类型的个数相同，本文输出值为钟差，故将输出节点设为1，利用神经网络建立x1,x2, ,xN与xn+1的映射关系(N为输入节点数)，采用滑动窗口的思想，在保证样本数目不变的前提下，不断使用新的预测数据来取代之前的已知数据，实现多历元的卫星钟差预报。利用MEA-BP模型进行卫星钟差预报的具体步骤如下：(1) 对卫星

21、钟差原始数据进行中位数探测，剔除数据中存在的粗差和钟跳。(2) 将钟差建模数据分为输入部分和输出部分，使用归一化方法将数据归化到区间1, 1，并根据输入值和输出值选择BP神经网络的拓扑结构。(3) 根据确定好的网络拓扑结构选择编码长度，定义迭代次数、种群大小、优胜子种群个数和临时子种群个数，求得子群体的大小。(4) 将待优化的权值和阈值根据步骤(3)设置的参数进行编码，并随机产生初始种群。(5) 根据式(10)和式(11)，先计算得到均方误差E，然后对其取倒数得到各种群和个体的得分函数score(10)(11)式中，为迭代后的输出值；yi为样本值；m为样本总数。(6) 将当前的得分与历史最优得

22、分比较，通过不断迭代，更新优胜种群及其最优得分。(7) 判断每次得到的最优个体分数是否发生变化或是否已达到最大迭代次数，若满足精度，输出最优个体，优化过程终止，否则重复步骤(5)和步骤(6)。(8) 根据步骤(7)获得的最优种群恢复BP神经网络需要的初始权值和阈值，代入神经网络进行训练和预报，对预报值进行反归一化运算，得到最终需要的卫星钟差值。算法流程如图 3所示。图 3MEA算法优化BP神经网络流程Fig. 3BP neural network model optimized by MEA algorithm图选项使用n维建模数据完成神经网络训练后，依据神经网络结构，输入节点数为N，输出节点

23、数为1，采用滑动窗口的思想，保证样本数目不变，以此类推，实现多步预报，具体见表 1。表 1神经网络多步预报Tab. 1Multi-step prediction of neural networkinputoutputxnN+1, ,xnxn+1xnN+2, ,xn+1xn+2xnN+m, ,xn+m1xn+m表选项得到一次差分序列的预报结果后，通过还原即可得到钟差预测值。3 算例分析为了验证本文算法的实用性，选用IGS站提供的多天GPS精密钟差产品数据进行试验分析。由于篇幅有限，本文只列出GPS系统week 1919第一天(对应时间为2016年10月16日)的卫星钟差数据，采样间隔为30 s

24、。本文选用5种不同原子钟类型的卫星PRN02、PRN10、PRN24、PRN25和PRN32，原子钟类型分别为：Block R Rb、Block A Cs、Block F Cs、Block F Rb和Block A Rb，使用前12 h的数据进行建模，预报接下来不同时刻的卫星钟差。将预报的钟差值和IGS发布的真实值进行对比，用公式(12)得到数据的均方根误差和误差的极差差值、均值，来评价分析预报结果的准确度(12)(13)(14)式中，为得到的预测钟差值；xi为IGS发布的真实钟差值；n为预测的钟差个数；xmax为误差的最大值；xmin为误差的最小值。建模采用12 h数据，卫星钟差数据采样间隔

25、为30 s，共1440个数据，所以设置MEA算法种群大小为1440，优胜子种群和临时子种群个数均为5，为使算法充分寻找最优个体，迭代次数设为200。文献17在预报不同时长的钟差时使用了不同的输入层和输出层，但实际进行钟差预报时选取较为复杂，并且网络模型的时效性受其结构的影响。本文将输入层设为1，由于隐含层的设置目前还未有可遵循的理论依据，只能依靠经验值进行确定，本文参考Kolmogorov定理29选择隐含层单元数(15)式中，O为隐含层个数；N为输入层个数。综上所述，本文将输入层节点设为1，输出层节点设为1，经计算得到隐含层节点为3，故最终确定的BP神经网络结构为1-3-1。3.1 试验1首先

26、验证一次差处理的优势。文献30使用一次差法对卫星数据进行处理，比较分析了一次差前后常用模型的预报精度，说明一次差后会造成预报过程中的误差累积，严重影响二次多项式模型的精度；而灰色模型本身就需要对数据进行累加和累减，所以一次差处理对该模型的影响难以给出定性结论。对神经网络而言，数据的非线性特性越明显，理论上预报精度越高。随机选取两颗原子钟不同的卫星进行预报。本文选取2号卫星(Rb)和10号卫星(Cs)举例说明，使用该天前12 h的钟差数据进行训练，预报接下来1 h的钟差，现使用未经优化的BP神经网络模型分别预报20次，对比其RMS值变化情况，如图 4所示。图 4一次差前后BP预报20次RMS的变

27、化Fig. 4Variation of RMS using BP to predict clock bias before and after once difference about 20 times图选项由图 4可知，对不作一次差处理的数据进行建模预报时，2号卫星和10号卫星预报20次的RMS值变化较大，且精度不高；而经过一次差处理后，两颗卫星的预报精度有了显著提升，并且多次预报的精度相当，有较强的稳定性。表明在同样使用BP神经网络模型进行卫星钟差预报的条件下，不管卫星原子钟是Rb钟还是Cs钟，对原始钟差进行一次差处理后，预报的精度和稳定性均优于一次差处理前。综上所述，一次差处理方法用于

28、BP神经网络，可大幅度提高卫星钟差预报的稳定性和精度，并且适用于本文使用的网络结构。接下来试验将在一次差基础上进行研究分析。3.2 试验2对比BP模型和MEA-BP模型对卫星钟差的预报性能。选取2号卫星(其他卫星均可)举例说明，为充分对比预报效果，使用该天前12 h的钟差数据分别独立预报10次接下来3、6和12 h的钟差值，两种模型预报的RMS值变化情况如图 5所示。图 5MEA-BP和BP预报10次RMS的变化Fig. 5Variation of RMS using MEA-BP and BP to predict SCB about 10 times图选项在本试验中，利用MEA-BP和BP

29、模型分别对2号卫星进行不同时段的预测，由图 5可知，MEA-BP和BP模型利用选取的网络结构，多次预报精度并无明显波动，证明了该网络结构的可行性；MEA-BP模型在预报时间增加时，比BP模型有更好的预报性能。在3个时段的钟差预测中，MEA-BP模型的预报精度均优于BP模型，说明通过MEA算法可以改善BP神经网络的初始权值和阈值，避免BP模型陷入局部最优，有效提高预报结果的精度。表明了MEA-BP模型用于钟差预报的可行性，是一种相对稳定的卫星钟差预报算法。3.3 试验3MEA-BP模型与几种常用模型的预报性能对比。对2号卫星、10号卫星、24号卫星、25号卫星和32号卫星一天中前12 h的数据进

30、行建模，分别预报2 h(240历元)、3 h(360历元)、6 h(720历元)和12 h(1440历元)的钟差。下面给出5颗卫星在以上4个时段的预报中，MEA-BP模型与二次多项式模型、二次多项式+周期项模型和灰色模型的具体精度对比情况，见表 2表 5。表 1给出了分步预报的情况，12 h的预报误差包含2、3及6 h，故给出预报12 h时5颗卫星均方根误差和误差的极差差值的条形统计图，及MEA-BP模型与3种常用模型的预测误差趋势走向图，如图 6图 7所示。表 25颗卫星2 h预报结果统计Tab. 2The statistics of prediction results of 2 hour

31、 for five satellitesnsPRNstatisticsQPSA-QPGM(1, 1)MEA-BPPRN02range1.3931.2951.0280.739mean-0.7670.398-0.2550.051RMS0.8230.5110.3240.150PRN10range0.3841.0800.5190.496mean0.680-0.363-0.448-0.217RMS0.6900.4760.4680.253PRN24range2.1825.0861.5811.700mean1.373-2.443-0.172-0.008RMS1.4602.8910.4090.358PRN25range0.5950.7420.3490.305mean0.137-0.633-0.81