零基线双差
伪距观测值可以表达为:
P = ρ + d t r − d t s + T + I + M P + O + ε P (1) P = \rho+dtr-dts+T+I+M_P+O+\varepsilon_P \tag{1} P=ρ+dtr−dts+T+I+MP+O+εP(1)
相位观测值可以表达为:
L
=
(
ρ
+
d
t
r
−
d
t
s
+
T
−
I
+
M
L
+
O
)
/
λ
+
N
+
ε
L
(2)
L = (\rho+dtr-dts+T-I+M_L+O)/\lambda + N +\varepsilon_L \tag{2}
L=(ρ+dtr−dts+T−I+ML+O)/λ+N+εL(2)
其中
ρ
\rho
ρ为地伪距,
d
t
r
dtr
dtr为接收机钟差,
d
t
s
dts
dts为卫星钟差,
T
T
T为对流层,
I
I
I为电离层,
M
M
M为多路径误差,
O
O
O为轨道误差,
ε
\varepsilon
ε为观测噪声。
零基线双差观测值
△
▽
P
r
,
b
i
,
j
=
△
▽
ρ
r
,
b
i
,
j
+
△
▽
d
t
r
r
,
b
i
,
j
−
△
▽
d
t
s
r
,
b
i
,
j
+
△
▽
T
r
,
b
i
,
j
+
△
▽
I
r
,
b
i
,
j
+
△
▽
M
P
r
,
b
i
,
j
+
△
▽
O
r
,
b
i
,
j
+
△
▽
ε
P
r
,
b
i
,
j
(3)
\triangle\triangledown P_{r,b}^{i,j} = \triangle\triangledown\rho_{r,b}^{i,j}+\triangle\triangledown {dtr}_{r,b}^{i,j}-\triangle\triangledown {dts}_{r,b}^{i,j}+ \triangle\triangledown T_{r,b}^{i,j}+\triangle\triangledown I_{r,b}^{i,j}+\triangle\triangledown {M_P}_{r,b}^{i,j}+\triangle\triangledown O_{r,b}^{i,j}+ \triangle\triangledown {\varepsilon_P}_{r,b}^{i,j}\tag{3}
△▽Pr,bi,j=△▽ρr,bi,j+△▽dtrr,bi,j−△▽dtsr,bi,j+△▽Tr,bi,j+△▽Ir,bi,j+△▽MPr,bi,j+△▽Or,bi,j+△▽εPr,bi,j(3)
r
,
b
r,b
r,b分别代表不同的两个接收机,
i
,
j
i,j
i,j分别代表双差的两颗卫星。
由于信号用一个天线公分出来,所以信号的传播路径完全一致,双差观测值中的大气误差可以认为被完全消除;两个接收机在非常相近的时刻内提取观测值(几ms内),卫星的轨道误差、卫星钟误差、接收机钟差可认为变化很小,可以被双差消掉;故最后双差后伪距观测量的数学表达可表示为:
△
▽
P
r
,
b
i
,
j
=
△
▽
ρ
r
,
b
i
,
j
+
△
▽
ε
P
r
,
b
i
,
j
(4)
\triangle\triangledown P_{r,b}^{i,j} = \triangle\triangledown\rho_{r,b}^{i,j}+\triangle\triangledown {\varepsilon_P}_{r,b}^{i,j}\tag{4}
△▽Pr,bi,j=△▽ρr,bi,j+△▽εPr,bi,j(4)
同理,零基线双差相位表达式为:
△
▽
L
r
,
b
i
,
j
=
△
▽
ρ
r
,
b
i
,
j
/
λ
+
N
r
,
b
i
,
j
+
△
▽
ε
L
r
,
b
i
,
j
(5)
\triangle\triangledown L_{r,b}^{i,j} = \triangle\triangledown\rho_{r,b}^{i,j}/\lambda+N_{r,b}^{i,j}+ \triangle\triangledown {\varepsilon_L}_{r,b}^{i,j}\tag{5}
△▽Lr,bi,j=△▽ρr,bi,j/λ+Nr,bi,j+△▽εLr,bi,j(5)
零基线双差评估接收机观测值质量
从(4)、(5)可以看出,双差后的观测值中还含有双差地伪距以及相位模糊度,剩下的就是双差后的接收机噪声了,只要消除了前面两项影响,那么就可以评估接收机观测值的噪声水平。
归化提取观测值时刻
△ ▽ ρ r , b i , j \triangle\triangledown\rho_{r,b}^{i,j} △▽ρr,bi,j的产生是由于两个接收机的实际提取观测值时刻不同,卫星的高速运动以及地球自转而产生的地伪距变化引起的,归根到底就是接收机的提取时刻没有严格一致,那么我们只需要把观测时刻归化到统一的提取时刻就可以消除该项的影响。
要归化统一的时刻,就必须知道接收机实际的提取时刻,因此,必须要估算出两个接收机的接收机钟差。
首先,对两个接收机分别进行单点定位,估算出各自的接收机钟差,求出接收机的实际提取时刻。然后通过多普勒观测值将观测值归算到同一时刻
设两接收机的实际提取时刻分别为 T r T_r Tr、 T b T_b Tb
d
t
=
T
r
−
T
b
(6)
dt = T_r - T_b\tag{6}
dt=Tr−Tb(6)
△
▽
ρ
r
,
b
i
,
j
=
d
t
∗
△
D
b
i
,
j
(7)
\triangle\triangledown\rho_{r,b}^{i,j} = dt * \triangle D_{b}^{i,j}\tag{7}
△▽ρr,bi,j=dt∗△Dbi,j(7)
其中
D
D
D为多普勒观测值。
对于伪距,扣除地伪距的影响后,即为伪距观测值的双差噪声;
对于相位观测值,还存在整周模糊度,这时候减去取整部分即可得到相位的观测噪声;