2. 全球导航卫星系统基础
RTK（Real-Time Kinematic，实时动态载波相位差分）定位技术的核心，是充分利用全球导航卫星系统（GNSS）提供的原始观测值进行高精度解算。要深入理解RTK算法，首先必须透彻掌握四大全球卫星导航系统（GPS、GLONASS、Galileo、BDS）的星座设计、信号体制以及三种基本观测值（伪距、载波相位、多普勒）的物理含义与数学表达。本章将系统梳理这些基础框架，为后续章节的差分模型、整周模糊度解算和误差处理奠定坚实的知识底座。
上图从宏观视角展示了本章的核心知识框架：左侧为四大系统的星座构型与信号频段，右侧为RTK算法所依赖的三种原始观测值。下面我们将逐一深入剖析每个模块的细节。
2.1 四大星座结构：从轨道设计到可用性
2.1.1 核心原理：轨道高度与几何构型
全球导航卫星系统的定位原理基于“距离交会”：接收机同时观测至少4颗卫星，通过解算卫星位置和信号传播时间来确定自身三维坐标与钟差。星座结构决定了卫星的几何分布，直接影响定位精度因子（DOP）和可用性。四大系统的轨道设计各有侧重：
- GPS（美国）：标称24颗工作卫星，分布于6个中地球轨道（MEO）平面，轨道倾角55°，轨道高度约20200 km，周期约11小时58分。这种设计保证了全球任何地点、任何时间至少可见4颗卫星，实际在轨运行通常超过30颗。
- GLONASS（俄罗斯）：24颗卫星均匀分布在3个轨道面内，轨道倾角64.8°，高度约19100 km。高倾角设计使其在高纬度地区（如俄罗斯本土）的覆盖性能优于GPS。
- Galileo（欧盟）：设计为24颗工作卫星加6颗备份，3个轨道面，倾角56°，高度约23222 km。其轨道高度最高，信号到达地面的路径损耗略大，但卫星运行更稳定。
- BDS（中国）：采用独特的“3+3”混合星座——3颗地球静止轨道（GEO，高度35786 km）、3颗倾斜地球同步轨道（IGSO，高度35786 km，倾角55°）以及24颗中地球轨道（MEO，高度21528 km）。GEO/IGSO卫星增强了亚太地区的服务性能。
知识结构表：四大系统星座参数对比
| 参数 | GPS | GLONASS | Galileo | BDS |
|---|---|---|---|---|
| 轨道面数 | 6 | 3 | 3 | 3类（GEO/IGSO/MEO） |
| 轨道倾角 | 55° | 64.8° | 56° | 55°（MEO/IGSO） |
| 轨道高度 | 20200 km | 19100 km | 23222 km | 21528 km（MEO） |
| 卫星数量 | 24+（实际30+） | 24 | 24+6 | 30+ |
| 多址方式 | CDMA | FDMA（L1/L2） | CDMA | CDMA |
| 坐标参考框架 | WGS-84 | PZ-90 | GTRF | CGCS2000 |
2.1.2 判断方法：如何评估星座可用性
在RTK作业中，我们通常使用“可见卫星数”和“PDOP值”两个指标来评估星座状态。PDOP（位置精度因子）反映了卫星几何构型对定位误差的放大程度，PDOP值越小越好（通常要求小于3）。实际操作中，可以通过接收机输出的NMEA语句中的GGA和GSV帧来实时监控：
// GGA语句示例（第7、8字段为可见卫星数和PDOP）
$GPGGA,123519,4807.038,N,01131.000,E,1,08,0.9,545.4,M,46.9,M,,*47
// 其中 08 表示可见卫星数，0.9 为PDOP值
当可见卫星数少于8颗或PDOP大于3时，RTK初始化时间会显著延长，甚至无法固定解。此时应检查周围是否有遮挡物（高楼、树木、山体），或考虑更换观测时段。
2.1.3 案例拆解：城市峡谷环境下的星座选择
场景：某测绘项目需要在上海陆家嘴金融区进行道路中线放样，两侧高楼密集，天空可见度不足40%。初始使用单GPS接收机，RTK固定率仅45%，频繁失锁。
分析过程：通过接收机日志发现，GPS卫星可见数在4~6颗之间波动，PDOP经常超过4。而BDS的GEO卫星（位于赤道上空固定位置）在上海地区仰角约45°，能够穿透部分楼宇间隙；Galileo卫星虽然数量少，但其轨道高度高，信号遮挡概率相对较低。
操作步骤：
- 将接收机设置为“GPS+BDS+Galileo”多系统联合解算模式。
- 在测区中央架设基准站，确保基准站天线视野开阔（选择楼顶平台）。
- 流动站作业时，保持天线竖直，避免身体和车辆遮挡。
结果对比：启用多系统后，可见卫星数提升至12~16颗，PDOP降至1.8~2.5，RTK固定率提升至92%，初始化时间从平均45秒缩短至8秒。
复盘经验：在复杂环境下，多系统融合是提升RTK可用性的最直接手段。但需注意，不同系统间存在系统间偏差（ISB），需要接收机固件或算法进行补偿。
2.1.4 常见误区与总结
误区1：认为卫星越多越好。实际上，卫星数量超过一定阈值后，对定位精度的提升趋于饱和，反而增加计算负担。关键在于卫星的几何分布是否均匀。
误区2：忽略GLONASS的FDMA特性。GLONASS L1/L2信号采用频分多址，每颗卫星使用不同的载波频率，导致其载波相位观测值的硬件延迟偏差与CDMA系统不同，在RTK解算中需要额外处理。
总结：星座结构决定了系统的“可用性天花板”。RTK算法设计必须考虑多系统融合时的系统间偏差、坐标框架统一和时差改正。
2.2 信号体制：从载波到测距码
2.2.1 核心原理：三层信号结构
GNSS信号由三层构成：载波（L波段正弦波）、测距码（伪随机噪声码，PRN）和导航电文（数据码）。测距码通过BPSK（二进制相移键控）或BOC（二进制偏移载波）调制在载波上。接收机通过复制本地码与接收码进行相关运算，得到信号传播时间，进而计算伪距。
- GPS：L1（1575.42 MHz）承载C/A码（粗捕获码，码率1.023 Mcps）和P(Y)码（精密码，10.23 Mcps）；L2（1227.60 MHz）主要承载P(Y)码；L5（1176.45 MHz）为新的民用信号，采用BPSK(10)调制，码率10.23 Mcps。
- GLONASS：L1（1602 MHz + k×0.5625 MHz，k为卫星频率编号）和L2（1246 MHz + k×0.4375 MHz），采用FDMA方式区分卫星。
- Galileo：E1（1575.42 MHz）采用CBOC（复合二进制偏移载波）调制，E5a/E5b（1176.45/1207.14 MHz）采用AltBOC（交替二进制偏移载波）调制，E6（1278.75 MHz）用于商业服务。
- BDS：B1I（1561.098 MHz）和B3I（1268.52 MHz）为老一代信号，B1C（1575.42 MHz）和B2a（1176.45 MHz）为新信号，采用BOC和QPSK调制。
经验提醒：RTK算法通常使用双频或三频观测值来消除电离层延迟。例如，GPS L1+L2双频消电离层组合，或L1+L2+L5三频组合。了解各系统的频率规划，是设计多频RTK解算策略的前提。
2.2.2 判断方法：信号质量评估指标
接收机输出的信噪比（SNR，单位dB-Hz）是衡量信号质量的核心指标。一般要求SNR > 35 dB-Hz才能保证稳定的载波相位跟踪。在RTK作业中，应关注以下阈值：
| SNR范围 (dB-Hz) | 信号状态 | RTK影响 |
|---|---|---|
| > 45 | 优秀 | 载波相位噪声极小，固定快 |
| 35 ~ 45 | 良好 | 正常作业，偶尔周跳 |
| 25 ~ 35 | 较差 | 周跳频繁，难以固定 |
| < 25 | 极差 | 无法跟踪载波，仅能使用伪距 |
2.2.3 实战场景：多路径效应对信号的影响
场景描述：在桥梁检测项目中，接收机架设在桥面护栏旁，下方为水面。观测数据中出现周期性SNR波动，载波相位观测值频繁发生周跳。
分析判断：水面反射GPS信号产生多路径效应，反射信号与直达信号叠加，导致SNR呈正弦波动。通过分析SNR残差序列，发现波动周期约10分钟，与桥梁固有振动周期无关，确认是水面反射所致。
解决方案：采用扼流圈天线（choke ring antenna）抑制低仰角反射信号，并将天线位置向桥面中央移动2米，避开水面反射区。处理后SNR波动幅度从±6 dB降至±1.5 dB，周跳数量减少90%。
2.2.4 常见误区与总结
误区：认为所有系统的同频段信号（如GPS L1和Galileo E1）可以完全等同处理。虽然频率相同，但码型、调制方式和电文格式不同，接收机需要分别跟踪，且系统间存在硬件延迟偏差。
总结：信号体制决定了观测值的精度上限。RTK算法需要针对不同信号设计对应的处理链路，包括码环、载波环和电文解码。
2.3 观测值类型：伪距、载波相位与多普勒
2.3.1 核心原理：三种观测值的物理模型
接收机输出的原始观测值包括三种，它们各自有不同的精度和应用场景：
- 伪距（P）：通过码相关测量得到，是卫星与接收机之间的几何距离加上各种误差项。精度为分米到米级，用于标准定位和RTK的初始解。
- 载波相位（φ）：通过载波跟踪环测量接收机复制的载波与接收到的卫星载波之间的相位差。精度为毫米级，但存在整周模糊度N（未知整数）。RTK的核心就是利用载波相位的高精度，并通过差分消除误差后固定N。
- 多普勒（f_d）：由于卫星与接收机之间的相对运动，接收到的载波频率发生偏移。多普勒观测值用于测速，也可用于检测周跳。
三种观测值的数学模型如下（以GPS L1为例）：
伪距: P = ρ + c·(δt_u - δt_s) + I + T + ε_P
载波相位: λ·φ = ρ + c·(δt_u - δt_s) - I + T + λ·N + ε_φ
多普勒: λ·f_d = v·e + c·(δf_u - δf_s) + ε_f
其中:
ρ = 卫星与接收机几何距离
δt_u = 接收机钟差, δt_s = 卫星钟差
I = 电离层延迟（伪距超前，载波滞后）
T = 对流层延迟
N = 整周模糊度（整数）
v = 卫星与接收机相对速度矢量
e = 视线方向单位矢量
观测值特性对比表
| 特性 | 伪距 | 载波相位 | 多普勒 |
|---|---|---|---|
| 测量精度 | 0.3~3 m | 0.5~5 mm | 0.1~1 mm/s |
| 模糊度 | 无 | 有（整数N） | 无 |
| 对电离层敏感度 | 延迟（+） | 超前（-） | 不敏感 |
| 主要用途 | 定位、授时 | 高精度定位（RTK/PPP） | 测速、周跳检测 |
| 周跳影响 | 无 | 导致N跳变 | 可检测周跳 |
2.3.2 判断方法：周跳检测与处理
周跳是载波相位观测值中整周计数的意外跳变，是RTK解算中最常见的故障源。检测周跳的常用方法包括：
- 多普勒积分法：将多普勒观测值积分，与载波相位变化量比较，差值超过阈值（如0.5周）则判定发生周跳。
- 电离层残差法：利用双频观测值计算电离层延迟变化率，若发生突变则可能存在周跳。
- M-W组合法：使用宽带组合（宽巷模糊度）和窄巷伪距组合，其线性组合值应保持稳定，若发生跳变则说明至少一个频率上发生周跳。
// M-W组合公式（无几何距离、无电离层影响）
MW = (f1·φ1 - f2·φ2)/(f1-f2) - (f1·P1 + f2·P2)/(f1+f2)
// 该值理论上为常数（宽巷模糊度），若发生跳变则存在周跳
2.3.3 案例拆解：RTK初始化中的模糊度固定
背景与目标：某变形监测项目要求水平精度±5mm，垂直精度±10mm。基准站与监测站距离约500m，使用双频GPS接收机。初始RTK解算时，模糊度固定成功率仅60%，固定后经常在几分钟内失锁。
关键数据：通过分析原始观测值，发现监测站周围有大型金属脚手架，导致多路径效应严重。载波相位观测值的SNR在35~40 dB-Hz之间波动，伪距噪声达到0.8m（正常应小于0.3m）。
分析判断过程：
- 首先检查卫星几何分布，PDOP为2.1，良好。
- 检查双频电离层残差，发现L1/L2电离层延迟变化率异常，波动幅度超过正常值3倍。
- 进一步检查M-W组合，发现宽巷模糊度在短时间内发生多次±1周的跳变。
实际操作步骤：
- 在监测站天线下方加装直径60cm的抑径板，吸收地面反射信号。
- 将接收机固件升级，启用多路径抑制算法（如Strobe Correlator）。
- 调整RTK解算策略，将高度截止角从10°提高到15°，减少低仰角多路径影响。
结果对比：处理后，SNR提升至42~48 dB-Hz，伪距噪声降至0.2m，模糊度固定成功率提升至95%，固定后连续观测24小时未发生失锁。
风险与代价：提高截止角会减少可见卫星数（从10颗降至8颗），但PDOP仅从2.1升至2.4，仍在可接受范围。抑径板增加了设备成本约800元。
复盘经验：周跳检测不能只看单一指标，应综合使用多普勒积分、电离层残差和M-W组合三种方法。在强多路径环境下，即使模糊度固定成功，也需要持续监控观测值质量。
2.3.4 常见误区与总结
误区1：将伪距和载波相位视为独立观测值。实际上，两者共享相同的几何路径和大部分误差源，RTK算法通过差分消除公共误差后，利用伪距辅助载波相位模糊度解算。
误区2：忽略多普勒观测值在RTK中的作用。多普勒不仅用于测速，更重要的是在动态RTK中辅助周跳检测和瞬时速度估计，提高滤波器的鲁棒性。
总结：三种观测值构成了RTK算法的“原材料”。伪距提供初始位置和模糊度浮点解，载波相位提供高精度测量值，多普勒提供速度约束和周跳检测手段。理解它们的物理特性和误差模型，是设计高效RTK解算器的关键。