什么是宽巷观测值？

目标：把“宽巷观测值、宽巷模糊度、MW 组合、PPP-AR 中的作用”讲清楚。

一、结论

**宽巷观测值（Wide-Lane, WL）**是把两个频点的载波相位做差分组合后得到的“虚拟载波观测值”。

其特点是：

• 频率变低，波长变长（$c=\lambda \cdot f$）；
• 同样的距离误差折算成更小的周数误差；
• 因此，宽巷整周模糊度比单频模糊度更容易固定。

宽巷不是让观测值本身更精密，而是让“判断整数周”这件事变得更容易。

二、先明确符号和单位

注意：

• ${\Phi }_i$ 是周单位相位；
• $L_i$ 是米单位相位；
• 二者关系为：

$$L_i={\lambda }_i\cdot {\Phi }_i$$

很多宽巷公式看起来不同，原因就是单位不同。

三、定义宽巷相位组合

1. 周单位形式

宽巷相位组合最基本的形式是：

$${\Phi }_{WL}={\Phi }_1-{\Phi }_2$$

对应的宽巷整周模糊度为：

$$N_{WL}=N_1-N_2$$

其中：

• ${\Phi }_{WL}$：宽巷相位组合，单位为周；
• ${\Phi }_1,{\Phi }_2$：两个频点的载波相位观测值，单位为周；
• $N_{WL}$：宽巷整周模糊度；
• $N_1,N_2$：两个频点各自的整周模糊度。

2. 频率

宽巷的频率为两个频率之差：

$$f_{WL}=f_1-f_2$$

其中：

• $f_{WL}$：宽巷等效频率；
• $f_1,f_2$：两个载波频率，通常假设 $f_1>f_2$。

如果不限定频率大小，也可写成：

$$f_{WL}=|f_1-f_2|$$

3. 波长

宽巷的波长为：

$${\lambda }_{WL}=\frac{c}{f_{WL}}=\frac{c}{f_1-f_2}$$

其中：

• ${\lambda }_{WL}$：宽巷等效波长；
• $c$：真空光速；
• $f_{WL}$：宽巷等效频率。

以 GPS L1/L2 为例：

$$f_1=1575.42\text{MHz},f_2=1227.60\text{MHz}$$

$$f_{WL}=347.82\text{MHz}$$

$${\lambda }_{WL}=\frac{299792458}{347.82\times 10^6}\approx 0.862\text{m}$$

对比原始载波波长：

因此，GPS L1/L2 宽巷波长明显长于任一单频载波波长。

计算波长（Python）

组合                f_WL (MHz)     波长 (m)
GPS L1-L2           347.8200     0.8619
GAL E1-E5a          398.9700     0.7514
BDS B1-B3           292.5780     1.0247

# 宽巷频率与波长计算
c = 299792458.0  # 光速 (m/s)

# 各频点频率 (Hz)  -- 来自提供的宏定义
freq = {
    'GPS_L1':   1.57542e9,
    'GPS_L2':   1.22760e9,
    'GAL_E1':   1.57542e9,
    'GAL_E5a':  1.17645e9,
    'BDS_B1':   1.561098e9,
    'BDS_B3':   1.26852e9,
}

# 组合定义：名称，f1键，f2键
combos = [
    ("GPS L1-L2", "GPS_L1", "GPS_L2"),
    ("GAL E1-E5a", "GAL_E1", "GAL_E5a"),
    ("BDS B1-B3", "BDS_B1", "BDS_B3"),
]

print(f"{'组合':<15} {'f_WL (MHz)':>12} {'波长 (m)':>10}")
for name, k1, k2 in combos:
    f_wl = freq[k1] - freq[k2]          # 宽巷频率 (Hz)
    wl_m = c / f_wl                     # 波长 (m)
    print(f"{name:<15} {f_wl/1e6:12.4f} {wl_m:10.4f}")

四、米单位宽巷相位组合：不是 MW 组合

如果载波相位已用米表示：

$$L_1={\lambda }_1\cdot {\Phi }_1,L_2={\lambda }_2\cdot {\Phi }_2$$

则宽巷相位组合的米单位形式为：

