CN109444935B - 一种低采样率的多普勒周跳探测和修复方法 - Google Patents
一种低采样率的多普勒周跳探测和修复方法 Download PDFInfo
- Publication number
- CN109444935B CN109444935B CN201811205995.XA CN201811205995A CN109444935B CN 109444935 B CN109444935 B CN 109444935B CN 201811205995 A CN201811205995 A CN 201811205995A CN 109444935 B CN109444935 B CN 109444935B
- Authority
- CN
- China
- Prior art keywords
- value
- doppler
- cycle slip
- formula
- carrier phase
- 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
Images
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/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
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
本发明公开了一种低采样率的多普勒周跳探测和修复方法,其特征在于,包括如下步骤:1)获取L1频段的观测数据;2)计算拉格朗日内插值;3)载波相位内插值进行差分运算;4)载波载波相位差值进行多项式拟合运算;5)求解多项式系数;6)计算载波相位外推值;7)获得多普勒外推值;8)获得多项式外推值多普勒周跳;9)获得拉格朗日内插多普勒周跳;10)多项式外推多普勒周跳与拉格朗日内插多普勒周跳做差;11)设置门限;12)判断周跳;13)周跳修复。这种方法不仅可以探测出小周跳,并且在采样率比较低的情况下能提高探测周跳的精度。
Description
技术领域
本发明应用于北斗未定导航定位领域,具体针对于滑坡形变监测高精度定位的改进周跳检测算法,尤其是一种低采样率的多普勒周跳探测和修复方法。
背景技术
目前,周跳问题广泛存在于卫星导航定位等工程测量中,周跳会破坏卫星至接收机的距离观测值的变化规律,对定位产生影响,并且在高精度定位中,电离层、对流层、伪距和多路径效应等误差也会对周跳的探测产生很大的影响。一直以来,提出有效且实用的周跳探测方法是载波相位数据处理中周跳探测的难点,只有有效地探测出周跳,才能保证载波相位周跳修复以及整周模糊度固定的精确性。在卫星导航定位的周跳探测与修复中,多普勒观测值免受周跳的影响,是一种非常稳定并独立于载波相位的观测值,可辅助载波相位数据进行周跳探测。但随着采样率的降低,历元间各项误差的相关性也降低,周跳探测能力显著减小。
发明内容
本发明的目的是针对现有技术的不足,而提供一种低采样率的多普勒周跳探测和修复方法。这种方法不仅可以探测出小周跳,并且在采样率比较低的情况下能提高探测周跳的精度。
实现本发明目的的技术方案是:
一种低采样率的多普勒周跳探测和修复方法,与现有技术部同在于,包括如下步骤:
2)计算拉格朗日内插值:采用拉格朗日插值方程将采样率为5s的载波相位观测值插值到1s采样间隔,获得载波相位内插值为公式(1):
式中,在选取的观测历元区间[ti,tj]中,t为插值历元,为观测历元ti对应的载波相位观测值,Ln为插值历元t对应的载波相位内插值,n为拉格朗日插值阶数,本文取值为7,同理,采用拉格朗日插值方程将采样率为5s的多普勒观测值插值到1s采样间隔,获得多普勒内插值为公式(2):
式中,di为观测历元ti对应的载波相位观测值,Dn为插值历元t对应的多普勒内插值;
3)载波相位内插值进行差分运算:对载波相位内插值的相邻两个观测历元按照公式(3)进行历元间的差分运算:
ΔLn(ti)=Ln(ti+1)-Ln(ti),(i=0,1,…,n) (3),
式中,ΔLn为历元ti对应的载波相位差值;
4)载波载波相位差值进行多项式拟合运算:将步骤3)中的载波相位差值ΔLn(t)进行多项式拟合运算,按照公式(4)将无周跳的m个内插值拟合成一个k阶多项式:
式中,ti为拟合历元的观测时间,t0为初始时间,a0,a1,…,ak为多项式系数,本文中m取值为14,k取值为4。
5)求解多项式系数:采用最小二乘法求解多项式系数a0,a1,…,ak,公式(4)的矩阵形式可表示为:
β=t·α (5),
其中:
则:
α=(tΤt)-1tΤβ (7),
式中,α为多项式系数a0,a1,…,ak的矩阵形式,β为载波相位差值ΔLn(t)的矩阵形式;
6)计算载波相位外推值:依据步骤5)获得的多项式系数外推出第m+1历元的载波相位外推值,如公式(8):
7)获得多普勒外推值:对步骤2)获得的多普勒内插值Dn进行多项式拟合,并外推第m+1历元的多普勒外推值,重复步骤4)获得多普勒外推值多项式拟合方程,如公式(9):
重复步骤5)和步骤6),获得多普勒外推值,其中多项式外推方程为公式(10):
8)获得多项式外推值多普勒周跳:将载波相位外推值和多普勒外推值按照公式(11)进行多普勒积分运算:
9)获得拉格朗日内插多普勒周跳:依据步骤2)中的载波相位内插值和多普勒内插值按照公式(12)进行多普勒积分运算:
式中,ΔNL为拉格朗日内插多普勒周跳;
10)多项式外推多普勒周跳与拉格朗日内插多普勒周跳按照公式(13)做差:
式中,ΔN为差后估计值;
11)设置门限:设置门限值δ;
步骤2)中所述的拉格朗日插值方程的阶数n为7阶,阶数太高或者太低都会引起插值精度的降低,并且阶数为7拉格朗日多项式在很多场合使用较多。
步骤3)中,载波相位差值历元间做差,运用了差分的思想,使得到的观测值中的电离层延迟和对流层延迟基本被消除,其他各种误差的变化值也非常小,了削弱观测误差。
步骤4)中k阶多项式中,m取值为14,k取值为4,多项式拟合的阶数选择中,由于星地距离对时间的四阶导数或者五阶导数一般已趋近于零了,其变化规律是随机的,无法再用多项式拟合了,因此阶数取4阶;拟合窗宽度越大,虽然外推值越准确,同时拟合后的中误差会很小,但导致计算量增大,但是当拟合窗宽度越小,外推的值会越粗糙。
步骤11)中门限值δ=0.5,门限值是根据情况而设定的,当δ较大时,表示只有当观测值偏离外推值很大的时候才认为它是异常的;当δ较小时,表示当观测值离外推值较小的时候就认为它是异常值了,一般δ取0.2到0.9之间的数。
本技术方案的优点在于:
本技术方案提出低采样率的多普勒周跳探测和修复方法,主要包括拉格朗日插值和多项式拟合两部分。利用拉格朗日插值法将采样率较大的多普勒观测值内插到1s的采样间隔中,并对内插值进行多项式拟合,外推出噪声较小的多普勒值,最后将多项式拟合外推多普勒值与内插多普勒值相减得到周跳值。与现有技术相比,步骤3)中载波相位差值历元间做差,运用了差分的思想,使得到的观测值中的电离层延迟和对流层延迟基本被消除,其他各种误差的变化值也非常小,了削弱观测误差。而对于采样频率的问题,利用拉格朗日插值法将低采样频率的载波相位观测值插值到较高的采样频率,有助于提高拟合精度。因此本技术方案不仅可以消除电离层与几何距离产生的噪声干扰,减小对周跳产生的影响,同时提高了多普勒观测值的采样率和周跳的检测精度。
这种方法不仅可以探测出小周跳,并且在采样率比较低的情况下能提高探测周跳的精度。
附图说明
图1为实施例的方法流程示意图。
具体实施方式
下面结合附图和实施例对本发明内容作进一步的阐述,但不是对本发明的限定。
实施例:
参照图1,一种低采样率的多普勒周跳探测和修复方法,包括如下步骤:
2)计算拉格朗日内插值:采用拉格朗日插值方程将采样率为5s的载波相位观测值插值到1s采样间隔,获得载波相位内插值为公式(1):
式中,在选取的观测历元区间[ti,tj]中,t为插值历元,为观测历元ti对应的载波相位观测值,Ln为插值历元t对应的载波相位内插值,n为拉格朗日插值阶数,同理,采用拉格朗日插值方程将采样率为5s的多普勒观测值插值到1s采样间隔,获得多普勒内插值为公式(2):
式中,di为观测历元ti对应的载波相位观测值,Dn为插值历元t对应的多普勒内插值;
3)载波相位内插值进行差分运算:对载波相位内插值的相邻两个观测历元按照公式(3)进行历元间的差分运算:
ΔLn(ti)=Ln(ti+1)-Ln(ti),(i=0,1,…,n) (3),
式中,ΔLn为历元ti对应的载波相位差值;
4)载波载波相位差值进行多项式拟合运算:将步骤3)中的载波相位差值ΔLn(t)进行多项式拟合运算,按照公式(4)将无周跳的m个内插值拟合成一个k阶多项式:
式中,ti为拟合历元的观测时间,t0为初始时间,a0,a1,…,ak为多项式系数;
5)求解多项式系数:采用最小二乘法求解多项式系数a0,a1,…,ak,公式(4)的矩阵形式可表示为:
β=t·α (5),
其中:
则:
α=(tΤt)-1tΤβ (7),
式中,α为多项式系数a0,a1,…,ak的矩阵形式,β为载波相位差值ΔLn(t)的矩阵形式;
6)计算载波相位外推值:依据步骤5)获得的多项式系数外推出第m+1历元的载波相位外推值,如公式(8):
7)获得多普勒外推值:对步骤2)获得的多普勒内插值Dn进行多项式拟合,并外推第m+1历元的多普勒外推值,重复步骤4)获得多普勒外推值多项式拟合方程,如公式(9):
重复步骤5)和步骤6),获得多普勒外推值,其中多项式外推方程为公式(10):
8)获得多项式外推值多普勒周跳:将载波相位外推值和多普勒外推值按照公式(11)进行多普勒积分运算:
9)获得拉格朗日内插多普勒周跳:依据步骤2)中的载波相位内插值和多普勒内插值按照公式(12)进行多普勒积分运算:
式中,ΔNL为拉格朗日内插多普勒周跳;
10)多项式外推多普勒周跳与拉格朗日内插多普勒周跳按照公式(13)做差:
式中,ΔN为差后估计值;
11)设置门限:设置门限值δ;
本例步骤2)中所述的拉格朗日插值方程的阶数n为7阶,阶数太高或者太低都会引起插值精度的降低。
步骤3)中,载波相位差值历元间做差,运用了差分的思想,使得到的观测值中的电离层延迟和对流层延迟基本被消除,其他各种误差的变化值也非常小,了削弱观测误差。
步骤4)中k阶多项式中,m取值为14,k取值为4,多项式拟合的阶数选择中,由于星地距离对时间的四阶导数或者五阶导数一般已趋近于零了,其变化规律是随机的,无法再用多项式拟合了,因此阶数取4阶;拟合窗宽度越大,虽然外推值越准确,同时拟合后的中误差会很小,但导致计算量增大,但是当拟合窗宽度越小,外推的值会越粗糙,本例中,通过取不同的值进行试验,在周跳探测时取m=14比较合适。
步骤11)中门限值δ=0.5,门限值是根据情况而设定的,当δ较大时,表示只有当观测值偏离外推值很大的时候才认为它是异常的;当δ较小时,表示当观测值离外推值较小的时候就认为它是异常值了,一般δ取0.2到0.9之间的数,本例中,经多次取值试验,δ=0.5可行。
Claims (4)
1.一种低采样率的多普勒周跳探测和修复方法,其特征在于,包括如下步骤:
2)计算拉格朗日内插值:采用拉格朗日插值方程将采样率为5s的载波相位观测值插值到1s采样间隔,获得载波相位内插值为公式(1):
式中,在选取的观测历元区间[ti,tj]中,t为插值历元,为观测历元ti对应的载波相位观测值,Ln为插值历元t对应的载波相位内插值,n为拉格朗日插值阶数,同理,采用拉格朗日插值方程将采样率为5s的多普勒观测值插值到1s采样间隔,获得多普勒内插值为公式(2):
式中,di为观测历元ti对应的载波相位观测值,Dn为插值历元t对应的多普勒内插值;
3)载波相位内插值进行差分运算:对载波相位内插值的相邻两个观测历元按照公式(3)进行历元间的差分运算:
ΔLn(ti)=Ln(ti+1)-Ln(ti),(i=0,1,…,n) (3),
式中,ΔLn为历元ti对应的载波相位差值;
4)载波相位差值进行多项式拟合运算:将步骤3)中的载波相位差值ΔLn(t)进行多项式拟合运算,按照公式(4)将无周跳的m个内插值拟合成一个k阶多项式:
式中,ti为拟合历元的观测时间,t0为初始时间,a0,a1,…,ak为多项式系数;
5)求解多项式系数:采用最小二乘法求解多项式系数a0,a1,…,ak,公式(4)的矩阵形式可表示为:
β=t·α (5),
其中:
则:
α=(tΤt)-1tΤβ (7),
式中,α为多项式系数a0,a1,…,ak的矩阵形式,β为载波相位差值ΔLn(t)的矩阵形式;
6)计算载波相位外推值:依据步骤5)获得的多项式系数外推出第m+1历元的载波相位外推值,如公式(8):
7)获得多普勒外推值:对步骤2)获得的多普勒内插值Dn进行多项式拟合,并外推第m+1历元的多普勒外推值,重复步骤4)获得多普勒外推值多项式拟合方程,如公式(9):
重复步骤5)和步骤6),获得多普勒外推值,其中多项式外推方程为公式(10):
8)获得多项式外推值多普勒周跳:将载波相位外推值和多普勒外推值按照公式(11)进行多普勒积分运算:
9)获得拉格朗日内插多普勒周跳:依据步骤2)中的载波相位内插值和多普勒内插值按照公式(12)进行多普勒积分运算:
式中,ΔNL为拉格朗日内插多普勒周跳;
10)多项式外推多普勒周跳与拉格朗日内插多普勒周跳按照公式(13)做差:
式中,ΔN为差后估计值;
11)设置门限:设置门限值δ;
2.根据权利要求1所述的低采样率的多普勒周跳探测和修复方法,其特征在于,步骤2)中所述的拉格朗日插值方程的阶数n为7阶。
3.根据权利要求1所述的低采样率的多普勒周跳探测和修复方法,其特征在于,步骤4)中k阶多项式中,m取值为14,k取值为4。
4.根据权利要求1所述的低采样率的多普勒周跳探测和修复方法,其特征在于,步骤11)中门限值δ=0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811205995.XA CN109444935B (zh) | 2018-10-17 | 2018-10-17 | 一种低采样率的多普勒周跳探测和修复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811205995.XA CN109444935B (zh) | 2018-10-17 | 2018-10-17 | 一种低采样率的多普勒周跳探测和修复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109444935A CN109444935A (zh) | 2019-03-08 |
CN109444935B true CN109444935B (zh) | 2022-10-21 |
Family
ID=65546647
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811205995.XA Active CN109444935B (zh) | 2018-10-17 | 2018-10-17 | 一种低采样率的多普勒周跳探测和修复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109444935B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824520B (zh) * | 2019-11-22 | 2021-11-09 | 中国科学院国家空间科学中心 | 一种实时计算电离层相位闪烁指数计算方法及系统 |
CN110967717A (zh) * | 2019-12-23 | 2020-04-07 | 合肥工业大学 | 一种基于小波变换法的周跳探测和修复方法 |
CN111323795B (zh) * | 2020-03-20 | 2022-03-22 | 湖南联智科技股份有限公司 | 一种北斗变形监测中多路径误差的削弱方法 |
CN111397645B (zh) * | 2020-04-06 | 2020-12-18 | 华中科技大学 | 一种相位差分解调方法及系统 |
CN112130806B (zh) * | 2020-09-25 | 2023-10-24 | 锐迪科创微电子(北京)有限公司 | 相位升采样方法及装置、计算机可读存储介质 |
CN113866808B (zh) * | 2021-09-29 | 2024-08-02 | 哈尔滨工程大学 | 一种gnss振荡器级实时周跳探测与修复方法 |
CN114363819B (zh) * | 2022-01-10 | 2024-01-09 | 中国人民解放军国防科技大学 | 一种地基高精度区域定位导航系统周跳修复方法 |
Citations (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5323322A (en) * | 1992-03-05 | 1994-06-21 | Trimble Navigation Limited | Networked differential GPS system |
JPH0777570A (ja) * | 1993-07-16 | 1995-03-20 | Caterpillar Inc | 衛星基地局航法システムにおけるサイクルスリップ検出方法及び装置 |
JPH07190769A (ja) * | 1993-12-27 | 1995-07-28 | Sokkia Co Ltd | Gps干渉測位方法 |
US5914685A (en) * | 1997-04-25 | 1999-06-22 | Magellan Corporation | Relative position measuring techniques using both GPS and GLONASS carrier phase measurements |
US6127968A (en) * | 1998-01-28 | 2000-10-03 | Trimble Navigation Limited | On-the-fly RTK positioning system with single frequency receiver |
WO2001022111A1 (en) * | 1999-09-24 | 2001-03-29 | Thales Geosolutions Group Limited | Global navigation satellite systems and methods |
CN1592855A (zh) * | 2001-10-22 | 2005-03-09 | 高通股份有限公司 | 开环跟踪gps信号的方法 |
EP1729145A1 (en) * | 2005-06-02 | 2006-12-06 | Gmv, S.A. | Method and system for providing GNSS navigation position solution with guaranteed integrity in non-controlled environments |
JP2008175730A (ja) * | 2007-01-19 | 2008-07-31 | Churyo Eng Kk | 移動局の速度計測装置 |
CN101449178A (zh) * | 2006-05-30 | 2009-06-03 | 诺基亚公司 | 支持相对定位 |
JP2010071686A (ja) * | 2008-09-16 | 2010-04-02 | Sumitomo Electric Ind Ltd | 測位装置、コンピュータプログラム及び測位方法 |
CN102035774A (zh) * | 2009-09-25 | 2011-04-27 | 上海伽利略导航有限公司 | 一种单频载波相位定位的方法及装置 |
CN102288978A (zh) * | 2011-07-20 | 2011-12-21 | 东南大学 | 一种cors基站周跳探测与修复方法 |
CN102565821A (zh) * | 2011-12-22 | 2012-07-11 | 浙江大学 | 多普勒频偏辅助的卫星导航信号载波周跳实时探测与修复的方法 |
WO2013016800A1 (en) * | 2011-07-29 | 2013-02-07 | Baseband Technologies Inc. | System, method, and computer program for a low power and low cost gnss receiver |
CN103576167A (zh) * | 2013-10-28 | 2014-02-12 | 中国科学院国家授时中心 | 基于hht和支持向量机的周跳探测与修复方法 |
JP2015068768A (ja) * | 2013-09-30 | 2015-04-13 | 日本電気株式会社 | 測位システムと装置と方法並びにプログラム |
CN104570013A (zh) * | 2014-12-30 | 2015-04-29 | 北京无线电计量测试研究所 | 一种用于频率驯服的实时gps载波相位周跳的探测方法 |
EP2966477A1 (en) * | 2014-07-09 | 2016-01-13 | ANavS GmbH | Method for determining the position and attitude of a moving object using low-cost receivers |
CN106168672A (zh) * | 2016-01-01 | 2016-11-30 | 广州中海达卫星导航技术股份有限公司 | 一种gnss多模单频rtk周跳探测方法及装置 |
CN106371115A (zh) * | 2016-06-27 | 2017-02-01 | 上海达华测绘有限公司 | 基于北斗短报文技术多模多频海上精密定位方法 |
CN106441174A (zh) * | 2016-09-09 | 2017-02-22 | 桂林电子科技大学 | 一种高边坡形变监测方法及系统 |
CN106526634A (zh) * | 2016-10-19 | 2017-03-22 | 闽江学院 | 一种基于自调整卡尔曼滤波的多普勒辅助载波相位平滑伪距方法 |
CN106772472A (zh) * | 2016-12-06 | 2017-05-31 | 中国地震局第监测中心 | 一种顾及gnss接收机钟跳的导出多普勒观测值构造方法 |
CN108317949A (zh) * | 2018-02-07 | 2018-07-24 | 桂林电子科技大学 | 一种rtk高精度差分定位形变监测系统及方法 |
CN108572377A (zh) * | 2018-04-13 | 2018-09-25 | 桂林电子科技大学 | 基于多普勒辅助的mw组合法探测与修复周跳改进的方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU639308B2 (en) * | 1990-03-09 | 1993-07-22 | Sigtec Navigation Pty Ltd | G p s receiver |
US7069019B2 (en) * | 2001-09-08 | 2006-06-27 | Sirf Technology, Inc. | System and method to estimate the location of a receiver |
US7498979B2 (en) * | 2006-04-17 | 2009-03-03 | Trimble Navigation Limited | Fast decimeter-level GNSS positioning |
US20060047413A1 (en) * | 2003-12-02 | 2006-03-02 | Lopez Nestor Z | GNSS navigation solution integrity in non-controlled environments |
US8131463B2 (en) * | 2003-12-02 | 2012-03-06 | Gmv Aerospace And Defence, S.A. | GNSS navigation solution integrity in non-controlled environments |
US7982667B2 (en) * | 2006-04-17 | 2011-07-19 | Trimble Navigation Limited | Post-processed accuracy prediction for GNSS positioning |
US9651667B2 (en) * | 2007-06-22 | 2017-05-16 | Trimble Inc. | Combined cycle slip indicators for regionally augmented GNSS |
US9958549B2 (en) * | 2012-09-20 | 2018-05-01 | Cornell University | Methods and apparatus for detecting spoofing of global navigation satellite system signals using carrier phase measurements and known antenna motions |
CN105652298B (zh) * | 2016-01-12 | 2018-05-11 | 广州市城市规划勘测设计研究院 | 一种bds三频伪距相位组合的周跳探测修复方法及装置 |
-
2018
- 2018-10-17 CN CN201811205995.XA patent/CN109444935B/zh active Active
Patent Citations (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5323322A (en) * | 1992-03-05 | 1994-06-21 | Trimble Navigation Limited | Networked differential GPS system |
JPH0777570A (ja) * | 1993-07-16 | 1995-03-20 | Caterpillar Inc | 衛星基地局航法システムにおけるサイクルスリップ検出方法及び装置 |
JPH07190769A (ja) * | 1993-12-27 | 1995-07-28 | Sokkia Co Ltd | Gps干渉測位方法 |
US5914685A (en) * | 1997-04-25 | 1999-06-22 | Magellan Corporation | Relative position measuring techniques using both GPS and GLONASS carrier phase measurements |
US6127968A (en) * | 1998-01-28 | 2000-10-03 | Trimble Navigation Limited | On-the-fly RTK positioning system with single frequency receiver |
WO2001022111A1 (en) * | 1999-09-24 | 2001-03-29 | Thales Geosolutions Group Limited | Global navigation satellite systems and methods |
CN1592855A (zh) * | 2001-10-22 | 2005-03-09 | 高通股份有限公司 | 开环跟踪gps信号的方法 |
EP1729145A1 (en) * | 2005-06-02 | 2006-12-06 | Gmv, S.A. | Method and system for providing GNSS navigation position solution with guaranteed integrity in non-controlled environments |
CN101449178A (zh) * | 2006-05-30 | 2009-06-03 | 诺基亚公司 | 支持相对定位 |
JP2008175730A (ja) * | 2007-01-19 | 2008-07-31 | Churyo Eng Kk | 移動局の速度計測装置 |
JP2010071686A (ja) * | 2008-09-16 | 2010-04-02 | Sumitomo Electric Ind Ltd | 測位装置、コンピュータプログラム及び測位方法 |
CN102035774A (zh) * | 2009-09-25 | 2011-04-27 | 上海伽利略导航有限公司 | 一种单频载波相位定位的方法及装置 |
CN102288978A (zh) * | 2011-07-20 | 2011-12-21 | 东南大学 | 一种cors基站周跳探测与修复方法 |
WO2013016800A1 (en) * | 2011-07-29 | 2013-02-07 | Baseband Technologies Inc. | System, method, and computer program for a low power and low cost gnss receiver |
CN102565821A (zh) * | 2011-12-22 | 2012-07-11 | 浙江大学 | 多普勒频偏辅助的卫星导航信号载波周跳实时探测与修复的方法 |
JP2015068768A (ja) * | 2013-09-30 | 2015-04-13 | 日本電気株式会社 | 測位システムと装置と方法並びにプログラム |
CN103576167A (zh) * | 2013-10-28 | 2014-02-12 | 中国科学院国家授时中心 | 基于hht和支持向量机的周跳探测与修复方法 |
EP2966477A1 (en) * | 2014-07-09 | 2016-01-13 | ANavS GmbH | Method for determining the position and attitude of a moving object using low-cost receivers |
CN104570013A (zh) * | 2014-12-30 | 2015-04-29 | 北京无线电计量测试研究所 | 一种用于频率驯服的实时gps载波相位周跳的探测方法 |
CN106168672A (zh) * | 2016-01-01 | 2016-11-30 | 广州中海达卫星导航技术股份有限公司 | 一种gnss多模单频rtk周跳探测方法及装置 |
CN106371115A (zh) * | 2016-06-27 | 2017-02-01 | 上海达华测绘有限公司 | 基于北斗短报文技术多模多频海上精密定位方法 |
CN106441174A (zh) * | 2016-09-09 | 2017-02-22 | 桂林电子科技大学 | 一种高边坡形变监测方法及系统 |
CN106526634A (zh) * | 2016-10-19 | 2017-03-22 | 闽江学院 | 一种基于自调整卡尔曼滤波的多普勒辅助载波相位平滑伪距方法 |
CN106772472A (zh) * | 2016-12-06 | 2017-05-31 | 中国地震局第监测中心 | 一种顾及gnss接收机钟跳的导出多普勒观测值构造方法 |
CN108317949A (zh) * | 2018-02-07 | 2018-07-24 | 桂林电子科技大学 | 一种rtk高精度差分定位形变监测系统及方法 |
CN108572377A (zh) * | 2018-04-13 | 2018-09-25 | 桂林电子科技大学 | 基于多普勒辅助的mw组合法探测与修复周跳改进的方法 |
Non-Patent Citations (16)
Title |
---|
"Cycle slip detection in single frequency GPS carrier observations using expected Doppler shift";Cederholm P;《Nordic journal of surveying and real estate research》;20140131;第1109-1114页 * |
"Detection and Processing of Carrier Phase Cycle Slips Using High Order Difference and Autocorrelation Function";CHEN Shiping;《Modern Navigation》;20180228;第104-108页 * |
"Improved cycle slip detection and repair method for GPS receiver based on TurboEdit algorithm";Duan Rong;《Chinese Journal of Scientific Instrument》;20151130;第2487-2494页 * |
"Improving TurboEdit Real-timing Cycle Slip Detection by the Construction of Threshold Model";ZHANG Xiaohong;《Geomatics and Information Science of Wuhan University》;20170330;第285-292页 * |
"Instantaneous Doppler-aided RTK positioning with single frequency receivers";M. Bahrami;《IEEE/ION Position, Location and Navigation Symposium》;20101231;第70-78页 * |
"MATLAB software for GPS cycle-slip processing";Dai Z;《GPS solutions》;20120228;第267-272页 * |
"Multipath sensitivity and carrier slip tolerance of an integrated Doppler DGPS navigation algorithm";J. Sennott;《IEEE Symposium on Position Location and Navigation. A Decade of Excellence in the Navigation Sciences》;19901231;第638-644页 * |
"一种改进型多普勒周跳探测方法";张炘;《航天控制》;20160530;第9-13页 * |
"利用多普勒观测检测周跳和粗差";常志巧;《测绘通报》;20080330;第28-30页 * |
"基于GNSS信号极化特性的多路径分析";纪元法;《第一届中国卫星导航学术年会论文集》;20101231;第209-214页 * |
"基于GPS载波相位的边坡形变监测方法";孙希延;《科学技术与工程》;20171231;第84-88页 * |
"基于星历计算的多普勒积分法的周跳探测与修复";董明;《大地测量与地球动力学》;20120330;第121-124页 * |
"多普勒观测量在GNSS数据处理中的应用与效果评估";李仁安;《工程勘察》;20170630;第54-58页 * |
"联合多普勒的载波相位精密测速方法";姚曜;《中国惯性技术学报》;20180330;第310-315页 * |
基于多项式拟合的周跳探测改进方法;张俊;《西部大开发(土地开发工程研究)》;20170920(第09期);第12-16页 * |
用多普勒积分检测周跳;申春明等;《城市勘测》;20060225(第01期);第24-26页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109444935A (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109444935B (zh) | 一种低采样率的多普勒周跳探测和修复方法 | |
CN109738917B (zh) | 一种北斗变形监测中的多路径误差削弱方法及装置 | |
CN108572377B (zh) | 基于多普勒辅助的mw组合法探测与修复周跳改进的方法 | |
CN103529462B (zh) | 一种用于全球导航卫星系统的动态周跳探测与修复方法 | |
EP3109670B1 (en) | Global navigation satellite system (gnss) spoofing detection with carrier phase and inertial sensors | |
CN108169774B (zh) | 支持rtppp和rtk的多模gnss单频周跳探测与修复方法 | |
JP6119838B2 (ja) | 状態検出方法、補正値処理装置、測位システム、および状態検出プログラム | |
CN107678050B (zh) | 基于粒子滤波的glonass相位频间偏差实时追踪和精密估计方法 | |
EP2348335A1 (en) | A receiver and method for authenticating satellite signals | |
JP5094344B2 (ja) | 異常衛星検知装置および測位装置 | |
CN108873029B (zh) | 一种实现导航接收机钟差建模的方法 | |
CN109799520B (zh) | 一种基于gnss三频的纯载波相位rtk定位方法 | |
US9810790B2 (en) | Method for evaluating a satellite signal in a global navigation satellite system with respect to a multipath error, receiver for a global navigation satellite system and motor vehicle | |
CN108196267B (zh) | 一种基于gnss cp技术的不间断时间传递方法 | |
JP2012073246A (ja) | 擬似的距離測定値の選定による位置推定の修正方法 | |
CN114935770B (zh) | 一种多历元加快精密单点定位收敛速度的方法及装置 | |
CN108181632A (zh) | 基于模糊度全微分的gnss单频数据周跳探测与修复方法 | |
CN115480279A (zh) | Gnss导航方法和终端、组合导航系统、存储介质 | |
US10024974B1 (en) | Method and apparatus applicable to positioning in NLOS environment | |
CN116009042A (zh) | 一种单站载波历元间差分实时探测相对形变的方法及系统 | |
CN115902963A (zh) | 一种单点定位数据处理方法、装置、电子设备及存储介质 | |
CN115267848A (zh) | 一种双频周跳探测及修复方法、系统、设备及介质 | |
CN105652298B (zh) | 一种bds三频伪距相位组合的周跳探测修复方法及装置 | |
Deo et al. | Cycle slip and clock jump repair with multi-frequency multi-constellation GNSS data for precise point positioning | |
CN106873003A (zh) | 一种基于多历元奇偶矢量积累的双模接收机raim方法 |
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 |