3. 惯性传感器原理与误差特性

在AUV组合导航系统中,惯性测量单元(IMU)是唯一能够自主、连续输出姿态与加速度信息的核心器件。它通常由三只正交安装的加速度计和三只正交安装的陀螺仪构成:加速度计测量比力(单位质量所受的非引力外力),陀螺仪测量角速度。二者经积分运算分别得到速度、位置与姿态,这一"航位推算"过程决定了惯导的短期精度高、长期误差随时间累积的基本特征。理解加速度计与陀螺仪的工作原理及其误差特性,是后续进行DVL、深度计组合与卡尔曼滤波设计的前提。

一、加速度计与陀螺仪的工作原理

加速度计的本质是一个质量-弹簧系统。以常见的石英挠性摆式加速度计为例,敏感质量在比力作用下产生位移,电容或电磁检测环节将位移转换为电信号,再经伺服回路反馈平衡,输出与比力成正比的电压或数字量。其输出模型可写为:

a_meas = (1 + S_a)·a_true + B_a + n_a

其中 S_a 为标度因数误差,B_a 为零偏,n_a 为随机噪声。需要特别强调:加速度计测量的是比力而非运动加速度,静止时其输出为重力加速度 g 在敏感轴上的投影。因此AUV在水平静止时,水平加速度计输出约为0,垂直加速度计输出约为 +g(取决于安装方向约定)。本文统一约定:载体坐标系为前-右-下(FRD),重力方向沿 +z 轴,静止时 z 轴加速度计输出 +g。

陀螺仪测量的是相对惯性空间的角速度。以MEMS陀螺为例,利用科里奥利效应:驱动质量沿某一方向振动,当载体绕敏感轴旋转时,科里奥利力使质量在正交方向产生位移,该位移经电容检测转换为角速度输出。其输出模型为:

ω_meas = (1 + S_g)·ω_true + B_g + n_g

陀螺仪输出经一次积分得到姿态角,姿态误差又通过重力投影耦合进加速度计通道,形成"姿态误差→比力投影误差→速度误差→位置误差"的传播链。这是惯导误差随时间发散的根本机理。

惯性传感器误差传播链(FRD坐标系,静止时z轴加速度计输出+g) 陀螺仪 ω_meas = ω + B_g + n_g 姿态解算 ∫ω dt → 姿态角 姿态误差 θ_err 零偏积分 + 随机游走 加速度计 a_meas = a + B_a + n_a 比力投影 f_b → f_n(含重力) 速度/位置误差 二次积分累积 重力投影误差 g·θ_err 主要误差源 • 零偏 B(常值) • 标度因数 S(比例) • 随机噪声 n(白噪声) • 零偏不稳定性(1/f) • 角度/速度随机游走 • 安装失准角 误差传播要点: 陀螺零偏 → 姿态角线性漂移 → 重力在水平轴投影产生虚假比力 → 速度误差随时间二次增长 → 位置误差三次增长。 加速度计零偏直接引入速度误差(一次增长)与位置误差(二次增长),是纯惯导发散最快的通道之一。

二、零偏与标度因数:可标定但会漂移的系统误差

零偏(Bias)是输入为零时的非零输出,单位分别为 °/h(陀螺)和 mg 或 m/s²(加速度计)。零偏分为常值零偏、零偏重复性和零偏稳定性:常值零偏可在实验室标定后补偿;零偏重复性指每次上电零偏的离散程度;零偏稳定性则描述一次通电过程中零偏的慢变波动,通常用Allan方差在1小时或10小时簇时间处的值表征。对AUV而言,陀螺零偏是航向误差的主要来源——若陀螺z轴零偏为 10 °/h,则纯惯导航向约 1 小时漂移 10°,若AUV以 1.5 m/s 航行,横向位置误差可达数百米量级。

标度因数(Scale Factor)是输出增量与输入增量的比值,理想值为1。标度因数误差表现为比例性偏差,例如标度因数误差 1000 ppm 意味着 100 °/s 的输入会产生 0.1 °/s 的附加输出。标度因数误差与零偏的区别在于:零偏与输入无关,标度因数误差与输入成正比,因此在大角速率机动或高动态条件下影响更显著。工程上常通过多位置转台试验或速率试验分别辨识零偏与标度因数,并用温度多项式建模补偿温漂。

三、随机噪声与Allan方差分析基础

随机噪声无法通过单次标定消除,只能用统计特性描述。Allan方差是识别惯性传感器随机误差类型的标准工具:将长时间静态采集的数据按簇时间 τ 分组,计算各组平均值的方差,再对 τ 作图(双对数坐标)。不同斜率对应不同噪声类型:

  • 斜率 -1/2:角度随机游走(ARW)或速度随机游走(VRW),由白噪声积分产生,对应曲线左段。
  • 斜率 0:零偏不稳定性,曲线在某一 τ 处出现"谷底",该处方差值即零偏稳定性。
  • 斜率 +1/2:速率随机游走,对应曲线右段上升部分。
  • 斜率 +1:速率斜坡,常由温漂或未建模的慢变趋势引起。

实操中,采集静止状态下至少 2~4 小时、采样率 100~200 Hz 的陀螺与加速度计原始数据,去除明显野值与开机瞬态,按对数间隔选取簇时间 τ(如 1 s 到 1000 s),计算Allan偏差 σ(τ),在双对数坐标上拟合各段斜率,读取 τ=1 s 处的 σ 值换算 ARW(单位 °/√h),读取谷底值作为零偏稳定性。这些参数直接作为卡尔曼滤波中陀螺与加速度计随机游走的过程噪声输入。

教学示例(合理假设,非实测):某MEMS陀螺静态数据Allan分析得 ARW ≈ 0.3 °/√h,零偏稳定性 ≈ 8 °/h,速率随机游走 ≈ 50 °/h/√h。若滤波中把零偏稳定性设得过小,滤波器会过度信任陀螺,导致航向修正不足;设得过大则姿态估计抖动加剧。参数应与Allan分析结果同量级。

四、工程上的处理顺序

面对一只新IMU,合理的处理顺序是:先做温度相关的零偏与标度因数标定并建立补偿表;再采集长静态数据做Allan分析,确定ARW、零偏稳定性与随机游走参数;最后将标定补偿后的数据与噪声参数一并送入组合导航滤波器。需要牢记:标定能消除系统性常值误差,Allan分析刻画的是不可消除的随机部分,二者共同构成IMU误差模型,缺一不可。