CN113514856B - 周跳检测方法、装置、电子设备及可读储存介质 - Google Patents
周跳检测方法、装置、电子设备及可读储存介质 Download PDFInfo
- Publication number
- CN113514856B CN113514856B CN202110388651.2A CN202110388651A CN113514856B CN 113514856 B CN113514856 B CN 113514856B CN 202110388651 A CN202110388651 A CN 202110388651A CN 113514856 B CN113514856 B CN 113514856B
- Authority
- CN
- China
- Prior art keywords
- satellites
- satellite
- cycle slip
- cycle
- determining
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 47
- 238000001914 filtration Methods 0.000 claims description 19
- 238000000034 method Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 11
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 239000005433 ionosphere Substances 0.000 description 16
- 230000000694 effects Effects 0.000 description 7
- 238000007796 conventional method Methods 0.000 description 4
- 235000008694 Humulus lupulus Nutrition 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
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/23—Testing, monitoring, correcting or calibrating of receiver elements
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
本发明提供了一种周跳检测方法、装置、电子设备及可读储存介质,其中,周跳检测方法,包括:获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,N为大于1的整数;依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,所述第二卫星为初始判定存在周跳的第一卫星;依据每一所述第二卫星对应的频点的数量,确定M个循环,M为大于1的整数;在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳。本发明能够有效降低周跳的虚假判断,避免了模糊度参数频繁的初始化,从而有助于提高卫星定位性能。
Description
技术领域
本发明涉及导航卫星技术领域,具体地,涉及一种周跳检测方法、装置、电子设备及可读储存介质。
背景技术
全球导航卫星系统(Global Navigation Satellite System,GNSS)能够提供全球性、全天候、三维的高精度定位、导航和授时服务。随着自动化和精细作业的需求,用户对定位精度的需求也与日俱增。
电离层误差作为GNSS最棘手的误差源之一,电离层通常距离地球表面48-965km,是日地空间环境的重要组成之一。电离层作为一种传播介质对无线电波产生反射、折射、散射和吸收作用,因而对GNSS的信号传播起着重要影响。
在电离层平静的状态下,电离层误差在时空域表现出稳定的变化特性,然而在电离层扰动状态下,电离层误差表现出不规则的变化特性,研究全球导航卫星系统在电离层扰动下的定位性对获得稳定可靠的定位结果能具有重要意义。
由于载波相位观测值比伪距观测值的精度高100倍左右,载波相位在卫星高精度定位中扮演着重要作用,然而载波相位受到周跳影响,当受到信号受到遮挡或者电离层的扰动情况下,存在不连续观测。常规判断周跳的方法是采用几何无关组合法,即电离层残差法来判断周跳。然而在电离层的扰动下,电离层残差较大,很容易错误的判断周跳的发生,错误的周跳判断会使得模糊度参数频繁的初始化,从而影响最后高精度定位的性能。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种周跳检测方法、装置、电子设备及可读储存介质。
根据本发明提供的周跳检测方法,包括:
获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,N为大于1的整数;
依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,所述第二卫星为初始判定存在周跳的第一卫星;
依据每一所述第二卫星对应的频点的数量,确定M个循环,M为大于1的整数;
在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳。
本发明实施例还提供了一种周跳检测装置,包括:
获取模块,用于获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,N为大于1的整数;
第一确定模块,用于依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,所述第二卫星为初始判定存在周跳的第一卫星;
第二确定模块,用于依据每一所述第二卫星对应的频点的数量,确定M个循环,M为大于1的整数;
检测模块,用于在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳。
本发明实施例还提供了一种电子设备,所述设备包括:处理器以及存储有计算机程序指令的存储器;
所述处理器执行所述计算机程序指令时实现上述的周跳检测方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现上述的周跳检测方法。
本发明实施例提供的周跳检测方法,获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,依据每一所述第二卫星对应的频点的数量,确定M个循环,在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳。本发明实施例中可以认为是对初步的周跳判定结果进行了进一步的验证,有效降低周跳的虚假判断,避免了模糊度参数频繁的初始化,从而有助于提高卫星定位性能。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明实施例提供的周跳检测方法的流程示意图;
图2为本发明实施例提供的周跳检测方法在一个具体应用例中的流程示意图;
图3为常规方法和本发明实施例所对应的迭代滤波探测周跳数目的对比图;
图4为常规方法的定位效果和本发明实施例对应的迭代滤波的定位效果的对比图;
图5为本发明实施例提供的周跳检测装置的结构示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
在本发明的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
如图1所示,本发明实施例提供的周跳检测方法,包括:
步骤101,获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,N为大于1的整数;
步骤102,依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,所述第二卫星为初始判定存在周跳的第一卫星;
步骤103,依据每一所述第二卫星对应的频点的数量,确定M个循环,M为大于1的整数;
步骤104,在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳。
本实施例中,目标历元可以对应的是某一个观测时刻,在该观测时刻下可以基于卫星信号,来获取到初始观测值。初始观测值中可以确定出可用卫星数,对应了上述的N个第一卫星。
容易理解的是,依据目标历元的初始观测值,基于常规的周跳检测方法,实际可以判断出是否存在周跳的情况,例如基于载波相位的几何无关组合(geometry-freecombination,GF)或者是宽巷相位减窄巷伪距组合(Melbourne-Wübbena,MW),来判断是否存在周跳;其中,GF又可以成为电离层残差组合。
然而,如上文中所示的,基于这些常规的周跳检测方法,在电离层的扰动下,电离层残差较大,很容易错误的判断周跳的发生,错误的周跳判断会使得模糊度参数频繁的初始化,从而影响最后高精度定位的性能。
因此,本实施例中,可以是通过这些常规的周跳检测方法初步检测周跳,确定出存在周跳的第一卫星,也就是上述的第二卫星。容易理解的是,本实施例中,可以是对初步的周跳检测结果进一步对周跳进行确认。
具体地,可以根据每一所述第二卫星对应的频点的数量,确定M个循环,并在M个循环中依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,进一步检测所述初始观测值是否存在周跳。
举例来说,对于一些卫星,可以会存在一个或多个频点,例如某一颗卫星可能会使用到L1频点与L2频点等;在一个可行的实施方式中,可以将各个第二卫星对应的频点数量相加,得到循环的数量,即上述的M的值。
而在M个循环中依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,进一步检测所述初始观测值是否存在周跳,可以基于例如下面的原理进行:首先,可用的初始观测值数量并不会因为假设不存在周跳而提出更多的卫星;这是因为如果不存在周跳,那么观测值应该是一致的。其次,载波相位的前后两次估计差别不应该太大,具体是表达式应该是差值的绝对值在基于中误差确定的某一数值范围以内。再者,是所有残差也应该不会变化太大。
当以上条件均满足,或者说满足其中的某几项时,一般可以认为不存在周跳的假设成立,反之,则可以认为确实存在周跳的情况。也就是说,基于以上条件的判断,可以检测出所述初始观测值是否存在周跳。
本发明实施例提供的周跳检测方法,获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,依据每一所述第二卫星对应的频点的数量,确定M个循环,在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳。本发明实施例中可以认为是对初步的周跳判定结果进行了进一步的验证,有效降低周跳的虚假判断,避免了模糊度参数频繁的初始化,从而有助于提高卫星定位性能。
可选地,上述步骤101,获取目标历元的初始观测值,包括:
获取目标历元的原始观测值;
对所述原始观测值滤波得到第一卫星的数量N,以及每一所述第一卫星对应的位置参数与残差。
本实施例中,原始观测值可以认为是原始的卫星信息,通过对原始观测值的滤波,得到滤波结果,也就是上述的第一卫星的数量N,以及每一所述第一卫星对应的位置参数与残差。具体可以是通过例如卡尔曼滤波等常用的滤波方式对原始观测值进行滤波的。
本实施例通过滤波来获得第一卫星的数量N,以及每一所述第一卫星对应的位置参数与残差,有助于提高后续周跳检测结果的准确度。
可选地,所述依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,包括:
依据所述初始观测值,采用几何无关组合,和/或,宽巷相位减窄巷伪距组合,从所述N个第一卫星中确定出第二卫星。
以下结合一个举例对基于几何无关组合GF确定跳变的过程进行说明:
常规的周跳方法是载波相位的几何无关组合GF(geometry-free combination),也称为电离层残差组合,载波观测值的无几何观测值如以下公式表示:
LGF=λ1Φ1-λ2Φ2=(1-γ2)I1+λ1N1-λ2N2+εGF
其中LGF是载波相位的几何无关组合,ΔLGF是前后历元间的差值变化,Φi在第i种频率上的载波相位,λi则是第i种频率上对应的波长,i的取值可以为1或2,Ni是第i种频率上的模糊度参数,t是当前历元,I1是在L1频点上的电离层误差,/>表示在当前历元下在L1频点上的电离层误差,εGF是观测值噪声和多路径等误差。
容易理解的是,常规GF方法的缺点有两点。第一,需要对ΔLGF选择合适的阈值来判断是否发生周跳,GF组合前提是假设电离层残差稳定变化,这样历元间差分会变化较小,如果出现较大残差,则认为存在周跳。然而电离层残差跟电离层的活跃程度密切相关,可见很难确定合适的阈值来判断是否发生了周跳。第二,一旦某颗卫星判断为存在周跳,该卫星的所有频点认为都存在周跳,所有频点的模糊度都需要初始化。然而研究表明不同L1频点相对L2频点在扰动情况下会表现的更加稳健。
而结合上文实施例中的描述,本发明实施例中,一方面,可以有效降低周跳的虚假判断,尤其在电离层扰动下,给的GF探测阈值比较保守的情况下。另一方面,由于在确定循环数量时,考虑了第二卫星的频点数量,可以有效区分不同频点是否发生周跳;另外,在一些可行的实施方式中,还可以进一步对每个频点进行标记,从而根究需求逐一进行初始化,这样可以使得某一频点数据对定位有更多的贡献。
相应地,在考虑第二卫星的频点数量的时,所述依据N与每一所述第二卫星对应的频点的数量,确定M个循环,可以包括:
在每一所述第二卫星对应的频点的数量均为P的情况下,将N与P的乘积确定为M的值,P为大于1的整数。
通常来说,第二卫星的频点的数量是一致的,例如,各第二卫星均可以包括第一频点与第二频点,设第二卫星对应的频点的数量为P,则可以将N和P的乘积来作为循环的数量。
可选地,上述在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳,包括:
在满足如下全部条件的情况下,确定初始观测值不存在周跳:
|N(i)-N(i-1)|≥α{Var_N(i)+Var_N(i-1)}
stdi/stdi-1≥β
其中,i为所述循环的序号,且1<i≤M;表示在第i个循环中假设在不存在周跳的情况下,确定出的卫星的数量;N(i)表示在第i个循环中依据所述位置参数确定的载波相位;Var_N(i)表示在第i个循环中所述载波相位的中误差;stdi表示在第i个循环中的残差;α与β均为预设的系数。
本实施例中,限定了同时满足以上三个条件的情况下,认为不存在周跳的假设成立,而在不满足至少一个条件的情况下,可以认为初始观测值中存在周跳。通过多个条件的限定,可以有效降低周跳的虚假判断。
在一个示例中,上述α与β可以分别取值为3与1.5;当然,两者的取值也可以进一步根据实际需要进行调整。
可选地,所述在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳之后,所述方法还包括:
在所述初始观测值存在周跳的情况下,基于所述目标历元的前一历元的观测值,确定所述目标历元的目标观测值;
在所述初始观测值不存在周跳的情况下,将所述目标历元的初始观测值,确定为所述目标历元的目标观测值。
换而言之,本实施例中,可以在认为不存在周跳的假设成立时,采用当前的滤波结果;不存在周跳的假设不成立时,采用上次滤波,也就是目标历元的前一历元的观测值的滤波结果。
可选地,所述依据每一所述第二卫星对应的频点的数量,确定M个循环,包括:
获取在目标历元之前的预设历元范围内出现周跳的历史数目;
在所述历史数目大于周跳数目阈值的情况下,依据每一所述第二卫星对应的频点的数量,确定M个循环。
容易理解的是,本发明实施例提供的周跳检测方法,相对于常规的周跳检测方法,需要额外的增加计算量,实际应用中,可以对定位效果和计算量之间的权衡。当周跳数目只有1-2个时,几何结构变化较小,对定位的最终影响较小,可以选择不用进行上述的循环迭代过程。但是一旦出现周跳数目较多的情况下,采用本发明实施例提供的周跳检测方法可以进一步判断周跳的发生与否,避免频繁的错误模糊度初始化,进而提高定位精度。
参见图2,以下结合一具体应用例对本发明实施例提供的周跳检测方法进行说明,在该应用例中,周跳检测方法包括:
步骤S11,获取当前观测值,记为t0;
步骤S12,基于当前观测值,使用常规GF和/或MW方法来判断周跳;
步骤S13,采用滤波方法来计算在现有周跳判断方法下,对应的滤波所有的位置参数、可用卫星数以及残差等信息,记为滤波结果solution0;
步骤S14,按照每颗可用卫星乘以对应的频点,进行固定次数的循环;
在可以卫星数的数量为n,且每一可用卫星对应的频点数量为nfreq的情况下,循环次数可以记为n*nfreq;
步骤S15,提出假设:假设没有周跳,将模糊度方差设置为0;
步骤S16,计算第i次循环中,因为假设不存在周跳而提出的卫星的数量、载波相位以及残差,记为solution_i;
步骤S17,对比solution_i与solution_i-1判断提出的卫星的数量、载波相位以及残差是否分别满足对应的预设条件;如满足,则认为假设成立,采用当前的滤波结果,若不满足,则认为假设不成立,采用上次滤波结果;
步骤S18,一直循环,直到循环次数结束。
图3示出了常规方法(对应conventional GF method of 2476)和本发明实施例所对应的迭代滤波探测(对应improved interative filter of 418)周跳数目的对比图;其中,number of cycle slips代表了周跳数目,UTC表示历元,单位小时;
图4示出了常规方法的定位效果和本发明实施例对应的迭代滤波的定位效果的对比图,其中,2D accuracy表示二维精度,height accuracy表示高度精度,3D accuracy表示三维精度,time为时间。
可见,本发明实施例提供的周跳检测方法可以有效降低周跳的虚假判断,尤其在电离层扰动下,给的GF探测阈值比较保守的情况下;可以有效区分不同频点是否发生周跳,并合理利用某一频点数据的贡献;在避免频繁的模糊度初始化的情况下,定位精度有明显的提升。
如图5所示,本发明实施例还提供了一种周跳检测装置,包括:
获取模块501,用于获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,N为大于1的整数;
第一确定模块502,用于依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,所述第二卫星为初始判定存在周跳的第一卫星;
第二确定模块503,用于依据每一所述第二卫星对应的频点的数量,确定M个循环,M为大于1的整数;
检测模块504,用于在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳。
可选地,获取模块501,包括:
第一获取单元,用于获取目标历元的原始观测值;
第二获取单元,用于对所述原始观测值滤波得到第一卫星的数量N,以及每一所述第一卫星对应的位置参数与残差。
可选地,第一确定模块502,可具体用于:
依据所述初始观测值,采用几何无关组合GF,和/或,宽巷相位减窄巷伪距组合MW,从所述N个第一卫星中确定出第二卫星。
可选地,第二确定模块503可具体用于:
在每一所述第二卫星对应的频点的数量均为P的情况下,将N与P的乘积确定为M的值,P为大于1的整数。
可选地,检测模块504,可以具体用于:
在满足如下全部条件的情况下,确定初始观测值不存在周跳:
|N(i)-N(i-1)|≥α{Var_N(i)+Var_N(i-1)}
其中,i为所述循环的序号,且1<i≤M;表示在第i个循环中假设在不存在周跳的情况下,确定出的卫星的数量;N(i)表示在第i个循环中依据所述位置参数确定的载波相位;Var_N(i)表示在第i个循环中所述载波相位的中误差;stdi表示在第i个循环中的残差;α与β均为预设的系数。
可选地,周跳检测装置还包括:
第三确定模块,用于在所述初始观测值存在周跳的情况下,基于所述目标历元的前一历元的观测值,确定所述目标历元的目标观测值;
第四确定模块,用于在所述初始观测值不存在周跳的情况下,将所述目标历元的初始观测值,确定为所述目标历元的目标观测值。
可选地,第二确定模块503,可以包括:
第三获取单元,用于获取在目标历元之前的预设历元范围内出现周跳的历史数目;
确定单元,用于在所述历史数目大于周跳数目阈值的情况下,依据每一所述第二卫星对应的频点的数量,确定M个循环。
需要说明的是,本发明实施例提供的周跳检测装置是能够执行上述周跳检测方法的装置,故上述周跳检测方法的所有实施例均适用于周跳检测装置,且均能达到相同或相似的有益效果。
本发明实施例还提供了一种电子设备,所述设备包括:处理器以及存储有计算机程序指令的存储器;
所述处理器执行所述计算机程序指令时实现上述的周跳检测方法。
本发明实施例还提供了一种计算机可读存储介质,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现上述的周跳检测方法。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (9)
1.一种周跳检测方法,其特征在于,包括:
获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,N为大于1的整数;
依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,所述第二卫星为初始判定存在周跳的第一卫星;
依据每一所述第二卫星对应的频点的数量,确定M个循环,M为大于1的整数;
在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳;
所述在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳,包括:
在满足如下全部条件的情况下,确定初始观测值存在周跳:
|N(i)-N(i-1)|≥α{Var_N(i)+Var_N(i-1)}
其中,i为所述循环的序号,且1<i≤M;表示在第i个循环中假设在不存在周跳的情况下,确定出的卫星的数量;N(i)表示在第i个循环中依据所述位置参数确定的载波相位;Var_N(i)表示在第i个循环中所述载波相位的中误差;stdi表示在第i个循环中的残差;α与β均为预设的系数。
2.根据权利要求1所述的方法,其特征在于,所述获取目标历元的初始观测值,包括:
获取目标历元的原始观测值;
对所述原始观测值滤波得到第一卫星的数量N,以及每一所述第一卫星对应的位置参数与残差。
3.根据权利要求1所述的方法,其特征在于,所述依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,包括:
依据所述初始观测值,采用几何无关组合GF,和/或,宽巷相位减窄巷伪距组合MW,从所述N个第一卫星中确定出第二卫星。
4.根据权利要求3所述的方法,其特征在于,依据N与每一所述第二卫星对应的频点的数量,确定M个循环,包括:
在每一所述第二卫星对应的频点的数量均为P的情况下,将N与P的乘积确定为M的值,P为大于1的整数。
5.根据权利要求1所述的方法,其特征在于,所述在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳之后,所述方法还包括:
在所述初始观测值存在周跳的情况下,基于所述目标历元的前一历元的观测值,确定所述目标历元的目标观测值;
在所述初始观测值不存在周跳的情况下,将所述目标历元的初始观测值,确定为所述目标历元的目标观测值。
6.根据权利要求1所述的方法,其特征在于,所述依据每一所述第二卫星对应的频点的数量,确定M个循环,包括:
获取在目标历元之前的预设历元范围内出现周跳的历史数目;
在所述历史数目大于周跳数目阈值的情况下,依据每一所述第二卫星对应的频点的数量,确定M个循环。
7.一种周跳检测装置,其特征在于,包括:
获取模块,用于获取目标历元的初始观测值,所述初始观测值包括N个第一卫星中每一所述第一卫星对应的位置参数与残差,N为大于1的整数;
第一确定模块,用于依据所述初始观测值,从所述N个第一卫星中确定出第二卫星,所述第二卫星为初始判定存在周跳的第一卫星;
第二确定模块,用于依据每一所述第二卫星对应的频点的数量,确定M个循环,M为大于1的整数;
检测模块,用于在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳;
所述在所述M个循环中,依据所述第二卫星的数量、以及与所述第二卫星对应的位置参数与残差,检测所述初始观测值是否存在周跳,包括:
在满足如下全部条件的情况下,确定初始观测值存在周跳:
|N(i)-N(i-1)|≥α{Var_N(i)+Var_N(i-1)}
其中,i为所述循环的序号,且1<i≤M;表示在第i个循环中假设在不存在周跳的情况下,确定出的卫星的数量;N(i)表示在第i个循环中依据所述位置参数确定的载波相位;Var_N(i)表示在第i个循环中所述载波相位的中误差;stdi表示在第i个循环中的残差;α与β均为预设的系数。
8.一种电子设备,其特征在于,所述设备包括:处理器以及存储有计算机程序指令的存储器;
所述处理器执行所述计算机程序指令时实现如权利要求1-6任意一项所述的周跳检测方法。
9.一种计算机存储介质,其特征在于,所述计算机存储介质上存储有计算机程序指令,所述计算机程序指令被处理器执行时实现如权利要求1-6任意一项所述的周跳检测方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110388651.2A CN113514856B (zh) | 2021-04-12 | 2021-04-12 | 周跳检测方法、装置、电子设备及可读储存介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110388651.2A CN113514856B (zh) | 2021-04-12 | 2021-04-12 | 周跳检测方法、装置、电子设备及可读储存介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113514856A CN113514856A (zh) | 2021-10-19 |
CN113514856B true CN113514856B (zh) | 2024-01-09 |
Family
ID=78062672
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110388651.2A Active CN113514856B (zh) | 2021-04-12 | 2021-04-12 | 周跳检测方法、装置、电子设备及可读储存介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113514856B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288978A (zh) * | 2011-07-20 | 2011-12-21 | 东南大学 | 一种cors基站周跳探测与修复方法 |
CN103529462A (zh) * | 2013-10-21 | 2014-01-22 | 西南交通大学 | 一种用于全球导航卫星系统的动态周跳探测与修复方法 |
CN106168672A (zh) * | 2016-01-01 | 2016-11-30 | 广州中海达卫星导航技术股份有限公司 | 一种gnss多模单频rtk周跳探测方法及装置 |
CN107132558A (zh) * | 2017-06-13 | 2017-09-05 | 武汉大学 | 惯性辅助的多频多模gnss周跳修复方法及系统 |
WO2018138904A1 (ja) * | 2017-01-30 | 2018-08-02 | 三菱電機株式会社 | 測位装置および測位方法 |
CN109884678A (zh) * | 2018-12-29 | 2019-06-14 | 北方信息控制研究院集团有限公司 | 实时动态新星模糊度快速固定的方法 |
CN111505694A (zh) * | 2020-05-07 | 2020-08-07 | 中航机载系统共性技术有限公司 | 一种面向机载的bds-3三天线阵多频点测姿方法 |
CN111781614A (zh) * | 2020-08-25 | 2020-10-16 | 中国石油大学(华东) | 一种周跳探测与修复方法、装置、设备及可存储介质 |
CN111947667A (zh) * | 2020-06-24 | 2020-11-17 | 火眼位置数智科技服务有限公司 | 一种基于运动学和动力学组合的低轨卫星实时高精度定轨方法 |
CN112505733A (zh) * | 2020-12-07 | 2021-03-16 | 北京理工雷科电子信息技术有限公司 | 一种适应于双天线动态定向的联合周跳探测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9651667B2 (en) * | 2007-06-22 | 2017-05-16 | Trimble Inc. | Combined cycle slip indicators for regionally augmented GNSS |
EP2333581A1 (en) * | 2009-12-09 | 2011-06-15 | Technische Universität München | Estimation of phase and code biases on multiple frequencies with a Kalman filter |
-
2021
- 2021-04-12 CN CN202110388651.2A patent/CN113514856B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288978A (zh) * | 2011-07-20 | 2011-12-21 | 东南大学 | 一种cors基站周跳探测与修复方法 |
CN103529462A (zh) * | 2013-10-21 | 2014-01-22 | 西南交通大学 | 一种用于全球导航卫星系统的动态周跳探测与修复方法 |
CN106168672A (zh) * | 2016-01-01 | 2016-11-30 | 广州中海达卫星导航技术股份有限公司 | 一种gnss多模单频rtk周跳探测方法及装置 |
WO2018138904A1 (ja) * | 2017-01-30 | 2018-08-02 | 三菱電機株式会社 | 測位装置および測位方法 |
CN107132558A (zh) * | 2017-06-13 | 2017-09-05 | 武汉大学 | 惯性辅助的多频多模gnss周跳修复方法及系统 |
CN109884678A (zh) * | 2018-12-29 | 2019-06-14 | 北方信息控制研究院集团有限公司 | 实时动态新星模糊度快速固定的方法 |
CN111505694A (zh) * | 2020-05-07 | 2020-08-07 | 中航机载系统共性技术有限公司 | 一种面向机载的bds-3三天线阵多频点测姿方法 |
CN111947667A (zh) * | 2020-06-24 | 2020-11-17 | 火眼位置数智科技服务有限公司 | 一种基于运动学和动力学组合的低轨卫星实时高精度定轨方法 |
CN111781614A (zh) * | 2020-08-25 | 2020-10-16 | 中国石油大学(华东) | 一种周跳探测与修复方法、装置、设备及可存储介质 |
CN112505733A (zh) * | 2020-12-07 | 2021-03-16 | 北京理工雷科电子信息技术有限公司 | 一种适应于双天线动态定向的联合周跳探测方法 |
Non-Patent Citations (7)
Title |
---|
GNSS周跳探测与修复融合算法研究;万军;《中国优秀硕士学位论文全文数据库 基础科学辑》;1-80 * |
GNSS整数模糊度估计与检验的理论和方法研究;张晶宇;《中国博士学位论文全文数据库 基础科学辑》;1-202 * |
Kalman-filter-based undifferenced cycle slip estimation in real-time precise point positioning;Pan Li, Xinyuan Jiang;《GPS Solutions 》;1-13 * |
Multi-GNSS triplefrequency differential code bias (DCB) determination with precise point positioning (PPP);Liu, T., Zhang, B., Yuan, Y., Li, Z., Wang, N.;《J. Geod》;765–784 * |
Zhang Xiaohong & Li Xingxing.《GPS Solutions》.2011,315–327 . * |
星载双频GPS载波相位周跳的实时探测与修复;白伟华 等;《数据采集与处理》;593-599 * |
综合利用双频观测值进行周跳的探测与修复;王威;隋立芬;符养;刘麦喜;《海洋测绘》;24-27 * |
Also Published As
Publication number | Publication date |
---|---|
CN113514856A (zh) | 2021-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9035826B2 (en) | Satellite differential positioning receiver using multiple base-rover antennas | |
US8704708B2 (en) | GNSS signal processing methods and apparatus with scaling of quality measure | |
JP6086901B2 (ja) | 複数のrtkエンジンを有するgnss測量受信機 | |
JP5421903B2 (ja) | 部分探索搬送波位相整数アンビギュイティ決定 | |
RU2451947C2 (ru) | Способ использования трех gps-частот для решения неопределенности фазы несущей | |
JP2013019893A (ja) | バイアスがかかった測定値に基づいた衛星航法システムの誤り検出 | |
EP0948748A1 (en) | Gps multipath detection method and system | |
WO2009014823A2 (en) | Methods and apparatus for geometry extra-redundant almost fixed solutions | |
JP2017111112A (ja) | 長周期波高測定用gnss基盤ブイにおける受信機の位相ワインドアップ補償のためのppp及びマグネチックコンパス統合システム{ppp and magnetic compass integration system for measuring of long period wave in buoy’s wind−up} | |
Ziedan | Improved multipath and NLOS signals identification in urban environments | |
CN105158778A (zh) | 多系统联合实施载波相位差分故障卫星剔除方法及其系统 | |
CN116819580A (zh) | 一种惯性辅助的双天线gnss海上船舶姿态确定方法 | |
CN113514856B (zh) | 周跳检测方法、装置、电子设备及可读储存介质 | |
CN114488235A (zh) | 一种双天线卫星定向方法、系统、存储介质和电子设备 | |
Zheng et al. | Multipath mitigation for improving GPS narrow-lane uncalibrated phase delay estimation and speeding up PPP ambiguity resolution | |
CN112462398A (zh) | 一种卫星信号失锁重定位方法及系统 | |
CN116953741B (zh) | 一种应用于全球导航卫星系统gnss的周跳探测与修复方法 | |
CN117607910B (zh) | 一种基于矢量跟踪新息矢量的欺骗检测方法及系统 | |
CN111830542B (zh) | 数据处理方法、装置、定位设备及存储介质 | |
CN113865592B (zh) | 适用于多频gnss精密导航定位中的多路径参数化方法及存储介质 | |
CN114114348A (zh) | 海上基线约束卡尔曼滤波精密单点定位方法 | |
CN117930290A (zh) | 一种周跳探测方法、控制装置、终端及存储介质 | |
CN111352135A (zh) | 伪距粗差的探测方法、装置及可读存储介质 | |
CN115755130A (zh) | 通过固定和确认载波距离模糊度来进行gnss定位 | |
CN115327585A (zh) | Gnss观测值定权方法、定位方法、装置、终端及介质 |
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 |