网络RTK定位学习总结

网络 RTK 定位流程总结

网络 RTK 的本质，是把"一个基准站 + 一条基线"的常规 RTK，改造成"一片参考站网 + 一个区域误差模型"的差分服务 ：在一个较大区域内均匀布设若干坐标已知的连续运行基准站（站间距 50~100 km），各站以 1 Hz 连续观测并把原始观测值经光纤/专网实时汇入数据处理中心；中心先做同步预处理（钟差、周跳、粗差、多路径），再在参考站之间解算并固定双差整周模糊度（BDS 用 B2/B3 超宽巷→宽巷→B1/B2/B3 三频逐级固定，GPS 用 L1/L2 宽巷→窄巷逐级固定，并以多频模糊度间的整数线性关系作为强约束），模糊度固定后即可从"观测值 − 已知站星几何距离 − 模型改正 − 已固定模糊度"中反算出每站每星的残余综合误差 （电离层、对流层、轨道等）说明：双差整周模糊度到残余综合误差的过程是通过假定所有站的参考卫星及某一个站的其他卫星为0，算出单星单站的整周模糊度，然后反算出的残余综合误差，这个综合误差不关心真实值只关注相对值，再把这些误差在网内做空间建模（三站线性内插、低阶曲面或 FKP 格网）这里是依据参考站的位置进行建模，按用户实时位置生成改正数或虚拟观测值（VRS / 主辅站 MAX / 区域改正数 FKP / 综合误差内插 CBI），经 RTCM 电文 + NTRIP/移动网络播发给流动站；流动站用概略坐标（伪距单点定位）接收对应改正信息若数据中心播放的是综合误差的梯度信息，则是单向通信，若生成VRS，则需要双向通信，改正后观测值的残余误差被压到厘米级、模糊度整数特性得以恢复，于是"解算流动站模糊度"就退化为一次超短基线的常规 RTK 固定（LAMBDA 搜索 + Ratio 检验），最终得到厘米级坐标并完成 CGCS2000/地方坐标与高程基准的转换。整条链路可以概括为"网---中心---模型---改正数---流动站"五段式：参考站网提供已知基准，数据中心把基准站间的模糊度固定下来，误差模型把离散的站上误差连续化，改正数把区域误差"搬运"到用户处，流动站在被"净化"过的观测值上完成最后一步整数固定。

主要依据

1. 李征航、黄劲松《GPS 测量与数据处理》（第 3 版），武汉大学出版社，2016 ------ §5.8.3 动态相对定位、§5.8.4 RTK、§5.9 网络 RTK 及 CORS、§5.10 差分 GPS。
2. 刘文建《北斗/GNSS 区域地基增强服务系统建立方法与实践》，武汉大学博士学位论文，2018 ------ 系统架构、基准站网建设、坐标框架维持、改正数生成与播发、用户服务。
3. 高猛《基于非差观测模型的 BDS/GPS 网络 RTK 算法研究》，辽宁工程技术大学博士学位论文，2019 ------ 参考站间多频模糊度固定、非差模糊度递推、非差区域误差改正、流动站多频非组合模糊度解算。

---

一、概述

网络 RTK 的本质，是把"一个基准站 + 一条基线"的常规 RTK，改造成"一片参考站网 + 一个区域误差模型"的差分服务 ：在一个较大区域内均匀布设若干坐标已知的连续运行基准站（站间距 50~100 km），各站以 1 Hz 连续观测并把原始观测值经光纤/专网实时汇入数据处理中心；中心先做同步预处理（钟差、周跳、粗差、多路径），再在参考站之间解算并固定双差整周模糊度（BDS 用 B2/B3 超宽巷→宽巷→B1/B2/B3 三频逐级固定，GPS 用 L1/L2 宽巷→窄巷逐级固定，并以多频模糊度间的整数线性关系作为强约束），模糊度固定后即可从"观测值 − 已知站星几何距离 − 模型改正 − 已固定模糊度"中反算出每站每星的残余综合误差 （电离层、对流层、轨道等），再把这些误差在网内做空间建模（三站线性内插、低阶曲面或 FKP 格网），按用户实时位置生成改正数或虚拟观测值（VRS / 主辅站 MAX / 区域改正数 FKP / 综合误差内插 CBI），经 RTCM 电文 + NTRIP/移动网络播发给流动站；流动站用概略坐标（伪距单点定位）接收对应改正信息，改正后观测值的残余误差被压到厘米级、模糊度整数特性得以恢复，于是"解算流动站模糊度"就退化为一次超短基线的常规 RTK 固定（LAMBDA 搜索 + Ratio 检验），最终得到厘米级坐标并完成 CGCS2000/地方坐标与高程基准的转换。整条链路可以概括为"网---中心---模型---改正数---流动站"五段式：参考站网提供已知基准，数据中心把基准站间的模糊度固定下来，误差模型把离散的站上误差连续化，改正数把区域误差"搬运"到用户处，流动站在被"净化"过的观测值上完成最后一步整数固定。

---

二、流程总览

阶段	名称	主要任务	输出
0	系统基础建设	基准站选址建站、坐标联测、框架维持	参考站网 + 精确 ITRF/CGCS2000 坐标
1	数据采集与传输	多频多系统观测、气象、1 Hz 实时回传	基准站原始观测数据流
2	预处理与质量控制	钟差/周跳/粗差/多路径处理、数据完整性	干净的单历元观测值
3	参考站间模糊度固定	超宽巷→宽巷→窄巷逐级固定（核心一）	各基线双差整周模糊度
4	双差→非差模糊度	用网络组合关系逐历元递推非差模糊度	非差整周模糊度（含基准）
5	区域误差建模	反算站上残余误差、空间内插/拟合（核心二）	误差模型 / 格网改正数
6	改正数生成与播发	VRS/MAX/FKP 生成 + RTCM 编码 + 播发	差分数据流（含用户位置）
7	流动站改正与固定	概略坐标→取改正→改正观测值→固定模糊度（核心三）	流动站固定解
8	坐标转换与输出	三维七参数、高程异常/大地高补偿	地方坐标与正常高成果
9	质量评估与保障	固定率、Ratio、残差、完整性监测	可用性与精度指标

---

三、详细展开

阶段 0　系统基础：参考站网与坐标基准

网络 RTK 的"绝对基准"来自参考站坐标，因此这一步决定服务精度的天花板。

• 站点数量与分布 ：基准站数量由覆盖范围、定位精度要求和外部环境（尤其是电离层延迟的空间相关性）共同决定，至少 3 个 ；站间距常规 RTK 限制在 15 km 内，而网络 RTK 站间距可放宽到 50~100 km（李征航 §5.9.1）。
• 站址条件：全波长双频（多频）接收机、数据传输设备、气象仪器；视野开阔、多路径弱、地质稳定；精确坐标已知（李征航 §5.9.1「系统的组成」）。
• 坐标框架 ：基准站需与 ITRF/CGCS2000 高精度联测，并建立坐标参考框架的动态更新与维持 机制（刘文建，第 3、6 章）。这是"区域地基增强"区别于普通 RTK 的根本：不仅给改正数，还给统一的坐标与高程基准。
• 系统构成：基准站网 + 数据处理与播发中心 + 数据通信链路 + 用户终端四部分（李征航 §5.9.1；刘文建的北斗/GNSS 地基增强系统即按此分层搭建，含网络 RTK 改正数子系统与广域实时精密定位子系统、播发子系统、用户终端）。

阶段 1　数据采集与实时传输

• 基准站连续跟踪全部可见卫星，采集测码伪距 + 载波相位（BDS B1/B2/B3、GPS L1/L2，可扩展三频/多系统），采样率典型 1 Hz，附带气象数据与设备状态。
• 通信分两类（李征航 §5.9.1）：固定站之间 （基准站---数据中心---播发中心）用光纤、专网等有线方式，强调稳定与低时延；播发中心---流动站用 GSM/GPRS/4G/5G 等无线移动通信。
• 数据以 RTCM 格式打包上传，中心需处理"传输时延---数据龄期"问题：观测值到达中心时已经"过时"，需要按历元严格对齐（这也是异步 RTK/低时延输出要单独研究的原因）。

阶段 2　预处理与质量控制

数据中心对多站数据做统一处理（李征航 §5.9.1；刘文建第 6 章「数据处理中心」）：

1. 格式解码与历元对齐（含不同采样率、不同系统时间基准 BDT/GPST 的归化）；
2. 钟差处理：站间接收机钟差在双差中消去，卫星钟差用广播星历/精密钟差改正；
3. 周跳探测与修复：用 MW 组合、无几何组合（GF）、多普勒/相位减伪距等检验量定位并修复周跳；
4. 粗差与多路径抑制：残差检验、抗差估计、必要时降权或剔除该卫星；
5. 模型化改正：天线相位中心（PCO/PCV）、相位缠绕、固体潮与海洋负荷潮、地球自转、相对论效应等，均按模型先扣掉（李征航 §5.7 系列改正公式）；
6. 完整性监测：卫星数、数据完整率、信噪比、多路径指标，为后续是否可用提供判据。

