2. IMU基础
在进入 Allan 方差与噪声分析之前，必须先建立对 IMU 本体的清晰认识。Allan 方差所处理的“数据”，本质上就是 IMU 内部加速度计与陀螺仪输出的角速度和比力序列；而坐标系定义、单位、采样率与数据格式，直接决定了后续噪声建模是否成立。本节从组成、原理、输出与坐标系四个层面，为后续章节打下物理与数据基础。
一、IMU 的组成与信号链
惯性测量单元（Inertial Measurement Unit，IMU）通常由三轴加速度计、三轴陀螺仪、信号调理电路、模数转换器（ADC）、温度传感器以及数字接口组成。部分器件还会集成磁力计，形成九轴 MARG 模块，但严格意义上的 IMU 只包含加速度计与陀螺仪。
其典型信号链为：敏感质量受惯性作用产生位移或应力 → 电容/压阻/谐振等检测方式转换为电信号 → 调理放大与滤波 → ADC 采样 → 数字滤波与温度补偿 → 通过 SPI/I²C/UART 输出。理解这条链路非常关键，因为 Allan 方差中观察到的量化噪声、角度随机游走、零偏不稳定性，分别对应链路中不同环节的物理与电气特性。
二、加速度计与陀螺仪工作原理
加速度计测量的是“比力”，即非引力外力产生的加速度。以最常见的 MEMS 电容式加速度计为例：敏感质量块通过弹性梁悬挂，当载体加速时，质量块因惯性产生位移，导致差分电容变化，电容变化量正比于加速度。静止水平放置时，Z 轴加速度计输出约为 1 g，这是重力反作用力的体现，而不是“重力加速度被直接测量”。
陀螺仪测量角速度。MEMS 陀螺仪多采用音叉式或振动环结构，利用科里奥利效应：驱动质量在平面内高频振动，当器件绕敏感轴旋转时，科里奥利力使质量块在正交方向产生位移，该位移被电容检测并解调为角速度。陀螺仪不直接输出角度，角度需通过对角速度积分获得，因此零偏和角度随机游走会随时间累积成角度误差。
关键区别：加速度计输出的是比力，包含重力分量；陀螺仪输出的是角速度，不含绝对角度。两者在 Allan 方差中的噪声项含义不同：加速度计常见速度随机游走，陀螺仪常见角度随机游走。
三、常见输出数据与单位
IMU 原始输出通常为三轴加速度和三轴角速度，部分器件还输出温度、时间戳和状态字。不同厂商的单位与量程不同，使用前必须统一。下表列出常见输出项及其典型单位与 Allan 方差对应关系。
| 输出项 | 物理含义 | 常见单位 | Allan 方差对应噪声 |
|---|---|---|---|
| 加速度计三轴 | 比力（含重力） | m/s²、g | 速度随机游走、零偏不稳定性 |
| 陀螺仪三轴 | 角速度 | °/s、rad/s | 角度随机游走、零偏不稳定性 |
| 温度 | 器件温度 | ℃ | 用于零偏温漂建模 |
| 时间戳 | 采样时刻 | s、μs | 决定采样间隔 τ₀ |
| 状态字 | 有效性与量程 | 位标志 | 数据质量筛查 |
采样率是 Allan 方差计算的关键参数。若 IMU 输出 200 Hz，则采样间隔 τ₀ = 0.005 s，Allan 方差的最小簇时间通常从 τ₀ 开始，逐倍递增。若采样率不稳定或存在丢包，Allan 方差曲线在高频段会出现异常抬升，容易被误判为量化噪声。
四、坐标系定义与安装关系
IMU 数据必须绑定坐标系才有物理意义。常见坐标系包括：
- 载体坐标系（b 系）：固连于 IMU 本体，通常 X 前、Y 右、Z 下，符合右手定则。
- 导航坐标系（n 系）：常用东北天（ENU）或北东地（NED），用于姿态解算与导航。
- 传感器坐标系（s 系）：各敏感轴实际指向，与 b 系之间存在安装误差和非正交误差。
在 Allan 方差分析中，通常直接使用传感器坐标系下的原始输出，不进行姿态旋转，因为噪声参数是器件固有属性。但若要做多 IMU 融合或与真值对比，则必须明确 b 系与 n 系之间的旋转矩阵。
实操建议：拿到 IMU 数据后，先确认三件事——单位是 °/s 还是 rad/s、加速度是 g 还是 m/s²、坐标系是 ENU 还是 NED。这三项不统一，后续 Allan 方差曲线的纵轴解读会完全错误。
五、案例：从数据手册到 Allan 方差输入
假设某 MEMS IMU 数据手册给出：陀螺仪量程 ±250 °/s，零偏不稳定性 8 °/h，角度随机游走 0.3 °/√h；加速度计量程 ±4 g，速度随机游走 0.05 m/s/√h。实际采集时，我们得到的是 200 Hz 的原始计数或物理量序列。
# 典型 IMU 输出数据列（CSV 示例）
timestamp, ax, ay, az, gx, gy, gz, temperature
0.000000, 0.012, -0.003, 9.807, 0.0012, -0.0008, 0.0005, 25.3
0.005000, 0.011, -0.002, 9.806, 0.0011, -0.0009, 0.0004, 25.3
0.010000, 0.013, -0.004, 9.808, 0.0013, -0.0007, 0.0006, 25.4
...
# 单位：加速度 m/s²，角速度 rad/s，温度 ℃
# 采样率：200 Hz，τ₀ = 0.005 s将上述数据按轴分别送入 Allan 方差计算流程，即可得到六条 Allan 偏差曲线。若某轴曲线在 τ = 1 s 附近出现明显“谷底”后抬升，通常对应零偏不稳定性；若在短 τ 段呈 -1/2 斜率，则对应角度随机游走或速度随机游走。这些判断都依赖于本节所讲的单位、坐标系和采样率前提。
六、常见误区与复盘
误区一：把加速度计输出直接当作“运动加速度”。静止时 Z 轴仍有 1 g 输出，这是重力反作用力，不是运动加速度。Allan 方差分析的是传感器噪声，不是运动轨迹。
误区二：忽略坐标系定义。同一组数据在 ENU 和 NED 下，Z 轴符号相反，若混用会导致零偏符号和融合结果错误。
误区三：采样率不稳定仍直接计算 Allan 方差。Allan 方差假设等间隔采样，丢包或抖动会污染高频段，建议先做时间戳检查与重采样。
误区四：把数据手册的噪声参数直接当作实测结果。手册值通常是典型值，实际器件受温度、安装应力和电路噪声影响，必须通过实测 Allan 方差重新标定。