2. GNSS基础

GNSS(全球导航卫星系统)是INS/GNSS紧耦合组合导航的“外部量测源”。要理解紧耦合(Tightly Coupled)为何能比松耦合(Loosely Coupled)在弱信号、高动态场景下表现更优,首先必须透彻理解四大卫星导航系统(GPS、GLONASS、北斗、Galileo)的星座设计、信号体制、观测值生成机理以及误差传播特性。紧耦合算法的本质,是在伪距/载波相位/多普勒层面直接融合INS预测值与GNSS原始观测量,因此本章内容将直接决定后续对“量测方程构建”“钟差估计”“周跳处理”等核心模块的理解深度。

GNSS四大系统观测值生成与误差传播链路 GPS (美国) L1/L2/L5, CDMA 31颗MEO, 24h周期 WGS-84坐标系 GLONASS (俄罗斯) G1/G2/G3, FDMA 24颗MEO, 11.25h周期 PZ-90.11坐标系 北斗 (中国) B1I/B3I/B1C/B2a, CDMA 30颗MEO+5GEO+3IGSO CGCS2000坐标系 Galileo (欧盟) E1/E5a/E5b/E6, CDMA 24颗MEO+6备份 GTRF坐标系 信号结构三要素 载波 (L1=1575.42MHz) 测距码 (C/A码, P码, B1I码) 导航电文 (星历, 钟差, 电离层参数) 三类观测值生成机理 伪距 ρ = r + c·(δt_u - δt_s) + I + T + ε 码相位测量, 精度~0.3-3m 载波相位 φ = λ⁻¹(r + c·δt + I - T) + N + ε 载波相位测量, 精度~0.5-5mm, 含整周模糊度N 多普勒 f_d = -f_r·(v·e)/c 载波频率变化率, 测速精度~0.01-0.1m/s 误差源四大类 卫星端 星历误差(1-2m), 钟差(1-5ns) 传播路径 电离层(2-10m), 对流层(0.2-2m) 接收机端 钟差, 噪声, 多路径(0.5-5m) 环境干扰 多路径, 电磁干扰, 遮挡

上图展示了GNSS紧耦合处理中必须关注的四大系统、信号结构、观测值生成与误差传播的完整链路。在后续章节中,我们将逐一深入每个环节,并重点讨论这些误差如何通过差分、滤波或约束在紧耦合框架中被消除或削弱。

2.1 四大卫星导航系统核心原理与对比

核心原理:所有GNSS系统均基于“时间测距”原理——卫星发射携带时间戳的无线电信号,接收机测量信号传播时间Δt,乘以光速c得到卫星与接收机之间的几何距离ρ。由于接收机时钟与卫星时钟存在偏差δt_u,实际测得的伪距ρ = r + c·(δt_u - δt_s) + 误差项。四大系统在星座设计、频率规划、多址方式上各有差异,但基本原理一致。

系统星座构成载波频率多址方式坐标框架时间基准
GPS31颗MEO (6轨道面, 55°倾角)L1:1575.42MHz, L2:1227.60MHz, L5:1176.45MHzCDMAWGS-84GPST (UTC-18s)
GLONASS24颗MEO (3轨道面, 64.8°倾角)G1:1602.0+k×0.5625MHz, G2:1246.0+k×0.4375MHzFDMA (k=-7~+6)PZ-90.11GLONASST (UTC+3h)
北斗30颗MEO+5GEO+3IGSOB1I:1561.098MHz, B3I:1268.52MHz, B1C:1575.42MHz, B2a:1176.45MHzCDMACGCS2000BDT (UTC-4s)
Galileo24颗MEO+6备份 (3轨道面, 56°倾角)E1:1575.42MHz, E5a:1176.45MHz, E5b:1207.14MHz, E6:1278.75MHzCDMAGTRFGST (UTC-13s)

判断方法:在紧耦合算法设计中,选择哪些系统参与融合,需考虑:①可见卫星数(城市峡谷中北斗GEO卫星具有高仰角优势);②频率数量(双频可消除电离层一阶项);③系统间时间偏差(需额外估计系统间钟差参数)。