这一步的目标只有一个：让残余误差足够小、足够"平滑"，以便下一步能把整周模糊度和误差分开。

阶段 3　参考站间整周模糊度固定（核心一）

这是整个网络 RTK 的前提条件 （高猛 §1.2、§4.1）。难点在于：参考站相距几十千米，双差后残余的电离层、对流层、轨道误差远大于半个波长，即使站坐标已知，直接取整也几乎必然失败。因此要"逐级固定 + 整数关系约束"。

BDS 三频方案（高猛 §4.2）

1. 超宽巷固定：用 B2/B3 的 MW 组合（宽巷相位减窄巷伪距），波长长、噪声小，配合严格的固定判据，数个历元内即可固定超宽巷整周模糊度（文中固定成功率 99.99%）；
2. 宽巷固定 ：把"已固定超宽巷与其他宽巷模糊度之间的整数线性关系"作为约束条件加入方程，同时估计宽巷模糊度、相对天顶对流层延迟 和电离层延迟参数，再搜索确定宽巷模糊度------用整数约束替代直接使用噪声较大的超宽巷组合，避免噪声污染参数估计；
3. 原始频率固定：用"宽巷模糊度与 B1/B2/B3 三频模糊度之间的整数线性关系"构造约束，代入载波相位参数估计模型，搜索得到三个原始频率的整周模糊度。

GPS 双频方案（高猛 §4.3） ：L1/L2 先固定宽巷（长波长 + 严格判据），再把宽巷与双频模糊度的线性关系作为约束，结合 L1/L2 无电离层组合确定 L1、L2 原始模糊度。

BDS/GPS 融合（§4.4）：先分别固定各系统的宽巷，再以 B1 与 L1 双差模糊度为对象联合估计，实现双系统统一固定（文中 B1 成功率 98.62%、L1 99.31%）。

长距离 / 低高度角处理（§4.5）：低高度角卫星大气延迟大、固定易错。做法是：先用高度角较好的卫星固定模糊度并提取其大气延迟信息，用这些信息约束低高度角卫星的宽巷模糊度，再以"宽巷---双频整数关系 + 无电离层组合"联合求解其双频模糊度。这样才能把 15° 截止高度角以上的卫星都用起来，改善网内几何构型和可用性。

关键点：单历元解算（只用当前历元的伪距 + 相位）可回避周跳处理，是实时服务的常用策略；而"整数线性关系约束"是把"逐级固定"的精度损失补回来的核心技巧。

阶段 4　双差 → 非差模糊度的转换

传统双差网络 RTK 的改正数是"相对主站"的，用户换主站就得换一套改正数，作业灵活性差、且双差观测值间强相关（高猛 §1.1.1、§5.1）。非差（无差）模型要摆脱这一点，就必须先把双差模糊度还原成非差模糊度。

• 常规做法是用转换矩阵，但参考站数量增加时矩阵维数急剧膨胀，运算困难。
• 高猛的替代方案（§4.6）：利用"双差模糊度 = 两站两星非差模糊度的组合"这一关系，从相邻基线的单个双差模糊度出发，配合每站每星的非差基准模糊度，逐历元线性递推出全网所有非差模糊度，实现单历元快速计算。
• 一个重要的性质：非差基准模糊度取值不同，只是把差异转移到流动站模糊度上，会被流动站自身的模糊度参数吸收，不影响最终定位结果（高猛摘要第 5 条）。

有了非差模糊度，全网所有参考站地位平等，不再区分主站/辅站，为下一步"非差误差改正"铺路。

阶段 5　区域误差建模与改正数生成（核心二）

第一步：反算站上残余误差。 对每个参考站、每颗卫星、每个频率，用已知坐标算几何距离，扣掉所有模型化改正与卫星钟差，再扣掉已固定的整周模糊度，剩下的就是残余综合误差（伪距改正数与载波相位改正数分别计算）：

ΔPri=Pri−ρri−c(δtr−δti)−Tri−Iri−Mri \Delta P^i_r = P^i_r - \rho^i_r - c(\delta t_r - \delta t^i) - T^i_r - I^i_r - M^i_r ΔPri=Pri−ρri−c(δtr−δti)−Tri−Iri−Mri