$$L_{WL}={\lambda }_{WL}\cdot ({\Phi }_1-{\Phi }_2)$$

代入

$$L=\lambda \cdot \Phi ,c=\lambda \cdot f,f_{WL}=f_1-f_2$$

可得：

$$\boxed{L_{WL}=\frac{f_1\cdot L_1-f_2\cdot L_2}{f_1-f_2}}$$

其中：

• $L_{WL}$：米单位宽巷相位组合；
• $L_1,L_2$：米单位载波相位观测值；
• $f_1,f_2$：两个频点频率。

需要特别注意：

$\frac{f_1\cdot L_1-f_2\cdot L_2}{f_1-f_2}$ 只是米单位宽巷相位组合，不是完整的 Melbourne-Wübbena（MW）组合。

MW 组合必须同时使用载波相位和伪距。

五、Melbourne-Wübbena 组合：相位 + 伪距

MW 组合的作用是更直接地估计宽巷模糊度。

1. 宽巷相位组合

$$L_{WL}=\frac{f_1{L}_1-f_2{L}_2}{f_1-f_2}$$

2. 窄巷伪距组合

$$P_{NL}=\frac{f_1{P}_1+f_2{P}_2}{f_1+f_2}$$

其中：

• $P_{NL}$：窄巷伪距组合，单位为米；
• $P_1,P_2$：两个频点的伪距观测值，单位为米。

3. MW 组合

常用 MW 组合可写为：

$$MW=\frac{L_{WL}-P_{NL}}{{\lambda }_{WL}}$$

理想情况下：

$$MW\approx N_{WL}+B_{WL}+{\epsilon }_{MW}$$

其中：

• $MW$：Melbourne-Wübbena 组合，单位为周；
• $N_{WL}$：宽巷整周模糊度；
• $B_{WL}$：宽巷相关硬件偏差，包含卫星端和接收机端偏差；
• ${\epsilon }_{MW}$：MW 组合噪声，主要来自伪距噪声和多路径。

MW 组合的优点是：几何距离、接收机钟差、卫星钟差、对流层延迟和一阶电离层项在组合中基本被消除或抵消，剩下的主要是：

$$\text{宽巷整数模糊度}+\text{宽巷硬件偏差}+\text{观测噪声}$$

因此，MW 组合常用于宽巷模糊度估计。

但要注意：

MW 组合不能无条件单历元直接取整。工程中通常需要平滑、周跳检测、粗差剔除，并改正或消除宽巷硬件偏差。

六、为什么宽巷更容易固定整数？

假设某个组合观测值存在距离误差 $e$，其折算到周单位的误差为：

$$\delta N=\frac{e}{\lambda }$$

其中：

• $\delta N$：误差折算成的周数；
• $e$：距离误差，单位为米；
• $\lambda$：观测值对应的波长，单位为米。

如果距离误差都是 $5\text{cm}$：

计算5cm对应的周数（Python）

# 根据宏定义计算波长与5cm误差对应的周数
c = 299792458.0          # 光速 m/s

# 系统频率 (Hz)
freq = {
    'GPS_L1': 1.57542e9,
    'GPS_L2': 1.22760e9,
    'GAL_E1': 1.57542e9,
    'GAL_E5a': 1.17645e9,
    'BDS_B1': 1.561098e9,
    'BDS_B3': 1.26852e9,
}

# 组合定义: (名称, f1键, f2键)
combos = [
    ("GPS", "GPS_L1", "GPS_L2"),
    ("GAL", "GAL_E1", "GAL_E5a"),
    ("BDS", "BDS_B1", "BDS_B3"),
]

err_m = 0.05  # 5 cm 距离误差

print(f"{'Sys':>4} {'Freq1 (MHz)':>13} {'Freq2 (MHz)':>13} "
      f"{'WL freq (MHz)':>13} {'λ1 (cm)':>8} {'λ2 (cm)':>8} {'λ_WL (cm)':>9} "
      f"{'5cm→L1(cyc)':>12} {'5cm→L2(cyc)':>12} {'5cm→WL(cyc)':>12}")
for sys_name, k1, k2 in combos:
    f1 = freq[k1]
    f2 = freq[k2]
    f_wl = f1 - f2               # 宽巷频率
    lam1 = c / f1 * 100          # cm
    lam2 = c / f2 * 100
    lam_wl = c / f_wl * 100
    cyc1 = err_m / (c / f1)      # 周数
    cyc2 = err_m / (c / f2)
    cyc_wl = err_m / (c / f_wl)
    print(f"{sys_name:>4} {f1/1e6:13.4f} {f2/1e6:13.4f} {f_wl/1e6:13.4f} "
          f"{lam1:8.4f} {lam2:8.4f} {lam_wl:9.4f} "
          f"{cyc1:12.4f} {cyc2:12.4f} {cyc_wl:12.4f}")

所以，宽巷更容易固定的根本原因是：

波长越长，同样的距离误差对应的周数误差越小，越容易判断最接近哪个整数。

七、宽巷组合的主要用途

1. 辅助整周模糊度固定

宽巷模糊度满足：

$$N_{WL}=N_1-N_2$$

工程流程通常是：

1. 用 MW 组合估计宽巷模糊度浮点值；
2. 对 MW 序列进行平滑和质量控制；
3. 改正或消除宽巷硬件偏差；
4. 固定 $N_{WL}$；
5. 利用 $N_{WL}=N_1-N_2$ 约束后续窄巷模糊度或IFLC模糊度。

宽巷固定后，模糊度搜索空间会明显缩小。

2. 周跳探测

宽巷/MW 组合常用于周跳探测。若某一频点发生周跳，宽巷模糊度变化为：

$$\Delta N_{WL}=\Delta N_1-\Delta N_2$$

其中：

• $\Delta N_1$：第一个频点的周跳量；
• $\Delta N_2$：第二个频点的周跳量；
• $\Delta N_{WL}$：宽巷组合中的周跳变化量。

如果 $\Delta N_1\ne \Delta N_2$，宽巷组合通常能反映出突变。

但如果两个频点发生相同周数的周跳：

$$\Delta N_1=\Delta N_2$$

则：

$$\Delta N_{WL}=0$$

这种情况仅靠宽巷组合无法识别，需要联合几何无关组合、多普勒预测、时间差分相位等方法。

3. PPP-AR 中的宽巷固定

PPP-AR 中，原始非差模糊度通常不能直接当作整数固定，因为它吸收了卫星端和接收机端的相位硬件偏差。

更严谨的流程是：

PPP-AR 依赖宽巷，但不是因为宽巷“不受偏差影响”，而是因为：

• 宽巷波长长，对随机噪声和残余误差更宽容；
• 宽巷硬件偏差可以通过 OSB、FCB、UPD 等产品改正，或通过差分方式消除；
• 宽巷先固定后，可以降低后续窄巷或消电离层组合模糊度的固定难度。

八、常见误区

误区 1：把米单位宽巷相位组合当成 MW 组合

错误理解：

$$\frac{f_1{L}_1-f_2{L}_2}{f_1-f_2}$$

就是 MW 组合。

正确理解：这个式子只是宽巷相位组合$L_{WL}$。MW 组合还要减去窄巷伪距组合$P_{NL}$：

$$MW=\frac{L_{WL}-P_{NL}}{{\lambda }_{WL}}$$

误区 2：宽巷模糊度可以直接取整

宽巷更容易固定，但仍要满足：

• 周跳已探测；
• 多路径和伪距噪声受控；
• 宽巷硬件偏差已处理；
• 浮点值足够接近整数；
• 固定结果通过质量检验。

误区 3：PPP-AR 中非差模糊度天然是整数

PPP 中的原始非差模糊度通常吸收了硬件偏差，不天然保持整数性。只有经过偏差改正、重参数化或差分消偏后，才具备整数固定条件。

误区 4：宽巷组合能发现所有周跳

宽巷组合对 $\Delta N_1\ne \Delta N_2$ 的周跳较敏感，但对 $\Delta N_1=\Delta N_2$ 的同幅周跳不敏感。因此不能单独依赖宽巷进行周跳检测。

九、小结

记住一句话：

宽巷的价值不是提高载波观测精度，而是通过变长等效波长，让整数模糊度更容易可靠固定；MW 组合则是工程上估计宽巷整数的重要工具。