CN115267847B - 一种单频周跳探测与修复方法、系统、设备及介质 - Google Patents
一种单频周跳探测与修复方法、系统、设备及介质 Download PDFInfo
- Publication number
- CN115267847B CN115267847B CN202210735243.4A CN202210735243A CN115267847B CN 115267847 B CN115267847 B CN 115267847B CN 202210735243 A CN202210735243 A CN 202210735243A CN 115267847 B CN115267847 B CN 115267847B
- Authority
- CN
- China
- Prior art keywords
- epoch
- target
- cycle slip
- slip detection
- value corresponding
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 148
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000008439 repair process Effects 0.000 title claims abstract description 47
- 238000001914 filtration Methods 0.000 claims abstract description 94
- 238000012545 processing Methods 0.000 claims abstract description 16
- 238000012549 training Methods 0.000 claims abstract description 12
- 238000000354 decomposition reaction Methods 0.000 claims description 35
- 230000006870 function Effects 0.000 claims description 21
- 238000001228 spectrum Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 9
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000013528 artificial neural network Methods 0.000 description 13
- 230000009286 beneficial effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种单频周跳探测与修复方法、系统、设备及介质,涉及卫星定位数据处理技术领域,所述方法包括:获取每个历元的伪距观测量和载波相位观测量,确定每个历元对应的单频周跳检测量;根据各个历元的单频周跳检测量,确定发生周跳的目标历元;对目标历元及其前N个历元的单频周跳检测量进行滤波处理,得到目标历元及其前N个历元中各个历元对应的Kalman滤波值;利用Kalman滤波值训练回归预测模型,通过训练后的回归预测模型对目标历元及其前一个历元对应的Kalman滤波值进行预测,得到目标历元及其前一个历元对应的模型预测值;根据目标历元及其前一个历元对应的Kalman滤波值和模型预测值,确定目标历元对应的组合预测值,根据组合预测值,完成周跳修复。
Description
技术领域
本发明涉及卫星定位数据处理技术领域,具体涉及一种单频周跳探测与修复方法、系统、设备及介质。
背景技术
北斗卫星导航系统是我国自行研制的全球卫星导航系统,可以为全球用户提供全天候、全天时、高精度的定位、导航和授时服务,是我国重要的时空基础设施。然而,由于北斗高精度定位所依赖的载波相位观测量中存在不可避免的周跳现象,卫星定位的精度也会因此受到影响。
周跳是指在全球卫星导航系统技术的载波相位测量中,由于卫星信号的失锁而导致的整周计数的跳变或中断。正确地探测并恢复周跳,是我们研究卫星定位精度的一个重要方向。北斗卫星导航接收机中,单频接收机由于价格低廉,因此应用范围广泛,单频接收机通过接收L1载波信号,测定载波相位观测值,进而实现定位。
目前,对于单频周跳的探测与修复,常用的方法有相位减伪距法、卡尔曼滤波法,但这些方法由于噪声的存在,导致卫星定位的精度不是很高。
发明内容
本发明所要解决的技术问题为:现有的用于单频周跳探测与修复的方法在进行周跳探测时,由于噪声的存在,对于小周跳修复精度不高,降低了卫星定位的精度。为解决该技术问题,本发明提供了一种单频周跳探测与修复方法、系统、设备及介质。
本发明解决上述技术问题的技术方案如下:
一种单频周跳探测与修复方法,包括:
步骤S1,获取每个历元的伪距观测量和载波相位观测量,对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量;
步骤S2,根据各个所述历元的单频周跳检测量,确定发生周跳的目标历元;
步骤S3,通过卡尔曼滤波算法对所述目标历元和所述目标历元的前N个历元的单频周跳检测量进行滤波处理,得到所述目标历元对应的Kalman滤波值和所述前N个历元中各个历元对应的Kalman滤波值,其中,所述N为常数,所述Kalman滤波值为对所述单频周跳检测量进行滤波后的单频周跳检测量;
步骤S4,通过所述前N个历元对应的Kalman滤波值训练预先构建的回归预测模型,得到训练后的回归预测模型;通过所述训练后的回归预测模型对所述目标历元的前一个历元对应的Kalman滤波值和所述目标历元对应的Kalman滤波值进行预测,得到所述前一个历元对应的模型预测值和所述目标历元对应的模型预测值,所述模型预测值为单频周跳预测值;
步骤S5,根据所述前一个历元对应的Kalman滤波值、所述前一个历元对应的模型预测值、所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值,确定所述目标历元对应的组合预测值;
步骤S6,根据所述组合预测值,完成周跳修复。
本发明的有益效果是:根据各个历元的观测数据(即伪距观测量和载波相位观测量)确定各个历元对应的单频周跳检测量,可消除观测数据中的部分噪声;根据单频周跳检测量,确定发生周跳的目标历元,利用卡尔曼滤波算法对目标历元及目标历元前的多个历元对应的单频周跳检测量进行滤波处理,可进一步减少单频周跳检测量中的噪声;再利用滤波处理得到的Kalman滤波值对回归预测模型进行在线学习训练,通过训练后的回归预测模型对目标历元及目标历元前一个历元对应的Kalman滤波值进行预测,得到模型预测值,根据Kalman滤波值和模型预测值,实现修复目标历元的周跳;本发明通过减少观测数据中的噪声,并联合卡尔曼滤波算法和回归预测模型,实现了对发生周跳的历元的探测与修复,提高了周跳修复精度和卫星定位的精度。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述根据各个所述历元的单频周跳检测量,确定发生周跳的目标历元,包括:
通过变分模态分解算法对每个所述单频周跳检测量进行信号分解,得到每个历元对应的多个本征模态函数分量;
对每个所述本征模态函数分量进行Hilbert变换,得到每个所述本征模态函数分量对应的Hilbert时频谱;
根据各个所述历元对应的Hilbert时频谱与预设阈值,确定发生周跳的目标历元。
采用上述进一步方案的有益效果是:利用变分模态分解算法对单频周跳检测量进行信号分解,便于后续确定发生周跳的历元;通过对本征模态函数分量进行Hilbert变换,再根据预设阈值确定发生周跳的历元,实现周跳的探测,为后续修复周跳奠定了基础。
进一步,所述根据所述组合预测值,完成周跳修复,包括:
将所述组合预测值替换所述目标历元的单频周跳检测量,得到目标历元对应的目标单频周跳检测量;
根据所述目标历元对应的目标单频周跳检测量、载波相位观测量和伪距观测量,以及所述目标历元的后一个历元的载波相位观测量和伪距观测量,确定所述目标历元对应的目标载波相位观测量;
将所述目标历元对应的目标载波相位观测量替换所述目标历元的载波相位观测量,完成周跳修复。
采用上述进一步方案的有益效果是:通过将组合预测值替换目标历元的单频周跳检测量,根据目标历元对应的目标单频周跳检测量、载波相位观测量和伪距观测量,以及目标历元的后一个历元的载波相位观测量和伪距观测量,确定目标历元对应的目标载波相位观测量,从而实现周跳的修复,提高了卫星定位的精度。
进一步,所述步骤S1中,所述获取每个历元的伪距观测量和载波相位观测量,具体包括:
通过单频接收机获取每个历元的伪距观测量和载波相位观测量,对于每个所述历元,该历元的伪距观测量和载波相位观测量分别为:
ρt=Rt+It+Tt+Δr,t-Δs,t+ερ,t
其中,ρt和分别为历元t对应的伪距观测量和载波相位观测量,Rt表示在历元t时单频接收机与卫星之间的几何距离,λ表示载波波长,Nt表示整周模糊度,It和Tt分别表示历元t对应的电离层延迟误差和对流层延迟误差,Δr,t和Δs,t分别表示历元t对应的接收机钟差和卫星钟差,ερ,t和/>分别表示历元t对应的所述伪距观测量的综合误差和所述载波相位观测量的综合误差,所述伪距观测量的综合误差包括多路径效应和伪距观测噪声,所述载波相位观测量的综合误差包括多路径效应和载波相位观测噪声;
所述步骤S1中,对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量,具体包括:
对于每个所述历元,将所述历元对应的载波相位观测量和所述伪距观测量作差,得到所述历元对应的观测量差量;
对于每个所述历元,根据所述历元对应的观测量差量、所述历元的下一个历元对应的观测量差量和载波波长,确定所述历元的单频周跳检测量。
采用上述进一步方案的有益效果是:根据获取的伪距观测量和载波相位观测量构造出单频周跳检测量,为后续探测单频周跳奠定了基础。
进一步,在通过变分模态分解算法对每个所述单频周跳检测量进行信号分解之前,所述方法还包括:
调整所述变分模态分解算法中的参数K和α,其中,参数K表示分解层数,参数α表示惩罚因子。
采用上述进一步方案的有益效果是:通过调整参数K和α,以使变分模态分解算法对单频周跳检测量的分解达到最优,再利用参数调整后的变分模态分解算法对单频周跳检测量进行信号分解,便于后续确定发生周跳的历元。
进一步,所述预设阈值为3倍标准差;
所述根据各个所述历元对应的Hilbert时频谱与预设阈值,确定发生周跳的目标历元,具体包括:
对于每个历元,确定该历元对应的Hilbert时频谱中的最大幅值;
当各所述最大幅值中存在大于3倍标准差的最大幅值时,将所述大于3倍标准差的最大幅值对应的历元确定为发生周跳的目标历元。
采用上述进一步方案的有益效果是:通过对本征模态函数分量进行Hilbert变换,再基于三倍标准差准则确定发生周跳的历元,为后续修复周跳奠定了基础。
进一步,所述步骤S5中,所述根据所述前一个历元对应的Kalman滤波值、所述前一个历元对应的模型预测值、所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值,确定所述目标历元对应的组合预测值,具体包括:
确定所述前一个历元的单频周跳检测量与Kalman滤波值之间的第一标准差σk-1,r;
确定所述前一个历元的单频周跳检测量与模型探测值之间的第二标准差σe-1,r;
确定所述目标历元对应的Kalman滤波值的第一权重wk和所述目标历元对应的模型预测值的第二权重we;
根据所述第一权重wk和所述第二权重we,对所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值进行加权求和,得到所述目标历元对应的组合预测值。
采用上述进一步方案的有益效果是:基于发生周跳的历元的前一个历元对应的Kalman滤波值、模型预测值设定权重,再根据设定的权重确定发生周跳的历元对应的组合预测值,实现对目标历元的Kalman滤波值和模型预测值进行校正,解决了单独使用回归预测模型进行周跳探测与修复存在的容易陷入局部最小值,以及单独使用卡尔曼滤波算法进行周跳探测与修复存在的欠滤波状态或过渡滤波状态,造成周跳修复精度不高的问题,提高了卫星定位精度。
为解决上述技术问题,本发明还提供了一种单频周跳探测与修复系统,包括:
数据预处理模块,用于获取每个历元的伪距观测量和载波相位观测量,对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量;
周跳探测模块,用于根据各个所述历元的单频周跳检测量,确定发生周跳的目标历元;
滤波处理模块,用于通过卡尔曼滤波算法对所述目标历元和所述目标历元的前N个历元的单频周跳检测量进行滤波处理,得到所述目标历元对应的Kalman滤波值和所述前N个历元中各个历元对应的Kalman滤波值,其中,所述N为常数,所述Kalman滤波值为对所述单频周跳检测量进行滤波后的单频周跳检测量;
模型处理模块,用于通过所述前N个历元对应的Kalman滤波值训练预先构建的回归预测模型,得到训练后的回归预测模型;通过所述训练后的回归预测模型对所述目标历元的前一个历元对应的Kalman滤波值和所述目标历元对应的Kalman滤波值进行预测,得到所述前一个历元对应的模型预测值和所述目标历元对应的模型预测值,所述模型预测值为单频周跳预测值;
周跳预测模块,用于根据所述前一个历元对应的Kalman滤波值、所述前一个历元对应的模型预测值、所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值,确定所述目标历元对应的组合预测值;
周跳修复模块,用于根据所述组合预测值,完成周跳修复。
为解决上述技术问题,本发明还提供了一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如前所述的单频周跳探测与修复方法。
为解决上述技术问题,本发明还提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如前所述的单频周跳探测与修复方法。
附图说明
图1为本发明的方法流程示意图;
图2为本发明实施例一中单频周跳检测量的示意图;
图3为本发明对图2的单频周跳检测量进行信号分解得到的本征模态函数分量的示意图;
图4为本发明的方法与现有方法对单频周跳进行修复的结果示意图;
图5为本发明的方法与现有方法对单频周跳进行修复的误差示意图;
图6为本发明的系统结构示意图;
图7为本发明的电子设备结构示意图。
具体实施方式
以下对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
实施例一
本实施例提供了一种单频周跳探测与修复方法,如图1所示,包括:
步骤S1,通过单频接收机获取每个历元的伪距观测量和载波相位观测量,对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量。
其中,所述步骤S1中,所述单频接收机只能接收L1载波信号,测定载波相位观测量进行定位,对于每个所述历元,该历元的伪距观测量和载波相位观测量分别为:
ρt=Rt+It+Tt+Δr,t-Δs,t+ερ,t
其中,ρt和分别为历元t对应的伪距观测量和载波相位观测量,Rt表示在历元t时单频接收机与卫星之间的几何距离,λ表示载波波长,Nt表示整周模糊度,It和Tt分别表示历元t对应的电离层延迟误差和对流层延迟误差,Δr,t和Δs,t分别表示历元t对应的接收机钟差和卫星钟差,ερ,t和/>分别表示历元t对应的所述伪距观测量的综合误差和所述载波相位观测量的综合误差,所述伪距观测量的综合误差包括多路径效应和伪距观测噪声,所述载波相位观测量的综合误差包括多路径效应和载波相位观测噪声。
对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量,具体包括:
对于每个所述历元,将所述历元对应的载波相位观测量和所述伪距观测量作差,以消除几何距离、接收机钟差、卫星钟差和对流层延迟误差,得到所述历元对应的观测量差量;
对于每个所述历元,根据所述历元对应的观测量差量、所述历元的下一个历元对应的观测量差量和载波波长,确定所述历元的单频周跳检测量。例如,对于历元t,其单频周跳检测量为:
其中,D(t)表示历元t的单频周跳观测量,分别为历元t+1、历元t对应的载波相位观测量,ρt+1、ρt分别为历元t+1、历元t对应的伪距观测量,λ表示载波波长。
步骤S2,根据各个所述历元的单频周跳检测量,确定发生周跳的目标历元。
其中,所述步骤S2具体包括:
通过变分模态分解算法对每个所述单频周跳检测量进行信号分解,得到每个历元对应的多个本征模态函数分量(简称IMF);
对每个所述本征模态函数分量进行Hilbert变换,得到每个所述本征模态函数分量对应的Hilbert时频谱;
根据各个所述历元对应的Hilbert时频谱与预设阈值,确定发生周跳的目标历元。
进一步,在通过变分模态分解算法对每个所述单频周跳检测量进行信号分解之前,所述方法还包括:
调整所述变分模态分解算法中的参数K和α,其中,参数K表示分解层数,参数α表示惩罚因子。通过调整参数K和α,以使变分模态分解算法对单频周跳检测量的分解达到最优。本实施例中,参数K的取值可以为6,参数α的取值可以为1000。
本实施例中,所述预设阈值为3倍标准差(3σ)。
所述根据各个所述历元对应的Hilbert时频谱与预设阈值,确定发生周跳的目标历元,具体包括:
对于每个历元,确定该历元对应的Hilbert时频谱中的最大幅值;
当各所述最大幅值中存在大于3倍标准差的最大幅值时,将所述大于3倍标准差的最大幅值对应的历元确定为发生周跳的目标历元。
步骤S3,通过卡尔曼滤波算法对所述目标历元和所述目标历元的前N个历元的单频周跳检测量进行滤波处理,得到所述目标历元对应的Kalman滤波值和所述前N个历元中各个历元对应的Kalman滤波值,其中,所述N为常数,所述Kalman滤波值为对所述单频周跳检测量进行滤波后的单频周跳检测量。本实施例中,N的取值为60。
进一步,在通过卡尔曼滤波算法对所述目标历元和所述目标历元的前N个历元的单频周跳检测量进行滤波处理之前,该方法还包括:
根据所述目标历元的前N个历元的单频周跳检测量调整卡尔曼滤波算法中的参数Q和R,其中,参数Q表示角度数据置信度,参数R表示角速度数据置信度。通过调整参数Q和R,以减少获取到的所述伪距观测量和所述载波相位观测量中的噪声,即减少所述单频周跳检测量中的噪声,提高对目标历元和目标历元的前N个历元的单频周跳检测量进行滤波处理的准确性,再利用目标历元的前N个历元对应的Kalman滤波值作为训练集输入回归预测模型,减少训练集噪声干扰,提高回归预测模型对待预测数据的预测准确率,以提高对周跳进行修复的准确度。
步骤S4,通过所述前N个历元对应的Kalman滤波值训练预先构建的回归预测模型,得到训练后的回归预测模型;通过所述训练后的回归预测模型对所述目标历元的前一个历元对应的Kalman滤波值和所述目标历元对应的Kalman滤波值进行预测,得到所述前一个历元对应的模型预测值和所述目标历元对应的模型预测值,所述模型预测值为单频周跳预测值。
本实施例中,所述回归预测模型选用对非线性时间序列建模效果良好的Elman神经网络,Elman神经网络结构中,输入层、隐含层、输出层的层数比值为2:6:1,利用所述前N个历元对应的Kalman滤波值对Elman神经网络进行在线学习,得到训练后的Elman神经网络,利用所述训练后的Elman神经网络对所述目标历元的前一个历元对应的Kalman滤波值和所述目标历元对应的Kalman滤波值进行预测,得到所述前一个历元对应的Elman神经网络预测值和所述目标历元对应的Elman神经网络预测值。
步骤S5,根据所述前一个历元对应的Kalman滤波值、所述前一个历元对应的模型预测值、所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值,确定所述目标历元对应的组合预测值。
其中,所述步骤S5具体包括:
确定所述前一个历元的单频周跳检测量与Kalman滤波值之间的第一标准差σk-1,r;
确定所述前一个历元的单频周跳检测量与模型探测值之间的第二标准差σe-1,r;
确定所述目标历元对应的Kalman滤波值的第一权重wk和所述目标历元对应的模型预测值的第二权重we;
根据所述第一权重wk和所述第二权重we,对所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值进行加权求和,得到所述目标历元对应的组合预测值。其中,所述目标历元对应的组合预测值为:
其中,Dr为所述目标历元对应的组合预测值,Dk为所述目标历元对应的Kalman滤波值,De为所述目标历元对应的模型预测值。
步骤S6,根据所述组合预测值,完成周跳修复。
其中,所述步骤S6具体包括:
将所述组合预测值替换所述目标历元的单频周跳检测量,得到目标历元对应的目标单频周跳检测量;
根据所述目标历元对应的目标单频周跳检测量、载波相位观测量和伪距观测量,以及所述目标历元的后一个历元的载波相位观测量和伪距观测量,确定所述目标历元对应的目标载波相位观测量;
将所述目标历元对应的目标载波相位观测量替换所述目标历元的载波相位观测量,完成周跳修复。例如,对于目标历元t,其目标单频周跳检测量为:
其中,D(t)'表示历元t的目标单频周跳观测量,分别为历元t+1、历元t对应的载波相位观测量、目标载波相位观测量,ρt+1、ρt分别为历元t+1、历元t对应的伪距观测量,λ表示载波波长。
为了进一步加深对本申请方案的理解,以下文所描述的仿真实验为例进行进一步的说明:
(1)实验条件
实验设备:北斗星通接收机
实验平台:Matlab R2020a软件
(2)实验内容
利用北斗星通接收机采集历元间隔为5s的北斗B1频段实测数据(实测数据包括伪距观测量和载波相位观测量),选取496个历元对应的北斗B1频段实测数据组成单频实验数据,采样频率为1Hz。
本方法首先根据获取的北斗B1频段实测数据确定每个历元的单频周跳检测量,其次利用变分模态分解算法对所述单频周跳检测量进行信号分解,得到6个本征模态函数分量,对每个所述本征模态函数分量进行Hilbert变换,从而确定发生周跳的目标历元,接着获取目标历元的前60个历元的载波相位观测量,得到单频周跳检测量集合,根据所述单频周跳检测量集合调整卡尔曼滤波算法中的参数Q和R,根据调整后得到的卡尔曼滤波算法对所述目标历元的前60个历元的载波相位观测量和所述目标历元的单频周跳检测量进行滤波处理,得到各个历元对应的Kalman滤波值,随后将所述目标历元的前60个历元的单频周跳检测量对应的Kalman滤波值作为训练集,利用所述训练集对回归预测模型进行在线学习,得到训练后的回归预测模型,利用所述训练后的回归预测模型对所述目标历元的前一历元对应的Kalman滤波值和所述目标历元对应的Kalman滤波值进行预测,得到所述目标历元的前一历元对应的模型预测值和所述目标对应的模型预测值。然后根据所述目标历元的前一历元对应的Kalman滤波值和模型预测值、所述目标历元对应的Kalman滤波值和模型预测值,确定所述目标历元对应的组合预测值,最后将所述组合预测值替换所述目标历元的单频周跳检测量,得到目标历元对应的目标单频周跳检测量,根据所述目标历元对应的目标单频周跳检测量、载波相位观测量和伪距观测量,以及所述目标历元的后一个历元的载波相位观测量和伪距观测量,确定所述目标历元对应的目标载波相位观测量,将所述目标历元对应的目标载波相位观测量替换所述目标历元的载波相位观测量,完成周跳修复。
(3)实验结果
以图2中的单频周跳检测量为例,利用变分模态分解算法对其进行信号分解,得到6个本征模态函数分量,对得到的6个本征模态函数分量进行Hilbert变换,变换结果如图3所示。根据3倍标准差准则可以得知,第3个本征模态函数分量(IMF3)和第5个本征模态函数分量(IMF5)中有发生周跳的历元。
图4中,横轴表示历元,纵轴表示差分观测量(即单频周跳检测量),原始值表示未进行周跳探测与修复前,根据实测数据得到的单频周跳检测量,“Elman神经网络”预测值表示单独使用Elman神经网络对周跳进行周跳探测及修复得到的用于替换单频周跳检测量的预测值,“Kalman滤波值”表示单独使用卡尔曼滤波算法进行周跳探测及修复得到的用于替换单频周跳检测量的预测值,“组合预测值”表示采用本方法进行周跳探测及修复得到的用于替换单频周跳检测量的预测值。
根据图4可以看出,使用单一预测模型(卡尔曼滤波算法或Elman神经网络)预测及修复周跳时,很容易得到超过1周的差分观测量,如卡尔曼滤波算法在历元64、65时的差分观测量的绝对值大于1,Elman神经网络在历元65时的差分观测量的绝对值也大于1,表明使用单一预测模型对周跳进行修复的效果并不理想,而通过本方法得到的所有差分观测量的绝对值均保持在0.5周以内,可见,本方法能够很好的探测并修复周跳。
图5中,横轴表示历元,纵轴表示误差,即通过不同方法进行周跳探测及修复得到的用于替换单频周跳检测量的预测值减去单频周跳检测量得到的差值,“Elman神经网络”表示单独使用Elman神经网络对周跳进行周跳探测及修复后得到的误差,“Kalman滤波”表示单独使用卡尔曼滤波算法进行周跳探测及修复后得到的误差,“KF-Elman组合”表示采用本方法进行周跳探测及修复后得到的误差,通过本方法进行周跳探测及修复得到的用于替换单频周跳检测量的预测值即组合预测值。
根据图5可以看出,通过本方法对周跳进行预测及修复所得到的误差的绝对值在0.5以内,而通过单一预测模型对周跳进行预测及修复所得到的误差的绝对值超过0.5,可见,本方法对周跳的修复效果优于单一预测模型对周跳的修复效果。
实施例二
基于与上述一种单频周跳探测与修复方法相同的原理,本实施例提供了一种单频周跳探测与修复系统,如图6所示,包括:
数据预处理模块,用于获取每个历元的伪距观测量和载波相位观测量,对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量。
周跳探测模块,用于根据各个所述历元的单频周跳检测量,确定发生周跳的目标历元。
所述周跳探测模块包括:
信号分解单元,用于通过变分模态分解算法对每个所述单频周跳检测量进行信号分解,得到每个历元对应的多个本征模态函数分量;
Hilbert变换单元,用于对每个所述本征模态函数分量进行Hilbert变换,得到每个所述本征模态函数分量对应的Hilbert时频谱;
周跳确定单元,用于根据各个所述历元对应的Hilbert时频谱与预设阈值,确定发生周跳的目标历元。
进一步,所述周跳探测模块还包括:
参数调整单元A,用于调整所述变分模态分解算法中的参数K和α,其中,参数K表示分解层数,参数α表示惩罚因子。
滤波处理模块,用于通过卡尔曼滤波算法对所述目标历元和所述目标历元的前N个历元的单频周跳检测量进行滤波处理,得到所述目标历元对应的Kalman滤波值和所述前N个历元中各个历元对应的Kalman滤波值,其中,所述N为常数,所述Kalman滤波值为对所述单频周跳检测量进行滤波后的单频周跳检测量。
进一步,所述滤波处理模块还包括:
参数调整单元B,用于根据所述目标历元的前N个历元的单频周跳检测量调整卡尔曼滤波算法中的参数Q和R,其中,参数Q表示角度数据置信度,参数R表示角速度数据置信度。
模型处理模块,用于通过所述前N个历元对应的Kalman滤波值训练预先构建的回归预测模型,得到训练后的回归预测模型;通过所述训练后的回归预测模型对所述目标历元的前一个历元对应的Kalman滤波值和所述目标历元对应的Kalman滤波值进行预测,得到所述前一个历元对应的模型预测值和所述目标历元对应的模型预测值,所述模型预测值为单频周跳预测值。
周跳预测模块,用于根据所述前一个历元对应的Kalman滤波值、所述前一个历元对应的模型预测值、所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值,确定所述目标历元对应的组合预测值。
所述周跳预测模块包括:
标准差确定单元,用于确定所述前一个历元的单频周跳检测量与Kalman滤波值之间的第一标准差σk-1,r,以及确定所述前一个历元的单频周跳检测量与模型探测值之间的第二标准差σe-1,r;
权重确定单元,用于确定所述目标历元对应的Kalman滤波值的第一权重wk和所述目标历元对应的模型预测值的第二权重we;
预测值确定单元,用于根据所述第一权重wk和所述第二权重we,对所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值进行加权求和,得到所述目标历元对应的组合预测值。
周跳修复模块,用于根据所述组合预测值,完成周跳修复。
实施例三
本实施例提供了一种电子设备,如图7所示,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如实施例一所述的单频周跳探测与修复方法。
实施例四
本实施例提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如实施例一所述的单频周跳探测与修复方法。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种单频周跳探测与修复方法,其特征在于,包括:
步骤S1,获取每个历元的伪距观测量和载波相位观测量,对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量;
步骤S2,根据各个所述历元的单频周跳检测量,确定发生周跳的目标历元;
步骤S3,通过卡尔曼滤波算法对所述目标历元和所述目标历元的前N个历元的单频周跳检测量进行滤波处理,得到所述目标历元对应的Kalman滤波值和所述前N个历元中各个历元对应的Kalman滤波值,其中,所述N为常数,所述Kalman滤波值为对所述单频周跳检测量进行滤波后的单频周跳检测量;
步骤S4,通过所述前N个历元对应的Kalman滤波值训练预先构建的回归预测模型,得到训练后的回归预测模型;通过所述训练后的回归预测模型对所述目标历元的前一个历元对应的Kalman滤波值和所述目标历元对应的Kalman滤波值进行预测,得到所述前一个历元对应的模型预测值和所述目标历元对应的模型预测值,所述模型预测值为单频周跳预测值;
步骤S5,根据所述前一个历元对应的Kalman滤波值、所述前一个历元对应的模型预测值、所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值,确定所述目标历元对应的组合预测值;
步骤S6,根据所述组合预测值,完成周跳修复。
2.根据权利要求1所述的方法,其特征在于,所述根据各个所述历元的单频周跳检测量,确定发生周跳的目标历元,包括:
通过变分模态分解算法对每个所述单频周跳检测量进行信号分解,得到每个历元对应的多个本征模态函数分量;
对每个所述本征模态函数分量进行Hilbert变换,得到每个所述本征模态函数分量对应的Hilbert时频谱;
根据各个所述历元对应的Hilbert时频谱与预设阈值,确定发生周跳的目标历元。
3.根据权利要求1所述的方法,其特征在于,所述根据所述组合预测值,完成周跳修复,包括:
将所述组合预测值替换所述目标历元的单频周跳检测量,得到目标历元对应的目标单频周跳检测量;
根据所述目标历元对应的目标单频周跳检测量、载波相位观测量和伪距观测量,以及所述目标历元的后一个历元的载波相位观测量和伪距观测量,确定所述目标历元对应的目标载波相位观测量;
将所述目标历元对应的目标载波相位观测量替换所述目标历元的载波相位观测量,完成周跳修复。
4.根据权利要求1至3中任一项所述的方法,其特征在于,所述步骤S1中,所述获取每个历元的伪距观测量和载波相位观测量,具体包括:
通过单频接收机获取每个历元的伪距观测量和载波相位观测量,对于每个所述历元,该历元的伪距观测量和载波相位观测量分别为:
ρt=Rt+It+Tt+Δr,t-Δs,t+ερ,t
其中,ρt和分别为历元t对应的伪距观测量和载波相位观测量,Rt表示在历元t时单频接收机与卫星之间的几何距离,λ表示载波波长,Nt表示整周模糊度,It和Tt分别表示历元t对应的电离层延迟误差和对流层延迟误差,Δr,t和Δs,t分别表示历元t对应的接收机钟差和卫星钟差,ερ,t和/>分别表示历元t对应的所述伪距观测量的综合误差和所述载波相位观测量的综合误差,所述伪距观测量的综合误差包括多路径效应和伪距观测噪声,所述载波相位观测量的综合误差包括多路径效应和载波相位观测噪声;
所述步骤S1中,对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量,具体包括:
对于每个所述历元,将所述历元对应的载波相位观测量和所述伪距观测量作差,得到所述历元对应的观测量差量;
对于每个所述历元,根据所述历元对应的观测量差量、所述历元的下一个历元对应的观测量差量和载波波长,确定所述历元的单频周跳检测量。
5.根据权利要求2所述的方法,其特征在于,在通过变分模态分解算法对每个所述单频周跳检测量进行信号分解之前,所述方法还包括:
调整所述变分模态分解算法中的参数K和α,其中,参数K表示分解层数,参数α表示惩罚因子。
6.根据权利要求2所述的方法,其特征在于,所述预设阈值为3倍标准差;
所述根据各个所述历元对应的Hilbert时频谱与预设阈值,确定发生周跳的目标历元,具体包括:
对于每个历元,确定该历元对应的Hilbert时频谱中的最大幅值;
当各所述最大幅值中存在大于3倍标准差的最大幅值时,将所述大于3倍标准差的最大幅值对应的历元确定为发生周跳的目标历元。
7.根据权利要求1至3中任一项所述的方法,其特征在于,所述步骤S5中,所述根据所述前一个历元对应的Kalman滤波值、所述前一个历元对应的模型预测值、所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值,确定所述目标历元对应的组合预测值,具体包括:
确定所述前一个历元的单频周跳检测量与Kalman滤波值之间的第一标准差σk-1,r;
确定所述前一个历元的单频周跳检测量与模型探测值之间的第二标准差σe-1,r;
确定所述目标历元对应的Kalman滤波值的第一权重wk和所述目标历元对应的模型预测值的第二权重we;
根据所述第一权重wk和所述第二权重we,对所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值进行加权求和,得到所述目标历元对应的组合预测值。
8.一种单频周跳探测与修复系统,其特征在于,包括:
数据预处理模块,用于获取每个历元的伪距观测量和载波相位观测量,对于每个所述历元,根据所述历元对应的伪距观测量和载波相位观测量,确定所述历元的单频周跳检测量;
周跳探测模块,用于根据各个所述历元的单频周跳检测量,确定发生周跳的目标历元;
滤波处理模块,用于通过卡尔曼滤波算法对所述目标历元和所述目标历元的前N个历元的单频周跳检测量进行滤波处理,得到所述目标历元对应的Kalman滤波值和所述前N个历元中各个历元对应的Kalman滤波值,其中,所述N为常数,所述Kalman滤波值为对所述单频周跳检测量进行滤波后的单频周跳检测量;
模型处理模块,用于通过所述前N个历元对应的Kalman滤波值训练预先构建的回归预测模型,得到训练后的回归预测模型;通过所述训练后的回归预测模型对所述目标历元的前一个历元对应的Kalman滤波值和所述目标历元对应的Kalman滤波值进行预测,得到所述前一个历元对应的模型预测值和所述目标历元对应的模型预测值,所述模型预测值为单频周跳预测值;
周跳预测模块,用于根据所述前一个历元对应的Kalman滤波值、所述前一个历元对应的模型预测值、所述目标历元对应的Kalman滤波值和所述目标历元对应的模型预测值,确定所述目标历元对应的组合预测值;
周跳修复模块,用于根据所述组合预测值,完成周跳修复。
9.一种电子设备,其特征在于,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如权利要求1至7中任一项所述的单频周跳探测与修复方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至7中任一项所述的单频周跳探测与修复方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210735243.4A CN115267847B (zh) | 2022-06-27 | 2022-06-27 | 一种单频周跳探测与修复方法、系统、设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210735243.4A CN115267847B (zh) | 2022-06-27 | 2022-06-27 | 一种单频周跳探测与修复方法、系统、设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115267847A CN115267847A (zh) | 2022-11-01 |
CN115267847B true CN115267847B (zh) | 2024-04-09 |
Family
ID=83763798
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210735243.4A Active CN115267847B (zh) | 2022-06-27 | 2022-06-27 | 一种单频周跳探测与修复方法、系统、设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115267847B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008039691A (ja) * | 2006-08-09 | 2008-02-21 | Toyota Motor Corp | 搬送波位相式測位装置 |
CN109407127A (zh) * | 2018-09-20 | 2019-03-01 | 昆明理工大学 | 一种北斗卫星导航系统载波相位周跳探测与修复方法 |
CN109581453A (zh) * | 2018-12-28 | 2019-04-05 | 北方信息控制研究院集团有限公司 | 基于周跳探测与修复的gnss分段平滑滤波方法 |
CN110208836A (zh) * | 2019-05-30 | 2019-09-06 | 东南大学 | 基于卡尔曼滤波的gnss高适应性周跳探测与修复方法 |
CN110531389A (zh) * | 2019-07-30 | 2019-12-03 | 湘潭大学 | Mw组合法与无几何相位法联合探测与修复周跳的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7961141B2 (en) * | 2008-12-09 | 2011-06-14 | Navcom Technology, Inc. | Methods and systems to increase accuracy in the navigation of single frequency receivers |
-
2022
- 2022-06-27 CN CN202210735243.4A patent/CN115267847B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008039691A (ja) * | 2006-08-09 | 2008-02-21 | Toyota Motor Corp | 搬送波位相式測位装置 |
CN109407127A (zh) * | 2018-09-20 | 2019-03-01 | 昆明理工大学 | 一种北斗卫星导航系统载波相位周跳探测与修复方法 |
CN109581453A (zh) * | 2018-12-28 | 2019-04-05 | 北方信息控制研究院集团有限公司 | 基于周跳探测与修复的gnss分段平滑滤波方法 |
CN110208836A (zh) * | 2019-05-30 | 2019-09-06 | 东南大学 | 基于卡尔曼滤波的gnss高适应性周跳探测与修复方法 |
CN110531389A (zh) * | 2019-07-30 | 2019-12-03 | 湘潭大学 | Mw组合法与无几何相位法联合探测与修复周跳的方法 |
Non-Patent Citations (3)
Title |
---|
Impacts of Ammonia Gas Plasma Surface Treatment on Polycrystalline-Silicon Junctionless Thin-Film Transistor;Yen-Wei Hsueh 等;IEEE Transactions on Plasma Science;20210131;第49卷(第1期);全文 * |
一种北斗三频实时周跳探测与修复新方法;满小三;孙付平;丁赫;刘帅;吴帅;叶险峰;;全球定位系统;20160215(第01期);全文 * |
多普勒积分重构与STPIR联合周跳探测与修复;蔡成林 等;测绘学报;20210215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115267847A (zh) | 2022-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
Zhang et al. | An improved robust adaptive Kalman filter for GNSS precise point positioning | |
CN108254773B (zh) | 一种多gnss的实时钟差解算方法 | |
CN108828640B (zh) | 一种卫星导航定位观测值定权方法及装置 | |
CN104714244A (zh) | 一种基于抗差自适应Kalman滤波的多系统动态PPP解算方法 | |
CN107678050A (zh) | 基于粒子滤波的glonass相位频间偏差实时追踪和精密估计方法 | |
CN110346823B (zh) | 可用于北斗精密单点定位的三频模糊度解算方法 | |
CN110687556A (zh) | 一种适用于laas的多径误差模型化方法 | |
CN110208836B (zh) | 基于卡尔曼滤波的gnss高适应性周跳探测与修复方法 | |
CN103728643A (zh) | 附有宽巷约束的北斗三频网络rtk模糊度单历元固定方法 | |
CN114355418A (zh) | 基于北斗地基增强系统的事后数据质量评估方法和系统 | |
CN111077550A (zh) | 一种应用于智能终端rtd定位的粗差探测方法及系统 | |
CN113358017B (zh) | 一种多站协同处理的gnss高精度变形监测方法 | |
CN107966722A (zh) | 一种gnss钟差解算方法 | |
CN116540303B (zh) | 一种历元间差分观测方程建立方法和地震位移计算方法 | |
CN115373005A (zh) | 卫星导航信号间高精度产品转化方法 | |
CN114859389A (zh) | 一种gnss多系统抗差自适应融合rtk解算方法 | |
CN112230254B (zh) | 一种gps载波相位多径误差的校正方法及装置 | |
CN104199054A (zh) | 一种用于北斗卫星导航系统共视数据的预处理方法 | |
CN111123315A (zh) | 非差非组合ppp模型的优化方法及装置、定位系统 | |
Yin et al. | A novel cycle slips detection model for the high precision positioning | |
CN115267847B (zh) | 一种单频周跳探测与修复方法、系统、设备及介质 | |
CN113031036A (zh) | 基于GNSS 30s采样频率数据的电离层相位闪烁因子构建方法 | |
Kuang et al. | Galileo real-time orbit determination with multi-frequency raw observations | |
CN115902968A (zh) | 基于北斗三号geo播发增强信息的ppp终端定位方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |