2. 惯性导航基础

在IMU预积分与因子图融合定位系统中,惯性测量单元(IMU)是整个状态估计链路的最前端传感器。理解IMU的物理测量原理、坐标系约定、姿态表示方法以及运动学方程,是后续推导预积分模型、构建因子图约束的数学前提。本章从传感器物理层出发,逐层建立从原始测量到姿态描述的完整认知框架。

一、IMU传感器原理与测量模型

典型的6轴IMU由三轴加速度计和三轴陀螺仪组成,分别测量载体在惯性空间中的比力(specific force)和角速度。加速度计的核心是MEMS质量块-弹簧结构,当载体存在加速度时,质量块因惯性产生位移,通过电容变化换算出比力值。陀螺仪则利用科里奥利力效应,测量绕三轴旋转的角速度。

需要特别强调的是,加速度计测量的并非纯粹的“运动加速度”,而是比力,即除去重力后的非引力外力加速度。其理想测量模型可表示为:

a_m = R_wb^T · (a_w - g_w) + b_a + n_a
ω_m = ω_b + b_g + n_g

其中 a_m、ω_m 为加速度计和陀螺仪的原始测量值,R_wb 为世界系到机体系的旋转矩阵,g_w 为世界系重力向量,b_a、b_g 为加速度计和陀螺仪的零偏,n_a、n_g 为高斯白噪声。零偏是IMU误差的主要来源,也是预积分中必须在线估计的关键状态量。

工程提示:消费级MEMS IMU的陀螺仪零偏通常在0.1~10 °/h量级,加速度计零偏在0.1~10 mg量级。零偏不补偿时,姿态会在数十秒内发散,因此因子图中通常将零偏建模为随机游走状态。

二、坐标系定义与变换关系

惯性导航涉及多个坐标系的定义与转换,工程中必须严格区分:

  • 惯性系(i系):原点在地心,坐标轴指向恒星,牛顿力学定律严格成立。
  • 世界系/导航系(w系或n系):原点在载体初始位置,东北天(ENU)或北东地(NED)指向,重力方向固定,是融合定位的参考基准。
  • 机体系(b系):原点在IMU中心,前-右-下(FRD)或前-左-上(FLU)指向,与传感器固连。
  • 相机系(c系):原点在相机光心,z轴沿光轴方向,用于视觉惯性对齐。

各坐标系之间的变换由旋转矩阵 R 和平移向量 t 描述。机体系到世界系的变换 T_wb = [R_wb | t_wb] 即为待估计的位姿状态。IMU预积分的本质,就是在世界系下对机体系测量进行积分,从而建立两个关键帧之间的相对运动约束。

三、姿态表示方法:欧拉角、旋转矩阵与四元数

姿态表示是惯性导航的核心数学工具,三种主流方法各有优劣:

表示方法参数个数优点缺点典型用途
欧拉角3直观易读,便于人机交互存在万向节死锁,非线性强姿态显示、航向输出
旋转矩阵9线性变换,无奇异点冗余参数,需正交约束坐标变换、误差状态推导
四元数4紧凑无奇异,插值平滑物理意义不直观,需归一化预积分、因子图优化

在因子图优化中,四元数是最常用的姿态表示。其乘法更新规则为:

q_new = q_old ⊗ δq
δq ≈ [1, 0.5·δθ]^T

其中 δθ 为姿态误差的小角度扰动。这种局部参数化方式避免了四元数模长约束,使优化问题转化为无约束最小二乘,是VINS-Mono、ORB-SLAM3等系统的标准做法。

四、惯性运动学方程与预积分思想

惯性导航的运动学方程描述了位姿、速度、零偏随时间的演化关系。在连续时间下,世界系中的运动学方程为:

ṗ_wb = v_wb
v̇_wb = R_wb · (a_m - b_a - n_a) + g_w
Ṙ_wb = R_wb · [ω_m - b_g - n_g]×
ḃ_a = n_ba
ḃ_g = n_bg

其中 [·]× 表示反对称矩阵。直接对上述方程积分会带来两个问题:一是积分结果依赖初始状态,每次优化迭代都需重新积分;二是计算量随关键帧数量增长。预积分技术通过将测量积分到相对运动量中,使积分结果仅依赖零偏,从而大幅降低重复计算。

预积分量的定义为:

Δp_ij = ∫∫ R_ib · (a_m - b_a) dt²
Δv_ij = ∫ R_ib · (a_m - b_a) dt
Δq_ij = ∫ q_ib ⊗ [0, 0.5·(ω_m - b_g)] dt

这三个预积分量构成了因子图中IMU约束的核心,与视觉观测约束共同参与位姿、速度、零偏的联合优化。

IMU惯性导航核心链路:从传感器到预积分因子 IMU传感器 加速度计 a_m 陀螺仪 ω_m 零偏 b_a, b_g + 噪声 坐标系定义 惯性系 i / 世界系 w 机体系 b / 相机系 c R_wb, t_wb 变换关系 ENU / NED 约定 姿态表示方法 欧拉角 旋转矩阵 四元数 万向节死锁 / 正交约束 / 无奇异 因子图优化首选四元数 惯性运动学方程 ṗ = v, v̇ = R·(a_m - b_a) + g Ṙ = R·[ω_m - b_g]× 连续时间积分 → 离散化递推 IMU预积分量:Δp_ij, Δv_ij, Δq_ij 仅依赖零偏 → 因子图约束 → 与视觉联合优化 主要误差源 • 零偏随机游走 • 白噪声积分漂移 • 尺度因子/轴间耦合 • 温度漂移 预积分核心价值 将高频IMU测量压缩为关键帧间相对约束 避免优化迭代中重复积分,提升实时性 零偏更新时通过一阶雅可比近似修正

五、工程案例:从原始IMU到预积分因子的数据流

以VINS-Mono系统为例,IMU数据以200Hz频率进入系统,视觉帧率为30Hz。在相邻两个视觉关键帧之间,通常有6~7个IMU测量区间。系统首先对IMU测量进行预积分,得到 Δp_ij、Δv_ij、Δq_ij 及其对零偏的雅可比矩阵,然后构建IMU残差因子:

r_imu = [
  r_p, r_v, r_q, r_ba, r_bg
]^T

该残差与视觉重投影残差、边缘化先验残差共同构成因子图的目标函数,通过Ceres或GTSAM求解。整个过程中,坐标系定义的一致性、四元数更新的正确性、零偏雅可比的准确性,直接决定融合定位的精度和鲁棒性。

常见误区:许多初学者将加速度计测量值直接当作运动加速度进行积分,忽略重力项和比力定义,导致轨迹在垂直方向迅速发散。正确做法是:先通过初始静止段估计重力方向和零偏,再在预积分中显式补偿重力。

六、复盘与要点归纳

本节从IMU物理原理出发,建立了“传感器测量→坐标系变换→姿态表示→运动学方程→预积分”的完整知识链路。核心要点包括:加速度计测量的是比力而非纯加速度;世界系、机体系、相机系的定义必须严格统一;四元数是因子图优化中最优的姿态表示;预积分将高频测量转化为关键帧间的相对约束,是IMU与视觉融合的数学桥梁。掌握这些基础,才能正确推导预积分残差及其雅可比,为后续因子图构建打下坚实基础。