一种基于改进TurboEdit的GPS周跳探测与修复方法

一种基于改进TurboEdit的GPS周跳探测与修复方法
党子越;夏凤雨
【摘要】基于TurboEdit算法中的MW组合,利用卡尔曼滤波分离宽巷模糊度以
提高模糊度精度.利用移动窗口的多项式拟合对原GF组合中的多项式拟合进行改进,并对周跳的浮点解进行计算.使用GPS数据检测方法的正确性,结果表明,该方法
能探测并修复小周跳、大周跳以及多个周跳.
【期刊名称】《大地测量与地球动力学》
【年(卷),期】2016(036)008
【总页数】4页(P723-726)
【关键词】GPS;周跳;TurboEdit算法;移动窗口;卡尔曼滤波
【作者】党子越;夏凤雨
【作者单位】长安大学地质工程与测绘学院,西安市雁塔路126号,710064;长安大
学地质工程与测绘学院,西安市雁塔路126号,710064
【正文语种】中文
【中图分类】P228
GPS周跳探测与修复方法[1-5]存在缺陷：高阶差分法不能探测5周以下的小周跳；多项式拟合法在拟合阶数和拟合长度的选择上容易引入人为误差，并且无法探测小周跳；电离层残差法中，当L1和L2频率上发生的周跳为L1和L2的频率之比时，则无法探测出该周跳；TurboEdit法中，窄巷伪距易受多路径效应影响[6]，且其
中用到的P1码一般无法获取，只能用C/A码代替，但由于C/A码精度较低，会
对周跳检测结果产生较大影响。

本文对TurboEdit方法进行改进，在MW组合模型中，利用卡尔曼滤波法，在解算多路径误差的同时将宽巷模糊度与多路径误差进行实时分离，得到较为准确的宽巷模糊度；再采用固定弧段长度的移动窗口多项式拟合法对GF组合进行改进，从而实现双频GPS的周跳探测与修复。

MW组合整周模糊度为[7]：
式(1)通过宽巷组合消除了电离层误差、接收机钟差以及卫星钟差，仅剩下模糊度项、多路径误差以及GPS测量噪声，因此该组合在消除多路径误差及测量噪声的前提下适用于周跳的探测与修复。

根据文献[7]，利用误差传播律计算出先验均方根中误差，然后采用递推公式递推出每一个历元的宽巷模糊度的均值及其均方根中误差。

通过(i+1)历元与i历元的宽巷模糊度之差来判断i历元观测值是否存在周跳：
根据载波相位、码伪距观测方程可导出无几何距离组合LGF[7-8]：
对LGF进行多项式拟合，拟合出的结果记为QGF，拟合的阶数为
min[(N/100+1),6]，N为历元个数。

由以下不等式来判断发生周跳的历元：
若式(4)中两个不等式同时成立，则判断i历元发生周跳。

式(5)为含伪距多路径与噪声的宽巷相位减窄巷伪距组合。

可以看出，含多路径误差的码观测值对确定宽巷模糊度具有非常大的影响：
运用宽巷相位减窄巷伪距组合模型求解宽巷模糊度时，需要L1相位上的P1码。

但是大多数情况下，P1码不易获取，常常采用C/A码代替P1码。

但是C/A码精度较低，不能满足正确求解模糊度的需求，且C/A码极易受伪距多路径误差的影响。

因此，在解算模糊度浮点解时，得到的结果精度往往较低。

要求解出正确的宽巷模糊度，需要将式(5)中的伪距多路径与宽巷模糊度进行实时分离。

根据伪距与时间间隔有一定相关性的特质，对多路径采用一阶高斯-马尔科夫过程建立如下模型[9]：
不发生周跳时，整周模糊度为常数。

因此，可以建立以宽巷模糊度和伪距多路径误差为参数的状态方程：
将式(7)作为离散系统的观测方程，结合式(5)可列出如下卡尔曼方程：
式中，和Qk分别为状态噪声及其协方差阵，vk和Rk分别为观测噪声及其方差阵。

通过滤波公式便可逐历元计算出整周模糊度的浮点解及其对应的方差阵，继而推算出整周模糊度的准确值。

在用TurboEdit方法探测周跳时，GF组合模型中采用多项式拟合法对电离层组合进行拟合。

而实际情况中，拟合阶数和拟合长度易受人为影响，若两者不符，则会产生较大的拟合误差，导致周跳探测过程出现问题。

