CN109765589B - 一种基于无电离层组合的三频gnss实时周跳固定技术 - Google Patents

一种基于无电离层组合的三频gnss实时周跳固定技术 Download PDF

Info

Publication number
CN109765589B
CN109765589B CN201910131040.2A CN201910131040A CN109765589B CN 109765589 B CN109765589 B CN 109765589B CN 201910131040 A CN201910131040 A CN 201910131040A CN 109765589 B CN109765589 B CN 109765589B
Authority
CN
China
Prior art keywords
cycle slip
combination
ewl
free
ionosphere
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
Application number
CN201910131040.2A
Other languages
English (en)
Other versions
CN109765589A (zh
Inventor
杨福鑫
李亮
赵琳
程建华
贾春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201910131040.2A priority Critical patent/CN109765589B/zh
Publication of CN109765589A publication Critical patent/CN109765589A/zh
Application granted granted Critical
Publication of CN109765589B publication Critical patent/CN109765589B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明属于卫星导航定位技术研究领域,具体涉及一种基于无电离层组合的三频GNSS实时周跳固定技术,包括以下步骤:接收广播星历电文、实时精密卫星修正电文和导航电文并解析;根据实时精密卫星轨道和钟差的修正量及主要误差的模型化修正,提高伪距和载波相位观测量的测距精度;在第k‑1和第k个历元间进行原始观测量差分;构建无电离层EWL组合;基于固定的超宽项周跳,构建无电离层WL组合;基于固定的超宽项和宽巷周跳,构建无电离层NL组合;本发明简化了传统多频组合周跳探测和修复的步骤,提高了计算效率以及周跳在电离层扰动情况下的固定成功率。

Description

一种基于无电离层组合的三频GNSS实时周跳固定技术
技术领域
本发明属于卫星导航定位技术研究领域,具体涉及一种基于无电离层组合的三频GNSS实时周跳固定技术。
背景技术
随着全球卫星导航系统(Global Navigation Satellite System,GNSS)现代化的发展,提供实时高精度的定位服务是GNSS未来发展的主要方向。载波相位观测量对于高精度GNSS定位至关重要,例如实时动态(Real-time Kinematic,RTK)差分定位技术和精密单点定位(Precise Point Positioning,PPP)技术,因为它的量测精度可达mm级。在实际应用中,由于信号失锁等原因而导致的载波相位观测量中出现整周模糊度不连续(即周跳)的现象,该现象可造成PPP定位精度以及RTK模糊度固定收敛时间的增加,进而降低高精度GNSS定位的性能。因此,周跳的探测与修复对于基于载波相位的高精度GNSS定位技术至关重要。
电离层延迟的处理是实现周跳探测与修复的关键因素之一,特别是在电离层闪烁或低采样率情况下。目前,在周跳探测与修复方法中电离层抑制的手段包括:(1)引入外部修正信息,但受限于电离层模型精度,导致周跳探测与修复的可靠性降低,同时引入了通信的负担。(2)将电离层延迟作为参数待估项,但该方法通常需要几分钟的时间才能获取到相对准确的电离层延迟误差,难以满足实时性的要求。(3)引入差分技术,但电离层闪烁将导致电离层延迟项的弱相关性,使得电离层修正后的残差较大。(4)随着GNSS现代化的发展,三频及以上信号的播发为多频观测量组合的电离层抑制周跳探测与修复方法提供了便利,但现有的多频观测量组合方法通常假设在高采样率下历元间的电离层延迟可以忽略,而随着历元间隔的增加,势必降低周跳探测与修复的可靠性。综上所述,设计一种新型可解决电离层扰动的三频周跳探测与修复技术具有相当的迫切性。
发明内容
为解决上述问题,本发明提供了一种基于无电离层组合的三频GNSS实时周跳固定方法,本方法同时采用实时精密卫星轨道及钟差修正信息、广播星历信息以及模型化修正的方法对原始伪距和载波相位观测量中的主要误差源进行修正。在此基础上将修正后的三频原始观测量进行历元差分,完成无电离层超宽巷(Extra-wide Lane,EWL)、宽巷(WideLane,WL)以及窄巷(Narrow Lane,NL)组合的构建,并通过EWL、WL和NL周跳的线性关系,最终实现原始观测量中周跳的确定。
一种基于无电离层组合的三频GNSS实时周跳固定技术,包括以下步骤:
(1)实时接收广播星历电文、实时精密卫星修正电文和导航电文并解析;
(2)根据实时精密卫星轨道和钟差的修正量及主要误差的模型化修正,提高伪距和载波相位观测量的测距精度;
(3)在第k-1和第k个历元间进行原始观测量差分;
(4)构建无电离层EWL组合;
(5)基于固定的超宽项周跳,构建无电离层WL组合;
(6)基于固定的超宽项和宽巷周跳,构建无电离层NL组合;
(7)如果步骤(5)或步骤(6)固定周跳失败,则等待到下一历元,进行多历元累积,直到步骤(5)和步骤(6)固定成功为止;
(8)利用EWL、WL和NL周跳的线性关系计算各原始观测量中的周跳值。
所述实时接收广播星历电文、实时精密卫星修正电文和导航电文并解析,包括:
接收卫星的个数为n,且接收卫星的个数满足如下条件:
n>3+sys
其中,sys表示系统个数。
所述在第k-1和第k个历元间进行原始观测量差分,包括:
历元差分的伪距及载波相位观测量为:
Figure BDA0001974609450000021
Figure BDA0001974609450000022
其中,Δ表示相邻历元的差值,P和L分别表示伪距和载波相位观测量,j表示卫星号,g表示频率,g的取值为1,2,3,
Figure BDA0001974609450000023
为修正后站星距离,I1表示在f1频段上的电离层延迟,
Figure BDA0001974609450000024
为电离层系数,λg表示在fg频段上的载波相位波长,εg和ξg分别表示在fg频段上的伪距和载波相位观测量噪声。
所述构建无电离层EWL组合,包括:
Figure BDA0001974609450000025
其中,
Figure BDA0001974609450000026
为无电离层超宽巷载波相位观测量组合;
通过就近取整方法实现超宽巷周跳的固定,超宽项周跳表示为:
Figure BDA0001974609450000027
EWL周跳的浮点解为:
Figure BDA0001974609450000031
其中,λEWL=c/(f2-f3)为超宽巷周跳组合的波长,c为光速;
采用就近取整的方式完成超宽巷周跳的固定:
Figure BDA0001974609450000032
其中,round表示就近取整。
所述基于固定的超宽项周跳,构建无电离层WL组合,包括:
Figure BDA0001974609450000033
其中,
Figure BDA0001974609450000034
为无电离层伪距观测量组合,
Figure BDA0001974609450000035
为无电离层宽巷载波相位组合;
由最小二乘输出的协方差矩阵和宽巷模糊度的浮点解,通过LAMBDA算法进行宽巷模糊度的固定,宽项周跳表示为:
Figure BDA0001974609450000036
结合超宽项周跳和无电离层EWL组合得:
Figure BDA0001974609450000037
其中,
Figure BDA0001974609450000038
Figure BDA0001974609450000039
分别表示在无电离层WL组合的伪距和载波相位观测量噪声;
联立n卫星的观测量,进行加权最小二乘求解,获得状态信息为:
Figure BDA00019746094500000310
及其协方差矩阵,其中(Δx,Δy,Δz)为三维位置变化量,Δdt为接收机钟差变化量,
Figure BDA00019746094500000311
为宽巷周跳浮点解;
通过宽巷周跳浮点解及其方差协方差阵采用LAMBDA的方法进行固定,采用比率检验的方法检验最优解。
所述基于固定的超宽项和宽巷周跳,构建无电离层NL组合,包括:
Figure BDA0001974609450000041
其中,
Figure BDA0001974609450000042
Figure BDA0001974609450000043
分别为f1、f2和f1、f3的无电离层窄巷载波相位组合。
所述利用EWL、WL和NL周跳的线性关系计算各原始观测量中的周跳值,包括:
Figure BDA0001974609450000044
本发明的有益效果在于:
针对传统周跳探测和修复方法在电离层扰动情况下成功率不高的问题,本发明充分利用多频观测量组合的优势,构建了无电离层EWL、WL以及NL组合,并利用不同观测量组合下的周跳间线性关系,逐次完成无电离层EWL、WL以及NL组合的周跳固定,即满足了在周跳固定过程中对电离层的抑制。本发明所提出的基于无电离层组合的三频GNSS实时周跳固定方法,简化了传统多频组合周跳探测和修复的步骤,提高了计算效率以及周跳在电离层扰动情况下的固定成功率。
附图说明
图1是本发明的工作流程图;
具体实施方式
下面结合附图对本发明做进一步描述。
实施例:
一种基于无电离层组合的三频GNSS实时周跳固定方法,具体步骤为:
步骤1,
广播星历电文、实时精密卫星修正电文和导航电文并解析;假设接收卫星个数为n(n>3+sys),其中,sys表示系统个数。实时精密卫星修正电文格式遵循RTCM 3.02标准中的状态空间表示(State Space Representation,SSR)。
步骤2,误差修正;
a.修正卫星轨道与钟差;由于SSR修正信息对星轨与星钟的校正是基于广播星历的,因而需要先利用广播星历解算出卫星位置和卫星钟差,然后根据SSR广播的卫星轨道及钟差改正量对卫星轨道及钟差误差进行校正。
b.模型化误差;如地球自旋转、天线相位缠绕以及地场位移效应等采用经验模型修正
步骤3,历元差分观测量构建;
在第k-1和第k个历元间进行GNSS原始观测量差分后,消除了随时间变化缓慢的误差,例如,码硬件延迟偏差,未校准相位延迟偏差以及对流层湿分量等。此时,历元差分的伪距及载波相位观测量如式1所示,包含了历元间的三维位置变化量、钟差变化量、电离层变化量以及噪声变化量。
步骤4,无电离层EWL组合构建;
a.通过式3的EWL组合消除了电离层误差和几何相关项,则EWL周跳的浮点解可表示为,
Figure BDA0001974609450000051
其中,λEWL=c/(f2-f3)为超宽巷周跳组合的波长,c为光速。
b.式6采用就近取整的方式完成超宽巷周跳的固定。
Figure BDA0001974609450000052
其中,round表示就近取整。
步骤5,无电离层WL组合构建;
a.结合式6求得的超宽项周跳和式4可得,
Figure BDA0001974609450000053
式中,
Figure BDA0001974609450000054
Figure BDA0001974609450000055
分别表示在无电离层WL组合的伪距和载波相位观测量噪声。
b.联立n卫星的观测量,进行加权最小二乘求解,获得状态信息为:
Figure BDA0001974609450000056
及其协方差矩阵,其中(Δx,Δy,Δz)为三维位置变化量,Δdt为接收机钟差变化量,
Figure BDA0001974609450000057
为宽巷周跳浮点解。
c.通过b中求解的宽巷周跳浮点解及其方差协方差阵采用LAMBDA的方法进行固定,采用比率检验(ratio-test)的方法检验最优解,检验次优解与最优解之间比率的门限值设为3。
d.如固定不成功,则进入下一历元,联立第k-1与第k个历元差分观测量和第k-1与第k+1个历元差分观测量,重复步骤1。反之,获得固定的宽巷周跳
Figure BDA0001974609450000058
进行一步。
步骤6,无电离层NL组合构建;
a.结合式6和式7求得的超宽项和宽巷周跳与式5得,
Figure BDA0001974609450000061
式中,
Figure BDA0001974609450000062
Figure BDA0001974609450000063
分别表示在无电离层NL组合的载波相位观测量噪声。
b.联立n卫星的观测量,进行加权最小二乘求解,获得状态信息为:
Figure BDA0001974609450000064
及其协方差矩阵,其中(Δx,Δy,Δz)为三维位置变化量,Δdt为接收机钟差变化量,
Figure BDA0001974609450000065
为宽巷周跳浮点解。
c.通过b中求解的窄巷周跳浮点解及其方差协方差阵采用LAMBDA的方法进行固定,采用比率检验的方法检验最优解,检验次优解与最优解之间比率的门限值设为3。
d.如固定不成功,则进入下一历元,联立第k-1与第k个历元差分观测量和第k-1与第k+1个历元差分观测量,重复步骤1。反之,获得固定的窄巷周跳
Figure BDA0001974609450000066
进行一步。
步骤7,利用EWL、WL和NL周跳的线性关系计算各原始观测量中的周跳值。
Figure BDA0001974609450000067
步骤8,根据固定的周跳对载波相位观测量进行修正。
整体执行流程图如图1所示。
随着全球卫星导航系统(Global Navigation Satellite System,GNSS)现代化的发展,提供实时高精度的定位服务是GNSS未来发展的主要方向。载波相位观测量对于高精度GNSS定位至关重要,例如实时动态(Real-time Kinematic,RTK)差分定位技术和精密单点定位(Precise Point Positioning,PPP)技术,因为它的量测精度可达mm级。在实际应用中,由于信号失锁等原因而导致的载波相位观测量中出现整周模糊度不连续(即周跳)的现象,该现象可造成PPP定位精度以及RTK模糊度固定收敛时间的增加,进而降低高精度GNSS定位的性能。因此,周跳的探测与修复对于基于载波相位的高精度GNSS定位技术至关重要。
电离层延迟的处理是实现周跳探测与修复的关键因素之一,特别是在电离层闪烁或低采样率情况下。目前,在周跳探测与修复方法中电离层抑制的手段包括:(1)引入外部修正信息,但受限于电离层模型精度,导致周跳探测与修复的可靠性降低,同时引入了通信的负担。(2)将电离层延迟作为参数待估项,但该方法通常需要几分钟的时间才能获取到相对准确的电离层延迟误差,难以满足实时性的要求。(3)引入差分技术,但电离层闪烁将导致电离层延迟项的弱相关性,使得电离层修正后的残差较大。(4)随着GNSS现代化的发展,三频及以上信号的播发为多频观测量组合的电离层抑制周跳探测与修复方法提供了便利,但现有的多频观测量组合方法通常假设在高采样率下历元间的电离层延迟可以忽略,而随着历元间隔的增加,势必降低周跳探测与修复的可靠性。综上所述,设计一种新型可解决电离层扰动的三频周跳探测与修复技术具有相当的迫切性。
本发明的一种基于无电离层组合的三频GNSS实时周跳固定方法包括:
步骤1,实时接收广播星历电文、实时精密卫星修正电文和导航电文并解析;
步骤2,根据实时精密卫星轨道和钟差的修正量及主要误差的模型化修正,提高伪距和载波相位观测量的测距精度;
步骤3,在第k-1和第k个历元间进行原始观测量差分,
Figure BDA0001974609450000071
Figure BDA0001974609450000072
式中,Δ表示相邻历元的差值,P和L分别表示伪距和载波相位观测量,j表示卫星号,g表示频率(g=1,2,3),
Figure BDA0001974609450000073
修正后站星距离,I1表示在f1频段上的电离层延迟,
Figure BDA0001974609450000074
为电离层系数,λg表示在fg频段上的载波相位波长,εg和ξg分别表示在fg频段上的伪距和载波相位观测量噪声。
步骤4,构建无电离层EWL组合,
Figure BDA0001974609450000075
式中,
Figure BDA0001974609450000076
为无电离层超宽巷载波相位观测量组合,并通过就近取整方法实现超宽巷周跳的固定,超宽项周跳表示为
Figure BDA0001974609450000077
步骤5,基于固定的超宽项周跳,构建无电离层WL组合,
Figure BDA0001974609450000078
式中,
Figure BDA0001974609450000079
为无电离层伪距观测量组合,
Figure BDA00019746094500000710
为无电离层宽巷载波相位组合。由最小二乘输出的协方差矩阵和宽巷模糊度的浮点解,通过LAMBDA(The Least-squaresAMBiguity Decorrelation Adjustment)算法进行宽巷模糊度的固定,宽项周跳表示为
Figure BDA0001974609450000081
步骤6,基于固定的超宽项和宽巷周跳,构建无电离层NL组合,
Figure BDA0001974609450000082
式中,
Figure BDA0001974609450000083
Figure BDA0001974609450000084
分别为f1、f2和f1、f3的无电离层窄巷载波相位组合。由最小二乘输出的协方差矩阵和宽巷模糊度的浮点解,通过LAMBDA算法进行宽巷模糊度的固定,窄巷周跳表示为
Figure BDA0001974609450000085
步骤7,根据步骤4、步骤5和步骤6得到的超宽巷、宽巷和窄巷周跳最终确定三频观测量上的周跳值。
步骤8,如步骤5或步骤6固定周跳失败,则等待到下一历元,进行多历元累积,直到步骤5和步骤6固定成功为止。
针对传统周跳探测和修复方法在电离层扰动情况下成功率不高的问题,本发明充分利用多频观测量组合的优势,构建了无电离层EWL、WL以及NL组合,并利用不同观测量组合下的周跳间线性关系,逐次完成无电离层EWL、WL以及NL组合的周跳固定,即满足了在周跳固定过程中对电离层的抑制。本发明所提出的基于无电离层组合的三频GNSS实时周跳固定方法,简化了传统多频组合周跳探测和修复的步骤,提高了计算效率以及周跳在电离层扰动情况下的固定成功率。
一种基于无电离层组合的三频GNSS实时周跳固定技术,包括:
步骤1,实时接收广播星历电文、实时精密卫星修正电文和导航电文并解析;
步骤2,根据实时精密卫星轨道和钟差的修正量及主要误差的模型化修正,提高伪距和载波相位观测量的测距精度;
步骤3,在第k-1和第k个历元间进行原始观测量差分,
Figure BDA0001974609450000086
Figure BDA0001974609450000087
式中,Δ表示相邻历元的差值,P和L分别表示伪距和载波相位观测量,j表示卫星号,g表示频率(g=1,2,3),
Figure BDA0001974609450000088
修正后站星距离,I1表示在f1频段上的电离层延迟,
Figure BDA0001974609450000091
为电离层系数,λg表示在fg频段上的载波相位波长,εg和ξg分别表示在fg频段上的伪距和载波相位观测量噪声。
步骤4,构建无电离层EWL组合,
Figure BDA0001974609450000092
式中,
Figure BDA0001974609450000093
为无电离层超宽巷载波相位观测量组合,并通过就近取整方法实现超宽巷周跳的固定,超宽项周跳表示为
Figure BDA0001974609450000094
步骤5,基于固定的超宽项周跳,构建无电离层WL组合,
Figure BDA0001974609450000095
式中,
Figure BDA0001974609450000096
为无电离层伪距观测量组合,
Figure BDA0001974609450000097
为无电离层宽巷载波相位组合。由最小二乘输出的协方差矩阵和宽巷模糊度的浮点解,通过LAMBDA算法进行宽巷模糊度的固定,宽项周跳表示为
Figure BDA0001974609450000098
步骤6,基于固定的超宽项和宽巷周跳,构建无电离层NL组合,
Figure BDA0001974609450000099
式中,
Figure BDA00019746094500000910
Figure BDA00019746094500000911
分别为f1、f2和f1、f3的无电离层窄巷载波相位组合。由最小二乘输出的协方差矩阵和宽巷模糊度的浮点解,通过LAMBDA算法进行宽巷模糊度的固定,窄巷周跳表示为
Figure BDA00019746094500000912
步骤7,根据步骤4、步骤5和步骤6得到的超宽巷、宽巷和窄巷周跳最终确定三频观测量上的周跳值。
步骤8,如步骤5或步骤6固定周跳失败,则等待到下一历元,进行多历元累积,直到步骤5和步骤6固定成功为止。
步骤4构建了无电离层EWL组合,如式3所示,周跳固定的方法不仅限于就近取整法。
步骤5、6中构建了无电离层WL和NL组合,且在构建无电离层WL和NL组合时,不仅限于式4和式5的组合形式,其它根据电离层频率特性组成的无电离层组合观测量进行周跳固定方法仍属于本专利权利范围。
步骤5、6中的周跳固定采用LAMBDA方法,但门限值不仅限于3。
步骤4、5、6中的周跳固定方法仅采用了历元差分。

Claims (7)

1.一种基于无电离层组合的三频GNSS实时周跳固定方法,其特征在于,包括以下步骤:
(1)实时接收广播星历电文、实时精密卫星修正电文和导航电文并解析;
(2)根据实时精密卫星轨道和钟差的修正量及主要误差的模型化修正,提高伪距和载波相位观测量的测距精度;
(3)在第k-1和第k个历元间进行原始观测量差分;
(4)构建无电离层EWL组合;
(5)基于固定的EWL周跳,构建无电离层WL组合;
(6)基于固定的EWL和WL周跳,构建无电离层NL组合;所述无电离层NL组合具体包括:
Figure FDA0003858508540000011
式中,
Figure FDA0003858508540000012
Figure FDA0003858508540000013
分别为f1、f2和f1、f3的历元间无电离层窄巷载波相位组合;f1、f2和f3分别代表1频、2频和3频对应的频率;
Figure FDA0003858508540000014
Figure FDA0003858508540000015
分别表示在无电离层NL组合的
Figure FDA0003858508540000016
Figure FDA0003858508540000017
载波相位观测量噪声;c为光速;j表示卫星号;
Figure FDA0003858508540000018
为修正后站星距离;
Figure FDA0003858508540000019
为卫星号为j的固定后的WL周跳;
Figure FDA00038585085400000110
为卫星号为j的固定后的EWL周跳;Δ表示相邻历元的差值;
(7)如果步骤(5)或步骤(6)固定周跳失败,则等待到下一历元,进行多历元累积,直到步骤(5)和步骤(6)固定成功为止;
(8)利用EWL、WL和NL周跳的线性关系计算各原始观测量中的周跳值。
2.根据权利要求1所述的一种基于无电离层组合的三频GNSS实时周跳固定方法,其特征在于,所述实时接收广播星历电文、实时精密卫星修正电文和导航电文并解析,包括:
接收卫星的个数为n,且接收卫星的个数满足如下条件:
n>3+sys
其中,sys表示系统个数。
3.根据权利要求1所述的一种基于无电离层组合的三频GNSS实时周跳固定方法,其特征在于,所述在第k-1和第k个历元间进行原始观测量差分,包括:
历元差分的伪距及载波相位观测量为:
Figure FDA0003858508540000021
Figure FDA0003858508540000022
其中,Δ表示相邻历元的差值,P和L分别表示伪距和载波相位观测量,j表示卫星号,g的取值为1,2,3,表示频段f的编号;
Figure FDA0003858508540000023
为修正后的站星距离,I1表示在f1频段上的电离层延迟,
Figure FDA0003858508540000024
为电离层系数,λg表示在fg频段上的载波相位波长,εg和ξg分别表示在fg频段上的伪距和载波相位观测量噪声。
4.根据权利要求3所述的一种基于无电离层组合的三频GNSS实时周跳固定方法,其特征在于,所述构建无电离层EWL组合,包括:
Figure FDA0003858508540000025
其中,Δ表示相邻历元的差值,
Figure FDA0003858508540000026
为无电离层EWL载波相位观测量组合;
通过就近取整方法实现EWL周跳的固定,EWL周跳表示为:
Figure FDA0003858508540000027
EWL周跳的浮点解为:
Figure FDA0003858508540000028
其中,λEWL=c/(f2-f3)为EWL周跳组合的波长,c为光速;
采用就近取整的方式完成EWL周跳的固定:
Figure FDA0003858508540000029
其中,round表示就近取整;Δ表示相邻历元的差值。
5.根据权利要求4所述的一种基于无电离层组合的三频GNSS实时周跳固定方法,其特征在于,所述基于固定的EWL周跳,构建无电离层WL组合,包括:
Figure FDA00038585085400000210
其中,Δ表示相邻历元的差值,
Figure FDA0003858508540000031
为无电离层伪距观测量组合,
Figure FDA0003858508540000032
为无电离层WL载波相位组合;
由最小二乘输出的协方差矩阵和WL模糊度的浮点解,通过LAMBDA算法进行WL模糊度的固定,WL周跳表示为:
Figure FDA0003858508540000033
结合EWL周跳和无电离层WL组合得:
Figure FDA0003858508540000034
联立n卫星的观测量,进行加权最小二乘求解,获得状态信息为:
Figure FDA0003858508540000035
及其协方差矩阵,其中(Δx,Δy,Δz)为三维位置变化量,Δdt为接收机钟差变化量,
Figure FDA0003858508540000036
为WL周跳浮点解;
通过WL周跳浮点解及其方差协方差阵采用LAMBDA的方法进行固定,采用比率检验的方法检验最优解;
其中,Δ表示相邻历元的差值;c为光速;
Figure FDA0003858508540000037
Figure FDA0003858508540000038
分别表示在无电离层WL组合的伪距和载波相位观测量噪声。
6.根据权利要求5所述的一种基于无电离层组合的三频GNSS实时周跳固定方法,其特征在于,构建无电离层NL组合时,包括:
Figure FDA0003858508540000039
其中,Δ表示相邻历元的差值;
Figure FDA00038585085400000310
Figure FDA00038585085400000311
分别为f1、f2和f1、f3的无电离层窄巷载波相位组合。
7.根据权利要求6所述的一种基于无电离层组合的三频GNSS实时周跳固定方法,其特征在于,所述利用EWL、WL和NL周跳的线性关系计算各原始观测量中的周跳值,包括:
Figure FDA0003858508540000041
CN201910131040.2A 2019-02-21 2019-02-21 一种基于无电离层组合的三频gnss实时周跳固定技术 Active CN109765589B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910131040.2A CN109765589B (zh) 2019-02-21 2019-02-21 一种基于无电离层组合的三频gnss实时周跳固定技术

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910131040.2A CN109765589B (zh) 2019-02-21 2019-02-21 一种基于无电离层组合的三频gnss实时周跳固定技术

Publications (2)

Publication Number Publication Date
CN109765589A CN109765589A (zh) 2019-05-17
CN109765589B true CN109765589B (zh) 2022-12-02

Family

ID=66456937

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910131040.2A Active CN109765589B (zh) 2019-02-21 2019-02-21 一种基于无电离层组合的三频gnss实时周跳固定技术

Country Status (1)

Country Link
CN (1) CN109765589B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208835B (zh) * 2019-05-21 2023-05-05 哈尔滨工程大学 一种基于消电离层组合的跨系统紧组合差分定位方法
CN110244332B (zh) * 2019-06-06 2021-09-14 中国科学院国家天文台 一种天线阵大气相位扰动修正方法及系统
CN110398762A (zh) * 2019-07-15 2019-11-01 广州中海达卫星导航技术股份有限公司 实时钟差估计中的模糊度固定方法、装置、设备及介质
CN110398758A (zh) * 2019-07-24 2019-11-01 广州中海达卫星导航技术股份有限公司 实时钟差估计中的粗差探测方法、装置、设备及存储介质
CN111045066A (zh) * 2019-12-30 2020-04-21 威海欧瑞亚信息科技有限公司 一种基于参数等价约化原理确定gnss位置变化的方法
CN111856525B (zh) * 2020-06-29 2023-01-03 哈尔滨工程大学 一种基于lstm神经网络的周跳探测和修复方法
CN112083462B (zh) * 2020-09-14 2022-12-13 哈尔滨工程大学 一种海洋精密单点定位方法
CN114675314B (zh) * 2022-05-30 2022-08-26 深圳市智联时空科技有限公司 一种重新收敛的精密单点定位方法
CN115166799B (zh) * 2022-07-11 2023-07-21 中国科学院精密测量科学与技术创新研究院 顾及硬件延迟时变特征的gnss精密单点定位方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006108227A1 (en) * 2005-04-11 2006-10-19 Yanming Feng Improved phase ambiguity resolution using three gnss signals
CN102650692A (zh) * 2011-02-25 2012-08-29 中国人民解放军61081部队 利用北斗三频观测量进行周跳探测与修复的方法
CN106772512A (zh) * 2016-11-11 2017-05-31 哈尔滨工程大学 一种基于消电离层‑噪声约束的三频模糊度解算方法
CN107132558A (zh) * 2017-06-13 2017-09-05 武汉大学 惯性辅助的多频多模gnss周跳修复方法及系统
CN107607966A (zh) * 2017-08-08 2018-01-19 北京大学 一种基于gnss三频载波相位的周跳检测与修复方法
CN109358350A (zh) * 2018-10-08 2019-02-19 中国人民解放军战略支援部队信息工程大学 一种北斗三频周跳探测方法与装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2849209B1 (fr) * 2002-12-19 2007-04-06 Agence Spatiale Europeenne Procede et systeme de navigation en temps reel a l'aide de signaux radioelectriques a trois porteuses emis par des satellites et de corrections ionospheriques

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006108227A1 (en) * 2005-04-11 2006-10-19 Yanming Feng Improved phase ambiguity resolution using three gnss signals
CN102650692A (zh) * 2011-02-25 2012-08-29 中国人民解放军61081部队 利用北斗三频观测量进行周跳探测与修复的方法
CN106772512A (zh) * 2016-11-11 2017-05-31 哈尔滨工程大学 一种基于消电离层‑噪声约束的三频模糊度解算方法
CN107132558A (zh) * 2017-06-13 2017-09-05 武汉大学 惯性辅助的多频多模gnss周跳修复方法及系统
CN107607966A (zh) * 2017-08-08 2018-01-19 北京大学 一种基于gnss三频载波相位的周跳检测与修复方法
CN109358350A (zh) * 2018-10-08 2019-02-19 中国人民解放军战略支援部队信息工程大学 一种北斗三频周跳探测方法与装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Ionosphere-Constrained Triple-Frequency Cycle Slip Fixing Method for the Rapid Re-Initialization of PPP;Fuxin Yang,et al;《sensors》;20181231;p1-16 *
Real-time detection and repair of cycle slips in triple-frequency GNSS measurements;Qile Zhao,et al;《GPS Solut》;20140725;p381-391 *
一种改进的基于BDS三频非差观测的周跳实时探测与修复模型;高杰等;《测绘工程》;20161231;第25-31页 *
基于贝叶斯压缩感知的周跳探测与修复方法;李慧等;《物理学报》;20161231;第1-8页 *

Also Published As

Publication number Publication date
CN109765589A (zh) 2019-05-17

Similar Documents

Publication Publication Date Title
CN109765589B (zh) 一种基于无电离层组合的三频gnss实时周跳固定技术
CN109709591B (zh) 一种面向智能终端的gnss高精度定位方法
CN109581452B (zh) 一种gnss参考站载波相位整周模糊度解算方法
CN108415049B (zh) 提高网络rtk双差宽巷模糊度固定正确率的方法
Loyer et al. Zero-difference GPS ambiguity resolution at CNES–CLS IGS Analysis Center
EP1678516B1 (en) Method for using three gps frequencies to resolve carrier-phase integer ambiguities
US6127968A (en) On-the-fly RTK positioning system with single frequency receiver
AU2008260578B2 (en) Distance dependant error mitigation in real-time kinematic (RTK) positioning
CN108828642B (zh) 一种ins辅助bds单频接收机的模糊度快速解算方法
CN109001781B (zh) 一种顾及电离层约束的bds三频模糊度解算方法
CN108445518A (zh) 一种基于双差模糊度固定解约束的gnss精密时间传递方法
CN111965673A (zh) 基于多gnss的单频精密单点定位算法的时间频率传递方法
CN109799520B (zh) 一种基于gnss三频的纯载波相位rtk定位方法
CN116299615B (zh) 一种实现单北斗实时ppp模糊固定的相位偏差估计方法
CN103543454B (zh) 一种嵌入在移动通讯网中的卫星定轨系统
CN105974440B (zh) 一种获取观测噪声的方法和装置
CN116882209B (zh) 一种leo实时定轨随机模型构建方法
CN112731496B (zh) 一种面向智能终端的gnss精密单点定位数据质量控制方法
Bisnath Relative Positioning and Real‐Time Kinematic (RTK)
CN115902968A (zh) 基于北斗三号geo播发增强信息的ppp终端定位方法
CN112230254B (zh) 一种gps载波相位多径误差的校正方法及装置
CN117421933A (zh) 一种多频多系统gnss卫星相位偏差估计方法
CN112444832A (zh) 一种一机多天线接收机的周跳修复方法
CN117665871A (zh) 一种载波距支持下的gnss非差固定解精密钟差确定方法
CN116125513A (zh) 一种gnss模糊度固定解检核方法及系统

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