Δφri=λφri−ρri−c(δtr−δti)−Tri+Iri−λNri−Mri \Delta \varphi^i_r = \lambda \varphi^i_r - \rho^i_r - c(\delta t_r - \delta t^i) - T^i_r + I^i_r - \lambda N^i_r - M^i_r Δφri=λφri−ρri−c(δtr−δti)−Tri+Iri−λNri−Mri

（双差形式下同理，见李征航式 5-124：∇Δφ⋅λ=∇Δρ−∇ΔN⋅λ+∇Δdorb+∇Δdion+∇Δdtrop+∑δi\nabla\Delta\varphi\cdot\lambda = \nabla\Delta\rho - \nabla\Delta N\cdot\lambda + \nabla\Delta d_{orb} + \nabla\Delta d_{ion} + \nabla\Delta d_{trop} + \sum\delta_i∇Δφ⋅λ=∇Δρ−∇ΔN⋅λ+∇Δdorb+∇Δdion+∇Δdtrop+∑δi，其中前三项之外的残余项正是网络 RTK 要建模的对象。）

第二步：误差的空间建模。 依据误差与距离的相关性分而治之：

• 电离层：强色散、空间相关性强，用参考站穿刺点 VTEC 建区域模型后内插（或直接用双频无几何组合提取）；
• 对流层 ：用相对天顶对流层延迟参数 + 投影函数建模（高猛在宽巷固定时同步估计该参数，避免与模糊度耦合）；
• 轨道误差：现势精密/预报星历已足够好，通常可视为零（李征航 §5.9.1 明确指出这一点）；
• 多路径与噪声：与距离无关、不可内插，只能靠站址环境、接收机性能、抗差处理限制。

第三步：内插到用户位置。 常用三种数学形式（高猛 §5.1、§5.2）：

1. 三站（或多站）平面线性内插（CBI，综合误差内插法）：以流动站周围 3 个基准站的同类误差为节点，按平面坐标线性内插；
2. 低阶多项式曲面拟合：误差视为坐标的二次/一次曲面，最小二乘拟合系数后代入用户坐标；
3. FKP 格网：中心把区域误差拟合成南北、东西方向的线性系数或格网改正数，直接播发参数。

高猛的创新在于：以上建模与改正全部在非差层面完成 ------ 改正数是"每站每星"的绝对量，不依附主站，因此可以跨 CORS 网融合，也不再受"必须与参考站同构双差"的约束（§5.2 引言、§5.4 小结）。

阶段 6　改正数生成与播发：VRS / MAX / FKP 三条技术路线

这是"同一个模型、三种交付方式"（李征航 §5.9.1「几种常用的方法」）：

（1）虚拟参考站 VRS（Virtual Reference Station） ------ 双向通信

1. 数据中心用各基准站观测资料估计残余系统误差，生成误差模型；
2. 用户先用测码伪距单点定位得到概略坐标，经 NMEA-GGA 上传给中心；
3. 中心把该点当作一个"虚拟基准站"，用模型估计该点的残余误差，按该点坐标与基准站坐标生成一组虚拟的载波相位观测值，播发给用户；
4. 用户把虚拟参考站的观测值与坐标当作真实基准站，按常规 RTK 解算自身位置。

优点：用户端与常规 RTK 完全一样，精度高；缺点：必须双向通信，用户容量受网络带宽与中心服务器负载限制。另外需要注意虚拟点的坐标与高程基准问题 ：中心常需对 VRS 点位做大地高补偿 ΔH\Delta HΔH 后再生成观测值，用户端再做三维七参数转换，才能得到地方坐标与正常高（刘文建关于 VRS-CORS 正常高测量的处理步骤：①计算 ΔH\Delta HΔH；②把 ΔH\Delta HΔH 加入 VRS 三维坐标；③用户定位成果即已含高程补偿，再做七参数转换）。

（2）主辅站 MAX（Master-Auxiliary Concept） ------ 可单向

中心选一个主站，其余为辅站；播发主站的完整改正数 + 各辅站相对于主站的改正数之差（以及坐标差）；用户收到后按自身位置在辅站之间内插出本站改正数。只做单向通信即可工作；若让中心代算内插（用户上传概略位置），则退化为双向。

（3）区域改正数 FKP（Flächen-Korrektur-Parameter） ------ 单向、用户数不受限

中心估计残余误差后，用南北方向和东西方向的误差系数描述区域误差，直接把参数播发出去；用户用自备的专用软件按位置计算本站残余误差。只需单向通信，用户数量不受限制，但对用户端软件有要求。

