2、卫星导航基础回顾：GNSS基本原理、伪距测量、载波相位、多普勒效应、定位解算流程
在深入北斗三号新信号B1C与B2a的接收特性之前，必须夯实卫星导航的公共理论基础。GNSS（全球导航卫星系统）的定位本质上是测距交会：通过测量用户接收机与多颗已知位置卫星之间的距离，解算用户的三维坐标和接收机钟差。本章将系统梳理GNSS的基本原理、四种核心测量量（伪距、载波相位、多普勒、载噪比）的物理意义与数学表达，并完整呈现从原始观测值到最终定位结果的解算流程。这些内容不仅是理解B1C/B2a新信号特性的前提，也是后续章节中分析信号捕获、跟踪、测量质量评估的基石。
上图展示了GNSS定位解算的完整链路：从卫星发射信号，经空间传播（含误差），到接收机前端处理、基带信号处理，最终生成四种核心测量量并参与定位解算。图中特别标注了B1C/B2A新信号在测量量生成环节的改进潜力，这是后续章节展开的重点。
2.1 GNSS基本原理：测距交会与时空基准
核心原理：GNSS定位的基本思想是空间测距交会。接收机测量其到至少四颗卫星的距离，以每颗卫星为球心、以测量距离为半径画球面，多个球面的交点即为接收机位置。由于接收机时钟与卫星时钟（基于原子时）存在偏差，实际测量的距离并非几何真距，而是包含钟差影响的"伪距"，因此需要至少4颗卫星同时解算4个未知数（三维坐标X/Y/Z和接收机钟差δtu）。
时空基准：所有GNSS系统都建立在精确的时空基准之上。时间基准采用各系统定义的导航系统时（如北斗的BDT，GPS的GPST），空间基准采用各自的地心地固坐标系（如北斗的CGCS2000，GPS的WGS-84）。卫星在轨位置由地面运控系统注入的星历参数计算，用户接收机通过解调电文获取这些参数。
ρ = r + c·(δt_u - δt_s) + I + T + ε其中，ρ为伪距观测值，r为卫星与接收机的几何距离，c为光速，δt_u为接收机钟差，δt_s为卫星钟差，I为电离层延迟（与频率相关），T为对流层延迟（与频率无关），ε为多路径和噪声等未模型化误差。
2.2 伪距测量：码相位测距原理
测量原理：伪距测量基于测距码（C/A码、P码、B1C码等）的相关特性。卫星和接收机各自产生相同的测距码副本，接收机通过滑动相关搜索，找到与接收信号相关峰值最大的本地码相位偏移量，该偏移量乘以码片宽度再乘以光速，即为信号传播时间对应的距离。由于码率有限（如B1C的码率为1.023Mcps，码片宽度约293m），伪距测量精度通常在分米到米级。
B1C/B2A的伪距改进：北斗三号B1C信号采用BOC(1,1)调制，B2A信号采用QPSK(10)调制，其码率更高（B2A码率10.23Mcps），码片宽度更窄（约29.3m），理论上伪距测量噪声更低。此外，B1C/B2A均设计了导频通道，导频通道不含电文数据，可以采用更长的相干积分时间，提升弱信号环境下的测距能力。
| 信号 | 调制方式 | 码率(Mcps) | 码片宽度(m) | 理论测距噪声(1σ) |
|---|---|---|---|---|
| B1I（旧） | BPSK(2) | 2.046 | 146.5 | ~0.3m |
| B1C（新） | BOC(1,1) | 1.023 | 293.0 | ~0.5m（但抗多径更好） |
| B2a（新） | QPSK(10) | 10.23 | 29.3 | ~0.05m |
2.3 载波相位测量：高精度但存在整周模糊度
测量原理：载波相位测量是跟踪接收信号载波与接收机本地载波之间的相位差。由于载波频率远高于码率（如B1C载波频率1575.42MHz，波长约19cm），相位测量精度可达毫米级。但载波相位观测值包含未知的整周数N（整周模糊度），且跟踪过程中可能发生周跳（信号失锁后重新锁定导致整周计数跳变）。
观测方程：载波相位观测值φ（以周为单位）可表示为：
φ = λ⁻¹·(r + c·(δt_u - δt_s) - I + T) + N + ε_φ注意电离层延迟项前的符号为负（与伪距相反），这是因为载波相位在电离层中传播速度加快（相速度超光速），而码信号传播速度减慢（群速度低于光速）。这一特性使得双频组合可以消除电离层延迟一阶项。
2.4 多普勒效应：测速与信号捕获的基础
物理原理：卫星高速运动（约3.87km/s）导致接收信号频率发生偏移。多普勒频移f_d与卫星-接收机之间的相对径向速度v_r成正比：f_d = f_T·v_r/c。对于B1C信号（1575.42MHz），最大多普勒频移约为±5kHz；对于B2a信号（1176.45MHz），最大约为±3.7kHz。多普勒频移不仅是接收机测速的依据，也是信号捕获阶段必须搜索的频率维度。
工程应用：在接收机基带处理中，多普勒频移的估计精度直接影响载波跟踪环路的初始频率设置。捕获阶段通常以500Hz或更小的步进搜索多普勒频率范围，而跟踪阶段则通过FLL（锁频环）和PLL（锁相环）精确跟踪多普勒变化。B1C/B2A新信号在捕获阶段需要同时搜索码相位和频率两个维度，其搜索空间大小与信号码长和多普勒范围直接相关。
2.5 定位解算流程：从观测值到PVT输出
标准流程：接收机完成信号跟踪和电文解调后，进入定位解算模块。典型流程如下：
- 卫星位置计算：根据电文中的星历参数（开普勒轨道参数或北斗的广播星历），计算各卫星在信号发射时刻的位置和速度。
- 误差修正：应用卫星钟差参数修正、电离层延迟修正（单频用户使用Klobuchar模型，双频用户使用无电离层组合）、对流层延迟修正（使用Saastamoinen或Hopfield模型）。
- 建立观测方程组：将修正后的伪距（或载波相位）观测值代入观测方程，形成线性化方程组。
- 迭代求解：采用最小二乘法或扩展卡尔曼滤波（EKF）迭代求解位置、速度、钟差和钟漂。通常需要2-3次迭代收敛。
- 质量检验：计算残差、DOP值（精度衰减因子），判断定位结果是否可用。
| 步骤 | 输入 | 输出 | 关键参数 |
|---|---|---|---|
| 1. 卫星位置 | 星历参数、信号发射时间 | 卫星ECEF坐标(X_s,Y_s,Z_s) | 轨道根数、摄动修正项 |
| 2. 误差修正 | 原始伪距、频率、气象参数 | 修正后伪距 | 电离层模型系数、对流层天顶延迟 |
| 3. 方程建立 | 修正后伪距、卫星坐标 | 线性化观测矩阵H | 方向余弦、几何矩阵 |
| 4. 迭代求解 | 观测矩阵、权阵 | 位置、钟差 | 收敛阈值（通常位置<1cm） |
| 5. 质量检验 | 残差向量、DOP | 定位精度指标 | PDOP<6为良好 |
背景：某接收机在电离层活跃期（总电子含量TEC=100TECU）进行单频B1I定位，垂直电离层延迟误差约15m。现改用B1C+B2a双频无电离层组合定位。
分析：双频组合伪距ρ_IF = (f1²·ρ1 - f2²·ρ2)/(f1² - f2²)，其中f1=1575.42MHz，f2=1176.45MHz。该组合消除了电离层一阶项，但放大了噪声约3倍（噪声放大因子 = √((f1²/(f1²-f2²))² + (f2²/(f1²-f2²))²) ≈ 2.98）。
结果：单频B1I伪距误差约15m（电离层）+0.3m（噪声）=15.3m；双频组合后电离层误差消除，噪声放大至约0.9m，总误差降至约1m。定位精度从约10m提升至约2m（水平）。
代价：双频组合需要同时跟踪B1C和B2a两个信号，接收机功耗和复杂度增加；且B2a信号在遮挡环境下可能比B1C更弱，导致组合观测值不可用。
2.6 常见误区与总结复盘
误区一：认为伪距精度等于定位精度。实际上，定位精度还取决于卫星几何分布（DOP值）和误差修正的完备性。即使伪距测量噪声为0.1m，若PDOP=10，定位误差仍可能达到1m。
误区二：混淆电离层延迟对伪距和载波相位的影响符号。伪距中电离层延迟为正（延迟），载波相位中为负（超前）。在双频组合时，若符号处理错误，不仅无法消除电离层误差，反而会加倍。
误区三：忽视多普勒对伪距测量的影响。在动态接收机中，多普勒频移会导致码环跟踪的额外延迟，需在接收机设计中考虑动态应力对环路带宽的要求。