CN111626406A - 一种基于一次差分的bds钟差预测方法 - Google Patents
一种基于一次差分的bds钟差预测方法 Download PDFInfo
- Publication number
- CN111626406A CN111626406A CN202010426241.8A CN202010426241A CN111626406A CN 111626406 A CN111626406 A CN 111626406A CN 202010426241 A CN202010426241 A CN 202010426241A CN 111626406 A CN111626406 A CN 111626406A
- Authority
- CN
- China
- Prior art keywords
- difference
- data
- phase data
- clock
- clock error
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18519—Operations control, administration or maintenance
Abstract
本发明公开了一种基于一次差分的BDS钟差预测方法,该方法为:获取BDS精密星历中t个历元的精密钟差数据(相位数据),将其转化为钟差频率数据,采用MAD方法探测钟差频率数据中的粗差所对应的历元位置,从而得到钟差相位数据中的粗差所对应的历元位置,对其进行剔除,然后用该历元时刻的前一个相位数据和后一个相位数据的平均值代替该历元时刻的相位数据,得到平滑的钟差相位数据序列;对平滑钟差相位数据进行一次差分处理,用差分数据作为建模数据,采用径向基函数神经网络预测差分数据,最后对差分数据进行还原得到预测的钟差相位数据。该方法使得波动曲线更加平稳,解决了常规钟差预测中出现的数据发散问题。
Description
技术领域
本发明属于卫星钟差预测的技术领域,尤其涉及一种基于一次差分的BDS钟差预测方法。
背景技术
北斗卫星导航系统(Beidou Navigation Satellite System,BDS)是由我国自主研发的导航定位系统,是国防战略中的重要发展焦点,也是国际地位和综合国力的象征。北斗导航系统的构建共分为三个阶段,基于北斗一号的北斗二号的建设基础,北斗三号将以全球覆盖,性能优化为主要目标,同时更好的实现全天候、全天时的导航定位、短报文通讯、海上搜救等服务,在国防、军事、经济等领域发挥着重大作用。
北斗卫星导航系统的基础是时间测量,精确位置的测量也是在获取精确时间的基础上完成的,但北斗精密钟差的获取具有较为严重的滞后性,难以实时获取,这在一定程度上限制了BDS的基础服务,因此选择合适的方法和模型对钟差进行预测具有重大意义和广阔的应用价值。
发明内容
基于以上现有技术的不足,本发明所解决的技术问题在于提供一种基于一次差分的BDS钟差预测方法,采用MAD探测法和一次差分法完成数据预处理,再利用径向基函数(Radial Basis Function,RBF)神经网络建模预测,最后进行数据还原处理得到预测钟差数据。
为了解决上述技术问题,本发明通过以下技术方案来实现:
本发明提供一种基于一次差分的BDS钟差预测方法,包括以下步骤:
步骤1:获取BDS卫星系统中ti时刻和ti+1时刻的精密钟差相位数据li和li+1,其中,t为任意时刻,i∈[1,n-1],n为历元个数;
步骤2:根据ti时刻和ti+1时刻所对应的相位数据以及采样间隔Δt确定与其对应的频率数据yi;
步骤3:采用MAD探测法检验频率数据序列{yi},对粗差值yp用其相邻的数值yp-1和yp+1的平均值代替,得到平滑频率数据序列{zi},然后通过{zi}的波动判定粗差剔除效果,其中p为粗差所在位置数,p∈[1,n-1];
MAD探测法的公式如下所示:
|yi|>a+b·MAD
其中,median为中位数运算符,a为频率数据序列的中位数,b为根据实际情况所确定的正整数,MAD为频率序列偏差绝对值的中位数;
步骤4:在得出平滑频率数据序列{zi}的前提下,记录粗差值yp所对应的相位数据lp,对lp用其相邻的数值lp-1和lp+1的平均值代替,得到平滑相位数据序列{xi};
步骤5:对平滑相位数据序列{xi}做一次差分得到差分数据序列{Δxi-1};
ΔXi=Xi+1-Xi
优选的,在步骤2中,确定频率数据yi的计算公式如下:
其中,Δt为采样间隔,l为相位数据,li和li+1是ti时刻和ti+1时刻所对应的相位数据,yi是ti时刻所对应的频率数据。
进一步的,在步骤4中,确定相位数据lp的计算公式如下所示:
其中,lp-1、lp、lp+1分别为不同时刻的钟差相位数据。
可选的,所述步骤6包括以下步骤:
步骤6.1:解算RBF基函数中心,在进行预测中RBF核函数选择高斯函数,表达式如下所示:
i=1,2,L,k
其中,h为网络隐含层节点数,||xi-ch||为欧式距离范数,ch为激活函数中心,σ2为方差,i为隐含层节点个数;
步骤6.2:解算隐含层与输出层之间的权值,计算公式如下所示:
其中,w为隐含层与输出层之间的权值,k为隐含层节点个数;
步骤6.3:解算RBF网络输出:
进一步的,所述步骤8包括以下步骤:
其中,RMSj为均方根误差,q为历元个数,q∈[1,k];
其中,MEj为残差均值。
由上,本发明的基于一次差分的BDS钟差预测方法采用MAD探测法和均值法进行数据预处理,然后对平滑相位数据进行一次差分,从而有效的减少噪声,利用RBF神经网络预测钟差的差分相位数据,最后做数据还原得到最终的预测钟差。本发明提高了钟差预测精度,解决了常规方法中可能导致误差过大的情况。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂,以下结合优选实施例,并配合附图,详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例的附图作简单地介绍。
图1为本发明的基于一次差分的BDS钟差预测方法的流程图;
图2是实验得到的基于一次差分的BDS钟差预测方法的残差图;
其中,(a)是实验一得到的基于一次差分的BDS钟差预测方法的残差图;
(b)是实验二得到的基于一次差分的BDS钟差预测方法的残差图;
(c)是实验三得到的基于一次差分的BDS钟差预测方法的残差图;
(d)是实验四得到的基于一次差分的BDS钟差预测方法的残差图;
(e)是实验五得到的基于一次差分的BDS钟差预测方法的残差图;
图3是实验得到的基于一次差分的BDS钟差预测方法的均方根误差图;
其中,(a)是实验一得到的基于一次差分的BDS钟差预测方法的均方根误差图;
(b)是实验二得到的基于一次差分的BDS钟差预测方法的均方根误差图;
(c)是实验三得到的基于一次差分的BDS钟差预测方法的均方根误差图;
(d)是实验四得到的基于一次差分的BDS钟差预测方法的均方根误差图;
(e)是实验五得到的基于一次差分的BDS钟差预测方法的均方根误差图;
图4是实验得到的基于一次差分的BDS钟差预测方法的残差均值图;
其中,(a)是实验一得到的基于一次差分的BDS钟差预测方法的残差均值图;
(b)是实验二得到的基于一次差分的BDS钟差预测方法的残差均值图;
(c)是实验三得到的基于一次差分的BDS钟差预测方法的残差均值图;
(d)是实验四得到的基于一次差分的BDS钟差预测方法的残差均值图;
(e)是实验五得到的基于一次差分的BDS钟差预测方法的残差均值图。
具体实施方式
下面结合附图详细说明本发明的具体实施方式,其作为本说明书的一部分,通过实施例来说明本发明的原理,本发明的其他方面、特征及其优点通过该详细说明将会变得一目了然。在所参照的附图中,不同的图中相同或相似的部件使用相同的附图标号来表示。
如图1所示,本发明的基于一次差分的BDS钟差预测方法包括以下步骤:
步骤1:获取BDS卫星系统中ti时刻和ti+1时刻的精密钟差相位数据li和li+1,其中,t为任意时刻,i∈[1,n-1],n为历元个数。
步骤2:根据ti时刻和ti+1时刻所对应的相位数据以及采样间隔Δt确定与其对应的频率数据yi。
本实施方式中,频率数据yi的计算公式如式(1)所示:
其中,li和li+1是ti时刻和ti+1时刻所对应的相位数据,yi是ti时刻所对应的频率数据,Δt是采样间隔。
步骤3:采用MAD(中位数)探测法检验频率数据序列{yi},对粗差值yp用其相邻的数值yp-1和yp+1的平均值代替,从而得到平滑频率数据序列{zi},然后通过{zi}的波动判定粗差剔除效果,其中p为粗差所在位置数,p∈[1,n-1];
本实施方式中,MAD探测公式如式(2)和(3)所示:
|yi|>a+b·MAD (2)
其中,median为中位数运算符,a为频率数据序列的中位数,b为根据实际情况所确定的正整数,MAD为频率序列偏差绝对值的中位数。
由公式(2)和(3)可知,将频率数据与序列中位数a加上MAD的若干倍相比较,当频率数据的绝对值小时,该频率数据为正常数据;当频率数据的绝对值大时,则该频率数据被定义为粗差。
步骤4:在得出平滑频率数据序列{zi}的前提下,记录粗差值yp所对应的相位数据lp,对lp用其相邻的数值lp-1和lp+1的平均值代替,从而得到平滑相位数据序列{xi}。
本实施方式中,确定相位数据lp的计算公式如式(4)所示:
其中,lp-1、lp、lp+1分别对应tp-1、tp、tp+1时刻的相位数据。
步骤5:对平滑相位数据序列{xi}做一次差分得到差分数据序列{Δxi-1};
ΔXi=Xi+1-Xi (5)
将式(5)改写为矩阵形式如式(6)所示:
步骤6.1:解算RBF基函数中心,在进行钟差预测中RBF核函数选择高斯函数,表达式如下所示:
其中,h为网络隐含层节点数,||xi-ch||为欧式距离范数,ch为激活函数中心,σ2为方差,i为隐含层节点个数;
步骤6.2:解算隐含层与输出层之间的权值,计算公式如下所示:
其中,w为隐含层与输出层之间的权值,k为隐含层节点个数。
步骤6.3:解算RBF网络输出:
其中,RMSj为均方根误差,q为历元个数,q∈[1,k]。
其中,MEj为残差均值,q为历元个数,q∈[1,k]。
本发明实施方式中,采用北斗精密卫星钟差作为实验数据,具体卫星号为C26、C27、C28、C29、C30,数据涉及时长24h,采样间隔为5min,以前12h的数据作为建模数据,后12h的数据作为预测数据的验证数据,共得到5组不同卫星的预测误差值,如图2~4所示。
本发明获取BDS精密星历中t个历元的精密钟差数据(相位数据),将其转化为钟差频率数据,采用MAD方法探测钟差频率数据中的粗差所对应的历元位置,从而得到钟差相位数据中的粗差所对应的历元位置,对其进行剔除,然后用该历元时刻的前一个相位数据和后一个相位数据的平均值代替该历元时刻的相位数据,得到平滑的钟差相位数据序列;对平滑钟差相位数据进行一次差分处理,用差分数据作为建模数据,采用径向基函数神经网络预测差分数据,最后对差分数据进行还原得到预测的钟差相位数据。该方法使得波动曲线更加平稳,解决了常规钟差预测中出现的数据发散问题。
以上所述是本发明的优选实施方式而已,当然不能以此来限定本发明之权利范围,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和变动,这些改进和变动也视为本发明的保护范围。
Claims (5)
1.一种基于一次差分的BDS钟差预测方法,其特征在于,包括以下步骤:
步骤1:获取BDS卫星系统中ti时刻和ti+1时刻的精密钟差相位数据li和li+1,其中,t为任意时刻,i∈[1,n-1],n为历元个数;
步骤2:根据ti时刻和ti+1时刻所对应的相位数据以及采样间隔Δt确定与其对应的频率数据yi;
步骤3:采用MAD探测法检验频率数据序列{yi},对粗差值yp用其相邻的数值yp-1和yp+1的平均值代替,得到平滑频率数据序列{zi},然后通过{zi}的波动判定粗差剔除效果,其中p为粗差所在位置数,p∈[1,n-1];
MAD探测法的公式如下所示:
|yi|>a+b·MAD
其中,median为中位数运算符,a为频率数据序列的中位数,b为根据实际情况所确定的正整数,MAD为频率序列偏差绝对值的中位数;
步骤4:在得出平滑频率数据序列{zi}的前提下,记录粗差值yp所对应的相位数据lp,对lp用其相邻的数值lp-1和lp+1的平均值代替,得到平滑相位数据序列{xi};
步骤5:对平滑相位数据序列{xi}做一次差分得到差分数据序列{Δxi-1};
ΔXi=Xi+1-Xi
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010426241.8A CN111626406A (zh) | 2020-05-19 | 2020-05-19 | 一种基于一次差分的bds钟差预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010426241.8A CN111626406A (zh) | 2020-05-19 | 2020-05-19 | 一种基于一次差分的bds钟差预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111626406A true CN111626406A (zh) | 2020-09-04 |
Family
ID=72258930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010426241.8A Pending CN111626406A (zh) | 2020-05-19 | 2020-05-19 | 一种基于一次差分的bds钟差预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111626406A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327340A (zh) * | 2021-01-06 | 2021-02-05 | 腾讯科技(深圳)有限公司 | 终端定位精度评估方法、装置、设备以及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044747A (zh) * | 2015-08-27 | 2015-11-11 | 交通信息通信技术研究发展中心 | 一种基于多星共视和滤波的时间同步装置及其方法 |
CN108490463A (zh) * | 2018-02-09 | 2018-09-04 | 东南大学 | 一种北斗卫星频间钟差偏差估计与建模方法 |
-
2020
- 2020-05-19 CN CN202010426241.8A patent/CN111626406A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105044747A (zh) * | 2015-08-27 | 2015-11-11 | 交通信息通信技术研究发展中心 | 一种基于多星共视和滤波的时间同步装置及其方法 |
CN108490463A (zh) * | 2018-02-09 | 2018-09-04 | 东南大学 | 一种北斗卫星频间钟差偏差估计与建模方法 |
Non-Patent Citations (8)
Title |
---|
张景元: ""基于神经网络的卫星钟差预报研究"", 《计算机工程与设计》 * |
杨嘉明: ""基于一阶差分LM-BP神经网络的卫星钟差预报模型"", 《中国卫星导航学术年会 2019》 * |
王宇: ""BDS星载原子钟钟差数据预处理与钟性能分析研究"", 《中国优秀硕士论文全文数据库》 * |
王宇谱: ""GNSS星载原子钟性能分析与卫星钟差建模预报研究"", 《中国优秀博士论文全文数据库》 * |
王宇谱: ""基于修正一次差的灰色模型及其在BDS卫星钟差预报中的应用"", 《中国卫星导航年会》 * |
许玉德: "《城市轨道交通轨道结构检测与修理》", 31 December 2018 * |
贾宇明: "《综合电子系统技术教育部重点实验室2009学术年会暨电子科技大学电子科学技术研究院第五届学术会议论文集》", 31 December 2009 * |
魏道坤: "卫星钟差预报模型研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327340A (zh) * | 2021-01-06 | 2021-02-05 | 腾讯科技(深圳)有限公司 | 终端定位精度评估方法、装置、设备以及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xiang et al. | A rainfall‐runoff model with LSTM‐based sequence‐to‐sequence learning | |
CN111651941A (zh) | 一种全球电离层电子总含量预测的算法 | |
CN110986747B (zh) | 一种滑坡位移组合预测方法及系统 | |
CN107092744B (zh) | 基于emd-svr的地表沉降量预测方法 | |
CN109685246B (zh) | 环境数据预估方法、装置及存储介质、服务器 | |
CN113379107B (zh) | 基于lstm和gcn的区域电离层tec预报方法 | |
Wen et al. | Short-Term PM2. 5 concentration prediction by combining GNSS and meteorological factors | |
CN111695473A (zh) | 基于长短时记忆网络模型的热带气旋强度客观监测方法 | |
Zhu et al. | A hybrid decomposition and Machine learning model for forecasting Chlorophyll-a and total nitrogen concentration in coastal waters | |
CN111626406A (zh) | 一种基于一次差分的bds钟差预测方法 | |
CN111365624A (zh) | 一种输卤管道泄漏检测的智能终端与方法 | |
Li et al. | A new PM2. 5 concentration forecasting system based on AdaBoost‐ensemble system with deep learning approach | |
CN116699731A (zh) | 一种热带气旋路径短期预报方法、系统及存储介质 | |
CN112946697B (zh) | 一种基于深度学习的卫星信号周跳探测与修复方法 | |
CN115755131A (zh) | 一种卫星定位的方法、装置及介质 | |
CN115408483A (zh) | 一种北斗导航服务定位性能智能预测方法及装置 | |
CN112700007B (zh) | 电离层电子含量预报模型的训练方法、预报方法及装置 | |
CN111814855B (zh) | 基于残差seq2seq神经网络的全球电离层电子总含量预测方法 | |
CN114037901A (zh) | 光伏发电预测导向的实时卫星近红外图像推算方法 | |
LU502828B1 (en) | Beidou navigation satellite system clock error prediction method based on first difference | |
Tang et al. | A review of machine learning for modeling air quality: Overlooked but important issues | |
Kibirige et al. | Using satellite data on remote transportation of air pollutants for PM2. 5 prediction in northern Taiwan | |
Nawab et al. | Comparison of ANN global horizontal irradiation predictions with satellite global horizontal irradiation using statistical evaluation | |
Li et al. | Improved receiver autonomous integrity monitoring algorithm | |
Yang et al. | Simulation System of Lake Eutrophication Evolution based on RS & GIS Technology—a Case Study in Wuhan East Lake |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |