武汉大学学报 信息科学版　　　　　　　　　　　　　　　2002年236

2　计算及结果分析

2.1　数据处理方法

精密单点定位计算过程主要分如下几个步骤:①为了方便计算,将精密星历拟合成多项式形式;②精密星历每15min给定一个卫星钟差值,这个间隔不能满足精密单点定位要求,又由于卫星钟差的变化较快而不能直接以30s为间隔对其进行线性内插,必须利用多个IGS跟踪站GPS观测数据与精密星历估计得到30s历元间隔的卫星钟差;③利用拟合的轨道多项式及卫星钟差与用户站观测数据一起进行精密单点定位计算。2.2　数据预处理

在精密单点定位中,相位平滑伪距等数据预处理工作2.2.1　,GPS观测数据

中的周跳是一项重要的工作。由于非差单点定位只有单站数据能利用,无法组成双差或三差观测值,一般消除周跳的方法如三差法、多项式拟合法并不适用。而Blewitt提出的利用双频双P码组合观测值修复周跳的方法很适合清除非差周跳[5]。

用于清除非差周跳的GPS观测值线性组合有以下几种。

1)Melbourne-Wübbena组合:

()(()())

L6i

=

f

1-f2

f1L1i

-f2L2i

-

f1

(11i)+f2P2(i))(5)(6)

b6(iL6((f1-f2)/c

Φ1,L22,L1=λ1 Φ2;i为历元序号;f1、f2分别为L1、L2的2 频率;P1、P2分别为P码伪距观测值;b6为宽巷整周模糊度。

2)Geometry-free组合:

(a)　X方向

(b)　Y方向

(c)　Z方向

图1　相位平滑伪距计算的结果与已知坐标在X、Y、Z方向上的差值

Fig.1　TheDifferencesinX,Y,ZComponentsbetweentheResultsComputedwith

SmoothedCodeObservationsandITRFCoordinates