第 6 章 MEMS 与惯性传感¶
IMU（惯性测量单元）给控制器"机体自己怎么动"的短时高频信息。这一章讲 MEMS 惯性器件的误差源体系（不是"分辨率"那类刻度，而是漂移与噪声那类误差）、Allan 方差如何把它们分门别类、以及从噪声密度到积分漂移的量级判断，并立它与滤波/融合的消费分工。器件的微纳制造工艺与封装不在覆盖范围，滤波与状态融合归 control-note 第 8 章，激光雷达/结构光的器件原理归属尚未定。具体型号的数字按厂商规格核实，未核实只给方法与量级。
误差源经 Allan 方差分开，再加权成漂移量级：
flowchart LR
src["误差源：偏置不稳定 / 角度随机游走 / 温漂"] --> allan["Allan 方差<br/>按 τ 斜率分类"]
allan --> drift["噪声密度 → 积分漂移量级"]
drift --> fuse["交滤波与融合（control 第 8 章）"]1. MEMS 机制（一句话边界）¶
- 加速度计：一个可动检验质量块 + 差分电容（或压阻）读其位移，牛顿第二定律下质量块惯性使电容失衡，反馈回路把它拉回并读出加速度。
- 陀螺仪（振动结构）：驱动质量做简谐振动，旋转时科里奥利力（Coriolis）在与驱动、自转轴正交方向感应出位移，读出即角速度。
工艺侧的梳齿、真空封装、应力隔离等不在本章展开。要点是：MEMS 惯性器件输出天生带漂移与噪声，它的价值在"短时高频相对稳"，不是"长时间不漂"——因此几乎总要被外部观测（视觉/激光/里程）修正，这就是下面误差源与 Allan 方差存在的理由。
量程与分辨率是一对反向指标：同一颗器件换大量程，最小可分辨增量随之变粗。选型先按任务的最大角速度/加速度乘安全系数定量程，再核这个量程下的噪声密度，顺序反了容易落在"量程够但噪声大到不可用"的角落。三轴封装还带来轴间对齐误差；焊锡冷却、灌封、螺钉紧固都会把机械应力传到芯片并改零偏，装到结构上再标定，比只看出厂标定值更贴近实际。
2. 误差源体系¶
IMU 的 datasheet 一堆指标，本质是这几类误差源，逐条对应不同的抑制与标定手段：
- 零偏（bias）：静止时的非零输出。零偏本身慢变（随时间、温度）；描述它的关键指标是零偏不稳定性（bias instability）——把 Allan 方差曲线的谷底值，是"平均能压到多低"的下限。
- 噪声密度 / 随机游走：白噪声速率（陀螺 \([^\circ/\mathrm{s}]/\sqrt{\mathrm{Hz}}\)，加计 \([\mathrm{m/s^2}]/\sqrt{\mathrm{Hz}}\)），积分后成角度随机游走 ARW/速度随机游走 VRW——短时误差主要来自它。
- 标度因数（scale factor）误差与交叉耦合：输出对输入的增益偏差与轴向串扰，随负载/频率出现，是与转动量成比例的误差，不能靠静止标定消。
- 温度效应：零偏温漂、标度温漂、热滞后——温度梯度还会诱发虚拟重力/对流应力，工程上靠恒温、温补系数或安装对称削。
- g 敏感（加速度对陀螺零偏的耦合）与振动整流：把宽带振动整流成低频零偏（在振源上尤其明显），是"车/腿一抖陀螺就漂"的器件成因。
这五类误差的可观测性不同：零偏与标度因数靠静态多位置与转台标定，噪声密度靠静止长采样，温度系数靠温箱扫温，g 敏感靠振动台或不同安装朝向的对照。选型时先问“误差预算里哪一项占大头”，再认定要看的指标——任务持续几十秒的，噪声密度决定成败；持续几分钟的，零偏不稳定性与温漂才是瓶颈。
3. Allan 方差：把误差源分开¶
单看一张 Allan 方差（Allan deviation, \(\sigma(\tau)\) vs 平均时间 \(\tau\)，log-log）图就能把上面几类误差按斜率分开，这是选型与建模读器件数据的核心工具：
| 曲线段（log-log 斜率） | 对应误差源 | 读法 |
|---|---|---|
| −1/2 | 白噪声（角度/速度随机游走 ARW/VRW） | 短时，随 \(\tau\) 下降 \(\propto1/\sqrt{\tau}\) |
| −1 | 量化噪声（一句名目） | 更短 \(\tau\) |
| 0（谷底） | 零偏不稳定性 | 谷底对应的 \(\tau\) 与 \(\sigma\) 是"最稳工作点" |
| +1/2 | 零偏随机游走（速率漂移 RRW） | 长时慢漂 |
谷底位置还提示在什么时间尺度上 IMU 最可信——超过它就得靠外部观测重置。Allan 方差的计算流程需要一段公开 IMU 数据或仿真噪声模型；原始文献与器件具体曲线未核实前，这里只给方法与形状，不写某型号的引文号与数字（标待核实）。
计算本身是“按簇平均”：把采样序列按簇长 \(\tau\) 分组，相邻簇均值作差，求差值的方差再除二，得到该 \(\tau\) 下的方差估计；逐倍增长 \(\tau\) 重复，就得到一条双对数曲线。两个常见陷阱：总时长不够时，长 \(\tau\) 端只有几对样本，曲线尾部的“上翘”可能只是统计噪声；把不同温度、不同振动条件下的数据混在一起算，会把温漂与振动整流误读成零偏随机游走。
读厂商给的曲线先核样本条件：采样率、总时长、是否恒温、单轴还是三轴同时采集。总时长至少要覆盖最长关注时间的十倍，否则谷底可能落在样本不足的区段；同一型号在不同量程与内建滤波配置下的曲线不能直接比较，换档就等于换了一颗器件。
4. 从噪声密度到积分漂移¶
短时最有用的一条换算：白噪声速率密度 \(N\) 积分 \(\Delta t\) 时间，随机游走的标准差按 \(\sqrt{\Delta t}\) 增长——
取 \(N=0.01\ ^\circ/\mathrm{s}/\sqrt{\mathrm{Hz}}\)，1 s 时约 0.010°，10 s 时约 0.032°，60 s 时约 0.078°——这就是"漂移随时间开方累积"的直观量级。
误差预算要把两项相加而不是只算一项：白噪声贡献按 \(\sqrt{\Delta t}\) 长，零偏不稳定性贡献随 \(\Delta t\) 近似线性长。短时以噪声项为主，时间拉长后线性项必然追上——两项量级接近的时间点，就是“换更高等级器件”或“加外部观测”的决策点。
它给出直觉：漂移随时间开方累积，60 s 已漂到近 0.08°（且还没算零偏不稳定性与温漂——那两项长期更凶）。所以纯 IMU 航位推算几分钟就不可信，必须被低频外部观测修正；这也是"IMU 精度看分辨率参数"这一误读的反例（分辨率是刻度，误差预算看这些漂移项）。
5. 消费侧分工¶
本章止于器件误差源：谁能压到多低、随什么漂、如何被 Allan 方差分类。把这些误差用进估计——互补滤波、卡尔曼/扩展卡尔曼、因子图、与视觉/激光/轮式里程的融合，以及四元数姿态传播——归 control-note 第 8 章（滤波）与第 9 章（里程计/观测器）。
接口与时间戳（IMU 话题、sensor_msgs/Imu、协方差字段）归 ROS2-note。器件侧只保证：给出的 \((N,\ \text{bias instability},\ \text{scale},\ \text{温度系数})\) 是滤波器输入噪声与建模的依据。
交接给滤波器的还有三件事：采样率与带宽（低通截止频率决定能看到多快的运动，也决定噪声被滤掉多少）、安装对齐（器件轴与机体轴的夹角、减震座引入的额外自由度）、时间戳精度（融合对时延极敏感，毫秒级偏差在高速旋转下就是可见的角度误差）。这三项不清楚，调滤波器只能靠试参。
多 IMU 并联时，误差源要相互独立才有收益：同型号两颗装在同一块板上、共用一路电源与同一处温度，零偏温漂高度相关，简单平均压不掉多少。真正独立的是拉开基线与供电的布置（振动与温度都不同），靠两路比对检出漂移、卡死与故障，再交给融合层做可信度加权。
6. 要点¶
- MEMS 惯性器件天生带漂移与噪声：加计靠惯性质量块-电容、陀螺靠科里奥利效应，价值在短时高频，不长时间保真。
- 误差源分五类：零偏（不稳定性=Allan 谷底）、噪声密度/随机游走、标度因数与交叉耦合、温度效应、g 敏感与振动整流——各有各的抑制手段。
- Allan 方差按 log-log 斜率（−1/2 白噪、谷底零偏不稳、+1/2 零偏游走）把误差分开，并指出最可信时间尺度。
- 积分漂移随 \(\sqrt{\Delta t}\) 累积，示例 60 s 已 ~0.08°（不含慢漂），纯 IMU 航位推算几分钟即失真。
- 分辨率 ≠ 精度：误差预算看漂移与噪声密度项，不看刻度位数。
- 器件误差源与滤波/融合消费分工：漂移与噪声归器件侧，
control-note第 8–9 章管估计，接口与话题归ROS2-note。