2、导航系统精度评估方法
精度评估，说白了就是回答一个问题：这个导航系统到底准不准？
我做了十几年导航试飞，见过太多理论上一套套、实际一飞就露馅的情况。所以这一章，咱们就聊聊精度评估的那些实打实的方法。
2.1 精度定义与度量
先搞清楚几个基本概念。精度不是简单的一个数，它有不同的度量方式。
2.1.1 圆概率误差（CEP）
CEP 是导航领域最经典的精度指标。它的定义是：以真实位置为圆心，包含 50% 定位点的圆半径。
什么意思呢？你想想看，如果 CEP 是 10 米，那就意味着有一半的定位点落在 10 米半径的圆内。另一半呢？在圆外。
我个人习惯用 CEP 做快速评估。为什么？因为它对异常点不敏感。你飞 100 个点，哪怕有 5 个点飘到天边去，CEP 可能也就涨个一两米。这在工程上其实挺实用的——我们更关心的是“大多数情况下”的表现。
CEP 计算公式：
CEP ≈ 0.589 × (σx + σy)其中 σx、σy 分别是东向和北向误差的标准差。
2.1.2 均方根误差（RMS）
RMS 就不一样了。它把所有误差都平方求和再开方，对大的误差特别敏感。
我记得有一次评估某型组合导航系统，CEP 显示 8 米，看起来不错。但 RMS 算出来 25 米。一查数据，原来有几次卫星失锁，误差瞬间飙到 100 多米。CEP 没怎么变，RMS 却暴露了问题。
所以我的建议是：CEP 和 RMS 要一起看。一个看常态，一个看极端。
2.1.3 2DRMS
2DRMS 是 RMS 的两倍，对应大约 95% 的置信区间。说白了，就是“95% 的情况下误差不超过这个值”。
| 指标 | 置信概率 | 特点 |
|---|---|---|
| CEP | 50% | 对异常点不敏感 |
| RMS | 63% | 对异常点敏感 |
| 2DRMS | 95% | 工程常用指标 |
我的经验：做系统验收时，甲方通常要求 2DRMS。但你自己做研发评估，建议三个指标都算一遍。能帮你发现很多隐藏问题。
2.2 静态精度测试方法
静态测试，就是把导航设备放在一个已知坐标的位置上，开机采集数据。
听起来简单吧？但坑不少。
2.2.1 测试条件
- 位置选择：开阔无遮挡，多路径效应小。我一般选楼顶或空旷操场。
- 时间长度：至少连续采集 24 小时。为什么？因为 GPS 的星座周期是 24 小时，你得覆盖所有卫星几何构型。
- 采样率：通常 1Hz 或 10Hz。别太低，否则会漏掉短时抖动。
2.2.2 数据处理
静态测试的数据处理流程，我总结为三步：
- 粗差剔除：把明显离谱的点去掉。比如定位点跑到 1 公里外去了，那肯定是多路径或失锁。
- 误差计算：每个采样点减去真实坐标，得到误差序列。
- 统计计算：算 CEP、RMS、2DRMS。
注意：我曾经遇到过一个案例，静态测试数据看起来很好，CEP 只有 3 米。但一上飞机动态测试，CEP 直接飙到 20 米。后来发现是接收机的动态跟踪环路有问题。所以静态测试只能反映一部分性能，千万别迷信。
2.3 动态精度测试方法
动态测试才是真功夫。导航系统最终是要在运动状态下用的。
2.3.1 测试方案
动态测试的核心是基准真值。你得知道“真实轨迹”是什么，才能算误差。
常用的基准手段：
- 差分 GPS：精度厘米级，够用。但要注意基准站和流动站的距离不能太远。
- 高精度惯导：比如光纤陀螺或激光陀螺的惯导系统。我习惯用这个做基准，因为它不依赖卫星信号。
- 全站仪/激光跟踪仪：适合小范围、高精度的测试场景。
2.3.2 典型动态轨迹
我个人建议至少包含以下几种运动模式：
| 运动模式 | 目的 | 参数建议 |
|---|---|---|
| 匀速直线 | 评估稳态精度 | 速度 50-100 km/h |
| 加减速 | 评估加速度影响 | 加速度 2-5 m/s² |
| 转弯 | 评估角速度影响 | 转弯半径 50-200 m |
| 爬升/俯冲 | 评估垂直通道 | 爬升率 5-10 m/s |
避坑指南：我曾经做过一次动态测试，跑完一看数据，误差曲线特别漂亮。后来才发现是基准系统和被测系统用了同一根天线！信号分路器导致两个系统接收到的卫星信号高度相关。从那以后，我坚持用独立天线做基准。
2.4 精度评估的数据处理流程
这部分我直接给流程，你照着做就行。
2.4.1 数据预处理
1. 时间对齐：基准数据和被测数据的时间戳要统一
2. 坐标转换：都转到同一个坐标系（比如 WGS-84 或当地水平坐标系）
3. 插值处理：如果采样率不一致，用线性插值对齐
4. 异常检测：3σ 准则或 IQR 方法剔除粗差2.4.2 误差计算
对于每个采样点 i：
ΔE_i = E_被测_i - E_基准_i
ΔN_i = N_被测_i - N_基准_i
ΔU_i = U_被测_i - U_基准_i
水平误差 = sqrt(ΔE_i² + ΔN_i²)
三维误差 = sqrt(ΔE_i² + ΔN_i² + ΔU_i²)2.4.3 统计计算
RMS = sqrt( mean(水平误差²) )
CEP = median(水平误差) × 1.1774 （假设误差服从瑞利分布）
2DRMS = 2 × RMS重要提醒：CEP 的系数 1.1774 只在误差服从瑞利分布时成立。如果你的数据明显偏态，建议直接用排序法算中位数。
2.4.4 结果呈现
我习惯用三种图来呈现结果：
- 误差散点图：看分布形态，有没有系统性偏差
- 误差时间序列：看有没有时变特性
- 累积分布函数（CDF）：看不同置信水平下的误差
嗯，这一章的内容就这些。精度评估说难不难，说简单也不简单。关键是要理解每个指标的物理意义，别光会算数。
下一章我们聊聊导航系统的完好性评估，那又是另一番天地了。