CN111812687A - 高频动态多星座gnss数据周跳低时延解算方法 - Google Patents
高频动态多星座gnss数据周跳低时延解算方法 Download PDFInfo
- Publication number
- CN111812687A CN111812687A CN202010778750.7A CN202010778750A CN111812687A CN 111812687 A CN111812687 A CN 111812687A CN 202010778750 A CN202010778750 A CN 202010778750A CN 111812687 A CN111812687 A CN 111812687A
- Authority
- CN
- China
- Prior art keywords
- cycle slip
- slip
- cycle
- frequency
- restoration
- 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
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/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
-
- 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)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种高频动态多星座GNSS数据周跳低时延解算方法,该方法包括逐颗卫星利用相位和伪距组合观测量对高频观测数据进行周跳初步修复,利用双频相位组合观测量对得到的周跳修复量初值进行小周跳修复,根据分别得到的周跳修复量计算最终的周跳。本发明基于高频观测数据的噪声特性,仅需利用单颗卫星的观测数据进行周跳探测与修复,且降低了周跳修复对伪距观测值的精度要求,具有复杂度低、计算效率高等优点,有利于算力资源受限的移动端超高频实时动态定位数据的在线处理。
Description
技术领域
本发明涉及周跳探测及修复技术领域,具体涉及一种高频动态多星座GNSS数据周跳低时延解算方法。
背景技术
高精度载波相位观测在精密定位和导航应用中起着重要的作用。但是地表观测GNSS信号相对较弱,易受干扰、遮挡等影响,载波相位观测值中会发生周跳。若周跳未能被准确的探测或修复,将导致整周模糊度参数失效,进而需要重新进行整周模糊度固定,影响GNSS导航定位的性能。在高动态RTK应用中,GNSS数据采集中周跳现象频繁发生,这直接影响动态定位在许多高速应用场景中的实用性,为此需要一种快速高效的低时延周跳修复方法。
目前,探测与修复周跳的方法主要有伪距相位组合法、电离层残差法、多项式拟合法、多普勒积分法、小波分析法、高次差法、卡尔曼滤波法等。试验证明这些方法都具有较好的探测及修复效果,但也存在不足之处。例如,伪距相位组合法,探测精度高、程序容易实现,但其探测能力主要取决于电离层延迟和多路径效应在历元间的变化以及伪距和载波相位观测噪声、载波波长大小等;电离层残差法是探测周跳很好的方法,但是传统的电离层残差法没有很好的解决周跳发生在哪个频率的问题,且对于组合周跳的问题一般需要一些辅助的方法联合求解;质量控制理论也用于周跳探测,通过整数最小二乘估计法可以获得更高的准确性,但由于需要进行搜索,周跳计算的算法复杂度高,周跳修复过程变得更加耗时。
自动驾驶、无人机等高速RTK应用场景需要实时处理超高频的GNSS数据,在这种情况下,低时延的周跳探测及修复至关重要。
发明内容
针对现有技术中的上述不足,本发明提供了一种高频动态多星座GNSS数据周跳低时延解算方法。
为了达到上述发明目的,本发明采用的技术方案为:
一种高频动态多星座GNSS数据周跳低时延解算方法,包括以下步骤:
S1、逐颗卫星利用相位和伪距组合观测量对高频观测数据进行周跳初步修复;
S2、利用双频相位组合观测量对步骤S1得到的周跳修复量初值进行小周跳修复;
S3、根据步骤S1和步骤S2分别得到的周跳修复量计算最终的周跳。
优选地,所述步骤S1具体为:
利用相位和伪距组合观测量计算每颗卫星单个频率上的周跳初始修复量。
优选地,所述周跳初始修复量的计算公式为:
优选地,所述步骤S1中设定周跳初始修复量满足条件:
ε(δN1,p)<4.5
ε(δN2,p)<3.5
其中,ε(δNi,p)表示周跳初始修复量与真实周跳值之间的偏差,i表示载波序号。
优选地,所述步骤S2具体为:
根据步骤S1得到的周跳修复量初值,利用双频相位组合观测量计算单个频率上的小周跳,得到未完全修复的周跳值。
优选地,所述小周跳的计算公式为:
δNi`=(e-(δO-MGF)d)βi
优选地,所述步骤S3中最终的周跳的计算公式为:
δNi=R(δNi,p)+(e-(δO-MGF)d)βi
其中,δNi表示最终修复的周跳。
本发明具有以下有益效果:
本发明基于高频观测数据的噪声特性,利用伪距与相位组合得到周跳初始修复量,再结合双频相位组合对初始修复量进行精确修复,仅需利用单颗卫星的观测数据进行周跳探测与修复,且降低了周跳修复对伪距观测值的精度要求,具有复杂度低、计算效率高等优点,有利于算力资源受限的移动端超高频实时动态定位数据的在线处理。
附图说明
图1为本发明的高频动态多星座GNSS数据周跳低时延解算方法流程示意图;
图2为本发明实施例中动态实验的运动轨迹示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,本发明基于高频数据误差特征和伪距与相位组合,提出了一种高频动态多星座GNSS数据周跳低时延解算方法O-MGF(Observation based MGF),包括以下步骤S1至S3:
S1、逐颗卫星利用相位和伪距组合观测量对高频观测数据进行周跳初步修复;
在本实施例中,本发明利用相位和伪距组合观测量计算每颗卫星单个频率上的周跳初始修复量。
计算单个频率上的周跳初始修复量的计算公式为:
本发明设定周跳初始修复量需要满足以下条件:
ε(δN1,p)<4.5
ε(δN2,p)<3.5
为了实现探测所有周跳,本发明设定伪距周跳探测结果与周跳真值差值需要满足以下条件:
即得到各卫星系统的伪距观测噪声需要满足以下条件:
S2、利用双频相位组合观测量对步骤S1得到的周跳修复量初值进行小周跳修复;
在本实施例中,本发明根据步骤S1得到的周跳修复量,利用双频相位组合观测量高效解算出双频相位数据中单个频率上的小周跳,得到未完全修复的周跳值。
本发明计算小周跳的计算公式为:
δNi`=(e-(δO-MGF)d)βi
S3、根据步骤S1和步骤S2分别得到的周跳修复量计算最终的周跳。
在本实施例中,本发明根据步骤S1得到的周跳初始修复量和步骤S2得到的未完全修复的周跳值,计算得到最终的周跳。
最终的周跳的计算公式为:
δNi=R(δNi,p)+(e-(δO-MGF)d)βi
其中,δNi表示最终修复的周跳。
由于e存在两个备选值{0,Sgn(δO-MGF)},因此最终的周跳δNi也存在两个备选值,分别对应为δNi(1),δNi(2)。基于周跳的整数性质,本发明选择较小的DNI(δNi)对应的δNi作为最终修复值,表示为:
本发明基于高频观测数据的噪声特性,利用伪距与相位组合得到周跳修复量,仅需利用单颗卫星的观测数据进行周跳探测与修复,且降低了周跳修复对伪距观测值的精度要求。
下面本发明采用多星座GNSS数据实验对本发明的探测及修复方法的性能和效果进行具体说明。
本发明采用2019年9月11日实测的动态数据和2020年1月13日实测的静态数据进行O-MGF周跳实验分析,数据均为板卡实测的双频多星座20HZ的高频数据,卫星系统包括GPS、BDS和GLONASS,动态数据和静态数据用不同品牌的接收机采集,时长分别为30分钟和40分钟。动态数据的运动轨迹如图2所示。
本发明通过对动态实测数据相应历元处加模拟周跳来验证本文周跳修复算法的准确性及高效性。模拟周跳如下:模拟周跳历元间隔为0.25s,模拟历元数为5832个。各卫星f1频率依次模拟周跳大小为:[-16,-14,-12,-10,-4,-3,-2,-1,1,2,3,4,6,8,10,12]周。f2频率依次模拟周跳大小为:[-12,-9,-6,-3,-2,-1,1,2,3,6,9,12]周。分别用本发明方法(O-MGF)和HMW-GF(Hatch-Melbourne-Wübbena和Geometry-free)组合方法进行周跳探测和修复,并对比分析修复结果。测试过程使用相同的软硬件环境,且仅统计两种方法周跳探测与修复模块的耗时,统计结果如表1。
表1模拟周跳修复结果统计
从表1可以看出,本发明方法在各个系统探测周跳成功率均高于HMW-GF组合。尤其在GLONASS系统中,本发明方法更是达到了百分百的修复率,GPS、BDS系统修复成功率次之。
经过统计,本发明方法探测耗时1684ms,HMW-GF组合耗时101310ms,两者周跳处理平均每个历元的时间延迟分别为0.06ms和3.47ms。解算效率方面O-MGF方法比HMW-GF组合方法快60倍。
本发明利用实测的20Hz三星座GNSS(GPS/BDS/GLONASS)动态数据进行试验分析,结果表明:O-MGF方法的成功率约为99.998%;与HMW-GF方法相比,O-MGF方法解算速度提高了60倍。O-MGF方法的算法复杂度低,计算效率高,有利于算力资源受限的移动端超高频实时动态定位数据的在线快速处理。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (7)
1.一种高频动态多星座GNSS数据周跳低时延解算方法,其特征在于,包括以下步骤:
S1、逐颗卫星利用相位和伪距组合观测量对高频观测数据进行周跳初步修复;
S2、利用双频相位组合观测量对步骤S1得到的周跳修复量初值进行小周跳修复;
S3、根据步骤S1和步骤S2分别得到的周跳修复量计算最终的周跳。
2.根据权利要求1所述的高频动态多星座GNSS数据周跳低时延解算方法,其特征在于,所述步骤S1具体为:
利用相位和伪距组合观测量计算每颗卫星单个频率上的周跳初始修复量。
4.根据权利要求3所述的高频动态多星座GNSS数据周跳低时延解算方法,其特征在于,所述步骤S1中设定周跳初始修复量满足条件:
ε(δN1,p)<4.5
ε(δN2,p)<3.5
其中,ε(δNi,p)表示周跳初始修复量与真实周跳值之间的偏差,i表示载波序号。
5.根据权利要求4所述的高频动态多星座GNSS数据周跳低时延解算方法,其特征在于,所述步骤S2具体为:
根据步骤S1得到的周跳修复量初值,利用双频相位组合观测量计算单个频率上的小周跳,得到未完全修复的周跳值。
7.根据权利要求6所述的高频动态多星座GNSS数据周跳低时延解算方法,其特征在于,所述步骤S3中最终的周跳的计算公式为:
δNi=R(δNi,p)+(e-(δO-MGF)d)βi
其中,δNi表示最终修复的周跳。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010778750.7A CN111812687B (zh) | 2020-08-05 | 2020-08-05 | 高频动态多星座gnss数据周跳低时延解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010778750.7A CN111812687B (zh) | 2020-08-05 | 2020-08-05 | 高频动态多星座gnss数据周跳低时延解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111812687A true CN111812687A (zh) | 2020-10-23 |
CN111812687B CN111812687B (zh) | 2023-08-01 |
Family
ID=72863659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010778750.7A Active CN111812687B (zh) | 2020-08-05 | 2020-08-05 | 高频动态多星座gnss数据周跳低时延解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111812687B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103529462A (zh) * | 2013-10-21 | 2014-01-22 | 西南交通大学 | 一种用于全球导航卫星系统的动态周跳探测与修复方法 |
CN104483688A (zh) * | 2014-11-03 | 2015-04-01 | 中国人民解放军63961部队 | 基于北斗卫星导航系统的高精度基线解算方法 |
CN107132558A (zh) * | 2017-06-13 | 2017-09-05 | 武汉大学 | 惯性辅助的多频多模gnss周跳修复方法及系统 |
CN109358350A (zh) * | 2018-10-08 | 2019-02-19 | 中国人民解放军战略支援部队信息工程大学 | 一种北斗三频周跳探测方法与装置 |
CN110727000A (zh) * | 2019-11-19 | 2020-01-24 | 西南交通大学 | 一种基于gnss高采样率数据的小周跳修复方法 |
-
2020
- 2020-08-05 CN CN202010778750.7A patent/CN111812687B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103529462A (zh) * | 2013-10-21 | 2014-01-22 | 西南交通大学 | 一种用于全球导航卫星系统的动态周跳探测与修复方法 |
CN104483688A (zh) * | 2014-11-03 | 2015-04-01 | 中国人民解放军63961部队 | 基于北斗卫星导航系统的高精度基线解算方法 |
CN107132558A (zh) * | 2017-06-13 | 2017-09-05 | 武汉大学 | 惯性辅助的多频多模gnss周跳修复方法及系统 |
CN109358350A (zh) * | 2018-10-08 | 2019-02-19 | 中国人民解放军战略支援部队信息工程大学 | 一种北斗三频周跳探测方法与装置 |
CN110727000A (zh) * | 2019-11-19 | 2020-01-24 | 西南交通大学 | 一种基于gnss高采样率数据的小周跳修复方法 |
Non-Patent Citations (7)
Title |
---|
冯威等: "一种改进的载波相位周跳探测与修复方法", 《测绘科学》 * |
冯威等: "一种改进的载波相位周跳探测与修复方法", 《测绘科学》, no. 06, 20 November 2010 (2010-11-20) * |
刘先冬: "动态GPS精密定位数据处理方法研究", 《中国知网 硕士电子期刊》 * |
刘先冬: "动态GPS精密定位数据处理方法研究", 《中国知网 硕士电子期刊》, 15 July 2012 (2012-07-15) * |
徐锐等: "一种改进的双频单P码周跳探测与修复方法", 《大地测量与地球动力学》 * |
徐锐等: "一种改进的双频单P码周跳探测与修复方法", 《大地测量与地球动力学》, no. 04, 15 August 2007 (2007-08-15) * |
谢兰天;赵乐文;叶世榕;: "北斗三频周跳探测与修复算法研究", 导航定位与授时, no. 05 * |
Also Published As
Publication number | Publication date |
---|---|
CN111812687B (zh) | 2023-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109459778B (zh) | 基于抗差方差分量估计的码伪距/多普勒联合测速方法及其应用 | |
CN108363084A (zh) | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 | |
CN103529462B (zh) | 一种用于全球导航卫星系统的动态周跳探测与修复方法 | |
CN110208836B (zh) | 基于卡尔曼滤波的gnss高适应性周跳探测与修复方法 | |
CN108254773B (zh) | 一种多gnss的实时钟差解算方法 | |
CN103197335B (zh) | 采用改进正则化方法抑制dgps整周模糊度病态性的方法 | |
DE102013206544A1 (de) | Fortschrittliche Positionierung mit globalen Navigationssatellitensystemen (GNSS) unter Verwendung von präziser Satelliteninformation | |
CN107678050A (zh) | 基于粒子滤波的glonass相位频间偏差实时追踪和精密估计方法 | |
CN102998681A (zh) | 一种卫星导航系统的高频钟差估计方法 | |
CN112462399A (zh) | 一种基于接收机的伪距偏差计算及改正的方法及装置 | |
CN109597105A (zh) | 一种顾及载波系统间偏差的gps/glonass紧组合定位方法 | |
CN114935770B (zh) | 一种多历元加快精密单点定位收敛速度的方法及装置 | |
CN110727000B (zh) | 一种基于gnss高采样率数据的小周跳修复方法 | |
CN115906496A (zh) | 一种卫星钟差估计方法 | |
CN113655504B (zh) | 一种卫星数据处理方法、装置、电子设备及存储介质 | |
CN116819585B (zh) | 基于非线性优化的gnss单点定位方法及导航方法 | |
CN111812687A (zh) | 高频动态多星座gnss数据周跳低时延解算方法 | |
CN115856973B (zh) | Gnss解算方法、装置、定位系统、电子设备及存储介质 | |
CN115902967B (zh) | 基于低轨导航增强卫星信号的导航定位方法、系统及飞行平台 | |
CN114152961A (zh) | 一种导航系统的周跳处理方法和装置 | |
CN105699991A (zh) | 一种基于跳变点阈值判决的自适应载波精度算法 | |
CN111999750B (zh) | 针对杆臂不准的实时单站周跳探测改进方法 | |
CN115826006B (zh) | 一种bds双频周跳探测组合方法 | |
CN115712130B (zh) | 一种伪距校正方法、装置、电子设备及存储介质 | |
Wang et al. | Geometry-free stochastic analysis of BDS triple frequency signals |
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 |