CN110727008A - 一种gnss整周模糊度的快速解算方法 - Google Patents

一种gnss整周模糊度的快速解算方法 Download PDF

Info

Publication number
CN110727008A
CN110727008A CN201911131222.6A CN201911131222A CN110727008A CN 110727008 A CN110727008 A CN 110727008A CN 201911131222 A CN201911131222 A CN 201911131222A CN 110727008 A CN110727008 A CN 110727008A
Authority
CN
China
Prior art keywords
integer ambiguity
integer
correction
value
gnss
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
Application number
CN201911131222.6A
Other languages
English (en)
Other versions
CN110727008B (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201911131222.6A priority Critical patent/CN110727008B/zh
Publication of CN110727008A publication Critical patent/CN110727008A/zh
Application granted granted Critical
Publication of CN110727008B publication Critical patent/CN110727008B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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

本发明公开了一种GNSS整周模糊度的快速解算方法，包括以下步骤：采集GNSS数据，并通过GNSS数据构建组合观测量
Figure DDA0002278348950000011
根据组合观测量，获取双差整周模糊度近似整数值的改正数候选值；对改正数候选值进行判断，确定改正数，并利用最终的改正数获取双差整周模糊度的终值，得到整周模糊度解算结果。本发明逐颗卫星进行整周模糊度解算，不需要模糊度搜索过程，流程简单，具有计算效率高的特点，且适用于动态模式下的整周模糊度快速解算。

Description

一种GNSS整周模糊度的快速解算方法
技术领域
本发明属于GNSS定位数据处理领域，具体涉及一种GNSS整周模糊度的快速解算方法。
背景技术
GNSS高精度定位依赖于毫米级精度的载波相位观测值，更是以整周模糊度的正确解算为基本前提，因此，整周模糊度解算方法对GNSS高精度动态定位至关重要。虽然目前已经存在多种GNSS载波相位整周模糊度解算方法，但现有的技术存在以下不足：
1、逐颗卫星进行模糊度解算的方法常常需要引入伪距观测量，由于伪距精度较低，进而导致此类方法直接解算相位整周模糊度的正确率较低，因此，此类一般只用于解算波长较长的相位宽巷组合量的模糊度；
2、以LAMBDA为代表的基于整数最小二乘的模糊度解算方法，虽然能够较好的解算出整周模糊度，但这种方法存在模糊度的搜索过程，一般需要较大的计算量，尤其是当观测卫星个数较多时方法的效率将降低，所以不能很好的适用于多星座的高动态定位场景；同时此类方法计算过程复杂，不利于定位终端的简化。
发明内容
针对现有技术中的上述不足，本发明提供的一种GNSS整周模糊度的快速解算方法解决了现有技术存在的问题。
为了达到上述发明目的，本发明采用的技术方案为：一种GNSS整周模糊度的快速解算方法，包括以下步骤：
S1、采集GNSS数据，并通过GNSS数据构建组合观测量
Figure BDA0002278348930000011
S2、根据组合观测量，获取双差整周模糊度近似整数值的改正数候选值；
S3、对改正数候选值进行判断，确定最终的改正数，并利用最终的改正数获取双差整周模糊度的终值，得到整周模糊度解算结果。
进一步地，所述步骤S1中组合观测量
Figure BDA0002278348930000021
为：
Figure BDA0002278348930000022
其中，
Figure BDA0002278348930000023
和
Figure BDA0002278348930000024
表示频率不同的以周为单位的载波相位观测值，λ1表示与
Figure BDA0002278348930000025
对应载波的波长，λ2表示与
Figure BDA0002278348930000026
对应载波的波长，表示卫星与接收机之间的双差计算，NF1表示与载波相位
Figure BDA0002278348930000028
对应的整周模糊度的近似整数值，NF2表示与载波相位
Figure BDA0002278348930000029
对应的整周模糊度的近似整数值。
进一步地，所述步骤S2中改正数候选值为
Figure BDA00022783489300000210
和
Figure BDA00022783489300000211
所述
Figure BDA00022783489300000212
和具体为：
其中，
Figure BDA00022783489300000215
表示组合观测量，βi表示与载波相位
Figure BDA00022783489300000216
对应的双差整周模糊度近似整数值的改正数计算系数，βi＝fi/(f1-f2)，f1表示载波相位
Figure BDA00022783489300000217
的频率，f2表示载波相位
Figure BDA00022783489300000218
的频率，i＝{1,2}，当i＝1时，i'＝2，当i＝2时，i'＝1；λi'表示与
Figure BDA00022783489300000219
对应载波的波长，Sgn(·)表示去符号运算，(·)d表示取小数运算。
进一步地，所述步骤S3中对改正数候选值进行判断，确定最终的改正数
Figure BDA00022783489300000224
的具体方法为：
其中，
Figure BDA00022783489300000221
和
Figure BDA00022783489300000222
表示改正数候选值，R(·)表示四舍五入运算，DNI(X)＝|R(X)-X|。
进一步地，所述步骤S3中双差整周模糊度的终值
Figure BDA00022783489300000223
为：
Figure BDA0002278348930000031
其中，
Figure BDA0002278348930000032
表示与载波相位
Figure BDA0002278348930000033
对应的双差整周模糊度的近似整数值，
Figure BDA0002278348930000034
表示最终的改正数。
本发明的有益效果为：
(1)本发明逐颗卫星进行整周模糊度解算，不需要模糊度搜索过程，流程简单，具有计算效率高的特点，且适用于多星座数据的快速整周模糊度解算。
(2)本发明无需接收机位置信息，适用于动态和静态模式下的整周模糊度解算，尤其适用于实时定位解算中的整周模糊度快速数据处理。
附图说明
图1为本发明提出的一种GNSS整周模糊度的快速解算方法流程图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
下面结合附图详细说明本发明的实施例。
如图1所示，一种GNSS整周模糊度的快速解算方法，包括以下步骤：
S1、采集GNSS数据，并通过GNSS数据构建组合观测量
Figure BDA0002278348930000035
S2、根据组合观测量，获取双差整周模糊度近似整数值的改正数候选值；
S3、对改正数候选值进行判断，确定最终的改正数，并利用最终的改正数获取双差整周模糊度的终值，得到整周模糊度解算结果。
在本实施例中，所述步骤1中GNSS数据包括载波相位观测值、整周模糊度近似整数值和载波的波长，忽略观测噪声，短基线双差后的载波相位观测值的观测方程为：
Figure BDA0002278348930000041
Figure BDA0002278348930000042
其中，
Figure BDA0002278348930000043
和
Figure BDA0002278348930000044
表示频率不同的以周为单位的载波相位观测值，λ1表示载波相位
Figure BDA0002278348930000045
的波长，λ2表示载波相位
Figure BDA0002278348930000046
的波长，
Figure BDA0002278348930000047
表示卫星与接收机之间的双差计算，ρ表示卫星与接收机之间的距离，N1表示载波相位
Figure BDA0002278348930000048
的整周模糊度，N2表示载波相位
Figure BDA0002278348930000049
的整周模糊度。
在本实施例中，所述步骤S1中构建组合观测量
Figure BDA00022783489300000410
的具体方法为：
A1、根据公式1和公式2整理得到：
Figure BDA00022783489300000411
A2、令
Figure BDA00022783489300000412
由公式3可以得到：
Figure BDA00022783489300000413
A3、根据公式4构建组合观测量，所述组合观测量
Figure BDA00022783489300000414
为：
其中，NFi表示与
Figure BDA00022783489300000416
对应的整周模糊度浮点解四舍五入取整值，
Figure BDA00022783489300000417
表示与
Figure BDA00022783489300000418
对应的整周模糊度的改正数，i＝{1,2}。
在本实施例中，所述步骤S2包括以下步骤：
S2.1、根据载波波长与频率反比的关系和组合观测量，可得到
Figure BDA00022783489300000419
和的关系为：
Figure BDA00022783489300000422
S2.2、将公式6和公式7进行取小数运算，得到公式8和公式9：
Figure BDA0002278348930000051
Figure BDA0002278348930000052
S2.3、令βi＝fi/(f1-f2)，得到公式8：
Figure BDA0002278348930000053
Figure BDA0002278348930000054
S2.4、当
Figure BDA0002278348930000055
和
Figure BDA0002278348930000056
时，可以得到整周模糊度的改正数表达式为：
Figure BDA0002278348930000057
Figure BDA0002278348930000058
其中，f1表示载波相位
Figure BDA0002278348930000059
的频率，f2表示载波相位
Figure BDA00022783489300000510
的频率，(·)d表示取小数运算，e为整数，且
Figure BDA00022783489300000511
Sgn(·)表示取符号运算，βi表示与
Figure BDA00022783489300000512
对应的整周模糊度的改正数计算系数，i＝{1,2}；
S2.5、根据整数值e和整周模糊度的改正数表达式，获取改正数候选值为
Figure BDA00022783489300000513
和
Figure BDA00022783489300000514
所述
Figure BDA00022783489300000515
和
Figure BDA00022783489300000516
具体为：
Figure BDA00022783489300000517
其中，i＝{1,2}，当i＝1时，i'＝2，当i＝2时，i'＝1；λi'表示与
Figure BDA00022783489300000518
对应载波的波长。
所述步骤S3中对改正数候选值进行判断，确定最终改正数的具体方法为：
其中，
Figure BDA00022783489300000521
表示与
Figure BDA00022783489300000522
对应的整周模糊度最终的改正数，R(·)表示四舍五入运算，DNI(X)＝|R(X)-X|。
所述步骤S3中双差整周模糊度的终值
Figure BDA0002278348930000061
为：
本发明逐颗卫星进行整周模糊度解算，不需要模糊度搜索过程，流程简单，具有计算效率高的特点，且适用于多星座数据的快速整周模糊度解算。本发明无需接收机位置信息，适用于动态和静态模式下的整周模糊度解算，尤其适用于实时定位解算中的整周模糊度快速数据处理。

Claims (5)

1.一种GNSS整周模糊度的快速解算方法，其特征在于，包括以下步骤：
S1、采集GNSS数据，并通过GNSS数据构建组合观测量
Figure FDA0002278348920000013
S2、根据组合观测量，获取双差整周模糊度近似整数值的改正数候选值；
S3、对改正数候选值进行判断，确定最终的改正数，并利用最终的改正数获取双差整周模糊度的终值，得到整周模糊度解算结果。
2.根据权利要求1所述的GNSS整周模糊度的快速解算方法，其特征在于，所述步骤S1中组合观测量
Figure FDA0002278348920000014
为：
Figure FDA0002278348920000011
其中，
Figure FDA0002278348920000015
和
Figure FDA0002278348920000016
表示频率不同的以周为单位的载波相位观测值，λ1表示与
Figure FDA0002278348920000017
对应载波的波长，λ2表示与
Figure FDA0002278348920000019
对应载波的波长，
Figure FDA0002278348920000018
表示卫星与接收机之间的双差计算，NF1表示与载波相位
Figure FDA00022783489200000110
对应的整周模糊度的近似整数值，NF2表示与载波相位
Figure FDA00022783489200000111
对应的整周模糊度的近似整数值。
3.根据权利要求1所述的GNSS整周模糊度的快速解算方法，其特征在于，所述步骤S2中改正数候选值为
Figure FDA00022783489200000112
和
Figure FDA00022783489200000113
所述和
Figure FDA00022783489200000115
具体为：
其中，
Figure FDA00022783489200000116
表示组合观测量，βi表示与载波相位对应的双差整周模糊度近似整数值的改正数计算系数，βi＝fi/(f1-f2)，f1表示载波相位的频率，f2表示载波相位
Figure FDA00022783489200000120
的频率，i＝{1,2}，当i＝1时，i'＝2，当i＝2时，i'＝1；λi'表示与
Figure FDA00022783489200000119
对应载波的波长，Sgn(·)表示去符号运算，(·)d表示取小数运算。
4.根据权利要求1所述的GNSS整周模糊度的快速解算方法，其特征在于，所述步骤S3中对改正数候选值进行判断，确定最终的改正数
Figure FDA00022783489200000121
的具体方法为：
Figure FDA0002278348920000021
其中，
Figure FDA0002278348920000022
和
Figure FDA0002278348920000023
表示改正数候选值，R(·)表示四舍五入运算，DNI(X)＝|R(X)-X|。
5.根据权利要求1所述的GNSS整周模糊度的快速解算方法，其特征在于，所述步骤S3中双差整周模糊度的终值为：
其中，
Figure FDA0002278348920000027
表示与载波相位
Figure FDA0002278348920000026
对应的双差整周模糊度的近似整数值，
Figure FDA0002278348920000028
表示最终的改正数。
CN201911131222.6A 2019-11-19 2019-11-19 一种gnss整周模糊度的快速解算方法 Active CN110727008B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911131222.6A CN110727008B (zh) 2019-11-19 2019-11-19 一种gnss整周模糊度的快速解算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911131222.6A CN110727008B (zh) 2019-11-19 2019-11-19 一种gnss整周模糊度的快速解算方法

Publications (2)

Publication Number Publication Date
CN110727008A true CN110727008A (zh) 2020-01-24
CN110727008B CN110727008B (zh) 2021-06-22

Family

ID=69225414

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911131222.6A Active CN110727008B (zh) 2019-11-19 2019-11-19 一种gnss整周模糊度的快速解算方法

Country Status (1)

Country Link
CN (1) CN110727008B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111522042A (zh) * 2020-05-08 2020-08-11 广东电网有限责任公司电力科学研究院 一种基于双卫星定位模块的定向装置及定向方法
CN111751855A (zh) * 2020-06-28 2020-10-09 北京建筑大学 Gnss单历元双差整周模糊度快速确定方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105204049A (zh) * 2015-10-12 2015-12-30 中国科学院光电研究院 一种基于三维转一维的载波相位整周模糊度搜索的定位方法
CN105223598A (zh) * 2015-09-28 2016-01-06 厦门理工学院 一种gnss载波相位整周模糊度单历元解算方法
CN105353393A (zh) * 2015-10-22 2016-02-24 北京理工大学 一种双频载波相位整周模糊度固定的方法
US20170269231A1 (en) * 2016-03-18 2017-09-21 Deere & Company Rapid recovery of precise position after temporary signal loss
US20190196022A1 (en) * 2017-12-21 2019-06-27 Valeo North America, Inc. Broadcast and utilization of precise gnss correction data
CN110456397A (zh) * 2019-07-29 2019-11-15 桂林电子科技大学 一种多天线超短基线定位监测方法、装置及存储介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105223598A (zh) * 2015-09-28 2016-01-06 厦门理工学院 一种gnss载波相位整周模糊度单历元解算方法
CN105204049A (zh) * 2015-10-12 2015-12-30 中国科学院光电研究院 一种基于三维转一维的载波相位整周模糊度搜索的定位方法
CN105353393A (zh) * 2015-10-22 2016-02-24 北京理工大学 一种双频载波相位整周模糊度固定的方法
US20170269231A1 (en) * 2016-03-18 2017-09-21 Deere & Company Rapid recovery of precise position after temporary signal loss
US20190196022A1 (en) * 2017-12-21 2019-06-27 Valeo North America, Inc. Broadcast and utilization of precise gnss correction data
CN110456397A (zh) * 2019-07-29 2019-11-15 桂林电子科技大学 一种多天线超短基线定位监测方法、装置及存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冯威: "GPS中央差分定位系统理论与应用研究", 《中国博士学位论文全文数据库基础科学辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111522042A (zh) * 2020-05-08 2020-08-11 广东电网有限责任公司电力科学研究院 一种基于双卫星定位模块的定向装置及定向方法
CN111751855A (zh) * 2020-06-28 2020-10-09 北京建筑大学 Gnss单历元双差整周模糊度快速确定方法
CN111751855B (zh) * 2020-06-28 2023-03-14 北京建筑大学 Gnss单历元双差整周模糊度快速确定方法

Also Published As

Publication number Publication date
CN110727008B (zh) 2021-06-22

Similar Documents

Publication Publication Date Title
CN114518586B (zh) 一种基于球谐展开的gnss精密单点定位方法
CN105676250B (zh) 一种基于gnss的单历元三频模糊度解算方法
CN104749594B (zh) 一种gps双频非差周跳探测与修复方法及装置
CN110058282B (zh) 一种基于双频gnss智能手机的ppp高精度定位方法
CN109765589B (zh) 一种基于无电离层组合的三频gnss实时周跳固定技术
CN110727008B (zh) 一种gnss整周模糊度的快速解算方法
CN111208541B (zh) 一种基于gpu的大规模gnss数据处理方法
CN109085628A (zh) 一种整周模糊度的固定方法及系统
CN110346823B (zh) 可用于北斗精密单点定位的三频模糊度解算方法
CN103197335A (zh) 采用改进正则化方法抑制dgps整周模糊度病态性的方法
CN113325446B (zh) 一种多模共频gnss载波相位时间传递方法及系统
CN114779301B (zh) 一种基于广播星历的卫星导航实时精密单点定位方法
CN110727000B (zh) 一种基于gnss高采样率数据的小周跳修复方法
CN112987059A (zh) 一种基于三频模糊度解算的整周模糊度计算方法
CN109884679A (zh) 一种单模gnss系统跨频点混合双差rtk解算方法
CN111562594A (zh) 基于站间单差和模糊度固定的北斗精密时间传递方法
CN112230254B (zh) 一种gps载波相位多径误差的校正方法及装置
CN105204049A (zh) 一种基于三维转一维的载波相位整周模糊度搜索的定位方法
CN115494535A (zh) 一种北斗三频非差相位模糊度的解算方法及装置
CN115308781B (zh) 基于bdgim辅助的相位平滑伪距高精度时间传递方法
CN116577815A (zh) 一种多频多gnss精密单点定位方法、装置及设备
CN115561793A (zh) 一种基于并行计算的实时北斗相位小数偏差快速估计方法
CN112987048B (zh) 一种适用于Andriod智能终端的高精度定位方法、智能终端以及存储介质
CN111766614B (zh) 一种守时实验室联合组网的亚纳秒实时时间比对方法
CN115826006B (zh) 一种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