黄德武 王建英

1 昆明理工大学城市学院，昆明市环城东路50号，650051 2 云南经济管理学院工程学院，云南省安宁市麒麟路17号，650106

周跳是影响GPS监测数据可靠性的因素之一，高精度的监测成果必须保证载波相位观测值中无周跳存在[1-2]。周跳探测方法主要有相位多次求差法[3]、电离层残差法[4-6]、非差法[7-8]、双频码相组合法[9]、拟合法[10-11]和TurboEdit方法[12-14]等。上述方法在实际应用中均存在一定局限性，如相位多次求差法对小周跳探测不敏感；双频码相组合法与电离层残差法不能确定周跳发生的频率位置及多值问题，并且对特殊组合周跳失效；拟合法对数据个数及拟合阶数要求较高，同时也会产生拟合误差和舍入误差；TurboEdit方法中MW组合观测值易受伪距观测值噪声的影响，部分周跳容易被噪声湮没。针对上述问题，本文组建相位方程ΔΦ=ΔN1-(λ2/λ1)ΔN2，考虑到周跳在ΔΦ、ΔN1、ΔN2序列上表现为异常值的特点，利用Grubbs准则探测ΔN1、ΔN2异常值，通过约束ΔN1、ΔN2的范围来探测不同组合周跳。

1 双频周跳探测模型

在双频条件下，设L1和L2上周跳大小分别为ΔN1、ΔN2，根据相位和伪距观测值可得:

式中，ΔΦ为电离层残差在历元间的变化值；λ为波长；Δion(t)为某时刻的电离层残差；Δε为总误差项，主要受电离层和多路径影响。

若采样间隔较小，电离层在历元间的变化Δion(t)很小，则式(1)可简化为：

(2)

从式(2)可以看出，针对2个未知数ΔN1、ΔN2，1个方程无法解算出周跳值，因此可通过计算附有限制条件ΔΦ=ΔN1-(λ2/λ1)ΔN2来确定ΔN1、ΔN2的取值区间。

2 双频周跳特征分析

3 Grubbs检验法

设存在观测序列Y=(Y1,Y2,…,Yn)，构建Grubbs检验统计量：