实时 PPP 观测数据预处理：原理、方法与矩阵推导-CSDN博客

适用对象：GNSS 数据处理初学者、PPP/PPP-AR 算法开发人员、实时定位系统工程师
本文主线：预处理边界 → 数据与产品检查 → 信号质量控制 → 钟跳与周跳 → 粗差与鲁棒估计 → 状态管理 → 工程实现

1. 结论

实时 PPP 预处理不只是周跳探测。

完整的预处理至少包括：

时间标签、数据完整性和单位检查；
卫星系统、频率和信号映射；
卫星健康、姿态和产品可用性检查；
高度角、C/N₀、锁定时间和 LLI 检查；
数据间断、接收机钟跳和系统间偏差跳变检测；
GF、MW、多普勒和创新联合周跳探测；
伪距、相位和多普勒粗差检查；
多路径、NLOS 和码观测质量检查；
观测协方差和鲁棒权构造；
相位弧段、模糊度和协方差重置；
PPP-AR 相位偏差基准连续性检查；
实时数据延迟、乱序、丢包和产品切换处理。

预处理的输出不应只有“有效”或“无效”，而应区分三类动作：

观测明显错误
→ 剔除

观测质量下降但仍有信息
→ 降权

相位连续性或整数基准被破坏
→ 新建弧段并重置模糊度

2. 为什么实时 PPP 必须严格预处理

PPP 使用非差伪距和载波相位。

一个未处理的异常可能被滤波器吸收到：

接收机钟差；
系统间偏差；
对流层延迟；
电离层延迟；
模糊度；
接收机坐标。

设正常创新为：

$νk0\boldsymbol{\nu}_k^{0}$

当前历元存在观测粗差：

$bk\mathbf b_k$

则实际创新为：

$νk=νk0+bk\boldsymbol{\nu}_k = \boldsymbol{\nu}_k^{0} + \mathbf b_k$

卡尔曼更新为：

$x^k+=x^k−+Kkνk\widehat{\mathbf x}_k^{+} = \widehat{\mathbf x}_k^{-} + \mathbf K_k \boldsymbol{\nu}_k$

粗差引起的状态污染为：

$Δxk=Kkbk\boxed{ \boldsymbol{\Delta x}_k = \mathbf K_k \mathbf b_k }$

式中：

k：历元编号；
x⁻：量测更新前的预测状态；
x⁺：量测更新后的状态；
K：卡尔曼增益；
ν：创新；
b：观测粗差；
Δx：粗差引起的状态偏移。

实时滤波具有记忆性。异常一旦进入状态，影响可能延续多个历元。

3. 预处理、误差改正和滤波质量控制的边界

环节	主要问题	典型内容
预处理	观测能否使用、怎样使用	时间、钟跳、周跳、粗差、筛选、降权
确定性改正	正常观测应减去哪些物理效应	轨道钟差、天线、潮汐、相位缠绕
滤波估计	怎样估计未知参数	坐标、钟差、大气、模糊度
后验质量控制	更新后的解是否可信	后验残差、状态跳变、位置一致性

例如：

相位缠绕没有改正
→ 观测模型错误
→ 不能只当成周跳删除

低C/N₀
→ 随机模型或筛选问题
→ 不一定直接剔除

检测到周跳
→ 相位连续性破坏
→ 必须处理模糊度状态

4. 推荐的总体流程

第一部分：观测模型与预处理接口

5. 非组合 PPP 观测模型

频率 f 上的伪距和载波相位可概念化为：

$Pfs=ρs+c(δtr−δts)+Ts+Ifs+bP,fs+εP,fsP_f^s = \rho^s + c \left( \delta t_r-\delta t_s \right) + T^s + I_f^s + b_{P,f}^s + \varepsilon_{P,f}^s$

$Lfs=ρs+c(δtr−δts)+Ts−Ifs+λfNfs+bL,fs+εL,fsL_f^s = \rho^s + c \left( \delta t_r-\delta t_s \right) + T^s - I_f^s + \lambda_f N_f^s + b_{L,f}^s + \varepsilon_{L,f}^s$

式中：

s：卫星；
f：频率或信号；
P：伪距，单位为米；
L：载波相位距离，单位为米；
ρ：几何距离；
c：真空光速；
δt_r、δt_s：接收机和卫星钟差；
T：对流层延迟；
I_f：频率 f 上的一阶电离层延迟；
λ_f：载波波长；
N_f：模糊度；
b_P、b_L：码和相位偏差；
ε_P、ε_L：噪声、多路径和未建模误差。

严格模型还需加入：

天线 PCO/PCV；
固体潮、海潮和极潮；
相位缠绕；
相对论效应；
卫星姿态和天线参考点。

