一种基于改进型bp神经网络的氢原子钟钟差预测方法
【技术领域】
[0001] 本发明涉及原子钟数据处理领域，特别涉及一种基于改进型BP(Back Propagation)神经网络的氢原子钟钟差预测方法。
【背景技术】
[0002] 钟差预测是指运用一种或多种算法通过对原子钟的历史钟差数据进行分析学习 后，算出下一刻或下一段时间内的钟差数据。原子钟钟差预测是原子钟权重分配、预测原子 钟的运行状态和对原子钟进行预调节、系统驾驭的重要依据。因此，精准的原子钟钟差预测 对提高原子时时标精度和原子钟系统驾驭精度有着重要作用。
[0003] 目前，常用的钟差预测算法有一次线性回归预测算法、基于支持向量机（SVM)预 测算法、小波神经网络预测算法、灰色模型和ARIMA模型的钟差预测算法、切比雪夫多项式 模型预测算法等。但是，针对钟差的预测算法主要运用在卫星钟（铷原子钟）当中，用于优 化导航电文的钟差参数、满足实时动态精密单点定位的需求和提供卫星自主导航所需的先 验信息。而针对地面氢原子钟钟差的预测算法较少，主要有一次线性回归算法和基于支持 向量机（SVM)预测算法。其中一次线性回归预测算法主要是拟合出钟差的变化趋势，算法 简单，容易实现，但预测精度比较低，SVM预测算法的预测精度较一次线性回归算法预测精 度要高。
【发明内容】
[0004] 为进一步提高对原子钟钟差数据的预测精度，为提高原子时时标精度和原子钟驾 驭精度提供依据，本文根据氢原子钟钟差特性、地面原子钟钟差数据与卫星钟钟差数据以 及与日长数据的相似性，在卫星钟差小波神经网络预测算法和小波神经网络对日长预报算 法中BP神经网络算法思想基础上，针对BP神经网络算法的不足，运用SVM算法中的预测惩 罚模型，提出了基于改进BP神经网络算法对氢原子钟钟差数据进行预测方法。
[0005] 改进型BP神经网络算法是在BP神经网络算法的基础之上，针对BP神经网络算法 的不足之处，运用SVM算法中的预测惩罚思想提出了改进方案，主要从以下两个方面进行 了算法的改进：
[0006] 1、）针对BP神经网络权值收敛到某个值，并不保证其为误差平面的全局最小值的 问题，本文通过增加惩罚性措施来提高预测精度，通过上一步预测误差S(n-l)即实测值 y(n-l)预测值与yjn-l)之间的差值，来修正下一步预测值。2、）针对预测数据中出现大 幅度的变化的问题，在本文当中，使用限制预测数据变化幅度的措施来提高预测数据的稳 定性和预测精度。
[0007] 本发明的特征如下：
[0008] 步骤1 :给系统参数WpVp0p〇^及学习率n赋初值并选定传递函数形式，一
【主权项】
1. 一种基于改进型BP神经网络的氢原子钟钟差预测方法，其特征在于：该方法包括下 述流程， 步骤1 :给系统参数Wp Vp Θ i、CDi及学习率η赋初值并选定传递函数形式，一般使 A=匕可简写为f，通常选_
其中Wiu为对应输入的连接权值，W %为对应输出单元的连接权值，f 输入单元到隐 层单元的传递函数，f2为隐层单元到输出单元的传递函数，Θ i为隐层单元的阈值即隐层单 元输入的最小氢钟钟差值，Oi为输出单元的阈值即预测单元输出的最小氢钟钟差值； 步骤2 :输入训练样本数据，根据如下公式求出各个神经元的实际输出： Γ = f (Wx- Θ ) (1) 其中Γ为每个神经元的实际输出的氢钟钟差值，f为神经元的激励函数，X为对应神经 元输入的氢钟钟差值，W为对应的连接权值，Θ为对应神经元的阈值； 步骤3 :根据给定的误差公式，求出系统性能误差：
(2) 其中outputs为输出钟差值的单元个数，tk为输出单元的钟差期望值，Γ k为输出单元 的钟差实际输出值；判断钟差的预测误差或训练次数是否达到预设值，如果达到退出训练， 否则进行下一步； 步骤4 :计算输出层每个输出单元的钟差误差项： 〇k= Γ k(l-Tk) (tk-Tk) (3) 步骤5 :计算隐层中每个单元的钟差误差项： '
(4) 步骤6 :根据设定的学习率η，调整各层的连接权值： ff(n+l) =ff(n)+nx〇k (5) 步骤7:返回第二步，重复计算，不断调整各连接权值，直到输出的钟差误差到达预设 值为止； 步骤8 :使用训练好的网络进行钟差预测； 步骤9:计算钟差预测误差 δ (n) = y (n) -yr (n) (6) 步骤10 :修正钟差预测值 yr(n) = yr(n) + y δ (n-1) (7) 其中μ为惩罚系数， 步骤11 :判断钟差预测值变化幅度，设置限制幅度ξ， 如果 yr(n)-yr(n-l) | ^ ξ (8) 则最终的钟差预测值^ (n)修改为：
(9) 否则，最终的钟差预测值h(n)保持不变； 步骤12 :从第二步开始重复以上步骤，直到完成预测； 通过以上步骤，改进型BP神经网络即完成了针对氢钟钟差数据的预测。
2.根据权利要求1所述的一种基于改进型BP神经网络的氢原子钟钟差预测方法，其特 征在于：首先，根据氢钟钟差数据，确定钟差数据正常波动范围，确定限制幅度值和预测惩 罚系数，并对改进型BP神经网络系统参数赋初值；然后，输入氢钟钟差数据的历史数据，对 钟差数据进行学习得到系统函数；最后，输入要预测点，根据系统函数输出预测值，根据限 制幅度值和预测惩罚系数对预测值进行修正，完成预测； 具体而言包括下述步骤， 步骤1，根据氢钟钟差特性分析钟差的正常波动范围，确定限制幅度值为〇. 00125ns和 预测惩罚系数0.4; 步骤2 :给系统参数Wp Vp Θ p (Di及权值学习率η赋初值并选定传递函数形式，一 般使A=丨2可简写为f，通常选择
步骤3 :输入训练样本数据，根据如下公式求出各个神经元的实际输出： Γ = f (Wx- Θ ) (1) 其中Γ为每个神经元的实际输出，f为神经元的激励函数，X为对应神经元的输入，W 为对应的连接权值，Θ为对应神经元的阈值； 步骤4 :根据给定的误差公式，求出系统性能误差：
其中outputs为输出单元的个数，tk为输出单元的期望输出值，Γ k为输出单元的实际 输出值；判断系统误差或训练次数是否达到预设值，如果达到退出训练，否则进行下一步； 步骤5 :计算输出层每个输出单元的误差项： 〇k= Γ k(l-Tk) (tk-Tk) (3) 步骤6 :计算隐层中每个单元的误差项：
步骤7:根据设定的学习 η，调整各层的连接权值： ff(n+l) =ff(n)+nx〇k (5) 步骤8 :返回第二步，重复计算，不断调整各连接权值，直到输出误差到达预设值为止； 步骤9 :使用训练好的网络进行预测； 步骤10:计算预测误差 δ (n) = y (n) -yr (η) (6) 步骤11 :修正预测值 yr(n) = yr(n) + y δ (η-1) (7) 其中μ =0.4为惩罚系数， 步骤12 :判断预测值变化幅度，设置限制幅度ξ = 0. 00125ns， 如果 yr(n)-yr(n-l) | ^ ξ (8) 则最终预测值(η)修改为：
否则，最终预测值^ (η)保持不变； 步骤13 :从第二步开始重复以上步骤，直到完成预测。
【专利摘要】一种基于改进型BP神经网络的氢原子钟钟差预测方法，本方法根据氢原子钟钟差特性、地面原子钟钟差数据与卫星钟钟差数据以及与日长数据的相似性，在卫星钟差小波神经网络预测算法和小波神经网络对日长预报算法中BP神经网络算法思想基础上，针对BP神经网络算法的不足，运用SVM算法中的预测惩罚模型，提出了基于改进BP神经网络算法对氢原子钟钟差数据进行预测方法。通过增加惩罚性措施来提高预测精度，来修正下一步预测值；使用限制预测数据变化幅度的措施来提高预测数据的稳定性和预测精度。该方法与现行的SVM预测算法和线性回归预测算法相比，其氢钟中差预测精度有明显提高，为提高氢原子钟的原子时标和驾驭精度提供了更高的依据。
【IPC分类】G06F19-00
【公开号】CN104679989
【申请号】CN201510041024
【发明人】朱江淼, 宋文峰
【申请人】北京工业大学
【公开日】2015年6月3日
【申请日】2015年1月27日