CN106199659B - 基于模糊数学的gnss单站双频观测数据周跳探测及处理方法 - Google Patents

基于模糊数学的gnss单站双频观测数据周跳探测及处理方法 Download PDF

Info

Publication number
CN106199659B
CN106199659B CN201610518199.6A CN201610518199A CN106199659B CN 106199659 B CN106199659 B CN 106199659B CN 201610518199 A CN201610518199 A CN 201610518199A CN 106199659 B CN106199659 B CN 106199659B
Authority
CN
China
Prior art keywords
epochs
value
observation
ionospheric delay
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.)
Active
Application number
CN201610518199.6A
Other languages
English (en)
Other versions
CN106199659A (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201610518199.6A priority Critical patent/CN106199659B/zh
Publication of CN106199659A publication Critical patent/CN106199659A/zh
Application granted granted Critical
Publication of CN106199659B publication Critical patent/CN106199659B/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/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

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (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单站双频观测数据周跳探测及处理方法，通过确定历元间电离层延迟变化量和历元间MW组合观测值变化量的隶属度判断观测值是否有周跳发生，并对难以确定是否有周跳发生的观测值进行降权处理；本发明可有效解决以往周跳探测数据处理过程中阈值难以可靠确定的问题，进而改善GNSS精密数据处理的精度和可靠性。

Description

基于模糊数学的GNSS单站双频观测数据周跳探测及处理方法
技术领域
本发明属于全球导航系统精密数据处理技术领域，涉及一种GNSS单站双频观测数据的周跳探测及处理方法，具体涉及一种基于模糊数学的GNSS单站双频观测数据周跳探测及处理方法。
背景技术
目前，GNSS(Global Navigation Satellite System)双频/多频数据的周跳探测方法主要有TurboEdit方法、电离层残差法、MW(Melbourne-Wübbena)组合结合电离层总电子含量变化率等，但这些方法都存在各自的局限性。电离层残差法基于电离层变化比较缓慢的假设进行周跳探测，如果电离层较为活跃或者信号失锁的时间较长，此时周跳探测的可靠性不高([文献1、2])。TurboEdit方法存在以下几点不足：①当观测值存在多路径误差的影响、电离层比较活跃或者卫星高度角较低时，导致伪距观测噪声放大，这时TurboEdit方法对1-2周的小周跳可能不灵敏，导致漏探([文献3、4])；②受伪距观测值精度的影响，会导致周跳的误判([文献4、5])。MW组合结合电离层总电子含量变化率探测周跳是一种单站双频的GNSS数据预处理方法，在电离层活跃期也能够很好的探测周跳，但这种方法对于观测数据的采样率要求较高，仅适用于高采样率观测数据([文献6])。
现有GNSS周跳探测方法都是在数据处理时先设定某一阈值，该阈值可能是不变的也可能是变化的，若由观测数据得到的检测量大于给定阈值可判定该观测值有周跳发生。通过设定阈值的方式进行周跳探测可能导致如下两类问题：①如给定的阈值过大，则对小周跳的探测不灵敏；②如给定的阈值过小，可能将正常观测值误判为有周跳发生的观测值。
[文献1]张亮，岳东杰.相位减伪距法与电离层残差法探测和修复周跳[J].测绘工程,2014,02:36-38.
[文献2]张晋升，匡翠林，谢荣安.三差法探测周跳的影响因素分析[J].测绘工程,2014,03:18-23.
[文献3]王振杰，聂志喜，欧吉坤.一种基于TurboEdit改进的GPS双频观测值周跳探测方法[J].武汉大学学报(信息科学版),2014,09:1017-1021.
[文献4]Cai Changsheng,Liu Zhizhao,Xia Pengfei,et al.Cycle SlipDetection and Repair for Undifferenced GPS Observations under HighIonospheric Activity[J].GPS Solutions,2012,doi:10.1007/s10291-012-0275-7.
[文献5]吴继忠，施闯，方荣新.TurboEdit单站GPS数据周跳探测方法的改进[J].武汉大学学报(信息科学版),2011,01:29-33.
[文献6]Liu Zhizhao.A New Automated Cycle Slip Detection and RepairMethod for a Single Dual-frequency GPS Receiver[J].J Geod,2011,85:171–183.
[文献7]Bansal,S,Banerjee,P.Smoothing of GPS-based timing data bymoving average technique utilizing fuzzy logic[C],International Conference onSignal Processing and Integrated Networks.2014:306-310.
发明内容
针对现有GNSS双频/多频观测数据周跳探测方法的不足，本发明基于模糊数学理论提出了一套单站双频周跳探测及处理方法，对难以确定是否有周跳发生的观测值进行降权处理。本发明可有效解决周跳探测阈值难以可靠确定的问题，进而改善GNSS精密数据处理的精度和可靠性。
本发明所采用的技术方案是：一种基于模糊数学的GNSS单站双频观测数据周跳探测及处理方法，其特征在于，包括以下步骤：
步骤1：获取历元间电离层延迟变化量的隶属度；
步骤2：根据历元间电离层延迟变化量的隶属度确定发生周跳的观测值；
步骤3：获取历元间MW组合观测值变化量的隶属度；
步骤4：根据历元间MW组合观测值变化量的隶属度确定发生周跳的观测值。
作为优选，步骤1的具体实现包括以下子步骤：
步骤1.1：计算所有可视卫星的历元间电离层延迟变化量(ΔLIF)i；
(ΔLIF)i＝LIF(t+1)-LIF(t)；
其中，λ1和λ2分别为历元时刻t的L1载波相位观测值、L2载波相位观测值、无距离组合观测值、L1载波相位波长和L2载波相位波长；和LIF(t+1)分别为历元时刻t+1的L1载波相位观测值、L2载波相位观测值和无距离组合观测值；
对于当前历元所有可视卫星的(ΔLIF)i，若前一历元不存在当前历元可视卫星信息，则往前遍历6个观测历元，如果仍获取不到对应的卫星信息，则该可视卫星失锁时间较长并将其标记为周跳；
对周跳检测量加权；历元间电离层延迟变化量采用高度角定权，即Pi＝sin(elei)；将(ΔLIF)i分别和权相乘得到加权后的历元间电离层延迟变化量；其中elei表示卫星高度角；
步骤1.2：获取历元间电离层延迟变化量的最大值max、最小值min；
步骤1.3：获取历元间电离层延迟变化量的模糊加权均值；
隶属函数选取如公式(1)所示，由于σ的值不会对求取模糊加权均值产生影响，故σ的初始值取σ＝max-min。公式(1)中的mean从最小值min开始遍历，循环至最大值结束，循环步长为0.001m。每次循环可得到各个(ΔLIF)i的隶属度，按公式(2)可计算出所有隶属度的总和使达到最大值的mean值便为模糊加权均值。
式中，mean、σ和μi分别为模糊均值、均方根误差和观测值的隶属度；xi为加权观测值，以(ΔLIF)i为例，xi＝(ΔLIF)i·Pi。
式中，为所有历元间电离层延迟变化量属度的总和，n为历元间电离层延迟变化量数量；
步骤1.4：获取历元间电离层延迟变化量的均方根误差σ；
其中，mean、σ和μi分别为模糊均值、均方根误差和观测值的隶属度；xi为加权观测值；n为观测值数量；
步骤1.5：获取历元间电离层延迟变化量的隶属度μi；
步骤1.1～1.4可得到模糊均值mean和均方根误差σ，将其代入至公式(1)便可得到历元间电离层延迟变化量的隶属度；
作为优选，步骤2中，当历元间电离层延迟变化量的隶属度小于指定的隶属度时，认为无距离组合观测值有周跳发生，否则对观测值按下式进行降权处理；
式中，(Pi)0和μi分别为历元间电离层延迟变化量原始权和隶属度；为周跳探测后获得的权。
作为优选，步骤3的具体实现包括以下子步骤：
步骤3.1：计算所有可视卫星的历元间MW组合观测值变化量(ΔMW)i；
其中，分别为历元时刻t的MW组合观测值、L1频率上的伪距观测值和L2频率上的伪距观测值，分别为历元时刻t+1的MW组合观测值、L1频率上的伪距观测值和L2频率上的伪距观测值，f1、f2分别为L1载波频率和L2载波频率。
步骤3.2：利用获取历元间电离层延迟变化量隶属度μi的原理，确定历元间MW组合观测值变化量的隶属度。
作为优选，步骤4中，利用MW组合观测值辅助探测周跳的阈值为：隶属度小于周跳探测的阈值0.5且|(ΔMW)i|＞3.0周。
本发明通过构造历元间电离层延迟变化量和历元间MW组合观测值变化量两个周跳检测量，并对周跳检测量进行加权，将得到的周跳检测量进行排序，剔除周跳较大的观测值，提高算法效率。通过获得周跳检测量的模糊加权均值及均方根误差，确定历元间电离层延迟变化量和历元间MW组合观测值变化量的隶属度。当历元间电离层延迟变化量的隶属度小于指定的隶属度时，认为观测值有周跳发生，否则对观测值进行降权处理；当历元间MW组合观测值变化量的隶属度小于0.5且|(ΔMW)i|＞3.0周时，判断观测值有周跳发生，MW组合观测值受伪距观测值噪声影响，不能利用其隶属度对载波相位观测值进行加权。通过对难以确定是否有周跳发生的观测值进行降权处理，可有效解决以往周跳探测阈值难以可靠确定的问题，进而改善GNSS精密数据处理的精度和可靠性。
附图说明
图1：本发明实施例的流程图。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。
请见图1，本发明提供的一种基于模糊数学的GNSS单站双频观测数据周跳探测及处理方法，通过确定历元间电离层延迟变化量和历元间MW组合观测值变化量的隶属度判断观测值是否有周跳发生，并对难以确定是否有周跳发生的观测值进行降权处理，本发明可有效解决以往周跳探测数据处理过程中阈值难以可靠确定的问题，进而改善GNSS精密数据处理的精度和可靠性。
其具体实现包括以下步骤：
步骤1：构造周跳检测量；
式(1)和式(2)分别为历元时刻t和t+1无距离组合观测值，式(2)-式(1)可得到历元间电离层延迟变化量(ΔLIF)i。
式(1)中λ1和λ2分别为历元时刻t的L1载波相位观测值、L2载波相位观测值、无距离组合观测值、L1载波相位波长和L2载波相位波长。式(2)中和LIF(t+1)分别为历元时刻t+1的L1载波相位观测值、L2载波相位观测值和无距离组合观测值。
式(3)和式(4)分别为历元时刻t和t+1的MW组合观测值，式(4)-式(3)可得到历元间MW组合观测值变化量(ΔMW)i。将以上得到的(ΔLIF)i和(ΔMW)i作为周跳检测量。
步骤2：对周跳检测量加权；
步骤1可得到当前历元所有可视卫星的(ΔLIF)i和(ΔMW)i，若前一历元不存在当前历元可视卫星信息，则往前遍历6个观测历元，如果仍获取不到对应的卫星信息，认为该可视卫星失锁时间较长并将其标记为周跳。历元间电离层延迟变化量采用高度角定权，即Pi＝sin(elei)。将(ΔLIF)i和(ΔMW)i分别和权相乘便可得到加权后的周跳检测量。
步骤3：确定历元间电离层延迟变化量的隶属度；
[文献7]中所提及的方法可从一系列离散点中剔除粗差较大的野值，本发明采用类似的方法确定(ΔLIF)i的隶属度，具体步骤如下：
步骤3.1：获取历元间电离层延迟变化量的最大值和最小值。获取最小值和最大值的主要目的是为了剔除较大的周跳值，提高算法的效率。对历元间电离层延迟变化量进行升序排序，从最小值开始遍历，当与其临近值绝对值之差小于0.1m时，该值即为最小值min；从最大值开始遍历，当与其临近值绝对值之差小于0.1m时，该值即为最大值max。
步骤3.2：获取历元间电离层延迟变化量的模糊加权均值。隶属函数选取如公式(5)所示，由于σ的值不会对求取模糊加权均值产生影响，故σ的初始值取σ＝max-min。公式(5)中的mean从最小值min开始遍历，循环至最大值结束，循环步长为0.001m。每次循环可得到各个(ΔLIF)i的隶属度，按公式(6)可计算出所有隶属度的总和使达到最大值的mean值便为模糊加权均值。
式中，mean、σ和μi分别为模糊均值、均方根误差和观测值的隶属度；xi为加权观测值，以(ΔLIF)i为例，xi＝(ΔLIF)i·Pi。
式中，为所有历元间电离层延迟变化量属度的总和，n为历元间电离层延迟变化量数量。
步骤3.3：获取历元间电离层延迟变化量的均方根误差。步骤3.2中虽可以得到比较稳健的模糊加权均值，但公式(5)中的σ是随意给定的一个初值，属于待估参数，故历元间电离层延迟变化量的隶属度仍属于未知。根据公式(7)可计算出均方根误差σ，将计算得到的σ代入到方程(5)中便可得到新的隶属度μi，如此反复迭代直至前后两次σ小于0.001m。
步骤3.4：获取历元间电离层延迟变化量的隶属度。步骤3.1～3.3可得到模糊均值mean和均方根误差，将其代入至公式(5)便可得到历元间电离层延迟变化量的隶属度。
步骤4：将历元间电离层延迟变化量作为周跳检测量确定发生周跳的观测值。模糊数学上认为隶属度为0.5时最模糊，故取0.5作为周跳探测的阈值。当历元间电离层延迟变化量的隶属度小于0.5时，认为无距离组合观测值有周跳发生，对隶属度大于等于0.5的观测值按公式(8)进行降权处理。
式中，(Pi)0和μi分别为历元间电离层延迟变化量原始权和隶属度；为周跳探测后获得的权。
λ1、λ2分别为L1、L2观测值的波长，当L1和L2上发生周跳值的比例近似等于λ2/λ1时，无距离组合观测值对该组周跳最不敏感，故必须利用MW组合观测值对该类周跳组合进行探测。
步骤5：确定历元间MW组合观测值变化量的隶属度，同步骤3。
步骤6：将历元间MW组合观测值变化量作为周跳检测量确定发生周跳的观测值；
由于MW组合观测值受伪距观测值噪声影响，既不能只靠隶属度这个阈值判断观测值是否有周跳发生，也不能利用其隶属度对载波相位观测值进行加权。本发明利用MW组合观测值辅助探测周跳的阈值为：隶属度小于0.5且|(ΔMW)i|＞3.0周。
应当理解的是，本说明书未详细阐述的部分均属于现有技术。
应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。

Claims (3)

1.一种基于模糊数学的GNSS单站双频观测数据周跳探测及处理方法，其特征在于，包括以下步骤：
步骤1：获取历元间电离层延迟变化量的隶属度；
具体实现包括以下子步骤：
步骤1.1：计算所有可视卫星的历元间电离层延迟变化量(ΔLIF)i；
(ΔLIF)i＝LIF(t+1)-LIF(t)；
其中，LIF(t)、λ1和λ2分别为历元时刻t的L1载波相位观测值、L2载波相位观测值、无距离组合观测值、L1载波相位波长和L2载波相位波长；和LIF(t+1)分别为历元时刻t+1的L1载波相位观测值、L2载波相位观测值和无距离组合观测值；
对于当前历元所有可视卫星的(ΔLIF)i，若前一历元不存在当前历元可视卫星信息，则往前遍历6个观测历元，如果仍获取不到对应的卫星信息，则该可视卫星失锁时间较长并将其标记为周跳；
对周跳检测量加权；历元间电离层延迟变化量采用高度角定权，即Pi＝sin(elei)；将(ΔLIF)i分别和权相乘得到加权后的历元间电离层延迟变化量；其中elei表示卫星高度角；
步骤1.2：获取历元间电离层延迟变化量的最大值max、最小值min；
步骤1.3：获取历元间电离层延迟变化量的模糊加权均值mean、历元间电离层延迟变化量的均方根误差σ、历元间电离层延迟变化量的隶属度μi；
式中，xi＝(ΔLIF)i·Pi；为所有历元间电离层延迟变化量隶属度的总和，n为历元间电离层延迟变化量数量；
由于σ的值不会对求取模糊加权均值产生影响，故σ的初始值取σ＝max-min；公式(1)中的mean从最小值min开始遍历，循环至最大值结束，循环步长为0.001m；每次循环可得到各个(ΔLIF)i的隶属度，按公式(2)可计算出所有隶属度的总和使达到最大值的mean值便为模糊加权均值；
将模糊均值mean和均方根误差σ代入至公式(1)便可得到历元间电离层延迟变化量的隶属度；
步骤2：根据历元间电离层延迟变化量的隶属度确定发生周跳的观测值；
步骤3：获取历元间MW组合观测值变化量的隶属度；
具体实现包括以下子步骤：
步骤3.1：计算所有可视卫星的历元间MW组合观测值变化量(ΔMW)i；
其中，P1(t)、P2(t)分别为历元时刻t的MW组合观测值、L1频率上的伪距观测值和L2频率上的伪距观测值，P1(t+1)、P2(t+1)分别为历元时刻t+1的MW组合观测值、L1频率上的伪距观测值和L2频率上的伪距观测值，f1、f2分别为L1载波频率和L2载波频率；
步骤3.2：利用获取历元间电离层延迟变化量隶属度μi的原理，确定历元间MW组合观测值变化量的隶属度；
步骤4：根据历元间MW组合观测值变化量的隶属度确定发生周跳的观测值。
2.根据权利要求1所述的基于模糊数学的GNSS单站双频观测数据周跳探测及处理方法，其特征在于：步骤2中，当历元间电离层延迟变化量的隶属度小于指定的隶属度时，认为无距离组合观测值有周跳发生，否则对观测值按下式进行降权处理；
式中，(Pi)0和μi分别为历元间电离层延迟变化量的原始权和隶属度；为周跳探测后获得的权。
3.根据权利要求1所述的基于模糊数学的GNSS单站双频观测数据周跳探测及处理方法，其特征在于：步骤4中，利用MW组合观测值辅助探测周跳的阈值为：隶属度小于周跳探测的阈值0.5且|(ΔMW)i|＞3.0周。
CN201610518199.6A 2016-07-04 2016-07-04 基于模糊数学的gnss单站双频观测数据周跳探测及处理方法 Active CN106199659B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610518199.6A CN106199659B (zh) 2016-07-04 2016-07-04 基于模糊数学的gnss单站双频观测数据周跳探测及处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610518199.6A CN106199659B (zh) 2016-07-04 2016-07-04 基于模糊数学的gnss单站双频观测数据周跳探测及处理方法

Publications (2)

Publication Number Publication Date
CN106199659A CN106199659A (zh) 2016-12-07
CN106199659B true CN106199659B (zh) 2018-08-21

Family

ID=57466035

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610518199.6A Active CN106199659B (zh) 2016-07-04 2016-07-04 基于模糊数学的gnss单站双频观测数据周跳探测及处理方法

Country Status (1)

Country Link
CN (1) CN106199659B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932793B (zh) * 2017-03-31 2019-07-23 武汉大学 一种北斗三频信号的实时周跳探测与修复方法
CN107561563B (zh) * 2017-08-30 2020-01-10 湖南航天电子科技有限公司 奇异点保留滤波降噪的周跳探测方法
CN110727000B (zh) * 2019-11-19 2022-09-13 西南交通大学 一种基于gnss高采样率数据的小周跳修复方法
EP4421534A1 (en) * 2021-12-30 2024-08-28 Huawei Technologies Co., Ltd. Data preprocessing method, data preprocessing apparatus, and chip
CN114814907A (zh) * 2022-04-07 2022-07-29 山东大学 一种基于弧段划分的北斗双频周跳探测与修复方法
CN114488228B (zh) * 2022-04-11 2022-07-01 南京北斗创新应用科技研究院有限公司 一种适用于动态载体平台的gnss多路径误差削弱方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105549046A (zh) * 2015-12-04 2016-05-04 航天恒星科技有限公司 Gnss接收机周跳探测与修复处理方法
CN105652298A (zh) * 2016-01-12 2016-06-08 广州市城市规划勘测设计研究院 一种bds三频伪距相位组合的周跳探测修复方法及装置
CN105676243A (zh) * 2016-01-11 2016-06-15 昆明理工大学 一种基于无几何相位和电离层残差法的北斗三频周跳探测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5553909B2 (ja) * 2010-12-13 2014-07-23 古野電気株式会社 サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるｇｎｓｓ受信装置、同ｇｎｓｓ受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105549046A (zh) * 2015-12-04 2016-05-04 航天恒星科技有限公司 Gnss接收机周跳探测与修复处理方法
CN105676243A (zh) * 2016-01-11 2016-06-15 昆明理工大学 一种基于无几何相位和电离层残差法的北斗三频周跳探测方法
CN105652298A (zh) * 2016-01-12 2016-06-08 广州市城市规划勘测设计研究院 一种bds三频伪距相位组合的周跳探测修复方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"不确定性数学理论在GNSS精确定位中的研究";康荣雷;《中国优秀硕士学位论文全文数据库基础科学辑》;20091115(第11期);A008-38 *
"智能GPS软件接收机载波跟踪环路设计";唐斌等;《北京航空航天大学学报》;20070731;第33卷(第7期);807-810 *

Also Published As

Publication number Publication date
CN106199659A (zh) 2016-12-07

Similar Documents

Publication Publication Date Title
CN106199659B (zh) 基于模糊数学的gnss单站双频观测数据周跳探测及处理方法
CN104749594B (zh) 一种gps双频非差周跳探测与修复方法及装置
CN110646820B (zh) Rtk定位数据的质检方法、装置、设备和存储介质
CN111077550A (zh) 一种应用于智能终端rtd定位的粗差探测方法及系统
CN108107455A (zh) 一种基于相位跳变的卫星钟差实时预报方法
CN107505642B (zh) 一种ins辅助的实时bds单频周跳探测方法
Feng et al. Integrity monitoring for carrier phase ambiguities
CN107796484A (zh) 一种基于北斗导航定位系统信噪比数据观测水位变化方法
JP2011247758A (ja) 位置算出方法及び位置算出装置
KR101433908B1 (ko) Gnss 데이터 품질 평가 방법 및 시스템
CN111505694A (zh) 一种面向机载的bds-3三天线阵多频点测姿方法
CN105652298B (zh) 一种bds三频伪距相位组合的周跳探测修复方法及装置
CN118226489A (zh) Gnss短时中断下基于广义最大相关熵的gnss/ins组合导航方法
CN111123315A (zh) 非差非组合ppp模型的优化方法及装置、定位系统
CN115728793B (zh) 一种基于dia理论的精密单点定位粗差探测与处理方法
Zhang et al. A Bayesian method of GNSS cycle slips detection based on ARMA model
Shi et al. A new real‐time cycle slip detection and repair approach based on BDS dual‐frequency carrier phase and Doppler observations
CN107462909A (zh) 一种北斗单频载波相位的周跳探测方法
CN110967706A (zh) 基于改进无迹卡尔曼粒子滤波的raim方法
CN115308781A (zh) 基于bdgim辅助的相位平滑伪距高精度时间传递方法
CN114252896A (zh) 一种单频实时精密单点定位方法
Closas et al. Integrity measures in direct-positioning
CN112558112B (zh) 一种gnss状态域慢变斜坡故障完好性监测方法
CN115390119A (zh) 一种应用于洪水环境下的北斗精密单点定位方法
CN115267845A (zh) 一种差分码偏差的获取方法、计算机设备和可读存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant