前言
GNSS定位中不少技术手段如PPP和长基线RTK需要将电离层延迟作为参数估计,电离层延迟的变化通常被描述为随机游走过程,而功率谱密度是描述随机游走的决定性参数。功率谱密度的平方乘以时间间隔可得到随机游走过程的方差,这在卡尔曼滤波中称过程噪声。
所以,本文题为自适应调整过程噪声,实际就是实时计算电离层随机游走过程的功率谱密度。
本文以非差非组合的PPP为例。
一、非差非组合PPP模型
已知原始观测方程如下:
电离层延迟参数的随机游走过程可被描述为:
已知过程噪声与功率谱密度的关系:
二、电离层功率谱密度计算
通过对GF组合进行历元间作差可得电离层延迟变化量(这里称为电离层观测值),以1、2频为例:
这里的公式有误,分母应该是
(
f
1
2
−
f
2
2
)
(f_1^2-f_2^2)
(f12−f22)。
该电离层观测值包含了纯净的电离层延迟变化量和观测噪声项:
上式中的观测噪声的标准差计算如下,公式来源于不同论文,请忽略符号差异:
其中
σ
ϵ
1
\sigma_{\epsilon 1}
σϵ1是载波相位观测值的标准差,这里假设不同频率的载波相位观测值精度相同;m和n为使用的频率。
电离层观测值的方差计算如下:
D
~
(
I
~
r
s
,
Δ
i
)
=
∑
i
=
1
i
=
n
{
I
~
r
s
(
i
,
i
+
Δ
i
)
}
2
/
n
(2.1)
\tilde{D}\left(\tilde{I}_r^s, \Delta i\right)=\sum_{i=1}^{i=n}\left\{ \tilde{I}_r^s(i, i+\Delta i)\right\}^2 / n \tag{2.1}
D~(I~rs,Δi)=i=1∑i=n{I~rs(i,i+Δi)}2/n(2.1)
这里的n为观测值个数。对该方差除以时间再开方即可得到电离层观测值的功率谱密度
q
~
(
I
r
s
,
Δ
i
)
\tilde{q}\left(I_r^s, \Delta i\right)
q~(Irs,Δi)。
进而可以计算纯净的电离层功率谱密度:
q
2
(
I
r
s
)
=
q
~
2
(
I
r
s
,
Δ
i
)
−
D
(
ε
ϕ
ϕ
m
)
/
Δ
T
(2.2)
q^2\left(I_r^s\right)=\tilde{q}^2\left(I_r^s, \Delta i\right)-D\left(\varepsilon_{\phi_{\phi_m}}\right) / \Delta T \tag{2.2}
q2(Irs)=q~2(Irs,Δi)−D(εϕϕm)/ΔT(2.2)
实际上,
D
(
ε
ϕ
ϕ
m
)
D\left(\varepsilon_{\phi_{\phi_m}}\right)
D(εϕϕm)无法获得。
三、具体实现方法
主要有三种具体实现方法,一是不平滑,二是三阶多项式平滑,三是Vondrak平滑。Vondrak平滑较复杂,不讨论,这里只讨论前两种方法。
3.1 不平滑
按上文方法计算每个历元的电离层观测值,进而计算一个滑动窗口内的方差。忽略观测噪声,将式(2.2)简化为:
q
2
(
I
r
s
)
=
q
~
2
(
I
r
s
,
Δ
i
)
q^2\left(I_r^s\right)=\tilde{q}^2\left(I_r^s, \Delta i\right)
q2(Irs)=q~2(Irs,Δi)
代入式(2.1)可得:
q
2
(
I
r
s
)
=
∑
i
=
1
i
=
n
{
I
~
r
s
(
i
,
i
+
Δ
i
)
}
2
/
(
n
⋅
Δ
T
)
q^2\left(I_r^s\right)=\sum_{i=1}^{i=n}\left\{ \tilde{I}_r^s(i, i+\Delta i)\right\}^2 / (n·\Delta T)
q2(Irs)=i=1∑i=n{I~rs(i,i+Δi)}2/(n⋅ΔT)
3.2 三阶多项式平滑
与不平滑方法不同的是计算方差时不是用直接计算得到的电离层观测值,而是用多项式平滑后的值。
三阶多项式拟合公式如下:
用最小二乘计算出拟合系数
a
1
a_1
a1、
a
2
a_2
a2、
a
3
a_3
a3、
a
4
a_4
a4,再带入历元数计算当前历元的电离层平滑值,进而计算电离层观测值的功率谱密度:
建议平滑窗口长度和方差计算窗口长度均取60s。
参考文献
- 徐宗秋2023 - 一种BDS非差非组合PPP中电离层功率谱密度估计方法
- 李军2023 - 电离层时变特性约束的BDS长距离RTK定位算法