三种方式的共同点：中心算模型、用户算位置；差别只在"在哪一端做内插、要不要双向链路、用户端是否需要专用软件"。

编码与链路 ：改正数按 RTCM 2.3/3.x 电文封装（测站坐标、天线描述、伪距/相位改正、FKP 参数、网络 RTK 改正消息等），经 NTRIP（互联网）或移动网络播发；播发频率要与观测采样率、模糊度更新周期匹配，并对卫星升降、周跳等事件做"及时重置"。

阶段 7　流动站改正与模糊度固定（核心三）

1. 概略定位：用伪距观测值做单点定位，得到米级概略坐标（同时作为 VRS 的虚拟点位置或 MAX/FKP 的内插依据）；
2. 获取改正数：按位置匹配网络（换网/换主站时自动切换，非差模式下无需"换主站"）；
3. 改正观测值 ：把内插得到的电离层、对流层、轨道等残余误差从伪距与相位中扣除，恢复载波相位的模糊度整数特性；
4. 参数估计与固定 ：
   • 组星间单差（非差模型下可任选基准星、无需考虑基准星切换带来的模糊度重排）；
   • 采用多频非组合 方式，把 BDS 的 B1/B2/B3 与 GPS 的 L1/L2 载波相位同时 作为参数估计（高猛 §6.2），用某一频率伪距作约束（BDS 选 B3，因其码长较长、多路径与噪声较小），按高度角定权（E<30°E<30°E<30° 时 P=sin⁡2EP=\sin^2EP=sin2E，否则取 1），伪距:相位权比约 1:100；
   • 多历元法方程叠加 → 最小二乘得位置改正数与模糊度浮点解及协方差 → LAMBDA 搜索 → Ratio 检验/固定成功率判据 → 固定解；
   • 输出坐标并给出精度评定。

为什么网络 RTK 能"快"？因为改正后残余误差远小于半波长，模糊度搜索空间极小，往往单历元即可固定；常规 RTK 则需要数分钟初始化。

阶段 8　坐标转换与成果输出

• 坐标系：广播星历解算得到的是 ITRF/CGCS2000 框架下的大地坐标，工程应用需经**三维七参数（或四参数 + 高程拟合）**转换到 2000 国家大地坐标系/地方坐标；
• 高程 ：GNSS 得到的是大地高 HHH，需用（似）大地水准面精化模型或高程异常 ζ\zetaζ 转换为正常高 h=H−ζh = H - \zetah=H−ζ；
• 基准一致性 ：区域增强系统必须保证"基准站坐标框架---改正数---用户成果"三者同基准，这也是刘文建论文把坐标参考框架更新与维持单列成章的原因。

阶段 9　质量评估与服务保障

• 精度：平面 cm 级、高程数 cm（网络 RTK 典型指标）；
• 可靠性：模糊度固定率/Ratio 值、验后残差、参考站间闭合差；
• 可用性：可用卫星数（含低高度角卫星的利用）、几何强度（PDOP）、数据完整率、时延；
• 完整性监测：对基准站坐标位移、周跳、异常观测做实时告警，避免"错误改正数"污染用户；
• 容量与扩展 ：VRS 受双向链路与服务器负载限制，FKP 用户数不受限；进一步的发展是 PPP-RTK（状态空间域 SSR 改正，用少得多的播发带宽服务无限用户）与多网融合。

---

四、双差模型与不差（非差）模型的分野

维度	双差（传统）网络 RTK	非差（无差）网络 RTK
主辅站	必须选主参考站，其余为辅站	全网站平等，不分主辅
改正数形式	相对主站的双差残差	每站每星的非差绝对改正数
相关性	双差观测值间强相关，忽略则不严密	无测站/卫星间相关性
换站/换网	需频繁切换主站，网被切成子网	可在网内自由漫游、跨 CORS 网融合
流动站组合	必须与参考站同构双差	可任意星间单差，无需考虑基准星切换
实现难点	改正数与流动站双差形式需严格统一	需先解出非差模糊度（基准问题）

高猛论文的全部工作，本质上就是沿着右列把"双差→非差"这条改造路线走通：参考站多频模糊度固定（§4）→ 非差模糊度递推（§4.6）→ 非差区域误差改正（§5）→ 流动站多频非组合固定（§6）。

---

五、关键公式索引