实战经验:在车载组合导航中,建议至少同时跟踪GPS+北斗双系统。北斗的GEO卫星在低纬度地区仰角稳定,能显著改善城市峡谷中的几何精度因子(GDOP)。但需注意,北斗GEO卫星的轨道高度约35786km,信号延迟比MEO卫星大,在紧耦合滤波器中需正确设置动态模型噪声。

2.2 信号结构与导航电文解码

核心原理:GNSS信号由三部分组成:载波(L-band正弦波)、测距码(伪随机噪声码PRN,用于码相位测量)和导航电文(包含卫星星历、钟差参数、电离层模型系数、卫星健康状态等)。以GPS L1信号为例,C/A码码率为1.023Mbps,码周期1ms,码长1023bit;P码码率10.23Mbps,周期7天。北斗B1I码码率为2.046Mbps,码周期1ms,码长2046bit。

操作流程:接收机信号处理流程为:①射频前端下变频至中频;②ADC采样(通常2-4bit量化);③载波环(Costas环或FLL)剥离载波;④码环(DLL)跟踪测距码;⑤帧同步后解码导航电文;⑥提取星历参数计算卫星位置和速度。

// 导航电文解码关键参数示例 (GPS L1 C/A)
// 星历参数: 开普勒轨道根数 + 摄动修正项
double sqrtA = 5153.65531;        // 半长轴平方根 (m^0.5)
double e = 0.002345678;           // 偏心率
double i0 = 0.959931;             // 轨道倾角 (rad)
double Omega0 = -2.567890;        // 升交点赤经 (rad)
double M0 = -1.234567;            // 平近点角 (rad)
double delta_n = 4.5e-9;          // 平均角速度修正 (rad/s)
double Cuc = 1.2e-6, Cus = 8.5e-6; // 纬度幅角修正
double Crc = 2.1e-6, Crs = 7.3e-6; // 轨道半径修正
double Cic = 1.5e-6, Cis = 6.2e-6; // 轨道倾角修正
// 卫星钟差参数: af0, af1, af2 (s, s/s, s/s^2)
double af0 = 2.3e-4, af1 = 1.2e-11, af2 = 0.0;

案例拆解:某车载测试中,接收机在隧道出口重新捕获信号时,发现GPS L1 C/A码跟踪环路失锁。分析日志发现,信号中断期间接收机本地振荡器频漂达±50Hz,导致重新捕获时多普勒搜索范围不足。解决方案:在紧耦合滤波器中,利用INS速度预测多普勒频移,将搜索范围从±5kHz缩小至±500Hz,捕获时间从8s缩短至1.2s。

常见误区:初学者常混淆“码相位测量”与“载波相位测量”的精度差异。码相位测量精度受码元宽度限制(C/A码码元宽度约293m,理论测距精度约2.9m),而载波相位测量精度可达毫米级(L1波长19.03cm,相位测量精度约0.5-2mm)。但载波相位存在整周模糊度问题,必须通过初始化或搜索算法求解。

2.3 观测值模型:伪距、载波相位与多普勒

核心原理:三种观测值分别对应不同的测量机理和应用场景:

  • 伪距观测值:ρ = r + c·(δt_u - δt_s) + I + T + ε_ρ,其中I为电离层延迟(与频率平方成反比),T为对流层延迟(与频率无关),ε_ρ为码噪声和多路径误差。伪距精度较低但无模糊度问题,适用于定位。
  • 载波相位观测值:φ = λ⁻¹·(r + c·(δt_u - δt_s) - I + T) + N + ε_φ,其中N为整周模糊度(整数),λ为载波波长。载波相位精度高但需要解算模糊度,适用于高精度相对定位和精密单点定位(PPP)。
  • 多普勒观测值:f_d = -f_r·(v·e)/c + δf_u - δf_s,其中v为接收机速度,e为卫星方向单位矢量,δf_u和δf_s分别为接收机和卫星钟漂。多普勒观测值直接反映速度信息,在紧耦合中用于约束INS速度误差。

判断方法:在紧耦合量测更新中,应根据应用场景选择观测值组合:①城市峡谷中,伪距+多普勒组合可保证鲁棒性;②开阔环境下,加入载波相位平滑伪距可提升精度;③高精度测绘场景,需全载波相位解算模糊度。

