第3期　　　　　　　　　　　　刘经南等:GPS非差相位精密单点定位技术探讨　　　　　　　　　　　　　　237

L4(i)=L1(i)-L2(i)=I+λ1n1-λ2n2

(7)

P4(i)=P1(i)-P2(i)

(8)

λ式中λ,1、L2的波长;n1、n2分别为2分别为L1、L1、L2的整周模糊度。

Melbourne-Wübbena组合消除了电离层、对

流层、钟差和计算的几何观测值的影响,而且具有较长的波长、较小的量测噪声等特点,因此适用于非差周跳的探测和修复。如果Melbourne-Wübbena的RMS小于0.5宽巷波长(43cm),利用它几乎可以确定所有的宽巷周跳。在实际计算中,采用递推的方法计算每一历元b6值及其残差误差σ:

(()<b6>

i

　　比较相邻历元b6值及其残差误差σ,可以判

断是否发生周跳。若发生周跳,则标记出发生周跳的历元,把此历元之前的数据作为一个数据弧段,并计算其b6均值及其残差误差σ,从下一个历元重新开始计算探测周跳,重复上述工作直到数据结束。弧段与弧段的周跳大小Δb6可以由两段之间的均值求得,并且Δb6与L1和L2周跳具有如下关系:

Δb6=Δn1-Δn2

(

11)

Δn2分别表示L1

和L2周跳。式中,Δn1

、

在完成利用Wübbena组合观测值确定所有宽巷周跳组合修复窄巷N个,再取周跳发生后的N,两个多项式在周跳发生历元时刻的差值可认为是窄巷周跳的大小,即可确定λ再利用1Δ

n1-λ2Δn2的大小。

=<b6>i-1+

i

b6i-<b6>i-1

22

σi=σi-1+

i

6ii-1

2

2

σi-(10)

(a)　X方向

(b)　Y方向

(c)　Z方向

图2　初始阶段非差相位观测值计算结果与已知坐标在X、Y、Z方向上的差值

Fig.2　PrecisePointPositioningSolutionConvergenceinX,Y,ZComponent