环节	公式要点	出处
双差观测方程	∇Δφλ=∇Δρ−∇ΔNλ+∇Δdorb+∇Δdion+∇Δdtrop+∑δi\nabla\Delta\varphi\lambda=\nabla\Delta\rho-\nabla\Delta N\lambda+\nabla\Delta d_{orb}+\nabla\Delta d_{ion}+\nabla\Delta d_{trop}+\sum\delta_i∇Δφλ=∇Δρ−∇ΔNλ+∇Δdorb+∇Δdion+∇Δdtrop+∑δi	李征航 式(5-124)
网络 RTK 的误差处理逻辑	用 3 个基准站解算站间残余误差，再按用户位置内插或生成虚拟观测值	李征航 §5.9.1
参考站误差反算	伪距/相位改正数 = 观测值 − 几何距离 − 钟差 − 模型改正 − 模糊度	高猛 §5.2.1、§5.2.2
双差↔非差	双差模糊度由两站两星非差模糊度组合，逐历元递推	高猛 式(4.51)、§4.6
流动站非组合模型	星间单差 + 各频率独立模糊度，伪距约束、高度角定权	高猛 式(6.14)--(6.16)
高程补偿	VRS 点位加入 ΔH\Delta HΔH，用户成果再做七参数转换	刘文建（VRS-CORS 正常高测量步骤）

---

六、工程实施要点与易错点

1. 模糊度固定是"必要条件"而非"可选步骤"：参考站间模糊度没固定，改正数就只是"含模糊度偏差的伪改正数"，用户永远只能得到浮点解。
2. 误差的空间相关性决定站间距：电离层活跃区（低纬、日间、太阳活动高年）应缩短站间距；测站高差大时对流层建模要加入高程方向，否则高程精度先崩。
3. 整数线性关系是长距离固定的"救命稻草"：单纯逐级用组合观测值固定会累积噪声，必须用整数约束把已固定的宽巷"传递"给原始频率。
4. 低高度角卫星要用起来：城市/遮挡环境下可用卫星数比精度更关键，长距离低高度角固定算法（高猛 §4.5）直接决定服务的可用性。
5. 通信时延与数据龄期：实时服务中改正数与观测值都不"新"，需按历元对齐并评估外推误差，否则会出现"改正数错位"导致的系统性偏差。
6. 基准一致性：坐标框架、高程基准、改正数模型三者必须同源，否则厘米级定位会被"基准差"吃掉。
7. VRS 的隐形成本：双向通信 + 中心并发计算，用户规模是瓶颈；大用户量场景优先考虑 FKP/SSR 类单向播发。

---

附录：3 基准站完整数值算例（双差 → 非差 → 改正数）

以高猛 §4.7 实验配置为例（云贵 CORS 网：AB=53 km，BC=51 km，CA=35 km，流动站距三站 31/27/25 km），取 B1 频率、卫星 p、q。

① 组网：3 站连成三角形，AB、BC 为 2 条独立基线（生成树），CA 为检核基线。

② 固定双差：在 AB、BC 上逐级固定，设 ∇ΔN_AB^pq = +7、∇ΔN_BC^pq = −3（CA 不固定、只检核）。

③ 补基准 + 递推（缺口 m+n−1 = 4 个）：取 N_A^p = N_A^q = N_B^q = N_C^q = 0，则

• N_B^p = ∇ΔN_AB^pq = +7
• N_C^p = ∇ΔN_AB^pq + ∇ΔN_BC^pq = +4
• 检核：CA 边双差应为 N_A^p − N_C^p − N_A^q + N_C^q = −4；独立固定 CA 对比，闭合差非零即报警。

站\星	p	q
A	0（基准）	0（基准）
B	+7	0（基准）
C	+4	0（基准）

④ 反算综合误差：对 A、B、C 每站每星做减法 ΔΦ_r^i = λφ_r^i − ρ_r\^i + c(δt_r−δt\^i) + T_model + λN_r\^i，得三组改正数时间序列 {Δ_A}、{Δ_B}、{Δ_C}。

⑤ 三站平面内插（CBI 显式形式，以 A 为参考、B/C 给梯度）：