这些属于误差建模，但其输入缺失或发生不连续时，应向预处理模块产生质量标志。

6. IF-PPP 与 UC-PPP 的预处理差异

6.1 双频无电离层组合 PPP

无电离层组合为：

$PIF=αP1+βP2P_{\mathrm{IF}} = \alpha P_1 + \beta P_2$

$LIF=αL1+βL2L_{\mathrm{IF}} = \alpha L_1 + \beta L_2$

其中：

$α=f12f12−f22\alpha = \frac{ f_1^2 }{ f_1^2-f_2^2 }$

$β=−f22f12−f22\beta = - \frac{ f_2^2 }{ f_1^2-f_2^2 }$

若任一组成频率缺失或发生周跳，该卫星的双频相位组合通常不能继续沿用旧 IF 模糊度。

6.2 非组合 PPP

非组合 PPP 可保留未受影响的频率，但必须：

为各频率维护独立模糊度；
正确传播频间协方差；
估计或改正电离层；
正确处理码、相位和频间偏差。

因此：

IF-PPP
→ 频率配对要求更强

UC-PPP
→ 状态更多，信号级管理要求更强

7. 线性化模型和观测筛选

线性化观测模型为：

$lk=Hkδxk+vk\boxed{ \mathbf l_k = \mathbf H_k \delta\mathbf x_k + \mathbf v_k }$

观测协方差为：

$Rk=Cov⁡(vk)\boxed{ \mathbf R_k = \operatorname{Cov} \left( \mathbf v_k \right) }$

设原始观测有 m 条，筛选后保留 q 条，定义选择矩阵：

$Sk∈Rq×m\mathbf S_k \in \mathbb R^{q\times m}$

筛选后的模型为：

$lkuse=Sklk\boxed{ \mathbf l_k^{\mathrm{use}} = \mathbf S_k \mathbf l_k }$

$Hkuse=SkHk\boxed{ \mathbf H_k^{\mathrm{use}} = \mathbf S_k \mathbf H_k }$

$Rkuse=SkRkSkT\boxed{ \mathbf R_k^{\mathrm{use}} = \mathbf S_k \mathbf R_k \mathbf S_k^{\mathrm T} }$

删除观测时，必须同步更新观测、设计矩阵、协方差和索引。

8. 组合观测的协方差传播

设原始观测为：

$y\mathbf y$

线性组合为：

$yc=Ty\boxed{ \mathbf y_c = \mathbf T \mathbf y }$

其中：

T：组合系数矩阵；
y_c：组合观测。

组合协方差为：

$Rc=TRyTT\boxed{ \mathbf R_c = \mathbf T \mathbf R_y \mathbf T^{\mathrm T} }$

因此，IF、GF、MW 和码多路径组合不能只计算组合值，还必须传播方差和相关性。

9. 实时改正不确定度怎样进入权阵

设已应用的实时改正向量为：

$c^=c+εc\widehat{\mathbf c} = \mathbf c + \boldsymbol{\varepsilon}_c$

改正后的观测模型为：

$l′=Hδx+v−Bcεc\mathbf l' = \mathbf H \delta\mathbf x + \mathbf v - \mathbf B_c \boldsymbol{\varepsilon}_c$

其中：

c：真实改正；
ĉ：服务端提供的改正；
ε_c：改正误差；
B_c：改正误差到观测域的映射矩阵；
v：接收机观测误差。

有效协方差为：

$Reff=Rv+BcQcBcT−CvcBcT−BcCvcT\boxed{ \mathbf R_{\mathrm{eff}} = \mathbf R_v + \mathbf B_c \mathbf Q_c \mathbf B_c^{\mathrm T} - \mathbf C_{vc} \mathbf B_c^{\mathrm T} - \mathbf B_c \mathbf C_{vc}^{\mathrm T} }$

其中：

R_v：原始观测协方差；
Q_c：实时改正误差协方差；
C_vc：观测误差与改正误差的互协方差。

若二者近似独立：

$Reff=Rv+BcQcBcT\boxed{ \mathbf R_{\mathrm{eff}} = \mathbf R_v + \mathbf B_c \mathbf Q_c \mathbf B_c^{\mathrm T} }$

所以实时改正不能只做“有效或无效”的硬判断。其质量指标、龄期和不确定度也应反映到观测权阵中。

10. 筛选后的可观测性检查

对单历元独立最小二乘，活动参数应满足：

$rank⁡(Hkuse)=nact\boxed{ \operatorname{rank} \left( \mathbf H_k^{\mathrm{use}} \right) = n_{\mathrm{act}} }$

其中：

n_act：当前单历元需要独立估计的活动参数数。

卡尔曼滤波还利用预测先验，其信息矩阵为：