观测值类型精度模糊度主要误差紧耦合应用
伪距0.3-3m无电离层、多路径位置约束
载波相位0.5-5mm有(整数)周跳、模糊度解算高精度位置约束
多普勒0.01-0.1m/s无钟漂、动态应力速度约束

案例拆解:某无人机紧耦合导航系统在飞行测试中,GPS信号被螺旋桨周期性遮挡,导致伪距出现周期性跳变。通过分析多普勒观测值,发现跳变时刻多普勒频率出现±3Hz的异常波动。进一步检查发现,遮挡导致载波环失锁,码环跟踪延迟增大。解决方案:在紧耦合滤波器中增加多普勒一致性检验,当伪距残差超过阈值且多普勒异常时,自动降低该卫星权重。

关键公式:双频电离层消除组合(以GPS L1/L2为例):
ρ_IF = (f1²·ρ1 - f2²·ρ2)/(f1² - f2²) = r + c·δt + T + ε_IF
该组合消除了电离层一阶项,但放大了码噪声约3倍。在紧耦合中,若使用双频接收机,建议在滤波器中直接使用消电离层组合,同时将观测噪声方差乘以3²=9倍。

2.4 误差源分类与抑制策略

核心原理:GNSS误差源可分为四大类:卫星端误差、传播路径误差、接收机端误差和环境干扰。理解各类误差的时空相关特性,是设计紧耦合滤波器状态向量和过程噪声矩阵的基础。

  1. 卫星端误差:星历误差(广播星历约1-2m,精密星历约2-5cm)、卫星钟差(广播钟差约1-5ns,精密钟差约0.1-0.3ns)。在紧耦合中,通常将卫星钟差视为已知(通过导航电文修正),剩余残差并入观测噪声。
  2. 传播路径误差:电离层延迟(白天可达10-30m,夜间2-10m,与频率平方成反比)、对流层延迟(天顶方向约2.3m,与频率无关,分为干分量和湿分量)。双频接收机可消除电离层一阶项,对流层湿分量需通过模型或参数估计。
  3. 接收机端误差:接收机钟差(需在滤波器中估计)、天线相位中心偏差(可标定)、码噪声和多路径(与环境和天线设计相关)。
  4. 环境干扰:多路径效应(城市峡谷中可达5-15m)、电磁干扰、信号遮挡。多路径误差具有非高斯、时间相关特性,是紧耦合中最难处理的误差源。

操作流程:误差抑制策略优先级:①双频消电离层组合(消除电离层一阶项);②对流层模型修正(Saastamoinen或UNB3模型)+ 湿分量参数估计;③载波相位平滑伪距(降低码噪声);④基于INS预测的多路径检测与剔除(利用INS短期高精度特性)。

风险提示:在紧耦合滤波器中,若将多路径误差简单地建模为白噪声,会导致滤波器过于自信(协方差过小),进而产生滤波发散。建议将多路径误差建模为一阶高斯马尔可夫过程,相关时间常数设为10-60s,并在城市峡谷场景中适当增大观测噪声方差。

案例拆解:某车载组合导航系统在高架桥下测试时,GPS伪距出现持续5s的+8m偏差。分析发现,高架桥的金属结构反射信号形成强多路径。通过对比INS预测位置与GPS伪距残差,发现残差呈现明显的非零均值特性。解决方案:在紧耦合滤波器中增加多路径检测模块,利用INS预测位置计算伪距残差,当残差超过3倍标准差时,将该卫星从量测更新中剔除。同时,将多路径误差建模为随机游走过程,在滤波器中实时估计。

总结复盘:GNSS误差处理的核心原则是“分类处理、模型匹配”。卫星端误差通过差分或精密产品消除;传播路径误差通过双频组合和模型修正;接收机端误差通过参数估计;环境干扰通过INS辅助检测和鲁棒滤波。在紧耦合框架中,INS提供了连续、短时高精度的位置/速度预测,这为GNSS误差的实时检测和剔除提供了独特的参考基准——这正是紧耦合相比松耦合的核心优势之一。