Δ(U)=ΔA+(∂Δ/∂x∂Δ/∂y)(xU−xAyU−yA),(∂Δ/∂x∂Δ/∂y)=(xB−xAyB−yAxC−xAyC−yA)−1(ΔB−ΔAΔC−ΔA)\Delta(U) = \Delta_A + \begin{pmatrix}\partial\Delta/\partial x & \partial\Delta/\partial y\end{pmatrix}\begin{pmatrix}x_U-x_A \\ y_U-y_A\end{pmatrix},\quad \begin{pmatrix}\partial\Delta/\partial x \\ \partial\Delta/\partial y\end{pmatrix} = \begin{pmatrix}x_B-x_A & y_B-y_A \\ x_C-x_A & y_C-y_A\end{pmatrix}^{-1}\begin{pmatrix}\Delta_B-\Delta_A \\ \Delta_C-\Delta_A\end{pmatrix}Δ(U)=ΔA+(∂Δ/∂x∂Δ/∂y)(xU−xAyU−yA),(∂Δ/∂x∂Δ/∂y)=(xB−xAxC−xAyB−yAyC−yA)−1(ΔB−ΔAΔC−ΔA)

内插结果 Δ(U) 既可直接作改正数播发（FKP/CBI），也可在 U 点生成虚拟观测值连坐标一起播发（VRS）。

3 站 = 3 条基线（2 独立 + 1 检核）、2 个独立双差、4 个基准模糊度、6 个非差模糊度、3 组改正数、1 个误差平面。 这也是网络 RTK"至少 3 个基准站"的根本原因：3 个点恰好撑起一个误差平面。

---

附录 B：转换矩阵法（∇ΔN=DN\nabla\Delta\mathbf{N}=\mathbf{D}\mathbf{N}∇ΔN=DN）

构造规则 ：以站 1 为基准站、卫星 1 为基准卫星，第 (i,j)(i,j)(i,j) 个双差在 D\mathbf{D}D 中占一行，仅 4 个非零元：Nij(+1)N_i^j(+1)Nij(+1)、Ni1(−1)N_i^1(-1)Ni1(−1)、N1j(−1)N_1^j(-1)N1j(−1)、N11(+1)N_1^1(+1)N11(+1)。

3 站 2 星实例 （未知数序 NAp,NAq,NBp,NBq,NCp,NCqN_A\^p,N_A\^q,N_B\^p,N_B\^q,N_C\^p,N_C\^qNAp,NAq,NBp,NBq,NCp,NCq）：

D=−1+1+1−10000−1+1+1−1,rank(D)=2,dim⁡ker⁡=6−2=4=m+n−1\mathbf{D}=\begin{bmatrix}-1&+1&+1&-1&0&0\\0&0&-1&+1&+1&-1\end{bmatrix},\quad \mathrm{rank}(\mathbf{D})=2,\quad \dim\ker = 6-2=4=m+n-1D=−10+10+1−1−1+10+10−1,rank(D)=2,dimker=6−2=4=m+n−1

解法一（基准约束，整数解） ：补 S\mathbf{S}S（选 NAp,NAq,NBq,NCq=0N_A^p,N_A^q,N_B^q,N_C^q=0NAp,NAq,NBq,NCq=0）拼成 6 阶方阵，

DSN=∇ΔN0⇒N=0,0,+7,0,+4,0T\begin{bmatrix}\mathbf{D}\\ \mathbf{S}\end{bmatrix}\mathbf{N}=\begin{bmatrix}\nabla\Delta\mathbf{N}\\ \mathbf{0}\end{bmatrix}\Rightarrow \mathbf{N}=0,0,+7,0,+4,0^{\mathsf T}DSN=∇ΔN0⇒N=0,0,+7,0,+4,0T

回代 DN=+7,−3\mathbf{D}\mathbf{N}=+7,-3DN=+7,−3 ✓；CA 闭合双差 =−4=-4=−4 ✓。整数进、整数出。

解法二（最小范数/伪逆） ：N=DT(DDT)−1∇ΔN=−1.833,+1.833,+1.667,−1.667,+0.167,−0.167T\mathbf{N}=\mathbf{D}^{\mathsf T}(\mathbf{D}\mathbf{D}^{\mathsf T})^{-1}\nabla\Delta\mathbf{N}=-1.833,+1.833,+1.667,-1.667,+0.167,-0.167^{\mathsf T}N=DT(DDT)−1∇ΔN=−1.833,+1.833,+1.667,−1.667,+0.167,−0.167T。回代同样满足，但与解法一相差一个零空间向量（1.833,−1.833,5.333,1.667,3.833,0.1671.833,-1.833,5.333,1.667,3.833,0.1671.833,−1.833,5.333,1.667,3.833,0.167），该差值即基准任意性，被流动站模糊度吸收；缺点是结果非整数，需再取整。

与递推法的关系 ：递推法等价于把 D\mathbf{D}D 按生成树拓扑重排后做稀疏三角分解 + 前代。同一方程组的稠密解（O((mn)3)O((mn)^3)O((mn)3)）与稀疏解（O(mn)O(mn)O(mn)），高猛选后者是为满足 1 Hz 实时性。

全组合的行冗余 ：全组合 (m2)(n2)\binom{m}{2}\binom{n}{2}(2m)(2n) 行，秩仍为 (m−1)(n−1)(m-1)(n-1)(m−1)(n−1)，多余行线性相关------这正是"多余基线只能做闭合检核"的矩阵解释。

---

附录 C：基准取 0 为何无害？改正数的"准确性"辨析

结论：反算出的残余综合误差绝对值含未知常数偏差（不是"真误差"），但其时变与空间相对部分准确；差分定位只用后者，故位置解不受影响。

机理 ：设真值 Nr∗ sN^{*\,s}_rNr∗s，递推值 Nrs=Nr∗ s+ar+bsN^s_r = N^{*\,s}_r + a_r + b^sNrs=Nr∗s+ar+bs（零空间结构），则相位改正数

Δ~rs=Δr∗ s−λar−λbs\tilde\Delta^s_r = \Delta^{*\,s}_r - \lambda a_r - \lambda b^sΔ~rs=Δr∗s−λar−λbs

流动站改正并做星间单差后：aaa 项（站类、全星相同）被消去/由接收机钟差参数逐历元吸收；bsb^sbs 项（星类、时间恒定）并入流动站模糊度估计值。流动站解出的模糊度 = "真模糊度 + 常数"，定位只需其恒定、不需其正确。类比：水准测量假定基准点高程为 0------高程带常数偏差，高差仍准。

会真正破坏精度的情形（常数变时变）：

1. 参考站周跳 → 改正数跳变一个波长（B1 约 19 cm），不可吸收，须实时探测重置；
2. 双差固定错误 → 错误常数混入真误差通道并经内插传给用户，三角形闭合检核负责拦截；
3. 内插误差（电离层/对流层空间结构偏离拟合平面）------真正的相对精度极限；
4. 多路径（站特异性、不可内插）。

补充：伪距改正数不含模糊度项，天然无基准问题。PPP-RTK（SSR）必须显式估计并播发 UPD，是因为 PPP 用户不做差分、无参数可吸收常数------非差网络 RTK 的"设 0 自由"是以用户端保留模糊度参数为前提的。

---

附录 D：误差的空间建模方法小结

问题本质：用每站每星历元的离散改正数采样，重建连续误差场 Δ(x,y,t)，在用户坐标取值。物理前提：电离层/对流层/轨道误差空间缓变（相关尺度≫站间距）；多路径与噪声空间尺度仅数米，不可内插，须建模前压制。

四类数学核：

1. 线性内插 LIM（Wanninger 1995，3 站）：3 站误差张成平面，两梯度由两站差分唯一确定------VRS 经典内核；
2. 低阶曲面 LCM（≥4 站）：二次曲面 + 最小二乘，冗余提可靠性；阶数须与站网密度匹配，过高会把多路径拟合进去；
3. FKP（Wübbena 1996）：中心把场拟合成线性区域多项式（参考面平行 WGS-84 椭球、站高程），播发南北/东西梯度参数；分散（电离层，按频率缩放）与非分散（对流层+轨道）两套系数；单向、用户数不限、坏站可重构平面；无气象设备时用标准大气，对流层精度受损（刘文建）；
4. CBI（高星伟 2002；改进版唐卫明 2006）：不做误差分离，直接内插双差综合误差（各误差均近似线性，其和亦然，避免分离引入的模型误差）；改进版加距离/高度角加权。非差改正数同属综合路线，且为绝对量，内插核任选。

内插截断误差 ≈ (D²/2)·|∂²Δ/∂x²|，50 km 站间距下典型厘米级以内------匹配模糊度固定需求。电离层活跃期二阶项增大 → 缩短站间距/升阶/分离建模；测站高差大时对流层需加高程项，平面模型失效。

建模 = 误差策略 × 数学核 × 执行位置：

方案	误差策略	数学核	执行位置	通信
VRS	综合/分离均可	LIM	中心代算	双向
MAX	分离（主辅站残差）	加权内插	用户/中心	单向或双向
FKP	分离（分散/非分散）	梯度参数	用户自算	单向
CBI/非差	综合打包	线性/加权内插	用户自算	单向

---

本文档依据上述三份文献的实际章节内容整理，公式为标准形式，符号含义以原书/原文为准。