CN111190200A - 一种动态环境下单频周跳检测与修复方法 - Google Patents
一种动态环境下单频周跳检测与修复方法 Download PDFInfo
- Publication number
- CN111190200A CN111190200A CN201911248588.1A CN201911248588A CN111190200A CN 111190200 A CN111190200 A CN 111190200A CN 201911248588 A CN201911248588 A CN 201911248588A CN 111190200 A CN111190200 A CN 111190200A
- Authority
- CN
- China
- Prior art keywords
- cycle slip
- matrix
- satellite
- positioning
- value
- 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/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/23—Testing, monitoring, correcting or calibrating of receiver elements
- G01S19/235—Calibration of receiver components
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
本发明涉及一种动态环境下单频周跳检测与修复方法,属于卫星导航高精度定位领域。本发明有一定的数理统计理论基础,与其他周跳检测修复方法需要双频点、不能同时修复大周跳和小周跳、需逐个剔除卫星后再次计算相比,该方法可以直接计算出各个卫星周跳数值大小及其对定位的影响,计算量大大减少,且通过试验表明该算法可以正确的检测和修复卫星周跳,是一种便捷快速有效的动态环境下单频周跳检测与修复方法。
Description
技术领域
本发明涉及一种动态环境下单频周跳检测与修复方法,属于卫星导航高精度定位领域。
背景技术
在使用载波相位进行高精度定位的时候,关键是实时整周模糊度的确定。由于接收机在求解出整周模糊度之后,在后续观测期间,可能受到多种干扰影响,导致接收机对卫星的载波观测值可能发生周跳。周跳会影响定位结果的可靠性,并且由于引入较大误差,若不修复会影响之后所有时刻的定位精度,因此持续的周跳检测与实时修复是高精度卫星导航定位系统中不可缺少的关键部分。
静态环境下,由于接收机相对静止,伪距载波等观测量变化较小,环境稳定,且位置无变化,因此易于周跳检测与修复。动态环境下,由于接收机载体运动复杂,受到信号阻挡、环境变化及大气扰动等影响,可能导致多周跳和连续周跳等现象出现,使得动态环境下周跳探测与修复变得更加困难。
目前周跳检测与修复的方法有很多,如多项式拟合法、高次差法、电离层残差法、伪距载波相位组合法等,但是这些方法都不适用于动态环境下的单频接收机。由于电离层残差法、伪距载波相位组合法需要对多个频点的载波进行组合,而单频接收机只能提供一个频点的载波观测值,因此不能使用电离层残差法、伪距载波相位组合法。多项式拟合法、高次差法可以用于单频接收机,但是需要接收机处于静止或低速运动的状态,因此也不适用于处于快速运动状态下的单频接收机。针对这一问题,本发明提出一种适合于快速运动状态下单频接收机的周跳检测与修复方法,克服现有技术的不足。
发明内容
本发明要解决的技术问题是:克服现有技术的不足,提出一种动态环境下单频周跳检测与修复方法,该方法,克服了现有周跳检测修复方法的缺点,电离层残差和伪距组合法利用双频且不能同时检测大周跳和小周跳,利用定位残差的方法计算量大需多次重复最小二乘计算,而本发明提供的方法既能同时检测和修复大周跳和小周跳,又可以仅利用单频数据,且不依赖于接收机的运动状态,在定位的同时即可进行周跳检测与修复,计算量小且速度快。
本发明采用了如下所述的技术方案:
一种动态环境下单频周跳检测与修复方法,包括如下步骤:
步骤1、建立接收机的初始定位矩阵为:
Y=Gα+ε
其中,Y为双差载波观测量减去双差载波整周模糊度,G为双差载波观测矩阵,α为待求解的差分位置三维坐标,ε为载波观测误差;
在进行双差计算时,选取某个卫星作为基准星,其他各个卫星为非基准卫星,其他各个卫星与基准星做双差计算,共有n个非基准卫星;
步骤2、在G矩阵中最后加入一列,数值为载波观测量对应的波长λ,将该增广矩阵记为X,相对应的,在待求解的差分位置三维坐标α最后加入一个未知参数,将该四维未知参数记为β,此时扩展定位矩阵为Y=Xβ+ε;
步骤3、对扩展定位矩阵Y=Xβ+ε进行最小二乘定位解算,然后计算残差平方和并根据得到的残差平方和计算双差载波观测量误差方差其中n是观测方程个数,n大于等于5,p为参数个数4,即X矩阵为n行、p列矩阵;
步骤5、计算扩展定位矩阵的投影矩阵H,如下:
H=X(XTX)-1XT
步骤6、计算X矩阵的每一行对定位结果的影响因子Di;
其中i=1,2,...,n;
hi为投影矩阵H的第i个对角元;
步骤7、将步骤6得到的影响因子Di与第二门限值进行比较,如果存在影响因子Di高于第二门限值,则表示第i个非基准卫星发生周跳,进入步骤8;如果所有影响因子Di均不高于第二门限值,则表示所有非基准卫星均未发生周跳,则使用步骤3中进行最小二乘定位解算时得到的作为定位结果进行输出,得到待求解的差分位置三维坐标α,进入步骤9;
步骤8、去掉Di最大的非基准卫星,进入步骤3,直至剩余的非基准卫星的个数小于5时,定位失败;
步骤9、将输出的定位结果代入到步骤8中去掉的非基准卫星对应观测方程中,得到载波观测误差ε,将载波观测误差ε除以波长λ并四舍五入取整数得到的结果作为非基准卫星周跳整数值的估计量,对非基准卫星进行修复。若没有去掉的非基准卫星,则不需要修复,进入步骤10;
步骤10、对定位结果中的最后一个元素进行t检验,若t检验结果中的Ρ值小于5%,则表示基准卫星的周跳值显著不为零,对进行四舍五入取整数后,作为基准卫星周跳整数的估计并加以修复,若t检验结果中的Ρ值不小于5%,则表示基准卫星未发生周跳,不对基准卫星进行修复。
有益效果
(1)本发明的优点在于克服了经典的周跳检测与修复方法的缺点,将动态场景下位置求解与周跳检测结合起来,而传统方法需要双频差分消除位置变化的影响,不能用于单频接收机,本发明将高精度的位置解算与周跳检测同时求解,因此只需单频数据就可以完成周跳检测与修复,从而可以在单频接收机上进行快速高精度的周跳检测与修复。同时,考虑到传统方法需要分两步先修复大周跳后再对小周跳进行修复的步骤,本发明可以直接对任意大小的周跳一次性检测与修复,提高了效率。
(2)在运算量的考虑上,其他周跳检测方法需在最小二乘定位之后,逐个剔除卫星后,对每一种卫星排列组合都需再次进行最小二乘定位解算,将所有排列组合对应的解算结果进行比较并判断是否有卫星发生了周跳,本发明所提供的方法不需要逐个组合进行定位解算,而是利用了定位影响因子D统计量,该统计量可衡量剔除单个观测量后对原方程的影响大小,比单纯的判断残差更能解释单个卫星对位置解算的影响大小,且由于该统计量利用了矩阵和统计学的知识进行推导和化简,从而避免了逐个排列组合进行矩阵求解的运算,可以直接计算出结果,大大减少了运算量,使得检测准确率高的同时,运算速度大大加快。
(3)本发明并不局限于单频动态场景下周跳检测与修复,也可以稍加修改用于静态、多频、非差环境下的周跳检测与修复。由于静态环境比动态环境更单纯,可以直接将本方法应用于静态环境下的周跳检测与修复,无需任何修改。对于多频点下的周跳检测与修复,只需在步骤2的G矩阵后再增加多列,每一列对应不同的波长,即可用后面相同的步骤进行周跳检测与修复。对于非差环境如PPP高精度定位等,只需要将步骤1中的定位方程改为非差定位观测方程即可。
(4)本发明实施了一种动态环境下单频周跳检测与修复方法,该方法包括以下步骤:求解出整周模糊度后构建定位方程,在定位方程中增加一列作为基准卫星是否发生周跳的估计量,进行最小二乘估计,计算出各个非基准卫星对应的最小二乘残差,根据残差平方和计算出载波误差方差计算出定位矩阵对应的投影矩阵的值,用前面三个步骤的结果计算D统计量,以此来判断各个非基准卫星是否发生周跳,将计算出的D统计量与预先设定的门限值进行比较,根据比较结果将发生周跳的卫星剔除后重新进行定位,并判断本次定位是否能够满足性能需求,成功定位之后将定位结果三维坐标代入到被剔除的卫星的观测方程中,计算周跳值并修复。对增加的一列的最小二乘回归系数做t检验,判断基准卫星是否发生周跳并修复。本发明有一定的数理统计理论基础,与其他周跳检测修复方法需要双频点、不能同时修复大周跳和小周跳、需逐个剔除卫星后再次计算相比,该方法可以直接计算出各个卫星周跳数值大小及其对定位的影响,计算量大大减少,且通过试验表明该算法可以正确的检测和修复卫星周跳,是一种便捷快速有效的动态环境下单频周跳检测与修复方法。
附图说明
图1为本发明动态环境下单频周跳检测与修复方法的流程图。
具体实施方式
一种动态环境下单频周跳检测与修复方法,包括以下步骤:
步骤1、建立接收机的初始定位矩阵为:
Y=Gα+ε
其中,Y为双差载波观测量减去双差载波整周模糊度,G为双差载波观测矩阵,α为待求解的差分位置三维坐标,ε为载波观测误差;
步骤2、在G矩阵中最后加入一列,数值为载波观测量对应的波长λ,将该增广矩阵记为X,相对应的,在待求解的差分位置三维坐标α最后加入一个未知参数,将该四维未知参数记为β,此时扩展定位矩阵为Y=Xβ+ε;
步骤3、对扩展定位矩阵Y=Xβ+ε进行最小二乘定位解算,然后计算残差平方和并根据得到的残差平方和计算双差载波观测量误差方差其中n是观测方程个数,n大于等于5,p为参数个数4,即X矩阵为n行、p列矩阵;
步骤5、计算扩展定位矩阵的投影矩阵H,如下:
H=X(XTX)-1XT
步骤6、计算X矩阵的每一行对定位结果的影响因子Di;
步骤7、将步骤6得到的影响因子Di与第二门限值进行比较,如果存在影响因子Di高于第二门限值,则表示第i个非基准卫星发生周跳,进入步骤8;如果所有影响因子Di均不高于第二门限值,则表示所有非基准卫星均未发生周跳,则使用步骤3中进行最小二乘定位解算时得到的作为定位结果进行输出,得到待求解的差分位置三维坐标α,进入步骤9;
步骤8、去掉Di最大的非基准卫星,进入步骤3,直至剩余的非基准卫星的个数小于5时,定位失败;
步骤9、将输出的定位结果代入到步骤8中去掉的非基准卫星对应观测方程中,得到载波观测误差ε,将载波观测误差ε除以波长λ并四舍五入取整数得到的结果作为非基准卫星周跳整数值的估计量,对非基准卫星进行修复。若没有去掉的非基准卫星,则不需要修复,进入步骤10;
步骤10、对定位结果中的最后一个元素进行t检验,若t检验结果中的Ρ值小于5%,则表示基准卫星的周跳值显著不为零,对进行四舍五入取整数后,作为基准卫星周跳整数的估计并加以修复,若t检验结果中的Ρ值不小于5%,则表示基准卫星未发生周跳,不对基准卫星进行修复;
其周跳检测与定位计算同时进行,所用的矩阵为双差几何矩阵G的增广矩阵X,在G矩阵中最后加入一列,数值为载波观测量对应的波长λ,将该增广矩阵记为X。与此对应的,待求解的位置三维向量最后增加一个整数待估参数,为基准卫星的周跳值。
投影矩阵H=X(XTX)-1XT,其中X是最小二乘定位时所用到的几何矩阵。
第i个卫星的D统计量Di的计算方法为:
D统计量的门限值为多次重复观测试验获得多个D统计量的95%分位点。
可以略加修改的应用于单频或双频、差分或非差、静态或动态环境下,对于大周跳、小周跳仅需一次就可检测并完成修复。
下面将就结合图1对本发明动态环境下单频周跳检测与修复方法作进一步的解释说明。
如图1所示,步骤1、求解出双差整周模糊度后,生成双差载波观测方程组,定位矩阵形式为Y=Gα+ε,其中Y为双差载波观测量减去已经求解出的双差载波整周模糊度,G为双差载波观测矩阵,α为待求解的差分位置三维坐标ε为载波观测误差向量。
在步骤1中,计算双差时,选取某个卫星作为基准星,其他各个卫星与基准星做双差计算。这里全部卫星为全部满足一定条件的可见星,即载噪比、仰角等大于预先设置好的门限值,且整周模糊度已经成功求解出的卫星。经过双差之后,观测方程数n比卫星数少一个。
步骤2、在G矩阵中最后加入一列,数值为载波观测量对应的波长λ,将该增广矩阵记为X,相对应的,在待求解的三维参数α最后加入一个未知参数,将该四维未知参数记为β。此时定位矩阵形式为Y=Xβ+ε。
步骤5、根据最小二乘定位矩阵计算对应的投影矩阵。
投影矩阵H=X(XTX)-1XT,其中X是最小二乘定位时所用到的几何矩阵。
在步骤5中,如果存在非基准卫星的周跳,那么根据步骤2中最小二乘定位最后一次使用的定位几何矩阵X计算出H矩阵的值。
步骤6、计算非基准卫星即X矩阵的每一行对定位结果的影响因子D。
步骤7、将步骤6得到的影响因子Di与第二门限值进行比较,如果存在影响因子Di高于第二门限值,则表示第i个非基准卫星发生周跳,进入步骤8;如果所有影响因子Di均不高于第二门限值,则表示所有非基准卫星均未发生周跳,则使用步骤3中进行最小二乘定位解算时得到的作为定位结果进行输出,得到待求解的差分位置三维坐标α,进入步骤9;
步骤7中计算出的D统计量为向量,每一颗非基准卫星对应一个Di值,将每一个卫星的Di值与门限值对比,作为判断卫星是否发生周跳的依据。D统计量的门限值可以根据多次重复观测试验获得多个D统计量的95%分位点。
步骤7与其他基于定位残差的方法的主要区别就是其他基于定位残差的方法需要将全部非基准卫星分成很多组,每组含有除了某一个或几个卫星之外的卫星,卫星的选取一般基于定位误差的大小和排序,然后对每组分别进行最小二乘解算,解算后计算相关的统计量,然后将这些统计量与门限值对比,作为是否发生周跳的依据,而本实施方案则利用D统计量,D统计量可以直接衡量剔除单个卫星后对定位结果的影响,由于D统计量其实本质上是经过矩阵运算推导出的结果,所以无需逐个卫星组合再次定位,简化了剔除卫星后再定位的运算,可以直接得到单个卫星对定位结果的影响大小,按影响大小依次剔除,没有冗余重复的运算,直到满足步骤4中的条件,因此在提高准确率的同时,大大减少了运算量。
步骤8、去掉Di最大的非基准卫星,进入步骤3,直至剩余的非基准卫星的个数小于5时,定位失败;
步骤9、将输出的定位结果代入到步骤8中去掉的非基准卫星对应观测方程中,得到载波观测误差ε,将载波观测误差ε除以波长λ并四舍五入取整数得到的结果作为非基准卫星周跳整数值的估计量,对非基准卫星进行修复。若没有去掉的非基准卫星,则不需要修复,进入步骤10;
每个定位周期执行卫星周跳检测与修复,下一定位周期重新进入步骤1时,重新使用全部卫星参与定位。
步骤10、判断定位结果中的最后一个元素为添加的最后一列λ波长的对应系数,是对基准卫星是否发生周跳的估计量,对进行t检验,若Ρ值小于5%,那么表示基准卫星的周跳估计显著不为零,对进行四舍五入取整数后,作为基准卫星周跳整数的估计并加以修复。
通过实验表明本发明针对动态环境下单频周跳检测与修复方法可以快速准确的识别出基准卫星与非基准卫星是否发生了周跳,是用户端接收机的一种简单高效快捷的高精度定位方法,为利用GPS、北斗及多模单频、多模多频联合高精度实时定位导航提供了可靠性保证。
实施例1
一种动态环境下单频周跳检测与修复方法,有十颗非基准卫星,包括如下步骤:
步骤1、建立接收机的初始定位矩阵为
Y=Gα+ε
有10颗卫星,卫星号从1到10,选取第10颗卫星作为基准星,经过差分后得到定位矩阵如下:
步骤2、在G矩阵中最后加入一列,数值为载波观测量对应的波长λ=0.19,将该增广矩阵记为X,相对应的,在待求解的差分位置三维坐标α最后加入一个未知参数,将该四维未知参数记为β,此时扩展定位矩阵为Y=Xβ+ε;
步骤5、计算扩展定位矩阵的投影矩阵H,如下:
H=X(XTX)-1XT
得到H的对角元为
(0.15,0.22,0.75,0.06,0.15,0.21,0.42,0.43,0.58);
步骤6、计算X矩阵的每一行对定位结果的影响因子Di;
其中i=1,2,...,9;
hi为投影矩阵H的第i个对角元;
步骤7、将步骤6得到的影响因子Di与第二门限值进行比较,第二门限值根据多次实验选取经验值1.0,可以发现D3>1.0,表示第3个非基准卫星发生周跳,进入步骤8;
步骤9、将输出的定位结果代入到步骤8中去掉的非基准卫星对应观测方程中,得到载波观测误差将载波观测误差ε除以波长λ=0.19,并四舍五入取整数得到的结果为3,非基准卫星3号星的周跳整数值为3,对非基准卫星3号星进行整周修复;
实施例2
一种动态环境下单频周跳检测与修复方法,有十颗非基准卫星,包括如下步骤:
步骤1、建立接收机的初始定位矩阵为
Y=Gα+ε
有10颗卫星,卫星号从1到10,选取第10颗卫星作为基准星,经过差分后得到定位矩阵如下:
步骤2、在G矩阵中最后加入一列,数值为载波观测量对应的波长λ=0.19,将该增广矩阵记为X,相对应的,在待求解的差分位置三维坐标α最后加入一个未知参数,将该四维未知参数记为β,此时扩展定位矩阵为Y=Xβ+ε;
步骤10、对定位结果中的最后一个元素进行t检验,定位结果的最小二乘估计为其协方差矩阵为对进行t检验,查标准t分布表得到P值为0.00000849,小于5%,表示基准卫星发生周跳,值为对基准卫星进行修复。
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (10)
1.一种动态环境下单频周跳检测与修复方法,其特征在于包括如下步骤:
步骤1、建立接收机的初始定位矩阵为:
Y=Gα+ε
其中,Y为双差载波观测量减去双差载波整周模糊度,G为双差载波观测矩阵,α为待求解的差分位置三维坐标,ε为载波观测误差;
步骤2、在G矩阵中最后加入一列,生成增广矩阵记为X,在待求解的差分位置三维坐标α最后加入一个未知参数,生成四维未知参数记为β,生成扩展定位矩阵为Y=Xβ+ε;
步骤5、计算扩展定位矩阵的投影矩阵H,如下:
H=X(XTX)-1XT
步骤6、计算X矩阵的每一行对定位结果的影响因子Di;
其中i=1,2,...,n;
hi为投影矩阵H的第i个对角元;
步骤7、将步骤6得到的影响因子Di与第二门限值进行比较,如果存在影响因子Di高于第二门限值,则表示第i个非基准卫星发生周跳,进入步骤8;如果所有影响因子Di均不高于第二门限值,则表示所有非基准卫星均未发生周跳,则使用步骤3中进行最小二乘定位解算时得到的作为定位结果进行输出,得到待求解的差分位置三维坐标α,进入步骤9;
步骤8、去掉Di最大的非基准卫星,进入步骤3,直至剩余的非基准卫星的个数小于5时,定位失败;
步骤9、将输出的定位结果代入到步骤8中去掉的非基准卫星对应观测方程中,得到载波观测误差ε,将载波观测误差ε除以波长λ并四舍五入取整数得到的结果作为非基准卫星周跳整数值的估计量,对非基准卫星进行修复;若没有去掉的非基准卫星,则不需要修复,进入步骤10;
2.根据权利要求1所述的一种动态环境下单频周跳检测与修复方法,其特征在于:所述的步骤1中,在进行双差计算时,选取某个卫星作为基准星,其他各个卫星为非基准卫星,其他各个卫星与基准星做双差计算。
3.根据权利要求2所述的一种动态环境下单频周跳检测与修复方法,其特征在于:共有n个非基准卫星。
4.根据权利要求1所述的一种动态环境下单频周跳检测与修复方法,其特征在于:所述的步骤2中,在G矩阵中最后加入的一列的数值为载波观测量对应的波长λ。
5.根据权利要求1所述的一种动态环境下单频周跳检测与修复方法,其特征在于:所述的步骤3中,n大于等于5。
6.根据权利要求1所述的一种动态环境下单频周跳检测与修复方法,其特征在于:所述的步骤3中,p=4。
7.根据权利要求1所述的一种动态环境下单频周跳检测与修复方法,其特征在于:所述的步骤3中,X矩阵为n行、p列矩阵。
8.根据权利要求1所述的一种动态环境下单频周跳检测与修复方法,其特征在于:所述的步骤4中,第一门限值为一设定值。
9.根据权利要求1所述的一种动态环境下单频周跳检测与修复方法,其特征在于:所述的步骤7中,第二门限值为一设定值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911248588.1A CN111190200B (zh) | 2019-12-09 | 2019-12-09 | 一种动态环境下单频周跳检测与修复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911248588.1A CN111190200B (zh) | 2019-12-09 | 2019-12-09 | 一种动态环境下单频周跳检测与修复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111190200A true CN111190200A (zh) | 2020-05-22 |
CN111190200B CN111190200B (zh) | 2022-05-13 |
Family
ID=70705764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911248588.1A Active CN111190200B (zh) | 2019-12-09 | 2019-12-09 | 一种动态环境下单频周跳检测与修复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111190200B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112558110A (zh) * | 2020-11-26 | 2021-03-26 | 广东星舆科技有限公司 | 一种基准站数据周跳检测方法、计算机介质及装置 |
CN113655504A (zh) * | 2021-07-19 | 2021-11-16 | 中铁第四勘察设计院集团有限公司 | 一种卫星数据处理方法、装置、电子设备及存储介质 |
Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000065593A (ja) * | 1998-08-26 | 2000-03-03 | Japan Aviation Electronics Industry Ltd | ハイブリッド形航法システム |
CN101403790A (zh) * | 2008-11-13 | 2009-04-08 | 浙江师范大学 | 单频gps接收机的精密单点定位方法 |
US20100141510A1 (en) * | 2008-12-09 | 2010-06-10 | Dai Liwen L | Methods and systems to increase accuracy in the navigation of single frequency receivers |
CN102116867A (zh) * | 2009-12-30 | 2011-07-06 | 中国科学院微电子研究所 | 一种在动态环境下探测并修复gps载波相位周跳的方法 |
CN103529462A (zh) * | 2013-10-21 | 2014-01-22 | 西南交通大学 | 一种用于全球导航卫星系统的动态周跳探测与修复方法 |
CN103698785A (zh) * | 2013-12-13 | 2014-04-02 | 合肥工业大学 | 基于贯序极限学习机的卫星信号周跳探测与修复方法 |
JP5553909B2 (ja) * | 2010-12-13 | 2014-07-23 | 古野電気株式会社 | サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム |
CN104570011A (zh) * | 2014-12-26 | 2015-04-29 | 中国电子科技集团公司第二十研究所 | 一种卫星导航相对定位装置及其载波相位周跳修复方法 |
US20150293233A1 (en) * | 2012-10-25 | 2015-10-15 | Fugron N.V. | Ppp-rtk method and system for gnss signal based position determination |
CN105137459A (zh) * | 2015-07-29 | 2015-12-09 | 昆明理工大学 | 一种北斗单频周跳探测方法 |
CN105467405A (zh) * | 2015-12-04 | 2016-04-06 | 航天恒星科技有限公司 | 星载gnss接收机的周跳探测与修复方法及装置 |
CN105607077A (zh) * | 2015-12-23 | 2016-05-25 | 北京时代民芯科技有限公司 | 一种星载双模四频gnss导航接收机 |
CN106168672A (zh) * | 2016-01-01 | 2016-11-30 | 广州中海达卫星导航技术股份有限公司 | 一种gnss多模单频rtk周跳探测方法及装置 |
KR101705882B1 (ko) * | 2015-10-05 | 2017-02-10 | 한국해양과학기술원 | 기준국 환경에서 수신기의 위치 정보를 이용한 사이클 슬립 검출 장치 및 방법 |
CN107462909A (zh) * | 2017-06-19 | 2017-12-12 | 昆明理工大学 | 一种北斗单频载波相位的周跳探测方法 |
CN107505642A (zh) * | 2017-10-23 | 2017-12-22 | 中国矿业大学 | 一种ins辅助的实时bds单频周跳探测方法 |
CN107783154A (zh) * | 2017-09-22 | 2018-03-09 | 北京时代民芯科技有限公司 | 一种接收机自主完好性故障检测及排除方法 |
CN108169774A (zh) * | 2017-12-26 | 2018-06-15 | 北方信息控制研究院集团有限公司 | 支持rtppp和rtk的多模gnss单频周跳探测与修复方法 |
CN108196281A (zh) * | 2017-11-22 | 2018-06-22 | 同济大学 | 一种基于位置域曲线约束的单频动态周跳探测与修复方法 |
CN108768393A (zh) * | 2017-12-19 | 2018-11-06 | 北京时代民芯科技有限公司 | 一种用于pll频率综合器的周跳抑制电路 |
CN109407127A (zh) * | 2018-09-20 | 2019-03-01 | 昆明理工大学 | 一种北斗卫星导航系统载波相位周跳探测与修复方法 |
KR101979184B1 (ko) * | 2018-12-28 | 2019-05-16 | 세종대학교산학협력단 | 전리층 변화율 추정을 이용한 다중 주파수 위성항법 반송파 측정치의 사이클 슬립 검출장치 및 그 방법 |
-
2019
- 2019-12-09 CN CN201911248588.1A patent/CN111190200B/zh active Active
Patent Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000065593A (ja) * | 1998-08-26 | 2000-03-03 | Japan Aviation Electronics Industry Ltd | ハイブリッド形航法システム |
CN101403790A (zh) * | 2008-11-13 | 2009-04-08 | 浙江师范大学 | 单频gps接收机的精密单点定位方法 |
US20100141510A1 (en) * | 2008-12-09 | 2010-06-10 | Dai Liwen L | Methods and systems to increase accuracy in the navigation of single frequency receivers |
CN102116867A (zh) * | 2009-12-30 | 2011-07-06 | 中国科学院微电子研究所 | 一种在动态环境下探测并修复gps载波相位周跳的方法 |
JP5553909B2 (ja) * | 2010-12-13 | 2014-07-23 | 古野電気株式会社 | サイクルスリップ検出装置、サイクルスリップ検出のための集積回路、同集積回路を備えるgnss受信装置、同gnss受信装置を備える情報端末機器、サイクルスリップ検出方法及びプログラム |
US20150293233A1 (en) * | 2012-10-25 | 2015-10-15 | Fugron N.V. | Ppp-rtk method and system for gnss signal based position determination |
CN103529462A (zh) * | 2013-10-21 | 2014-01-22 | 西南交通大学 | 一种用于全球导航卫星系统的动态周跳探测与修复方法 |
CN103698785A (zh) * | 2013-12-13 | 2014-04-02 | 合肥工业大学 | 基于贯序极限学习机的卫星信号周跳探测与修复方法 |
CN104570011A (zh) * | 2014-12-26 | 2015-04-29 | 中国电子科技集团公司第二十研究所 | 一种卫星导航相对定位装置及其载波相位周跳修复方法 |
CN105137459A (zh) * | 2015-07-29 | 2015-12-09 | 昆明理工大学 | 一种北斗单频周跳探测方法 |
KR101705882B1 (ko) * | 2015-10-05 | 2017-02-10 | 한국해양과학기술원 | 기준국 환경에서 수신기의 위치 정보를 이용한 사이클 슬립 검출 장치 및 방법 |
CN105467405A (zh) * | 2015-12-04 | 2016-04-06 | 航天恒星科技有限公司 | 星载gnss接收机的周跳探测与修复方法及装置 |
CN105607077A (zh) * | 2015-12-23 | 2016-05-25 | 北京时代民芯科技有限公司 | 一种星载双模四频gnss导航接收机 |
CN106168672A (zh) * | 2016-01-01 | 2016-11-30 | 广州中海达卫星导航技术股份有限公司 | 一种gnss多模单频rtk周跳探测方法及装置 |
CN107462909A (zh) * | 2017-06-19 | 2017-12-12 | 昆明理工大学 | 一种北斗单频载波相位的周跳探测方法 |
CN107783154A (zh) * | 2017-09-22 | 2018-03-09 | 北京时代民芯科技有限公司 | 一种接收机自主完好性故障检测及排除方法 |
CN107505642A (zh) * | 2017-10-23 | 2017-12-22 | 中国矿业大学 | 一种ins辅助的实时bds单频周跳探测方法 |
CN108196281A (zh) * | 2017-11-22 | 2018-06-22 | 同济大学 | 一种基于位置域曲线约束的单频动态周跳探测与修复方法 |
CN108768393A (zh) * | 2017-12-19 | 2018-11-06 | 北京时代民芯科技有限公司 | 一种用于pll频率综合器的周跳抑制电路 |
CN108169774A (zh) * | 2017-12-26 | 2018-06-15 | 北方信息控制研究院集团有限公司 | 支持rtppp和rtk的多模gnss单频周跳探测与修复方法 |
CN109407127A (zh) * | 2018-09-20 | 2019-03-01 | 昆明理工大学 | 一种北斗卫星导航系统载波相位周跳探测与修复方法 |
KR101979184B1 (ko) * | 2018-12-28 | 2019-05-16 | 세종대학교산학협력단 | 전리층 변화율 추정을 이용한 다중 주파수 위성항법 반송파 측정치의 사이클 슬립 검출장치 및 그 방법 |
Non-Patent Citations (4)
Title |
---|
F. ZANGENEH-NEJAD: "Cycle slip detection and repair of undifferenced single-frequency GPS carrier phase observations", 《GPS SOLUT》 * |
何海波等: "长基线双频GPS动态测量中的周跳修复算法", 《测绘科学技术学报》 * |
王凌轩: "INS辅助的实时动态GNSS单频周跳探测", 《武汉大学学报.信息科学版》 * |
由立华: "基于载噪比的BD2/GPS多模加权最小二乘定位算法设计", 《中国卫星导航系统管理办公室学术交流中心.第六届中国卫星导航学术年会论文集—S07北斗/GNSS用户终端技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112558110A (zh) * | 2020-11-26 | 2021-03-26 | 广东星舆科技有限公司 | 一种基准站数据周跳检测方法、计算机介质及装置 |
CN113655504A (zh) * | 2021-07-19 | 2021-11-16 | 中铁第四勘察设计院集团有限公司 | 一种卫星数据处理方法、装置、电子设备及存储介质 |
CN113655504B (zh) * | 2021-07-19 | 2023-11-07 | 中铁第四勘察设计院集团有限公司 | 一种卫星数据处理方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN111190200B (zh) | 2022-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103529462B (zh) | 一种用于全球导航卫星系统的动态周跳探测与修复方法 | |
CN110531390B (zh) | 多普勒积分法联合电离层残差法检测与修复周跳 | |
CN106168672B (zh) | 一种gnss多模单频rtk周跳探测方法及装置 | |
CN103076618B (zh) | 一种星载测量型gnss接收机地面验证和性能评测方法 | |
CN111190200B (zh) | 一种动态环境下单频周跳检测与修复方法 | |
CN102288978A (zh) | 一种cors基站周跳探测与修复方法 | |
CN106066485B (zh) | 一种周跳探测方法及装置 | |
CN107607966A (zh) | 一种基于gnss三频载波相位的周跳检测与修复方法 | |
CN107219537A (zh) | 一种融合选星与完好性检测的多系统兼容定位方法 | |
CN106842236B (zh) | Gnss接收机周跳探测与修复处理装置 | |
CN102998681A (zh) | 一种卫星导航系统的高频钟差估计方法 | |
CN111965668A (zh) | 一种面向卫星多故障的raim方法 | |
CN106932793A (zh) | 一种北斗三频信号的实时周跳探测与修复方法 | |
CN106093991A (zh) | 一种用于gnss定位的模糊度快速恢复方法及系统 | |
WO2019024895A1 (zh) | Gnss卫星钟差实时估计质量控制方法 | |
CN113126129B (zh) | 基于空间信号质量加权估计的gbas完好性监测方法 | |
CN110032709A (zh) | 一种用于地理坐标转换中异常点的定位与估值方法 | |
CN103886600B (zh) | 一种连续相位板加工误差识别方法 | |
CN105549046B (zh) | Gnss接收机周跳探测与修复处理方法 | |
CN107783154A (zh) | 一种接收机自主完好性故障检测及排除方法 | |
CN105652298B (zh) | 一种bds三频伪距相位组合的周跳探测修复方法及装置 | |
CN109581430A (zh) | 一种基于伪卫星的gbas电离层空间梯度的监测方法 | |
CN110727000A (zh) | 一种基于gnss高采样率数据的小周跳修复方法 | |
CN112083463A (zh) | 模糊度是否固定正确的检测方法及装置、定位终端 | |
CN104199054A (zh) | 一种用于北斗卫星导航系统共视数据的预处理方法 |
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 |