$Nk=(Pk−)−1+(Hkuse)T(Rkuse)−1Hkuse\boxed{ \mathbf N_k = \left( \mathbf P_k^{-} \right)^{-1} + \left( \mathbf H_k^{\mathrm{use}} \right)^{\mathrm T} \left( \mathbf R_k^{\mathrm{use}} \right)^{-1} \mathbf H_k^{\mathrm{use}} }$

应检查：

卫星几何；
各系统是否有观测支撑；
新增模糊度和电离层状态是否可初始化；
信息矩阵条件数；
解是否过度依赖先验。

11. 周跳后的状态和协方差重置

将状态分为保留状态和需要重置的模糊度：

$x=[xKNj]\mathbf x = \begin{bmatrix} \mathbf x_K \\ N_j \end{bmatrix}$

若新弧段模糊度独立初始化，可写为：

$x+=[x^K−N^j,0]\boxed{ \mathbf x^{+} = \begin{bmatrix} \widehat{\mathbf x}_K^{-} \\ \widehat N_{j,0} \end{bmatrix} }$

$P+=[PKK−00σN,02]\boxed{ \mathbf P^{+} = \begin{bmatrix} \mathbf P_{KK}^{-} & \mathbf0 \\ \mathbf0 & \sigma_{N,0}^{2} \end{bmatrix} }$

式中：

x_K：继续保留的状态；
N_j：发生周跳的模糊度；
N_j,0：新弧段初值；
σ_N,0²：新模糊度方差。

若新模糊度由当前观测联合初始化，交叉协方差不一定为零，应按初始化模型传播。

不能只把模糊度数值改成零，同时保留旧弧段的交叉协方差。

第二部分：数据流、时间和元数据

12. 历元装配、乱序和迟到数据

设相邻历元时间为：

$t_{k-1}$

$t_k$

实际采样间隔为：

$Δtk=tk−tk−1\boxed{ \Delta t_k = t_k-t_{k-1} }$

标称间隔为：

$Δt0\Delta t_0$

时间残差为：

$rt,k=Δtk−Δt0\boxed{ r_{t,k} = \Delta t_k-\Delta t_0 }$

超限可能表示：

数据丢失；
重复或乱序；
时间标签跳变；
采样率变化；
接收机重启。

实时流应设置有限历元装配窗口：

窗口内迟到
→ 允许补齐和排序

超过截止时间
→ 不回写已经完成的滤波历元

13. 单位、缺失值和信号映射

必须检查：

伪距单位；
相位是周还是米；
多普勒单位和符号；
C/N₀ 是否为 dB-Hz；
时间系统和周数；
卫星系统、频段和跟踪属性；
GLONASS FDMA 频率通道；
缺失值编码。

RINEX 观测中的缺失值可以写为空字段或 0.0。实时 RTCM 和厂商二进制流则应按各自协议处理。

相位弧段应按：

卫星 + 具体信号或跟踪模式

维护，而不是只按卫星或频段管理。

14. LLI、半周和信号强度字段

RINEX 4.02 中，LLI 仅用于载波相位观测。

LLI 位	含义
Bit 0	上一历元到当前历元可能失锁
Bit 1	当前历元可能存在半周模糊度
Bit 2	MBOC 信号采用特定 BOC 跟踪，噪声可能增大

0 或空白表示“正常或状态未知”，不等于绝对没有周跳。

SSI 是 1～9 的压缩信号强度指标；当原始强度以 dB-Hz 给出时，RINEX 规定了标准映射。

Sxx 观测则可保存更详细的信号强度，并由头文件说明单位。

因此：

SSI
→ 适合紧凑质量标志

Sxx或接收机原始C/N₀
→ 更适合定量权模型

15. 卫星健康、姿态和配置事件

应拒绝或标记：

不健康卫星；
无有效星历或产品的卫星；
产品定义与信号不匹配的卫星；
卫星机动、食或姿态异常；
接收机重启；
天线、射频链路或频标切换；
固件、跟踪模式或产品源切换。

这些事件可能造成：

码或相位偏差变化；
相位缠绕不连续；
系统间偏差变化；
模糊度整数基准变化。

16. 实时产品的三类时间指标

16.1 网络传输延迟

$τnet=trecv−tgen\boxed{ \tau_{\mathrm{net}} = t_{\mathrm{recv}} - t_{\mathrm{gen}} }$

16.2 改正外推时间

$τeval=tobs−t0\boxed{ \tau_{\mathrm{eval}} = t_{\mathrm{obs}} - t_0 }$

16.3 数据流陈旧时间

$τstale=tnow−tlast\boxed{ \tau_{\mathrm{stale}} = t_{\mathrm{now}} - t_{\mathrm{last}} }$

式中：

