CN110531389A - Mw组合法与无几何相位法联合探测与修复周跳的方法 - Google Patents
Mw组合法与无几何相位法联合探测与修复周跳的方法 Download PDFInfo
- Publication number
- CN110531389A CN110531389A CN201910694656.0A CN201910694656A CN110531389A CN 110531389 A CN110531389 A CN 110531389A CN 201910694656 A CN201910694656 A CN 201910694656A CN 110531389 A CN110531389 A CN 110531389A
- Authority
- CN
- China
- Prior art keywords
- combination
- frequency
- frequency point
- combined
- cycle slip
- 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.)
- Granted
Links
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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/40—Correcting position, velocity or attitude
- G01S19/41—Differential correction, e.g. DGPS [differential GPS]
-
- 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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种MW组合法与无几何相位法联合探测与修复周跳的方法,包括:步骤S1:选择北斗导航定位系统的三频数据(B1、B2、B3),获取伪距观测值和载波相位观测值;步骤S2:构建三频MW组合观测模型,建立MW方程,进行历元间差分得到周跳检测量△N12和△N13;步骤S3:构建无几何相位观测模型;步骤S4:根据组合观测数据的电离层系数和观测噪声同时相对小的原则选取其中的两组参数,再根据△N1=△N2=△N3时计算两组参数三个频点的周跳,选取满足实际情况的B组合系数;步骤S5:将△N12、△N13和B组合联立得到方程组,计算三频频点上的周跳,对周跳进行修复。本发明解决了无几何法无法探测不敏周跳和MW法无法探测相等周跳的缺陷。
Description
技术领域
本发明属于全球导航定位系统GNSS数据处理技术领域,具体为一种MW组合法与无几何相位法联合探测与修复周跳的方法。
背景技术
随着北斗产业的快速发展和应用领域的不断扩大,人们对于北斗卫星导航系统定位精度的要求也越来也高,北斗卫星导航系统实现高精度定位的前提在于获取高质量的载波相位测量值,利用载波相位观测量进行定位的精度能达到厘米级甚至亚毫米级,可满足更高精度的定位要求。然而载波相位观测量中特有的周跳现象将会对卫星定位结果产生很大的影响。因此在利用载波相位观测量进行高精度定位时,必须首先要对北斗卫星导航系统的原始载波相位观测数据进行预处理,以便得到质量较好的观测量,而数据预处理中的关键步骤就是对周跳进行探测与修复。
发明内容
针对现有技术存在的上述问题,本发明提供了一种MW组合法与无几何相位法联合探测与修复周跳的方法,解决了无几何法无法探测不敏周跳和MW法无法探测相等周跳的缺陷;利用三频数据,可以得到三频数据组合的优良特性,比单、双频组合具有更多优良特性的线性组合。
为了实现上述目的,本发明所采用的技术方案是:
一种MW组合法与无几何相位法联合探测与修复周跳的方法,包括如下步骤:步骤S1:选择北斗导航定位系统的三频数据(B1、B2、B3),获取观测数据伪距观测值和载波相位观测值,对观测的数据进行预处理,剔除粗差;步骤S2:用步骤S1获取的观测数据构建三频MW组合观测模型,建立不同频点组合观测量的MW方程,并进行历元间差分得到周跳检测量ΔN12和ΔN13;步骤S3:用步骤S1获取的观测数据构建无几何相位观测模型;步骤S4:选取多组组合系数,进行解算,根据组合观测数据的电离层系数和观测噪声同时相对小的原则选取其中的两组参数,再根据ΔN1=ΔN2=ΔN3时计算两组参数时的周跳,选取满足实际情况的B组合系数;步骤S5:将步骤S2中的ΔN12、ΔN13和B组合联立得到方程组,计算三频频点上的周跳,并对周跳进行修复。
步骤S2中,B1频点和B2频点的MW组合为L12,L12满足下列公式:
其中:f1、f2:分别为B1频点和B2频点的载波频率;
L1、L2:分别为B1频点和B2频点以米(m)为单位的载波相位观测量,其中
P1、P2:分别为B1频点和B2频点的伪距观测量;
λ1、λ2:别为B1频点和B2频点的波长;
分别是B1频点和B2频点的载波相位观测值。
MW组合的周跳检测量为:
将方程在相邻历元间求差,得到:
ΔN12=ΔN1-ΔN2
其中:ΔN1、ΔN2:分别是B1、B2频点的周跳;
N1、N2:分别为B1频点和B2频点的整周模糊度;
λ12:B1频点和B2频点的MW组合后的波长。
ΔN13的原理与ΔN12相同,且ΔN13=ΔN1-ΔN3。
步骤S3包括如下两步:
步骤S31:对步骤S1获取的观测数据建立三频载波相位观测方程;
根据载波频率与电离层延迟的关系可以得到:p1=(f1/f2)2,p2=(f1/f3)2,其中:
L1、L2、L3分别为三种频率信号的载波相位观测值,所述载波相位观测值以米为单位;
λ1、λ2、λ3:分别为三频信号的波长;
f1、f2、f3:分别为三频信号的频率;
N1、N2、N3:分别为三频的整周模糊度;
分别表示三频的测量误差;
r:接收机到卫星的几何距离;
T:对流层延迟;
I1:电离层延迟;
δu:接收钟差距离;
δs:卫星钟差距离;
B3频点的载波相位观测值;
p1、p2:分别是B2频点电离层系数和B3频点的电离层系数;
步骤S32:对步骤S31的三频载波相位观测方程进行线性组合,得到三频载波相位组合的观测方程:
其中,
式中:α、β、γ:分别表示组合的系数;
N:组合的整周模糊度;
组合的电离层系数;
组合的测量误差;
λ:组合的波长;
组合的载波相位值;
当α+β+γ=1时,三频载波相位组合的几何距离保持不变;当α+β+γ=0时,三频载波相位组合的几何距离被消除,此时可得到三频无几何相位组合。
步骤S4包括如下具体步骤:
步骤S41:选取多组系数,选取的多组系数均满足条件:α+β+γ=0,对多组系数进行解算;
步骤S42:根据解算结果,根据组合后的电离层系数和观测噪声同时小的原则选取其中的两组参数B、C,再根据ΔN1=ΔN2=ΔN3计算两组参数的周跳检测量,选取符合周跳数值范围的B组合,组合B的周跳检测量为:
ΔN=a2λ1ΔN1+b2λ2ΔN2+c2λ3ΔN3
式中,a2、b2、c2:分别表示B组合的α、β、γ的值,且其中a2+b2+c2=0。
步骤S5中包括以下步骤:
步骤S51:将ΔN12和ΔN13组合与B组合联立得到方程组
步骤S52:将每个组合求得的周跳值代入步骤S51中的方程组中进行求解,对所得的解进行取整,便可得到相应频点的载波相位观测量中的周跳值。
与现有技术相比,本发明的有益效果是:联合B1、B2频点的MW法,B1、B3频点的MW法,以及一组最佳组合系数的三频无几何法进行周跳的解算,解决了无几何法无法探测不敏周跳和MW法无法探测相等周跳的缺陷;利用三频数据,可以得到三频数据组合的优良特性,比单、双频组合具有更多优良特性的线性组合。
附图说明
图1为本发明一个实施例的算法流程图;
图2为本发明一个实施例的干净数据时B1频点和B2频点的MW组合探测结果;
图3为本发明一个实施例的干净数据时三频无几何相位组合法的探测结果;
图4为本发明一个实施例的原始数据B1B2频率的MW组合探测结果;
图5为本发明一个实施例的原始数据B1B3频率的MW组合探测结果;
图6为本发明一个实施例的原始数据三频无几何相位组合探测结果。
具体实施方式
下面结合实施例对本发明提供的MW组合法与无几何相位法联合探测与修复周跳的方法作进一步详细、完整地说明。下面描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
一种MW组合法与无几何相位法联合探测与修复周跳的方法,包括如下步骤:
步骤S1:选择北斗导航定位系统的三频数据(B1、B2、B3),获取观测数据伪距观测值和载波相位观测值,对观测的数据进行预处理,剔除粗差;
步骤S2:用步骤S1获取的观测数据构建三频MW组合观测模型,建立不同频点组合观测量的MW方程,并进行历元间差分得到周跳检测量ΔN12和ΔN13;
步骤S3:用步骤S1获取的观测数据构建无几何相位观测模型;
步骤S4:选取多组组合系数,进行解算,根据组合观测数据的电离层系数和观测噪声同时相对小的原则选取其中的两组参数,再根据ΔN1=ΔN2=ΔN3时计算两组参数时的周跳量,选取满足实际情况的B组合系数;
步骤S5:将步骤S2中的ΔN12、ΔN13和B组合联立得到方程组,计算三频频点上的周跳,并对周跳进行修复。
所述定位系统是北斗卫星导航系统,相应的,所述B1、B2、B3三个频点为北斗导航定位系统的三频数据。
步骤S2中所述的MW(Melbourne-Wubbena)组合具有较长的波长,消除了站星间的几何距离以及电离层项,减弱了多路径效应和接收机噪声影响,理论计算结果中只有宽巷模糊度。
步骤S2中,B1频点和B2频点的MW组合为L12,L12满足下列公式:
其中:
f1、f2:分别为B1频点和B2频点的载波频率;
L1、L2:分别为B1频点和B2频点以米(m)为单位的载波相位观测量,其中
P1、P2:分别为B1频点和B2频点的伪距观测量;
λ1、λ2:别为B1频点和B2频点的波长;
分别是B1频点和B2频点的载波相位观测值。
MW组合的周跳检测量为:
将方程在相邻历元间求差,得到:
ΔN12=ΔN1-ΔN2 (3)
其中:
ΔN1、ΔN2:分别是B1、B2频点的周跳;
N1、N2:分别为B1频点和B2频点的整周模糊度。
由上可知,MW组合法无法检测出B1频点和B2频点在同一历元发生相同大小周跳的情况。
ΔN13的原理与ΔN12相同,用上述同样的方法获取B1频点和B3频点的MW组合并计算ΔN13,得ΔN13=ΔN1-ΔN3,其中,ΔN3是B3频点的周跳,求取过程和ΔN12相同,在此不再赘述。
步骤S3包括如下两步:
步骤S31:对步骤S1获取的观测数据建立三频载波相位观测方程;
步骤S32:对步骤S31的三频载波相位观测方程进行线性组合,得到三频载波相位组合的观测方程。
基于北斗卫星导航系统同时播发的三个频点的卫星信号,从中可以选择比单、双频组合具有更多优良特性的线性组合,对周跳探测与修复的准确性更高。
步骤S31中,三频线性组合方程的基本原理如下:
根据北斗载波频率与电离层延迟的关系可以得到:p1=(f1/f2)2,p2=(f1/f3)2,其中:
L为以米为单位的载波相位观测值,下标分别表示北斗三种频率信号,即L1、L2、L3分别为三种频率信号的载波相位观测值,所述载波相位观测值以米为单位;
λ1、λ3、λ3:分别为三频信号的波长;
f1、f2、f3:分别为三频信号的频率;
N1、N2、N3:分别为三频的整周模糊度;
分别表示三频的测量误差;
r:接收机到卫星的几何距离;
T:对流层延迟;
I1:电离层延迟;
δu:接收钟差距离;
δs:卫星钟差距离;
p1、p2:分别为是B2频点电离层系数和B3频点的电离层系数。
步骤S32时,将以上北斗三个频率的载波相位观测方程进行线性组合,可以得到北斗三频载波相位组合的观测方程:
其中,
式中:
α、β、γ:分别表示组合的系数;
N:组合的整周模糊度;
N3:B3频点的整周模糊度;
组合的电离层系数;
组合的测量误差;
λ:组合的波长;
组合的载波相位值。
由以上分析可知,当α+β+γ=1时,三频载波相位组合的几何距离保持不变;当α+β+γ=0时,三频载波相位组合的几何距离被消除,此时可得到三频无几何相位组合:
将方程(9)在相邻历元间求差后,便可得到三频相位无几何距离组合的周跳检测量
ΔN=αλ1ΔN1+βλ2ΔN2+γλ3ΔN3
三频无几何相位组合探测周跳的标准为:
|αλ1ΔN1+βλ2ΔN2+γλ3ΔN3|≥lσΔN (10)
其中:
l为常数,取l=3,lσΔN为探测阈值。根据误差传播定律,取周,则无几何相位法周跳检测量的标准差可表示为:
为了能够达到较高的探测灵敏度,三频无几何相位系数应满足:
其中,在三频无几何相位组合法进行探测周跳的过程中,探测不出满足αλ1ΔN1+βλ2ΔN2+γλ3ΔN3=0的周跳。
三频无几何相位法消去了卫星到测站接收机之间的站星距离,其周跳检测量主要受组合观测量的电离层延迟误差以及观测噪声的影响。
步骤S4包括如下两步:
步骤S41:选取多组系数,选取的多组系数均满足条件:α+β+γ=0,对多组系数进行解算。
本实施例选取了A、B、C、D四组系数。
步骤S41中,在选取组合系数时,事先对选取的四组系数进行解算,解算的结果如表1所示:
表1:四组系数解算结果表
步骤S42:根据解算结果,根据组合后的电离层系数和观测噪声同时小的原则选取其中的两组参数B、C,再根据ΔN1=ΔN2=ΔN3条件计算两组参数的周跳检测量,选取满足实际情况的B组合系数,组合B的周跳检测量为:
ΔN=a2λ1ΔN1+b2λ2ΔN2+c2λ3ΔN3
式中,a2、b2、c2:分别表示B组合的α、β、γ的值,且其中a2+b2+c2=0
具体的,对于本实施例,步骤S42包括如下三步:
步骤S421:对步骤S41的解算结果进行评估,根据组合后的电离层系数和观测噪声同时小的原则选取其中的两组参数B、C。
对解算结果进行评估的依据为电离层系数和观测噪声σΔN的值,由表1可知,组合B的电离层系数和观测噪声σΔN均比较小,组合C有最小的电离层系数但观测噪声σΔN较大,两种组合各有优势,但由于MW组合法无法检测出ΔN1=ΔN2=ΔN3的情况,若此时三频无几何相位法的周跳检测量也在阈值范围内,即|ΔN|<3σΔN,则两种方法均会失效。
步骤S422:对步骤S421选取的两组系数B和C,分别计算在ΔN1=ΔN2=ΔN3时各自的周跳检测量。
对于组合C,周跳检测量为:
ΔN=a1λ1ΔN1+b1λ2ΔN2+c1λ3ΔN3 (13)
其中:a1+b1+c1=0,a1、b1、c1:分别表示C组合的α、β、γ的值,
将ΔN1=ΔN2=ΔN3代入上式(13),若要满足则有:
-6.260<ΔN1=ΔN2=ΔN3<6.260 (14)
由上可知,当北斗三个频点的载波相位观测量中同时发生(-6,6)周范围内的周跳时,组合C对于周跳的探测也会失效。
对于组合B,周跳检测量为:
ΔN=a2λ1ΔN1+b2λ2ΔN2+c2λ3ΔN3 (15)
其中:a2+b2+c2=0。将ΔN1=ΔN2=ΔN3代入上式,若要满足则有:
-0.743<ΔN1=ΔN2=ΔN3<0.743 (16)
由上可知,北斗三个频点的载波相位观测量中发生的周跳值应该都相等且绝对值小于0.743。
本领域技术人员认同一周是最小的周跳,即周跳是整数。由上述计算可知,对于B组合,在ΔN1=ΔN2=ΔN3,且周跳检测量在上述阈值范围内时,计算的三个频点的周跳绝对值小于1,达不到最小的周跳,即这种情况不存在,不满足现实情况。换句话说,对于B组合,在ΔN1=ΔN2=ΔN3,其周跳检测量不在上述阈值范围内,此时,B组合对应的检测方法有效。
步骤S423:针对步骤S422的计算结果,选取满足实际情况的B组合系数,作为方程组中的联立方程,进行周跳的探测与修复。
步骤S5中包括以下步骤:
步骤S51:将ΔN12和ΔN13组合与B组合联立得到方程组,
步骤S52:将每个组合求得的周跳值代入步骤S51中的方程组中进行求解,对所得的解进行取整,便可得到相应频点的载波相位观测量中的周跳值。
最后,对上述方法进行验证,为了验证方法的有效性,包括如下三步:
第一步:采用一组采样率为1s的载波相位观测数据和伪距观测数据对提出的方法进行有效性分析,本数据是不含周跳的干净数据;
第二步:检测另一组采样率为1s的原始载波相位观测数据和伪距观测数据;
第三步:将第二步中求得的周跳值代入步骤S51中的方程组进行求解,对所得的解进行取整便可得到相应频点的载波相位观测量中的周跳值。
本实施例中,第一步中用MW法和三频无几何方法分别进行周跳探测,图2和图3是实验的结果。第二步的实验结果如图4,图5,图6所示。由图4,图5,图6的探测结果可以看出,在第100历元处的小周跳和250历元处的大周跳均能被两种方法检测出来;对于第150历元处的相等周跳,MW的两个组合均探测失效,此时三频无几何相位法可以探测出;在第250历元处的不敏周跳,三频无几何相位法探测失效,此时MW的两个组合均可以探测出。将以上各个历元中三个组合求得的周跳值代入步骤S51中的方程组进行求解,对所得的解进行取整便可得到相应频点的载波相位观测量中的周跳值。
本实施例联合B1、B2频点的MW法,B1、B3频点的MW法,以及一组最佳组合系数的三频无几何法进行周跳的解算,计算上简洁、快速、有效,且周跳探测效果明显,可以很好解决无几何法无法探测不敏周跳和MW法无法探测相等周跳的缺陷。
最后有必要在此说明的是:以上实施例只用于对本发明的技术方案作进一步详细地说明,不能理解为对本发明保护范围的限制,本领域的技术人员根据本发明的上述内容做出的一些非本质的改进和调整均属于本发明的保护范围。
Claims (5)
1.一种MW组合法与无几何相位法联合探测与修复周跳的方法,其特征在于,包括如下步骤:
步骤S1:选择北斗导航定位系统的三频数据(B1、B2、B3),获取观测数据伪距观测值和载波相位观测值,对观测的数据进行预处理,剔除粗差;
步骤S2:用步骤S1获取的观测数据构建三频MW组合观测模型,建立不同频点组合观测量的MW方程,并进行历元间差分得到周跳检测量△N12和△N13;
步骤S3:用步骤S1获取的观测数据构建无几何相位观测模型;
步骤S4:选取多组组合系数,进行解算,根据组合观测数据的电离层系数和观测噪声同时相对小的原则选取其中的两组参数,再根据△N1=△N2=△N3时计算两组参数三个频点的周跳,选取满足实际情况的B组合系数;
步骤S5:将步骤S2中的△N12、△N13和B组合联立得到方程组,计算三频频点上的周跳,并对周跳进行修复。
2.根据权利要求1所述的MW组合法与无几何相位法联合探测与修复周跳的方法,其特征在于:步骤S2中,B1频点和B2频点的MW组合为L12,L12满足下列公式:
其中:f1、f2:分别为B1频点和B2频点的载波频率;
L1、L2:分别为B1频点和B2频点以米(m)为单位的载波相位观测量,其中
P1、P2:分别为B1频点和B2频点的伪距观测量;
λ1、λ2:别为B1频点和B2频点的波长;
分别是B1频点和B2频点的载波相位观测值;
MW组合的周跳检测量为:
将方程在相邻历元间求差,得到:
△N12=△N1-△N2
其中:△N1、△N2:分别是B1、B2频点的周跳;
N1、N2:分别为B1频点和B2频点的整周模糊度;
λ12:B1频点和B2频点的MW组合后的波长;
△N13的原理与△N12相同,且△N13=△N1-△N3。
3.根据权利要求2所述的MW组合法与无几何相位法联合探测与修复周跳的方法,其特征在于:步骤S3包括如下两步:
步骤S31:对步骤S1获取的观测数据建立三频载波相位观测方程;
根据载波频率与电离层延迟的关系可以得到:p1=(f1/f2)2,p2=(f1/f3)2,其中:
L1、L2、L3分别为三种频率信号的载波相位观测值,所述载波相位观测值以米为单位;
λ1、λ2、λ3:分别为三频信号的波长;
f1、f2、f3:分别为三频信号的频率;
N1、N2、N3:分别为三频的整周模糊度;
分别表示三频的测量误差;
r:接收机到卫星的几何距离;
T:对流层延迟;
I1:电离层延迟;
δu:接收钟差距离;
δs:卫星钟差距离;
B3频点的载波相位观测值;
p1、p2:分别是B2频点电离层系数和B3频点的电离层系数;
步骤S32:对步骤S31的三频载波相位观测方程进行线性组合,得到三频载波相位组合的观测方程:
其中,
式中:α、β、γ:分别表示组合的系数;
N:组合的整周模糊度;
组合的电离层系数;
组合的测量误差;
λ:组合的波长;
组合的载波相位值;
当α+β+γ=1时,三频载波相位组合的几何距离保持不变;当α+β+γ=0时,三频载波相位组合的几何距离被消除,此时可得到三频无几何相位组合。
4.根据权利要求3所述的MW组合法与无几何相位法联合探测与修复周跳的方法,其特征在于:步骤S4包括如下具体步骤:
步骤S41:选取多组系数,选取的多组系数均满足条件:α+β+γ=0,对多组系数进行解算;
步骤S42:根据解算结果,根据组合后的电离层系数和观测噪声同时小的原则选取其中至少两组参数,设为B组合和C组合,再根据△N1=△N2=△N3计算两组参数的周跳,选取符合周跳值范围的组合,设组合B符合周跳值范围,则组合B的周跳检测量为:
△N=a2λ1△N1+b2λ2△N2+c2λ3△N3
式中,a2、b2、c2:分别表示B组合的α、β、γ的值,且其中a2+b2+c2=0。
5.根据权利要求4所述的MW组合法与无几何相位法联合探测与修复周跳的方法,其特征在于:步骤S5中包括以下步骤:
步骤S51:将△N12和△N13组合与B组合联立得到方程组
步骤S52:将每个组合求得的周跳值代入步骤S51中的方程组中进行求解,对所得的解进行取整,便可得到相应频点的载波相位观测量中的周跳值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910694656.0A CN110531389B (zh) | 2019-07-30 | 2019-07-30 | Mw组合法与无几何相位法联合探测与修复周跳的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910694656.0A CN110531389B (zh) | 2019-07-30 | 2019-07-30 | Mw组合法与无几何相位法联合探测与修复周跳的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110531389A true CN110531389A (zh) | 2019-12-03 |
CN110531389B CN110531389B (zh) | 2022-11-29 |
Family
ID=68661084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910694656.0A Active CN110531389B (zh) | 2019-07-30 | 2019-07-30 | Mw组合法与无几何相位法联合探测与修复周跳的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110531389B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111239779A (zh) * | 2020-03-03 | 2020-06-05 | 东南大学 | 一种无盲点的gnss三频组合周跳探测与修复方法 |
CN111505667A (zh) * | 2020-03-21 | 2020-08-07 | 哈尔滨工程大学 | 一种基于动对动平台多径及观测噪声异常完好性检测方法 |
CN113671545A (zh) * | 2021-08-27 | 2021-11-19 | 中国科学院国家授时中心 | 一种星地精密时间同步及载波相位周跳探测方法 |
CN114814907A (zh) * | 2022-04-07 | 2022-07-29 | 山东大学 | 一种基于弧段划分的北斗双频周跳探测与修复方法 |
CN115267847A (zh) * | 2022-06-27 | 2022-11-01 | 湘潭大学 | 一种单频周跳探测与修复方法、系统、设备及介质 |
CN115932922A (zh) * | 2022-12-28 | 2023-04-07 | 辽宁工程技术大学 | 一种基于bds四频数据的周跳探测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5935194A (en) * | 1994-02-18 | 1999-08-10 | Trimble Navigation Limited | Method for using external constraints to improve the speed and reliability of phase ambiguity resolution in real-time kinematic initialization |
US20140002299A1 (en) * | 2007-06-22 | 2014-01-02 | Trimble Navigation Limited | Combined cycle slip indicators for regionally augmented gnss |
CN105467412A (zh) * | 2015-12-04 | 2016-04-06 | 昆明理工大学 | 一种北斗三频周跳探测与修复方法 |
-
2019
- 2019-07-30 CN CN201910694656.0A patent/CN110531389B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5935194A (en) * | 1994-02-18 | 1999-08-10 | Trimble Navigation Limited | Method for using external constraints to improve the speed and reliability of phase ambiguity resolution in real-time kinematic initialization |
US20140002299A1 (en) * | 2007-06-22 | 2014-01-02 | Trimble Navigation Limited | Combined cycle slip indicators for regionally augmented gnss |
CN105467412A (zh) * | 2015-12-04 | 2016-04-06 | 昆明理工大学 | 一种北斗三频周跳探测与修复方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111239779A (zh) * | 2020-03-03 | 2020-06-05 | 东南大学 | 一种无盲点的gnss三频组合周跳探测与修复方法 |
CN111505667A (zh) * | 2020-03-21 | 2020-08-07 | 哈尔滨工程大学 | 一种基于动对动平台多径及观测噪声异常完好性检测方法 |
CN111505667B (zh) * | 2020-03-21 | 2023-09-29 | 哈尔滨工程大学 | 一种基于动对动平台多径及观测噪声异常完好性检测方法 |
CN113671545A (zh) * | 2021-08-27 | 2021-11-19 | 中国科学院国家授时中心 | 一种星地精密时间同步及载波相位周跳探测方法 |
CN114814907A (zh) * | 2022-04-07 | 2022-07-29 | 山东大学 | 一种基于弧段划分的北斗双频周跳探测与修复方法 |
CN115267847A (zh) * | 2022-06-27 | 2022-11-01 | 湘潭大学 | 一种单频周跳探测与修复方法、系统、设备及介质 |
CN115267847B (zh) * | 2022-06-27 | 2024-04-09 | 湘潭大学 | 一种单频周跳探测与修复方法、系统、设备及介质 |
CN115932922A (zh) * | 2022-12-28 | 2023-04-07 | 辽宁工程技术大学 | 一种基于bds四频数据的周跳探测方法 |
CN115932922B (zh) * | 2022-12-28 | 2024-04-09 | 辽宁工程技术大学 | 一种基于bds四频数据的周跳探测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110531389B (zh) | 2022-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110531389A (zh) | Mw组合法与无几何相位法联合探测与修复周跳的方法 | |
CN106168672B (zh) | 一种gnss多模单频rtk周跳探测方法及装置 | |
Geng et al. | Improving the estimation of fractional-cycle biases for ambiguity resolution in precise point positioning | |
CN105467412B (zh) | 一种北斗三频周跳探测与修复方法 | |
CN101419275B (zh) | 基于多接收机的局域机场监视方法和系统 | |
CN106066485B (zh) | 一种周跳探测方法及装置 | |
CN103941272B (zh) | Gps、glonass和bds联合解算的定位方法 | |
CN110687556B (zh) | 一种适用于laas的多径误差模型化方法 | |
CN109001781B (zh) | 一种顾及电离层约束的bds三频模糊度解算方法 | |
CN105842719B (zh) | 一种顾及对流层影响的cors基准站网基线模糊度解算方法 | |
CN109738917A (zh) | 一种北斗变形监测中的多路径误差削弱方法及装置 | |
CN104614741B (zh) | 一种不受glonass码频间偏差影响的实时精密卫星钟差估计方法 | |
CN106125113B (zh) | 一种利用多系统gnss观测值的高精度基线解算方法 | |
CN104459722B (zh) | 一种基于多余观测分量的整周模糊度可靠性检验方法 | |
CN105738926A (zh) | 一种glonass系统接收机间相位频间偏差标定方法 | |
CN106842236B (zh) | Gnss接收机周跳探测与修复处理装置 | |
CN110018507A (zh) | 一种基于星座间作差的组合精密单点定位方法及系统 | |
CN109782306A (zh) | 一种高精度卫星导航接收机性能分析方法 | |
CN107728181A (zh) | 一种实时周跳探测修复方法 | |
Paziewski et al. | High-precision multi-constellation GNSS. Methods, selected applications and challenges | |
CN106292265A (zh) | 一种基于导航卫星的多地时间同步方法 | |
CN105549046B (zh) | Gnss接收机周跳探测与修复处理方法 | |
CN104898144A (zh) | 一种北斗系统参考站整周模糊度单历元确定方法 | |
CN109143289B (zh) | 一种gnss单站位移监测方法 | |
CN110646823A (zh) | 一种基于Helmet后验定权法的GPS\BDS紧组合精密单点定位方法 |
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 |