ISB 与ISPB
1.ISB
1.1 估计法(code phase)
接收机
r
r
r 和 卫星
s
s
s 在频率
j
j
j 上的载波相位观测方程可
ϕ
r
,
j
s
=
ρ
r
s
+
d
t
r
−
d
t
s
+
α
r
,
j
s
+
λ
j
(
δ
r
,
j
s
+
φ
r
,
j
−
δ
j
s
−
φ
j
s
+
N
r
,
j
s
)
+
ϵ
r
,
j
s
p
r
,
j
s
=
ρ
r
s
+
d
t
r
−
d
t
s
+
a
r
,
j
s
+
d
r
,
j
G
−
d
j
s
+
e
r
,
j
s
\begin{aligned} \phi_{r, j}^{s}=& \rho_{r}^{s}+d t_{r}-d t^{s}+\alpha_{r, j}^{s} \\ &+\lambda_{j}\left(\delta_{r, j}^{s}+\varphi_{r, j}-\delta_{j}^{s}-\varphi_{j}^{s}+N_{r, j}^{s}\right)+\epsilon_{r, j}^{s} \end{aligned}\\ p_{r, j}^{s}=\rho_{r}^{s}+d t_{r}-d t^{s}+a_{r, j}^{s}+d_{r, j}^{G}-d_{j}^{s}+e_{r, j}^{s}
ϕr,js=ρrs+dtr−dts+αr,js+λj(δr,js+φr,j−δjs−φjs+Nr,js)+ϵr,jspr,js=ρrs+dtr−dts+ar,js+dr,jG−djs+er,js
ϕ
r
,
j
s
\phi_{r, j}^{s}
ϕr,js为载波相位观测值,
ρ
r
s
\rho_{r}^{s}
ρrs 为接收机到卫星的距离;
d
t
r
、
d
t
s
d t_r、d t^{s}
dtr、dts 分别为接收机端和卫星端钟差,
α
r
,
j
s
\alpha_{r, j}^{s}
αr,js 为与频率相关的大气层延迟;
λ
j
\lambda_{j}
λj 为波长,
δ
r
,
j
s
\delta_{r, j}^{s}
δr,js 表示接收机相位硬件延迟,
φ
r
,
j
\varphi_{r, j}
φr,j 接收机初始相位,
δ
j
s
\delta_{j}^{s}
δjs 卫星端相位硬件延迟,
φ
j
s
\varphi_{j}^{s}
φjs 卫星端初始相位;
d
r
,
j
s
d_{r, j}^{s}
dr,js接收端伪距的硬件延迟,
d
j
s
d_{j}^{s}
djs卫星端伪距硬件延迟;
N
r
,
j
s
N_{r, j}^{s}
Nr,js 为整数模糊度;
ϵ
r
,
j
s
\epsilon_{r, j}^{s}
ϵr,js为包含噪声在内的其他误差。
同一卫星,流动站(下标2)与参考站(下标1)的观测值单差,得:
ϕ
12
,
j
s
=
ϕ
2
,
j
s
−
ϕ
1
,
j
s
=
ρ
12
s
+
d
t
12
+
α
12
,
j
s
+
λ
j
(
δ
12
,
j
B
2
+
φ
12
,
j
+
N
12
,
j
s
)
+
ϵ
12
,
j
s
p
12
,
j
s
=
p
2
,
j
s
−
p
1
,
j
s
=
ρ
12
s
+
d
t
12
+
a
12
,
j
s
+
d
12
,
j
B
2
+
e
12
,
j
s
\begin{aligned} \phi_{12, j}^{s}=& \phi_{2, j}^{s}-\phi_{1, j}^{s}=\rho_{12}^{s}+d t_{12}+\alpha_{12, j}^{s} \\ &+\lambda_{j}\left(\delta_{12, j}^{B2}+\varphi_{12, j}+N_{12, j}^{s}\right)+\epsilon_{12, j}^{s} \\ p_{12, j}^{s}=& p_{2, j}^{s}-p_{1, j}^{s}=\rho_{12}^{s}+d t_{12}+a_{12, j}^{s}+d_{12, j}^{B2}+e_{12, j}^{s} \end{aligned}
ϕ12,js=p12,js=ϕ2,js−ϕ1,js=ρ12s+dt12+α12,js+λj(δ12,jB2+φ12,j+N12,js)+ϵ12,jsp2,js−p1,js=ρ12s+dt12+a12,js+d12,jB2+e12,js
j
=
1
,
…
,
f
j=1, \ldots, f
j=1,…,f ,表示信号频率标号,
s
=
1
B
2
,
…
,
m
B
2
,
s=1_{B2}, \ldots, m_{B2},
s=1B2,…,mB2,表示追踪到的BDS-2的卫星标号,。 同理对BDS-3有:
ϕ
12
,
j
q
=
ϕ
2
,
j
q
−
ϕ
1
,
j
q
=
ρ
12
q
+
d
t
12
+
α
12
,
j
q
+
λ
j
(
δ
12
,
j
B
3
+
φ
12
,
j
+
N
12
,
j
q
)
+
ϵ
12
,
j
q
p
12
,
j
q
=
p
2
,
j
q
−
p
1
,
j
q
=
ρ
12
q
+
d
t
12
+
a
12
,
j
q
+
d
12
,
j
B
3
+
e
12
,
j
q
\begin{aligned} \phi_{12, j}^{q}=& \phi_{2, j}^{q}-\phi_{1, j}^{q}=\rho_{12}^{q}+d t_{12}+\alpha_{12, j}^{q} \\ &+\lambda_{j}\left(\delta_{12, j}^{B3}+\varphi_{12, j}+N_{12, j}^{q}\right)+\epsilon_{12, j}^{q} \\ p_{12, j}^{q}=& p_{2, j}^{q}-p_{1, j}^{q}=\rho_{12}^{q}+d t_{12}+a_{12, j}^{q}+d_{12, j}^{B3}+e_{12, j}^{q} \end{aligned}
ϕ12,jq=p12,jq=ϕ2,jq−ϕ1,jq=ρ12q+dt12+α12,jq+λj(δ12,jB3+φ12,j+N12,jq)+ϵ12,jqp2,jq−p1,jq=ρ12q+dt12+a12,jq+d12,jB3+e12,jq
q
=
1
B
3
,
…
,
m
B
3
q=1_{B3}, \ldots, m_{B3}
q=1B3,…,mB3表示观测到的BDS-3的卫星标号。对于零基线或者短基线而言,我们可以认为两个站对同一卫星同一频率信号的大气层延迟是一样的。如果我们选择一个参考卫星,
s
=
1
B
2
s=1_{B2}
s=1B2 ,
q
=
1
B
3
q=1_{B3}
q=1B3,则有:
ϕ
12
,
j
1
B
2
s
=
ϕ
12
,
j
s
−
ϕ
12
,
j
1
B
2
=
ρ
12
1
B
2
s
+
λ
j
N
12
,
j
1
B
2
s
+
ϵ
12
,
j
1
B
2
s
p
12
,
j
1
B
2
s
=
p
12
,
j
s
−
p
12
,
j
1
B
2
=
ρ
12
1
B
2
s
+
e
12
,
j
1
B
2
s
ϕ
12
,
j
1
B
3
q
=
ϕ
12
,
j
q
−
ϕ
12
,
j
1
B
3
=
ρ
12
1
B
3
q
+
λ
j
N
12
,
j
1
B
3
q
+
ϵ
12
,
j
1
B
3
q
p
12
,
j
1
B
3
q
=
p
12
,
j
q
−
p
12
,
j
1
B
3
=
ρ
12
1
B
3
q
+
e
12
,
j
1
B
3
q
\begin{array}{l} \phi_{12, j}^{1_{B2} s}=\phi_{12, j}^{s}-\phi_{12, j}^{1_{B2}}=\rho_{12}^{1_{B2} s}+\lambda_{j} N_{12, j}^{1_{B2} s}+\epsilon_{12, j}^{1_{B2} s} \\ p_{12, j}^{1_{B2} s}=p_{12, j}^{s}-p_{12, j}^{1_{B2}}=\rho_{12}^{1_{B2} s}+e_{12, j}^{1_{B2} s} \end{array}\\ \begin{aligned} \phi_{12, j}^{1_{B3} q} &=\phi_{12, j}^{q}-\phi_{12, j}^{1_{B3}}=\rho_{12}^{1_{B3} q}+\lambda_{j} N_{12, j}^{1_{B3} q}+\epsilon_{12, j}^{1_{B3} q} \\ p_{12, j}^{1_{B3} q} &=p_{12, j}^{q}-p_{12, j}^{1_{B3}}=\rho_{12}^{1_{B3} q}+e_{12, j}^{1_{B3} q} \end{aligned}
ϕ12,j1B2s=ϕ12,js−ϕ12,j1B2=ρ121B2s+λjN12,j1B2s+ϵ12,j1B2sp12,j1B2s=p12,js−p12,j1B2=ρ121B2s+e12,j1B2sϕ12,j1B3qp12,j1B3q=ϕ12,jq−ϕ12,j1B3=ρ121B3q+λjN12,j1B3q+ϵ12,j1B3q=p12,jq−p12,j1B3=ρ121B3q+e12,j1B3q
对于BDS-2和BDS-3的重叠频率(B1I、B3I),我们可以对不同系统的单差间求差,
ϕ
12
,
j
1
B
2
q
=
ϕ
12
,
j
q
−
ϕ
12
,
j
1
B
2
=
ρ
12
1
B
2
q
+
λ
j
(
δ
12
,
j
B
2
B
3
+
N
12
,
j
1
1
q
)
+
ϵ
12
,
j
1
B
2
q
q
p
12
,
j
1
B
2
q
=
p
12
,
j
q
−
p
12
,
j
1
B
2
=
ρ
12
1
B
2
q
+
d
12
,
j
B
2
B
3
+
e
12
,
j
1
G
q
\begin{array}{l} \phi_{12, j}^{1_{B2 q}}=\phi_{12, j}^{q}-\phi_{12, j}^{1_{B2}}=\rho_{12}^{1_{B2} q}+\lambda_{j}\left(\delta_{12, j}^{B2 B3}+N_{12, j}^{1_{1} q}\right)+\epsilon_{12, j}^{1_{B2 q} q} \\ p_{12, j}^{1_{B2} q}=p_{12, j}^{q}-p_{12, j}^{1_{B2}}=\rho_{12}^{1_{B2} q}+d_{12, j}^{B2 B3}+e_{12, j}^{1_{G} q} \end{array}
ϕ12,j1B2q=ϕ12,jq−ϕ12,j1B2=ρ121B2q+λj(δ12,jB2B3+N12,j11q)+ϵ12,j1B2qqp12,j1B2q=p12,jq−p12,j1B2=ρ121B2q+d12,jB2B3+e12,j1Gq
根据系统间双差结果,
I
S
B
=
δ
12
,
j
B
2
B
3
=
δ
2
,
j
B
2
B
3
−
δ
1
,
j
B
2
B
3
,
I
S
P
B
=
d
12
,
j
B
2
B
3
=
d
2
,
j
B
2
B
3
−
d
1
,
j
B
2
B
3
ISB = \delta_{12, j}^{B2 B3}=\delta_{2, j}^{B2 B3}-\delta_{1, j}^{B2 B3} ,ISPB = d_{12, j}^{B2 B3}=d_{2, j}^{B2 B3}-d_{1, j}^{B2 B3}
ISB=δ12,jB2B3=δ2,jB2B3−δ1,jB2B3,ISPB=d12,jB2B3=d2,jB2B3−d1,jB2B3, 分别表示BDS-2和BDS-3的相位和伪距的系统间偏差(BDS-2和BDS-3特定频率信号在同一个接收机上不同的硬件延迟),
r
=
1
,
2
r=1,2
r=1,2。
Characterization of between-receiver GPS-Galileo inter-system biases and their effect on mixed ambiguity resolution Dennis Odijk • Peter J. G. Teunissen
1.2 ISB的随机模型
五系统 GNSS(GPS + GLONASS + BDS + Galileo + QZSS)非差非组合模型(引入ISB)
{
p
r
,
j
s
,
G
=
u
r
s
,
G
⋅
x
+
c
d
t
ˉ
r
G
+
M
w
r
s
,
G
⋅
Z
W
D
r
+
γ
j
G
⋅
I
ˉ
r
,
1
s
,
G
+
ε
r
,
j
s
,
G
p
r
,
j
s
,
R
=
u
r
s
,
R
⋅
x
+
c
d
t
ˉ
r
G
+
c
ISB
r
R
+
c
Θ
r
,
j
s
,
R
+
M
w
r
s
,
R
⋅
Z
W
D
r
+
γ
j
R
⋅
I
ˉ
r
,
1
s
,
R
+
ε
r
,
j
s
,
R
p
r
,
j
s
,
C
=
u
r
s
,
C
⋅
x
+
c
d
t
ˉ
r
G
+
c
I
S
B
r
C
+
M
w
r
s
,
C
⋅
Z
W
D
r
+
γ
j
c
⋅
I
ˉ
r
,
1
s
,
C
+
ε
r
,
j
s
,
C
p
r
,
j
s
,
E
=
u
r
s
,
E
⋅
x
+
c
d
t
ˉ
r
G
+
c
I
S
B
r
E
+
M
w
r
s
,
E
⋅
Z
W
D
r
+
γ
j
E
⋅
I
ˉ
r
,
l
s
,
E
+
ε
r
,
j
s
,
E
p
r
,
j
s
,
J
=
u
r
s
,
J
⋅
x
+
c
d
t
ˉ
r
G
+
c
I
S
B
r
J
+
M
w
r
s
,
J
⋅
Z
W
D
r
+
γ
j
J
⋅
I
ˉ
r
,
1
s
,
J
+
ε
r
,
j
s
,
J
{
l
r
,
j
s
,
G
=
u
r
s
,
G
⋅
x
+
c
d
t
ˉ
r
G
+
M
w
r
s
,
G
⋅
Z
W
D
r
−
γ
j
G
⋅
I
ˉ
r
,
1
s
,
G
+
λ
j
s
,
G
⋅
N
ˉ
r
,
j
s
,
G
+
ξ
r
,
j
s
,
G
l
r
,
j
s
,
R
=
u
r
s
,
R
⋅
x
+
c
d
t
ˉ
r
G
+
c
I
S
B
r
R
+
M
w
r
s
,
R
⋅
Z
W
D
r
−
γ
j
R
⋅
I
ˉ
r
,
1
s
,
R
+
λ
j
s
,
R
⋅
N
ˉ
r
,
j
s
,
R
+
ξ
r
,
j
s
,
R
l
r
,
j
s
,
C
=
u
r
s
,
C
⋅
x
+
c
d
t
ˉ
r
G
+
c
I
S
B
r
C
+
M
w
r
s
,
C
⋅
Z
W
D
r
−
γ
j
c
⋅
I
ˉ
r
,
1
s
,
C
+
λ
j
s
,
C
⋅
N
ˉ
r
,
j
s
,
C
+
ξ
r
,
j
s
,
C
(
4.4
)
l
r
,
j
s
,
E
=
u
r
s
,
E
⋅
x
+
c
d
t
ˉ
r
G
+
c
I
S
B
r
E
+
M
w
r
s
,
E
⋅
Z
W
D
r
−
γ
j
E
⋅
I
ˉ
r
,
1
s
,
E
+
λ
j
s
,
E
⋅
N
ˉ
r
,
j
s
,
E
+
ξ
r
,
j
s
,
E
l
r
,
j
s
,
J
=
u
r
s
,
J
⋅
x
+
c
d
t
ˉ
r
G
+
c
I
S
B
r
J
+
M
w
r
s
,
J
⋅
Z
W
D
r
−
γ
j
J
⋅
I
ˉ
r
,
1
s
,
J
+
λ
j
s
,
J
⋅
N
ˉ
r
,
j
s
,
J
+
ξ
r
,
j
s
,
J
\begin{aligned} &\left\{\begin{array}{l} p_{r, j}^{s, G}=\mathbf{u}_{r}^{s, G} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+\mathrm{Mw}_{r}^{s, G} \cdot \mathrm{ZWD}_{r}+\gamma_{j}^{G} \cdot \bar{I}_{r, 1}^{s, G}+\varepsilon_{r, j}^{s, G} \\ p_{r, j}^{s, R}=\mathbf{u}_{r}^{s, R} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+c \operatorname{ISB}_{r}^{R}+c \Theta_{r, j}^{s, R}+\mathrm{Mw}_{r}^{s, R} \cdot \mathrm{ZWD}_{r}+\gamma_{j}^{R} \cdot \bar{I}_{r, 1}^{s, R}+\varepsilon_{r, j}^{s, R} \\ p_{r, j}^{s, C}=\mathbf{u}_{r}^{s, C} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+c \mathrm{ISB}_{r}^{C}+\mathrm{Mw}_{r}^{s, C} \cdot \mathrm{ZWD}_{r}+\gamma_{j}^{c} \cdot \bar{I}_{r, 1}^{s, C}+\varepsilon_{r, j}^{s, C} \\ p_{r, j}^{s, E}=\mathbf{u}_{r}^{s, E} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+c \mathrm{ISB}_{r}^{E}+\mathrm{Mw}_{r}^{s, E} \cdot \mathrm{ZWD}_{r}+\gamma_{j}^{E} \cdot \bar{I}_{r, \mathrm{l}}^{s, E}+\varepsilon_{r, j}^{s, E} \\ p_{r, j}^{s, J}=\mathbf{u}_{r}^{s, J} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+c \mathrm{ISB}_{r}^{J}+\mathrm{Mw}_{r}^{s, J} \cdot \mathrm{ZWD}_{r}+\gamma_{j}^{J} \cdot \bar{I}_{r, 1}^{s, J}+\varepsilon_{r, j}^{s, J} \end{array}\right.\\ &\left\{\begin{array}{l} l_{r, j}^{s, G}=\mathbf{u}_{r}^{s, G} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+\mathrm{Mw}_{r}^{s, G} \cdot \mathrm{ZWD}_{r}-\gamma_{j}^{G} \cdot \bar{I}_{r, 1}^{s, G}+\lambda_{j}^{s, G} \cdot \bar{N}_{r, j}^{s, G}+\xi_{r, j}^{s, G} \\ l_{r, j}^{s, R}=\mathbf{u}_{r}^{s, R} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+c \mathrm{ISB}_{r}^{R}+\mathrm{Mw}_{r}^{s, R} \cdot \mathrm{ZWD}_{r}-\gamma_{j}^{R} \cdot \bar{I}_{r, 1}^{s, R}+\lambda_{j}^{s, R} \cdot \bar{N}_{r, j}^{s, R}+\xi_{r, j}^{s, R} \\ l_{r, j}^{s, C}=\mathbf{u}_{r}^{s, C} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+c \mathrm{ISB}_{r}^{C}+\mathrm{Mw}_{r}^{s, C} \cdot \mathrm{ZWD}_{r}-\gamma_{j}^{c} \cdot \bar{I}_{r, 1}^{s, C}+\lambda_{j}^{s, C} \cdot \bar{N}_{r, j}^{s, C}+\xi_{r, j}^{s, C} \quad(4.4) \\ l_{r, j}^{s, E}=\mathbf{u}_{r}^{s, E} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+c \mathrm{ISB}_{r}^{E}+\mathrm{M} \mathrm{w}_{r}^{s, E} \cdot \mathrm{ZWD}_{r}-\gamma_{j}^{E} \cdot \bar{I}_{r, 1}^{s, E}+\lambda_{j}^{s, E} \cdot \bar{N}_{r, j}^{s, E}+\xi_{r, j}^{s, E} \\ l_{r, j}^{s, J}=\mathbf{u}_{r}^{s, J} \cdot \mathbf{x}+c d \bar{t}_{r}^{G}+c \mathrm{ISB}_{r}^{J}+\mathrm{Mw}_{r}^{s, J} \cdot \mathrm{ZWD}_{r}-\gamma_{j}^{J} \cdot \bar{I}_{r, 1}^{s, J}+\lambda_{j}^{s, J} \cdot \bar{N}_{r, j}^{s, J}+\xi_{r, j}^{s, J} \end{array}\right. \end{aligned}
⎩⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎧pr,js,G=urs,G⋅x+cdtˉrG+Mwrs,G⋅ZWDr+γjG⋅Iˉr,1s,G+εr,js,Gpr,js,R=urs,R⋅x+cdtˉrG+cISBrR+cΘr,js,R+Mwrs,R⋅ZWDr+γjR⋅Iˉr,1s,R+εr,js,Rpr,js,C=urs,C⋅x+cdtˉrG+cISBrC+Mwrs,C⋅ZWDr+γjc⋅Iˉr,1s,C+εr,js,Cpr,js,E=urs,E⋅x+cdtˉrG+cISBrE+Mwrs,E⋅ZWDr+γjE⋅Iˉr,ls,E+εr,js,Epr,js,J=urs,J⋅x+cdtˉrG+cISBrJ+Mwrs,J⋅ZWDr+γjJ⋅Iˉr,1s,J+εr,js,J⎩⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎧lr,js,G=urs,G⋅x+cdtˉrG+Mwrs,G⋅ZWDr−γjG⋅Iˉr,1s,G+λjs,G⋅Nˉr,js,G+ξr,js,Glr,js,R=urs,R⋅x+cdtˉrG+cISBrR+Mwrs,R⋅ZWDr−γjR⋅Iˉr,1s,R+λjs,R⋅Nˉr,js,R+ξr,js,Rlr,js,C=urs,C⋅x+cdtˉrG+cISBrC+Mwrs,C⋅ZWDr−γjc⋅Iˉr,1s,C+λjs,C⋅Nˉr,js,C+ξr,js,C(4.4)lr,js,E=urs,E⋅x+cdtˉrG+cISBrE+Mwrs,E⋅ZWDr−γjE⋅Iˉr,1s,E+λjs,E⋅Nˉr,js,E+ξr,js,Elr,js,J=urs,J⋅x+cdtˉrG+cISBrJ+Mwrs,J⋅ZWDr−γjJ⋅Iˉr,1s,J+λjs,J⋅Nˉr,js,J+ξr,js,J
其中,
ISB
r
Q
=
(
d
r
,
I
F
12
Q
−
d
r
,
I
F
2
G
)
+
(
d
D
Q
−
d
D
G
)
(
Q
≠
G
)
\operatorname{ISB}_{r}^{Q}=\left(d_{r, \mathrm{IF}_{12}}^{Q}-d_{r, \mathrm{IF}_{2}}^{G}\right)+\left(d D^{Q}-d D^{G}\right) \quad(Q \neq G)
ISBrQ=(dr,IF12Q−dr,IF2G)+(dDQ−dDG)(Q=G)
式中,
I
S
B
r
Q
\mathrm{ISB}_{r}^{Q}
ISBrQ 表示
Q
Q
Q 系统的 ISB 参数。可以看出,ISB 不仅源于不同 GNSS 系统
对应的接收机硬件延迟差异
(
d
r
,
1
F
2
Q
−
d
r
,
I
F
2
G
,
\left(d_{r, 1 \mathrm{~F}_{2}}^{Q}-d_{r, \mathrm{IF}_{2}}^{G},\right.
(dr,1 F2Q−dr,IF2G, 与接收机有关
)
,
),
), 而且还源于不同
GNSS 钟差产品相应的不同钟差基准约束引入的时间差异
(
d
D
Q
−
d
D
G
\left(d D^{Q}-d D^{G}\right.
(dDQ−dDG 与接收机
无关
)
)
) 。该模型的待估参数向量 X 为:
X
=
[
x
,
c
d
t
ˉ
r
G
,
c
ISB
r
Q
,
Z
W
D
r
,
I
‾
r
,
1
G
,
I
‾
r
,
1
Q
,
N
‾
r
,
j
G
,
N
‾
r
,
j
Q
]
T
\mathbf{X}=\left[\mathbf{x}, c d \bar{t}_{r}^{G}, c \operatorname{ISB}_{r}^{Q}, \mathbf{Z W D}_{r}, \overline{\mathbf{I}}_{r, 1}^{G}, \overline{\mathbf{I}}_{r, 1}^{Q}, \overline{\mathbf{N}}_{r, j}^{G}, \overline{\mathbf{N}}_{r, j}^{Q}\right]^{T}
X=[x,cdtˉrG,cISBrQ,ZWDr,Ir,1G,Ir,1Q,Nr,jG,Nr,jQ]T
模型化 ISB 参数的三个动态模型为即时间常数、随机游走过程和白噪声过程。
周锋. 多系统GNSS非差非组合精密单点定位相关理论和方法研究 [D]; 华东师范大学, 2018.
1.2.1常数估计
ISB 参数可作为时间常数来估计:
ISB
r
,
0
s
(
k
)
=
ISB
r
s
(
k
−
1
)
\operatorname{ISB}_{r, 0}^{s}(k)=\operatorname{ISB}_{r}^{s}(k-1)
ISBr,0s(k)=ISBrs(k−1)
其中,
σ
I
S
B
,
.
0
Q
(
k
)
2
=
σ
I
S
B
2
\sigma_{\mathrm{ISB}_{, .0}^{Q}(k)}^{2}=\sigma_{\mathrm{ISB}}^{2}
σISB,.0Q(k)2=σISB2
式中,
k
k
k 表示历元号
;
ISB
r
,
0
s
(
k
)
; \operatorname{ISB}_{r, 0}^{s}(k)
;ISBr,0s(k) 表示第
k
k
k 历元 ISB 的初值;
ISB
r
s
(
k
−
1
)
\operatorname{ISB}_{r}^{s}(k-1)
ISBrs(k−1) 是第
k
−
1
k-1
k−1
历元 ISB 估计值
;
σ
I
S
B
,
0
2
(
k
)
2
; \sigma_{\mathrm{ISB}_{, 0}^{2}(k)}^{2}
;σISB,02(k)2 是第
k
k
k 历元 ISB 参数的先验方差,而
σ
I
S
B
s
(
k
−
1
)
2
\sigma_{\mathrm{ISB}^{s}(k-1)}^{2}
σISBs(k−1)2 是第
k
−
1
k-1
k−1
历元 ISB 参数更新后的方差
1.2.2随机游走过程
ISB 参数可看作随机游走过程来估计:
ISB
r
,
0
s
(
k
)
=
ISB
r
Q
(
k
−
1
)
+
ω
I
S
B
r
,
0
e
(
k
)
,
ω
I
S
B
r
,
0
s
(
k
)
∼
N
(
0
,
σ
ω
I
S
B
r
,
0
2
(
k
)
2
)
\operatorname{ISB}_{r, 0}^{s}(k)=\operatorname{ISB}_{r}^{Q}(k-1)+\omega_{\mathrm{ISB}_{r, 0}^{e}(k)}, \omega_{\mathrm{ISB}_{r, 0}^{s}(k)} \sim N\left(0, \sigma_{\omega_{\mathrm{ISB}_{r, 0}^{2}(k)}^{2}}\right)
ISBr,0s(k)=ISBrQ(k−1)+ωISBr,0e(k),ωISBr,0s(k)∼N(0,σωISBr,02(k)2)
其中,
σ
I
S
B
r
,
0
s
(
k
)
2
=
σ
I
S
B
r
s
(
k
−
1
)
2
+
σ
ω
I
S
B
r
,
0
s
(
k
)
2
\sigma_{\mathrm{ISB}_{r, 0}^{s}(k)}^{2}=\sigma_{\mathrm{ISB}_{r}^{s}(k-1)}^{2}+\sigma_{\omega_{\mathrm{ISB}_{r, 0}^{s}(k)}^{2}}
σISBr,0s(k)2=σISBrs(k−1)2+σωISBr,0s(k)2
式中,ISB 参数的变化部分
(
ω
I
S
B
r
,
0
e
(
k
)
)
\left(\omega_{\mathrm{ISB}_{r, 0}^{e}(k)}\right)
(ωISBr,0e(k)) 对应的方差为
ω
I
S
B
r
,
0
s
(
k
)
2
\omega_{\mathrm{ISB}_{r, 0}^{s}(k)}^{2}
ωISBr,0s(k)2
1.3.3.白噪声过程
ISB 参数可看作白噪声过程来估计:
ISB
r
,
0
Q
(
k
)
∼
N
(
ISB
r
,
s
p
p
Q
(
k
)
,
σ
I
S
B
r
,
p
r
i
e
(
k
)
2
)
\operatorname{ISB}_{r, 0}^{Q}(k) \sim N\left(\operatorname{ISB}_{r, \mathrm{spp}}^{Q}(k), \sigma_{\mathrm{ISB}_{r, \mathrm{pri}}^{e}(k)}^{2}\right)
ISBr,0Q(k)∼N(ISBr,sppQ(k),σISBr,prie(k)2)
式中,
ISB
r
,
spp
Q
(
k
)
\operatorname{ISB}_{r, \text { spp }}^{Q}(k)
ISBr, spp Q(k) 表示伪距单点定位
(
S
P
P
)
(\mathrm{SPP})
(SPP) 计算的
I
S
B
\mathrm{ISB}
ISB 值
;
σ
I
S
B
r
,
p
r
i
s
(
k
)
2
; \sigma_{\mathrm{ISB}_{r, \mathrm{pri}}^{s}(k)}^{2}
;σISBr,pris(k)2 是 ISB 参数
的先验方差。
2.ISPB
It has been known that it is the ISPB originating in stations that can inhibit inter-GNSS double-difference AR (e.g. Odijk and Teunissen 2013 ). Regarding PPP, the hurdle is that station UPDs can no longer be removed through single differencing between satellites belonging to different GNSS. In contrast to Eq. 5, we have inter-GNSS single-difference ambiguities
N
^
i
,
m
k
q
=
L
i
,
m
k
λ
G
,
m
−
L
i
,
m
q
λ
C
,
m
=
N
˘
i
,
m
k
q
+
b
ˉ
i
,
m
G
C
−
b
ˉ
m
k
q
λ
G
,
n
N
^
i
,
1
k
−
λ
C
,
n
N
^
i
,
1
q
=
λ
G
,
1
N
^
i
,
3
k
−
λ
C
,
1
N
^
i
,
3
q
−
λ
G
,
m
g
G
+
1
N
~
i
,
m
k
+
λ
C
,
m
g
C
+
1
N
~
i
,
m
q
\begin{array}{l} \hat{N}_{i, \mathrm{~m}}^{k q}=\frac{L_{i, \mathrm{~m}}^{k}}{\lambda_{\mathrm{G}, \mathrm{m}}}-\frac{L_{i, \mathrm{~m}}^{q}}{\lambda_{\mathrm{C}, \mathrm{m}}}=\breve{N}_{i, \mathrm{~m}}^{k q}+\bar{b}_{i, \mathrm{~m}}^{\mathrm{GC}}-\bar{b}_{\mathrm{m}}^{k q} \\ \lambda_{\mathrm{G}, \mathrm{n}} \hat{N}_{i, 1}^{k}-\lambda_{\mathrm{C}, \mathrm{n}} \hat{N}_{i, 1}^{q}=\lambda_{\mathrm{G}, 1} \hat{N}_{i, 3}^{k}-\lambda_{\mathrm{C}, 1} \hat{N}_{i, 3}^{q} \\ \quad-\frac{\lambda \mathrm{G}, \mathrm{m}}{g_{\mathrm{G}}+1} \widetilde{N}_{i, \mathrm{~m}}^{k}+\frac{\lambda \mathrm{C}, \mathrm{m}}{g \mathrm{C}+1} \widetilde{N}_{i, \mathrm{~m}}^{q} \end{array}
N^i, mkq=λG,mLi, mk−λC,mLi, mq=N˘i, mkq+bˉi, mGC−bˉmkqλG,nN^i,1k−λC,nN^i,1q=λG,1N^i,3k−λC,1N^i,3q−gG+1λG,mN
i, mk+gC+1λC,mN
i, mq
where satellite
k
k
k is from GPS and
q
q
q is from BeiDou. Equation
8
b
8 \mathrm{~b}
8 b can be further developed into
λ
G
,
n
N
^
i
,
1
k
q
+
(
λ
G
,
n
−
λ
C
,
n
)
N
^
i
,
1
q
=
λ
G
,
1
N
^
i
,
3
k
−
λ
C
,
1
N
^
i
,
3
q
−
λ
G
,
m
g
G
+
1
N
~
i
,
m
k
q
−
(
λ
G
,
m
g
G
+
1
−
λ
C
,
m
g
C
+
1
)
N
~
i
,
m
q
\begin{array}{r} \lambda_{\mathrm{G}, \mathrm{n}} \hat{N}_{i, 1}^{k q}+\left(\lambda_{\mathrm{G}, \mathrm{n}}-\lambda_{\mathrm{C}, \mathrm{n}}\right) \hat{N}_{i, 1}^{q}=\lambda_{\mathrm{G}, 1} \hat{N}_{i, 3}^{k}-\lambda_{\mathrm{C}, 1} \hat{N}_{i, 3}^{q} \\ -\frac{\lambda \mathrm{G}, \mathrm{m}}{g_{\mathrm{G}}+1} \widetilde{N}_{i, \mathrm{~m}}^{k q}-\left(\frac{\lambda \mathrm{G}, \mathrm{m}}{g_{\mathrm{G}}+1}-\frac{\lambda \mathrm{C}, \mathrm{m}}{g_{\mathrm{C}}+1}\right) \widetilde{N}_{i, \mathrm{~m}}^{q} \end{array}
λG,nN^i,1kq+(λG,n−λC,n)N^i,1q=λG,1N^i,3k−λC,1N^i,3q−gG+1λG,mN
i, mkq−(gG+1λG,m−gC+1λC,m)N
i, mq
and afterward
λ
G
,
n
N
^
i
,
1
k
q
=
λ
G
,
1
N
^
i
,
3
k
−
λ
C
,
1
N
^
i
,
3
q
−
(
λ
G
,
n
−
λ
C
,
n
)
N
^
i
,
1
q
−
λ
G
,
m
g
G
+
1
N
~
i
,
m
k
q
−
(
λ
G
,
m
g
G
+
1
−
λ
C
,
m
g
C
+
1
)
N
~
i
,
m
q
=
λ
G
,
n
(
N
ˇ
i
,
1
k
q
+
b
ˉ
i
,
n
G
C
−
b
ˉ
n
k
q
)
\begin{aligned} \lambda_{\mathrm{G}, \mathrm{n}} \hat{N}_{i, 1}^{k q}=& \lambda_{\mathrm{G}, 1} \hat{N}_{i, 3}^{k}-\lambda_{\mathrm{C}, 1} \hat{N}_{i, 3}^{q}-\left(\lambda_{\mathrm{G}, \mathrm{n}}-\lambda \mathrm{C}, \mathrm{n}\right) \hat{N}_{i, 1}^{q} \\ &-\frac{\lambda \mathrm{G}, \mathrm{m}}{g_{\mathrm{G}}+1} \widetilde{N}_{i, \mathrm{~m}}^{k q}-\left(\frac{\lambda \mathrm{G}, \mathrm{m}}{g_{\mathrm{G}}+1}-\frac{\lambda \mathrm{C}, \mathrm{m}}{g \mathrm{C}+1}\right) \widetilde{N}_{i, \mathrm{~m}}^{q} \\ =& \lambda_{\mathrm{G}, \mathrm{n}}\left(\check{N}_{i, 1}^{k q}+\bar{b}_{i, \mathrm{n}}^{\mathrm{GC}}-\bar{b}_{\mathrm{n}}^{k q}\right) \end{aligned}
λG,nN^i,1kq==λG,1N^i,3k−λC,1N^i,3q−(λG,n−λC,n)N^i,1q−gG+1λG,mN
i, mkq−(gG+1λG,m−gC+1λC,m)N
i, mqλG,n(Nˇi,1kq+bˉi,nGC−bˉnkq)
We note that
{
b
ˉ
i
,
m
G
C
=
b
ˉ
i
,
m
G
−
b
ˉ
i
,
m
C
≠
0
b
ˉ
i
,
n
G
C
=
b
ˉ
i
,
n
G
−
b
ˉ
i
,
n
C
≠
0
\left\{\begin{array}{l} \bar{b}_{i, \mathrm{~m}}^{\mathrm{GC}}=\bar{b}_{i, \mathrm{~m}}^{\mathrm{G}}-\bar{b}_{i, \mathrm{~m}}^{\mathrm{C}} \neq 0 \\ \bar{b}_{i, \mathrm{n}}^{\mathrm{GC}}=\bar{b}_{i, \mathrm{n}}^{\mathrm{G}}-\bar{b}_{i, \mathrm{n}}^{\mathrm{C}} \neq 0 \end{array}\right.
{bˉi, mGC=bˉi, mG−bˉi, mC=0bˉi,nGC=bˉi,nG−bˉi,nC=0
which are the wide-lane and narrow-lane ISPB for station
i
i
i between GPS and BeiDou, respectively. Equation 9 is used to compute inter-GPS/BeiDou single-difference narrow-lane ambiguities
(
N
^
i
,
1
k
q
)
.
\left(\hat{N}_{i, 1}^{k q}\right) .
(N^i,1kq). On the right-hand side of this equation, we then need to know
N
^
i
,
3
k
\hat{N}_{i, 3}^{k}
N^i,3k and
N
^
i
,
3
q
\hat{N}_{i, 3}^{q}
N^i,3q which are directly PPP output, and also
N
~
i
,
m
k
q
\tilde{N}_{i, \mathrm{~m}}^{k q}
N~i, mkq which is provided by resolving Eq.
8
a
8 \mathrm{a}
8a. However,
N
~
i
,
m
q
\tilde{N}_{i, \mathrm{~m}}^{q}
N~i, mq cannot be derived since we are unable to resolve undifferenced wide-lane ambiguities with Eq. 8 a;
N
^
i
,
1
q
\hat{N}_{i, 1}^{q}
N^i,1q cannot be obtained either due to its dependence on
N
~
i
,
m
q
\tilde{N}_{i, \mathrm{~m}}^{q}
N~i, mq (see Eq. 3b).
Despite this difficulty, we note that the absolute values of the coefficients for
N
~
i
,
m
q
\tilde{N}_{i, \mathrm{~m}}^{q}
N~i, mq and
N
^
i
,
1
q
\hat{N}_{i, 1}^{q}
N^i,1q in Eq. 9 satisfy that
{
∣
λ
G
,
m
g
G
+
1
−
λ
C
,
m
g
C
+
1
∣
<
0.82
c
m
λ
G
,
n
−
λ
C
,
n
∣
<
0.14
c
m
\left\{\begin{array}{l} \left|\frac{\lambda_{\mathrm{G}, \mathrm{m}}}{g_{\mathrm{G}}+1}-\frac{\lambda_{\mathrm{C}, \mathrm{m}}}{g_{\mathrm{C}}+1}\right|<0.82 \mathrm{~cm} \\ \lambda_{\mathrm{G}, \mathrm{n}}-\lambda_{\mathrm{C}, \mathrm{n}} \mid<0.14 \mathrm{~cm} \end{array}\right.
{∣∣∣gG+1λG,m−gC+1λC,m∣∣∣<0.82 cmλG,n−λC,n∣<0.14 cm
Both are far smaller than
λ
G
,
n
=
10.695
c
m
,
\lambda_{\mathrm{G}, \mathrm{n}}=10.695 \mathrm{~cm},
λG,n=10.695 cm, hence implying that minor errors within
N
~
i
,
m
q
\tilde{N}_{i, \mathrm{~m}}^{q}
N~i, mq and
N
^
i
,
1
q
,
\hat{N}_{i, 1}^{q},
N^i,1q, if present, will have limited impact on
N
^
i
,
1
k
q
.
\hat{N}_{i, 1}^{k q} .
N^i,1kq. We therefore introduce the following approximations
{
N
~
i
,
m
q
≈
N
^
i
,
m
q
=
L
i
,
m
q
λ
C
,
m
N
^
i
,
1
q
≈
g
C
+
1
g
C
N
^
i
,
3
q
−
1
g
C
−
1
N
^
i
,
m
q
\left\{\begin{aligned} \tilde{N}_{i, \mathrm{~m}}^{q} & \approx \hat{N}_{i, \mathrm{~m}}^{q}=\frac{L_{i, \mathrm{~m}}^{q}}{\lambda_{\mathrm{C}, \mathrm{m}}} \\ \hat{N}_{i, 1}^{q} & \approx \frac{g \mathrm{C}+1}{g_{\mathrm{C}}} \hat{N}_{i, 3}^{q}-\frac{1}{g_{\mathrm{C}}-1} \hat{N}_{i, \mathrm{~m}}^{q} \end{aligned}\right.
⎩⎪⎪⎪⎨⎪⎪⎪⎧N~i, mqN^i,1q≈N^i, mq=λC,mLi, mq≈gCgC+1N^i,3q−gC−11N^i, mq
which refer to Eq. 3 . Contrasting Equations
3
a
3 \mathrm{a}
3a and
4
,
4,
4, we find that the approximation errors of Eq. 12 consist in the
difference between
(
b
ˉ
i
,
m
C
−
b
ˉ
m
q
)
\left(\bar{b}_{i, \mathrm{~m}}^{\mathrm{C}}-\bar{b}_{\mathrm{m}}^{q}\right)
(bˉi, mC−bˉmq) and
(
b
ˉ
^
i
,
m
C
−
b
ˉ
^
m
q
)
\left(\hat{\bar{b}}_{i, \mathrm{~m}}^{\mathrm{C}}-\hat{\bar{b}}_{\mathrm{m}}^{q}\right)
(bˉ^i, mC−bˉ^mq) where the
former contributes to the estimation of the latter. The impact of this approximation on
N
^
i
,
1
k
q
\hat{N}_{i, 1}^{k q}
N^i,1kq will be inspected in Sect. 2.4 . With Eqs.
8
a
8 \mathrm{a}
8a and 9 along with the substitution of Eq. 12 , we can calculate
N
^
i
,
m
k
q
\hat{N}_{i, \mathrm{~m}}^{k q}
N^i, mkq and
N
^
i
,
1
k
q
,
\hat{N}_{i, 1}^{k q},
N^i,1kq, similar to Eq.
5
,
5,
5, aiming at estimating station ISPBs.
We first presume that ISPBs
b
ˉ
i
,
m
G
C
\bar{b}_{i, \mathrm{~m}}^{\mathrm{GC}}
bˉi, mGC and
b
ˉ
i
,
n
G
C
\bar{b}_{i, \mathrm{n}}^{\mathrm{GC}}
bˉi,nGC should remain the same for all pairs of GPS/BeiDou satellites observed at station
i
i
i. From Eqs. 8 a and
9
,
9,
9, we can form double-difference ambiguities between stations
i
i
i and
j
j
j such that satellite FCBs
b
ˉ
m
k
q
\bar{b}_{\mathrm{m}}^{k q}
bˉmkq and
b
ˉ
n
k
q
\bar{b}_{\mathrm{n}}^{\mathrm{k} q}
bˉnkq cancel completely, leading to
{
N
^
i
j
,
m
k
q
=
N
ˇ
i
j
,
m
k
q
+
b
ˉ
i
j
,
m
G
C
N
^
i
j
,
1
k
q
=
N
ˇ
i
j
,
1
k
q
+
b
ˉ
i
j
,
n
G
C
\left\{\begin{array}{l} \hat{N}_{i j, \mathrm{~m}}^{k q}=\check{N}_{i j, \mathrm{~m}}^{k q}+\bar{b}_{i j, \mathrm{~m}}^{\mathrm{GC}} \\ \hat{N}_{i j, 1}^{k q}=\check{N}_{i j, 1}^{k q}+\bar{b}_{i j, \mathrm{n}}^{\mathrm{GC}} \end{array}\right.
{N^ij, mkq=Nˇij, mkq+bˉij, mGCN^ij,1kq=Nˇij,1kq+bˉij,nGC
where
{
N
~
i
j
,
∗
k
q
=
N
~
i
,
∗
k
q
−
N
~
j
,
∗
k
q
N
^
i
j
,
∗
k
q
=
N
^
i
,
∗
k
q
−
N
^
j
,
∗
k
q
b
ˉ
i
j
,
∗
G
C
=
b
ˉ
i
,
∗
G
C
−
b
ˉ
j
,
∗
G
C
\left\{\begin{array}{l} \tilde{N}_{i j, *}^{k q}=\tilde{N}_{i, *}^{k q}-\tilde{N}_{j, *}^{k q} \\ \hat{N}_{i j, *}^{k q}=\hat{N}_{i, *}^{k q}-\hat{N}_{j, *}^{k q} \\ \bar{b}_{i j, *}^{G C}=\bar{b}_{i, *}^{\mathrm{GC}}-\bar{b}_{j, *}^{\mathrm{GC}} \end{array}\right.
⎩⎪⎨⎪⎧N~ij,∗kq=N~i,∗kq−N~j,∗kqN^ij,∗kq=N^i,∗kq−N^j,∗kqbˉij,∗GC=bˉi,∗GC−bˉj,∗GC
Note that for intra-system double-difference ambiguities, the ISPB terms
b
ˉ
i
j
,
m
G
C
\bar{b}_{i j, \mathrm{~m}}^{\mathrm{GC}}
bˉij, mGC and
b
ˉ
i
j
,
n
G
C
\bar{b}_{i j, \mathrm{n}}^{\mathrm{GC}}
bˉij,nGC will be absent, and
N
^
i
j
,
m
k
q
\hat{N}_{i j, \mathrm{~m}}^{k q}
N^ij, mkq and
N
^
i
j
,
1
k
q
\hat{N}_{i j, 1}^{k q}
N^ij,1kq
have the integer property naturally. In case of inter-system AR, however, we have to correct for
b
ˉ
i
j
,
m
G
C
\bar{b}_{i j, \mathrm{~m}}^{\mathrm{GC}}
bˉij, mGC and
b
ˉ
i
j
,
n
G
C
\bar{b}_{i j, \mathrm{n}}^{\mathrm{GC}}
bˉij,nGC before being able to resolve
N
^
i
j
,
m
k
q
\hat{N}_{i j, \mathrm{~m}}^{k q}
N^ij, mkq and
N
^
i
j
,
1
k
q
\hat{N}_{i j, 1}^{k q}
N^ij,1kq
Similar to Eq. 7 , 7, 7, for all stations collecting GPS and BeiDou data simultaneously, we can calculate differential ISPBs
between stations
i
i
i and
j
j
j using
{
b
ˉ
^
i
j
,
m
G
C
=
N
^
i
j
,
m
k
q
−
[
N
^
i
j
,
m
k
q
]
b
^
i
j
,
n
G
C
=
N
^
i
j
,
1
k
q
−
[
N
^
i
j
,
1
k
q
]
\left\{\begin{array}{l} \hat{\bar{b}}_{i j, \mathrm{~m}}^{\mathrm{GC}}=\hat{N}_{i j, \mathrm{~m}}^{k q}-\left[\hat{N}_{i j, \mathrm{~m}}^{k q}\right] \\ \hat{b}_{i j, \mathrm{n}}^{\mathrm{GC}}=\hat{N}_{i j, 1}^{k q}-\left[\hat{N}_{i j, 1}^{k q}\right] \end{array}\right.
⎩⎨⎧bˉ^ij, mGC=N^ij, mkq−[N^ij, mkq]b^ij,nGC=N^ij,1kq−[N^ij,1kq]
Here, we have presumed that
b
ˉ
i
j
,
m
G
C
\bar{b}_{i j, \mathrm{~m}}^{\mathrm{GC}}
bˉij, mGC and
b
ˉ
i
j
,
n
G
C
\bar{b}_{i j, \mathrm{n}}^{\mathrm{GC}}
bˉij,nGC are both merely the fractional parts of ISPB and their integer parts have been assimilated into
N
˘
i
j
,
m
k
q
\breve{N}_{i j, \mathrm{~m}}^{k q}
N˘ij, mkq and
N
˘
i
j
,
1
k
q
,
\breve{N}_{i j, 1}^{k q},
N˘ij,1kq, respectively. For any given baseline, we will have a good number of interGPS/BeiDou double-difference ambiguities, and differential ISPB between the two relevant stations can thus be computed by averaging all estimates of
b
ˉ
i
j
,
m
G
C
\bar{b}_{i j, \mathrm{~m}}^{\mathrm{GC}}
bˉij, mGC and
b
ˉ
i
j
,
n
G
C
\bar{b}_{i j, \mathrm{n}}^{\mathrm{GC}}
bˉij,nGC. The resultant ISPBs can afterwards be converted into undifferenced values (e.g.,
b
ˉ
i
,
m
G
C
\bar{b}_{i, \mathrm{~m}}^{\mathrm{GC}}
bˉi, mGC and
b
ˉ
i
,
n
G
C
\bar{b}_{i, \mathrm{n}}^{\mathrm{GC}}
bˉi,nGC ) by choosing a reference station ISPB.