t_gen：产品生成时刻；
t_recv：用户接收时刻；
t_0：改正参考时刻；
t_obs：观测时刻；
t_last：最后一次有效产品到达时刻；
t_now：当前时刻。

三者不能混成一个“产品龄期”。

还应检查：

SSR Provider ID 和 Solution ID；
IOD SSR；
轨道改正对应的广播星历 IOD；
参考点是 APC 还是 CoM；
信号和跟踪模式；
相位偏差不连续计数器；
整数属性和 MW 一致性标志。

IOD 或产品版本变化本身不等于模糊度基准一定跳变。只有发现产品不兼容、相位偏差不连续或模型基准变化时，才需要转换或重置状态。

第三部分：高度角、C/N₀ 和随机模型

17. 高度角筛选

硬截止条件为：

$Ei≥Emin⁡\boxed{ E_i \ge E_{\min} }$

其中：

E_i：卫星高度角；
E_min：截止高度角。

建议：

极低高度角
→ 剔除

中低高度角
→ 降权

高高度角
→ 正常使用

截止角过高会减少卫星并恶化几何。

18. C/N₀、SNR 和 SSI 的区别

指标	常见形式	含义
C/N₀	dB-Hz	载波功率与噪声功率谱密度之比
SNR	dB 或线性比	指定带宽内的信号与噪声功率比
SSI	`1～9`	压缩信号强度指标

算法使用前必须确认字段的物理含义和单位。

高 C/N₀ 不等于一定没有多路径；低 C/N₀ 也不等于观测完全无用。

19. C/N₀ 三段式质量控制

设硬剔除和软降权阈值分别为：

$γdrop\gamma_{\mathrm{drop}}$

$γsoft\gamma_{\mathrm{soft}}$

处理规则为：