本文采用一种基于移动窗口的多项式拟合法[10]，对LGF(i)进行拟合。

取移动窗口大小为m、拟合阶数为n，则：选取当前历元前m个LGF进行多项式拟合。

根据式(6)判断当前历元是否发生周跳，若发生周跳，则先对该历元的载波相位观测值进行修复，将新的值赋予LGF(i)，窗口向前移动，直到最后一个历元。

本文实验数据采用大小为20的窗口，拟合阶数为3，拟合的情况最佳。

探测出周跳后，采用移动窗口的多项式拟合法对L1和L2相位观测值分别进行拟
合[11]：
若已探测出i历元发生周跳，则利用i历元前固定长度的窗口对i历元的L1和L2
相位观测值分别进行拟合，通过式(11)计算L1和L2相位分别发生的周跳值：
式中，mL1、mL2分别为L1和L2相位发生的周跳值，2(i)分别为i历元时L1和
L2相位分别拟合出的载波相位值，L1(i)、L2(i)分别为i历元时L1和L2相位的载
波相位观测值。

通过对周跳值的计算，可修复该历元发生的周跳，直至遍历最后一个历元。

依此类推，可修复已探测出的周跳。

通过对拟合阶数和拟合长度的实验得出，拟合阶数为4、拟合长度为16时，解算出的周跳浮点解取整误差最小。

采用2010年年积日为60 d的AMC站观测数据amc20600.10°。

该数据包含L1
相位、L2相位、C/A码、P1码、P2码，数据采样率为30 s。

选取PRN32号卫星，截取第1 h共240个历元作为观测数据。

通过图1的MW组合图形和GF组合图形可看出，该数据不含周跳，图像的不规则波动是由噪声和伪距多路径效应引起的。

通过以下3种方案检验本文方法的正确性。

方案一：在第50个历元分别给(L1，L2)相位加入(9,7)周跳；在第100个历元分别给(L1，L2)相位加入(2,2)周跳；
方案二：在第100个历元分别给(L1，L2)相位加入(-60,-60)周跳，在第200个历元分别给(L1，L2)相位加入(77,60)周跳；
方案三：在第50个历元分别给(L1，L2)相位模拟加入(2,-1)周跳，在第100个历
元分别给(L1，L2)相位加入(30,30)周跳；在第150个历元分别给(L1，L2)相位加
入(77,60)周跳。

比较图1和图2可明显看出，相比P1码数据，C/A码数据噪声更大。

从图3可
看出，经过卡尔曼滤波分离后的宽巷模糊度消除了伪距多路径效应和噪声等误差，固定出了一个较为精确的宽巷模糊度。

从图4～6可看出，MW组合对形如(9,7)和(77,60)的特殊比例周跳和(2，-1)这样的小周跳具有探测能力，却对(2,2)、(30,30)和(-60，-60)的等周周跳不敏感；GF 组合可探测出(2,2)、(30,30)和(-60，-60)的等周周跳，却对周跳(9,7)和(77,60)不敏感。

且无论是MW组合还是GF组合，都对多个周跳具有探测能力。

因此，将MW组合和GF组合联合起来探测周跳，可以得到准确的探测结果。

周跳修复时，采用移动窗口的多项式拟合法对上述方案加入的周跳值进行计算，计算结果见表1。

从表1可看出，用移动窗口的多项式拟合法计算周跳时，对于模拟加入的小周跳、大周跳、多个周跳、等值周跳和特殊比例周跳都能够计算出准确的周跳值，得到的
周跳浮点解取整后与周跳值完全相等。

因此，该方法可以用来计算可靠的周跳值。

本文在TurboEdit方法的基础上，提出一种改进的周跳探测与修复方法。

实验证明，本文提出的针对MW组合的卡尔曼滤波法能够提高计算出的整周模糊度的精度，为后续周跳探测打下坚实基础；改进的GF组合中的移动窗口多项式拟合法有效减少了人为选取拟合阶数和拟合长度所产生的误差，较好地探测出周跳发生的位置；周跳修复所采用的移动窗口的多项式拟合法能够精确地计算出周跳的大小。

本文所提出的改进方法对小周跳、大周跳以及多个周跳均能较好地探测并修复。

About the first author:DANG Ziyue,majors in surveying engineering, E-mail: **********************.。