${Ci<γdrop剔除γdrop≤Ci<γsoft降权Ci≥γsoft正常使用\boxed{ \begin{cases} C_i < \gamma_{\mathrm{drop}} & \text{剔除} \\ \gamma_{\mathrm{drop}} \le C_i < \gamma_{\mathrm{soft}} & \text{降权} \\ C_i \ge \gamma_{\mathrm{soft}} & \text{正常使用} \end{cases} }$

式中：

C_i：当前观测的 C/N₀；
γ_drop：硬剔除阈值；
γ_soft：软降权阈值。

不存在适用于所有接收机和信号的统一阈值。

应使用：

历史残差统计；
接收机和天线类型；
信号频率；
静态或动态环境；

进行标定，并采用滑动平均或滞回避免频繁切换。

20. 高度角与 C/N₀ 权模型

一种经验模型为：

$σi2=a2+b2sin⁡pEi10max⁡(Cref−Ci,0)10\boxed{ \sigma_i^2 = a^2 + \frac{ b^2 }{ \sin^p E_i } 10^{ \frac{ \max \left( C_{\mathrm{ref}}-C_i, 0 \right) }{ 10 } } }$

式中：

σ_i²：观测方差；
E_i：高度角；
C_i：C/N₀；
C_ref：参考 C/N₀；
a、b、p：经验参数。

注意：

该式不是统一标准；
码、相位和多普勒应分别标定；
不同系统和信号可使用不同参数；
高度角和 C/N₀ 相关，可能重复计权；
应设置方差上下限。

21. 低成本和手机接收机

低成本设备常见：

相位中断频繁；
C/N₀ 波动大；
码噪声和多路径大；
时钟不连续；
半周状态变化；
占空比运行；
天线方向快速变化。

因此不能直接套用测地型接收机的固定阈值和随机模型。

第四部分：接收机钟跳与公共异常

22. 接收机钟跳的类型

钟跳可能表现为：

时间标签跳变；
多颗卫星伪距公共跳变；
多颗卫星相位公共跳变；
码和相位跳变量不一致；
整毫秒或其他离散调整。

时间标签跳变还会影响：

卫星位置计算时刻；
地球自转改正；
产品插值；
多普勒积分。

因此应先确认时间标签，再估计观测域公共跳变量。

23. 单系统公共模式估计

定义第 i 颗卫星的历元间异常量：

$di=ΔPi−Δρ^i−ΔT^i−ΔI^id_i = \Delta P_i - \Delta\widehat{\rho}_i - \Delta\widehat T_i - \Delta\widehat I_i$

若存在接收机公共钟跳：

$d=1Jc+e\boxed{ \mathbf d = \mathbf1 J_c + \mathbf e }$

加权估计为：

$J^c=1TWd1TW1\boxed{ \widehat J_c = \frac{ \mathbf1^{\mathrm T} \mathbf W \mathbf d }{ \mathbf1^{\mathrm T} \mathbf W \mathbf1 } }$

残差为：

$rc=d−1J^c\boxed{ \mathbf r_c = \mathbf d - \mathbf1\widehat J_c }$

其中：

J_c：以米表示的公共钟跳；
W：权阵；
e：剩余误差。

24. 多系统公共钟跳与 ISB 跳变

以 GPS 为参考系统，定义：

$j=[JcomJRJEJC]T\mathbf j = \begin{bmatrix} J_{\mathrm{com}} & J_R & J_E & J_C \end{bmatrix}^{\mathrm T}$

各系统对应的设计矩阵行可写为：

$aG=[1000]\mathbf a_G = \begin{bmatrix} 1&0&0&0 \end{bmatrix}$

$aR=[1100]\mathbf a_R = \begin{bmatrix} 1&1&0&0 \end{bmatrix}$

$aE=[1010]\mathbf a_E = \begin{bmatrix} 1&0&1&0 \end{bmatrix}$

$aC=[1001]\mathbf a_C = \begin{bmatrix} 1&0&0&1 \end{bmatrix}$

组合成：

$d=Acj+e\boxed{ \mathbf d = \mathbf A_c \mathbf j + \mathbf e }$

加权估计为：

$j^=(AcTWAc)−1AcTWd\boxed{ \widehat{\mathbf j} = \left( \mathbf A_c^{\mathrm T} \mathbf W \mathbf A_c \right)^{-1} \mathbf A_c^{\mathrm T} \mathbf W \mathbf d }$

式中：

J_com：各系统共同的接收机钟跳；
J_R、J_E、J_C：相对 GPS 的附加系统跳变；
A_c：系统归属设计矩阵。

只有被观测的系统才能估计对应跳变；系统缺失时应移除相应参数。

25. 接收机钟跳与产品跳变的区别

接收机异常通常呈现同一接收机的多卫星公共模式。

产品异常可能表现为：

单颗卫星跳变；
单个星座整体跳变；
产品源切换时多个卫星同时变化；
相位偏差不连续计数器改变；
Provider ID 或 Solution ID 改变。

因此应同时检查观测、产品元数据和其他接收机结果。

第五部分：周跳和相位弧段

26. 周跳的本质

连续跟踪时：

$N_{f,k}^{s} = N_{f,k-1}^{s}$

发生周跳后：

$Nf,ks=Nf,k−1s+ΔNf,ks\boxed{ N_{f,k}^{s} = N_{f,k-1}^{s} + \Delta N_{f,k}^{s} }$

周跳后必须：

终止旧弧段；
建立新弧段；
重置对应模糊度；
重置 GF、MW 和平滑统计；
暂停该弧段的整数固定。

27. 几何无关组合 GF

双频相位以米为单位：

$LGF=L1−L2\boxed{ L_{\mathrm{GF}} = L_1-L_2 }$

构造四维观测向量：

$yGF,k=[L1,kL2,kL1,k−1L2,k−1]T\mathbf y_{\mathrm{GF},k} = \begin{bmatrix} L_{1,k} & L_{2,k} & L_{1,k-1} & L_{2,k-1} \end{bmatrix}^{\mathrm T}$

系数向量为：

$tGF=[1−1−11]T\mathbf t_{\mathrm{GF}} = \begin{bmatrix} 1&-1&-1&1 \end{bmatrix}^{\mathrm T}$

历元差分为：

$ΔLGF,k=tGFTyGF,k\boxed{ \Delta L_{\mathrm{GF},k} = \mathbf t_{\mathrm{GF}}^{\mathrm T} \mathbf y_{\mathrm{GF},k} }$

其一般方差为：

$σΔGF,k2=tGFTQGF,ktGF\boxed{ \sigma_{\Delta\mathrm{GF},k}^{2} = \mathbf t_{\mathrm{GF}}^{\mathrm T} \mathbf Q_{\mathrm{GF},k} \mathbf t_{\mathrm{GF}} }$

其中：

Q_GF,k：包含频间和跨历元相关性的协方差。

判据可写为：

$∣ΔLGF,k−ΔI^GF,k∣>kGFσΔGF,k\boxed{ \left\vert{} \Delta L_{\mathrm{GF},k} - \widehat{\Delta I}_{\mathrm{GF},k} \right\vert{} > k_{\mathrm{GF}} \sigma_{\Delta\mathrm{GF},k} }$

GF 对电离层快速变化和特定周跳组合存在盲区。

28. Melbourne–Wübbena 组合

假设：

$f_1>f_2$

且 L₁、L₂、P₁、P₂ 均以米为单位。

宽巷相位组合为：

$LW=f1L1−f2L2f1−f2L_W = \frac{ f_1L_1-f_2L_2 }{ f_1-f_2 }$

窄巷码组合为：

$PN=f1P1+f2P2f1+f2P_N = \frac{ f_1P_1+f_2P_2 }{ f_1+f_2 }$

MW 组合为：

$MW=LW−PNλW\boxed{ MW = \frac{ L_W-P_N }{ \lambda_W } }$

其中：

$λW=cf1−f2\lambda_W = \frac{ c }{ f_1-f_2 }$

将 MW 写成线性形式：

$\mathbf t_{\mathrm{MW}}^{\mathrm T} \mathbf y_{\mathrm{MW}}$

则方差为：

$σMW2=tMWTRMWtMW\boxed{ \sigma_{\mathrm{MW}}^{2} = \mathbf t_{\mathrm{MW}}^{\mathrm T} \mathbf R_{\mathrm{MW}} \mathbf t_{\mathrm{MW}} }$

MW 受伪距噪声和多路径影响较大，应采用滑动统计和动态阈值。

29. 多普勒辅助相位连续性

设载波相位以周为单位，多普勒为 Hz。

定义：

$yD,k=[ΦkΦk−1DkDk−1]T\mathbf y_{D,k} = \begin{bmatrix} \Phi_k & \Phi_{k-1} & D_k & D_{k-1} \end{bmatrix}^{\mathrm T}$

系数向量为：

$tD=[1−1−sDΔtk2−sDΔtk2]T\mathbf t_D = \begin{bmatrix} 1 & -1 & -\dfrac{s_D\Delta t_k}{2} & -\dfrac{s_D\Delta t_k}{2} \end{bmatrix}^{\mathrm T}$

相位连续性残差为：

$rD,k=tDTyD,k\boxed{ r_{D,k} = \mathbf t_D^{\mathrm T} \mathbf y_{D,k} }$

方差为：

$σrD,k2=tDTQD,ktD+σmodel2\boxed{ \sigma_{r_D,k}^{2} = \mathbf t_D^{\mathrm T} \mathbf Q_{D,k} \mathbf t_D + \sigma_{\mathrm{model}}^{2} }$

式中：

s_D：接收机多普勒符号系数；
Δt_k：历元间隔；
Q_D,k：相位和多普勒的联合协方差；
σ_model²：动态和积分模型误差。

多普勒符号必须用实际设备数据验证。

30. 码多路径组合

双频码多路径组合可写为：

$MP1=P1−f12+f22f12−f22L1+2f22f12−f22L2\boxed{ MP_1 = P_1 - \frac{ f_1^2+f_2^2 }{ f_1^2-f_2^2 } L_1 + \frac{ 2f_2^2 }{ f_1^2-f_2^2 } L_2 }$

$MP2=P2−2f12f12−f22L1+f12+f22f12−f22L2\boxed{ MP_2 = P_2 - \frac{ 2f_1^2 }{ f_1^2-f_2^2 } L_1 + \frac{ f_1^2+f_2^2 }{ f_1^2-f_2^2 } L_2 }$

式中全部观测以米为单位。

组合中仍包含常数模糊度和硬件偏差，因此通常分析：

弧段内去均值序列；
历元差分；
滑动 RMS；
周跳前后的突变。

它适合评估码多路径，不应直接当成纯多路径真值。

31. 联合周跳判决

推荐联合使用：

LLI和锁定时间
+
GF
+
MW
+
多普勒
+
滤波创新

单一方法存在盲区，多个指标也并非统计独立。

因此应采用：

分层判决；
状态机；
联合协方差；
真实数据标定；

而不是简单地“任一超限即周跳”或“多数投票即正确”。

32. 周跳探测与修复

实时 PPP 的默认安全策略是：

高置信检测到周跳
→ 新建弧段
→ 重置模糊度

只有在下列条件满足时才考虑周跳修复：

多频冗余充分；
整数周跳解唯一；
通过独立整数验证；
修复后的创新和残差正常；
错误修复风险可接受。

第六部分：粗差、创新和鲁棒估计

33. 预更新创新检验

预测创新为：

$νk=zk−h(x^k−)\boxed{ \boldsymbol{\nu}_k = \mathbf z_k - \mathbf h \left( \widehat{\mathbf x}_k^{-} \right) }$

创新协方差为：

$Sk=HkPk−HkT+Rk\boxed{ \mathbf S_k = \mathbf H_k \mathbf P_k^{-} \mathbf H_k^{\mathrm T} + \mathbf R_k }$

粗差检验应优先在量测更新之前执行。

34. 全局卡方检验

对创新协方差作 Cholesky 分解：

$Sk=LkLkT\mathbf S_k = \mathbf L_k \mathbf L_k^{\mathrm T}$

白化创新为：

$ek=Lk−1νk\mathbf e_k = \mathbf L_k^{-1} \boldsymbol{\nu}_k$

全局统计量为：

$Tk=νkTSk−1νk\boxed{ T_k = \boldsymbol{\nu}_k^{\mathrm T} \mathbf S_k^{-1} \boldsymbol{\nu}_k }$

在近似高斯、模型正确且 S_k 非奇异时：

$Tk∼χm2T_k \sim \chi_m^2$

接受条件为：

$Tk≤χm,1−α2\boxed{ T_k \le \chi_{m,1-\alpha}^{2} }$

其中：

m：创新维数；
α：显著性水平。

35. 局部异常检验

假设第 i 条观测存在附加粗差，定义单位向量：

$ci\mathbf c_i$

局部统计量为：

$wi=ciTSk−1νkciTSk−1ci\boxed{ w_i = \frac{ \mathbf c_i^{\mathrm T} \mathbf S_k^{-1} \boldsymbol{\nu}_k }{ \sqrt{ \mathbf c_i^{\mathrm T} \mathbf S_k^{-1} \mathbf c_i } } }$

在零假设下：

$wi∼N(0,1)w_i \sim \mathcal N \left( 0,1 \right)$

同时检验多条观测时，应考虑多重检验和观测相关性。

36. 全协方差条件下的鲁棒降权

设：

$R=LRLRT\mathbf R = \mathbf L_R \mathbf L_R^{\mathrm T}$

对观测白化后，定义鲁棒因子矩阵：

$A=diag⁡(α1,…,αm)\mathbf A = \operatorname{diag} \left( \alpha_1, \ldots, \alpha_m \right)$

其中：

$\alpha_i \le 1$

Huber 型因子为：

$αi={1,∣ui∣≤cHcH∣ui∣,∣ui∣>cH\boxed{ \alpha_i = \begin{cases} 1, & \left\vert{} u_i \right\vert{} \le c_H \\ \dfrac{ c_H }{ \left\vert{} u_i \right\vert{} }, & \left\vert{} u_i \right\vert{} > c_H \end{cases} }$

等效协方差为：

$R~=LRA−1LRT\boxed{ \widetilde{\mathbf R} = \mathbf L_R \mathbf A^{-1} \mathbf L_R^{\mathrm T} }$

鲁棒降权不能代替：

周跳后的状态重置；
产品基准跳变处理；
失锁后的新弧段。

37. 剔除、降权和重置的边界

异常	推荐动作
极低 C/N₀	剔除或停用
一般低高度角或低 C/N₀	降权
单历元伪距粗差	剔除或鲁棒降权
相位周跳	重置模糊度和弧段
接收机钟跳	修复公共跳变或调整钟差状态
产品超龄	拒绝、降级或增大产品方差
相位偏差不连续	转换或重置相关模糊度
数据间断	无法证明连续时重建弧段

第七部分：实时状态管理与 PPP-AR

38. PPP-AR 的额外检查

PPP-AR 还要检查：

相位偏差是否覆盖当前信号；
Signal Integer Indicator；
Wide-Lane Integer Indicator；
Signal Discontinuity Counter；
MW Consistency Indicator；
Provider ID、Solution ID 和产品源；
码、相位和钟差基准是否一致；
GLONASS FDMA 接收机相关偏差；
产品切换是否改变整数基准。

IGS SSR 允许同一信号的卫星偏差具有公共时变偏移。用户不能仅凭偏差数值变化就判断周跳，应结合不连续计数器和产品基准处理。

39. 实时处理的因果性和状态机

实时历元 k 只能使用：

$Dk={z1,…,zk}\boxed{ \mathcal D_k = \left\{ \mathbf z_1, \ldots, \mathbf z_k \right\} }$

不能使用未来观测。

状态机应按“卫星 + 具体信号”维护。

40. 后验检查和位置域影响

后验残差为：

$rk+=zk−h(x^k+)\boxed{ \mathbf r_k^{+} = \mathbf z_k - \mathbf h \left( \widehat{\mathbf x}_k^{+} \right) }$

后验残差适合：

长期质量统计；
多路径识别；
产品异常监测；
固定解释放判断。

但不能替代预更新创新检验，因为异常可能已被状态吸收。

还应监测：

坐标跳变；
钟差和 ISB 跳变；
对流层突变；
模糊度重置数量；
固定解与浮点解差异；
独立传感器一致性。

第八部分：工程实施

41. 推荐处理顺序

顺序很重要：

先排查时间和产品公共异常
再判断多颗卫星是否同时周跳

42. 阈值标定原则

动态阈值可概念化为：

$τi,k=kτσi,k2+σmodel,i,k2\boxed{ \tau_{i,k} = k_{\tau} \sqrt{ \sigma_{i,k}^{2} + \sigma_{\mathrm{model},i,k}^{2} } }$

式中：

τ_i,k：当前阈值；
k_τ：检验倍数；
σ_i,k²：观测噪声方差；
σ_model²：模型误差方差。

阈值应根据：

采样率；
接收机和天线；
信号类型；
静态或动态场景；
电离层活动；
历史误报率和漏检率；

进行标定。

43. 推荐检查清单

43.1 历元级

时间连续性；
数据乱序、重复和缺失；
可用系统和卫星数；
接收机公共钟跳；
产品延迟和陈旧状态；
全局创新检验；
输出延迟。

43.2 卫星级

健康和姿态；
轨道钟差和 IOD；
高度角和遮挡；
产品连续性；
单星粗差；
卫星退出和重新进入。

43.3 信号级

单位和观测类型；
C/N₀、SSI、LLI 和锁定时间；
半周状态；
GF、MW 和多普勒指标；
码多路径；
弧段编号；
模糊度重置标志。

44. 监控与日志

建议记录：

原始观测和时间；
卫星、信号、高度角和 C/N₀；
LLI、半周和锁定时间；
GF、MW、多普勒和创新统计量；
原始方差、产品方差和鲁棒因子；
Provider ID、Solution ID、IOD 和产品源；
弧段编号和重置原因；
剔除、降权和降级原因；
状态和坐标跳变量。

45. 常见概念混淆

45.1 预处理就是周跳探测

错误。周跳只是预处理的一部分。

45.2 低 C/N₀ 必须一律删除

错误。极低信号可删除，一般低质量更适合降权。

45.3 高 C/N₀ 就没有多路径

错误。强反射信号也可能具有高 C/N₀。

45.4 SSI 等于精确 C/N₀

错误。SSI 是压缩等级；定量模型应优先使用 Sxx 或原始 C/N₀。

45.5 LLI 为零就证明没有周跳

错误。零值也可能表示状态未知。

45.6 数据间断一定发生周跳

不一定。但无法证明连续性时，应建立新弧段。

45.7 删除周跳历元后可继续使用旧模糊度

错误。相位连续性已破坏，必须处理状态和协方差。

45.8 IOD 变化就一定要重置模糊度

错误。应检查产品兼容性和相位偏差不连续信息。

45.9 鲁棒估计可以替代周跳处理

错误。降权不能恢复模糊度连续性。

45.10 观测越多越好

错误。错误观测、强相关观测和失配权模型可能降低可靠性。

45.11 后验残差小就证明没有粗差

错误。粗差可能已经被状态吸收。

45.12 实时改正可以当作无误差常数

错误。改正的不确定度和相关性也应进入随机模型。

45.13 Hatch 平滑是 PPP 必需步骤

错误。它是可选辅助工具，不是所有 PPP 的标准步骤。

46. 总结

实时 PPP 预处理的本质是：

在观测进入具有记忆性的滤波器之前，保证数据、产品、时间、信号、整数基准和随机模型相互一致。

除了周跳探测，还应处理：

数据装配、乱序和时间异常；
单位、频率和信号映射；
卫星健康、姿态和配置事件；
实时轨道、钟差、码偏差和相位偏差；
产品延迟、外推时间和陈旧状态；
C/N₀、SSI、LLI、锁定时间和半周状态；
接收机钟跳、多系统 ISB 和产品公共跳变；
GF、MW、多普勒和码多路径组合；
创新全局检验和局部检验；
全协方差条件下的鲁棒降权；
状态可观测性检查；
相位弧段、模糊度和协方差重置；
PPP-AR 整数属性和相位偏差连续性；
后验质量和位置域监控。

最重要的处理原则是：

质量下降但仍有信息
→ 降权

明显错误
→ 剔除

相位连续性或整数基准破坏
→ 重置状态

多卫星公共异常
→ 优先检查时间 钟跳和产品

单星单信号局部异常
→ 优先检查周跳 粗差和遮挡

参考资料

IGS/RTCM RINEX Committee：RINEX Version 4.02。
International GNSS Service：IGS SSR Format Version 1.00。
International GNSS Service：Real-Time Service Formats and Products。
G. Blewitt：An Automatic Editing Algorithm for GPS Data。
F. Guo 等：An Efficient Method to Compensate Receiver Clock Jumps in Real-Time Precise Point Positioning。
N. Jiang 等：Detecting and Repairing Inter-System Bias Jumps with Satellite Clock Preprocessing。
B. Li、L. Lou、Y. Shen：GNSS Elevation-Dependent Stochastic Modeling and Its Impacts on Statistic Testing。
P. D. Groves、Z. Jiang：GNSS NLOS and Multipath Detection Using C/N₀ and Consistency Checking。
J. Kouba：A Guide to Using International GNSS Service Products。
P. J. G. Teunissen、O. Montenbruck：Springer Handbook of Global Navigation Satellite Systems。