简介imu-utils 是一套面向机器人导航、虚拟现实、飞行控制及移动设备开发者的 IMU 标定工具包用于消除加速度计、陀螺仪与磁力计的系统误差和随机误差提升惯性测量单元输出精度。资源包共 235 个文件约 4.32MB以 txt 说明文档、sample 样例数据、h/hpp 头文件、cpp 源码、yaml 与 launch 配置、m 脚本等为主并包含 code_utils 与 imu_utils 两大模块兼顾核心标定算法与数据处理、接口设计等辅助功能。内容预览可见多项式拟合、Allan 方差分析、加速度计与陀螺仪标定等实现代码中已修正错误引用稳定性更有保障。目前已有 438 人学习下载适合需要将 IMU 集成到项目中、希望减少标定技术障碍的工程师与研究者参考也可按自身硬件环境定制标定流程。1. imu-utils 到底解决什么问题从一次姿态漂移翻车说起如果你用过消费级 IMU 做姿态估计大概率遇到过这个场景设备静止放在桌上上位机里的横滚角却以每分钟零点几度的速度缓慢爬升十分钟后偏了好几度。你以为是滤波算法没调好换了三套 Madgwick 参数结果漂移依旧。真正的问题往往不在算法而在传感器本身的零偏、标度因数误差和轴间非正交——这些出厂校准没做干净后面怎么补都是治标。imu-utils 就是冲着这件事来的它是一套面向常见六轴/九轴 IMU 的标定工具集核心做三件事——零偏估计、标度因数与轴间正交校正、以及把标定结果导出成后续算法能直接吃的参数文件。它不绑定某一款芯片常见做法是配合 MPU 系列、ICM 系列、BMI 系列这类通过 I2C/SPI 读原始数据的模块使用。适合谁做机器人底盘、云台、可穿戴姿态追踪、以及任何需要「先标定再上算法」的嵌入式开发者。如果你现在还在拿未标定的原始数据直接喂给 EKF这篇值得往下看。2. 标定前必须搞懂的三个误差模型2.1 零偏、标度因数、非正交到底怎么影响输出IMU 的原始输出可以粗略写成这样一个模型测量值 真实值 × 标度因数 零偏 轴间耦合 噪声。零偏是传感器静止时输出不为零的那部分它随温度变化是姿态漂移的头号元凶。标度因数是「输入 1°/s 实际输出多少 LSB」的斜率误差表现为大角速度时比例失真。轴间非正交则是 X 轴转动时 Y、Z 轴也跟着有输出本质是三个敏感轴没有严格两两垂直。这三者里零偏最容易标也最容易被忽略标度因数和正交误差需要转台或至少精确的角度激励才能标得准。imu-utils 的设计思路是把这三项分开处理零偏用静止采样标度因数和正交用多位置法或六面法。理解这个分层后面看参数才不会懵。2.2 为什么六面法比单点静止标定更可靠单点静止只能估零偏估不了标度因数和正交。六面法把 IMU 依次以六个面朝下静置每个位置采集一段数据利用重力在各轴上的投影作为已知激励反解出完整的 3×3 校正矩阵加零偏向量。它的好处是激励覆盖了正负方向能同时约束标度因数和轴间耦合。代价是操作麻烦需要保证每个面真的水平。常见做法是用一个加工过的立方体夹具把 IMU 固定上去翻面时靠夹具保证垂直度。没有夹具的话靠手摆误差会直接进标定结果这也是很多人标完发现没改善的原因。2.3 温度漂移被低估的第二个变量零偏随温度变化消费级 IMU 的零偏温漂可以到每摄氏度零点几度每秒。imu-utils 支持在标定过程中记录温度做一阶线性补偿bias(T) bias0 k × (T - T0)。如果你的应用场景温度变化大比如户外机器人、车载只标常温零偏是不够的。实操上我一般会做一次升温标定让设备从冷启动开始边采集边等它自然升温记录零偏随温度的变化曲线拟合出 k。这一步不需要额外硬件只要耐心等十几分钟。3. 用 imu-utils 跑通一次完整标定从采集到导出3.1 环境准备与数据采集脚本imu-utils 通常以 Python 包或 C 库的形式提供Python 侧负责数据采集和参数拟合C 侧负责把参数烧进固件或运行时加载。先装依赖再写一个最小采集脚本。# collect_imu.py # 采集 IMU 原始数据用于后续标定 import serial import csv import time # 根据实际串口修改波特率常见 115200 或 921600 PORT /dev/ttyUSB0 BAUD 115200 DURATION 10.0 # 每个位置采集 10 秒 OUTPUT imu_raw.csv def read_line(ser): # 假设固件输出格式: ax,ay,az,gx,gy,gz,temp line ser.readline().decode(ascii, errorsignore).strip() if not line or line.startswith(#): return None parts line.split(,) if len(parts) 7: return None try: return [float(p) for p in parts[:7]] except ValueError: return None def main(): ser serial.Serial(PORT, BAUD, timeout1) time.sleep(2) # 等待串口稳定 with open(OUTPUT, w, newline) as f: writer csv.writer(f) writer.writerow([ax, ay, az, gx, gy, gz, temp]) t0 time.time() while time.time() - t0 DURATION: row read_line(ser) if row: writer.writerow(row) ser.close() print(f采集完成写入 {OUTPUT}) if __name__ __main__: main()这段脚本做的是最朴素的采集按行读串口解析七列数据写 CSV。关键参数是DURATION静止零偏采集建议至少 10 秒六面法每个面也建议 10 秒以上采样率 100Hz 的话每个位置有 1000 个点足够做平均。PORT和BAUD必须和你的硬件一致读不到数据先查这两个。3.2 零偏标定静止采样的参数与判据零偏标定最简单但判据要卡严。把 IMU 静止放好采集一段数据取均值就是零偏。但「静止」不是嘴上说的要看数据的标准差。# bias_calib.py import numpy as np import csv def load_csv(path): data [] with open(path) as f: reader csv.DictReader(f) for row in reader: data.append([float(row[k]) for k in [ax,ay,az,gx,gy,gz]]) return np.array(data) def estimate_bias(data, gyro_std_thresh0.5, accel_std_thresh0.05): # data: N x 6, 前3列加速度后3列角速度 accel data[:, :3] gyro data[:, 3:] a_std accel.std(axis0) g_std gyro.std(axis0) # 判据标准差过大说明没放稳标定无效 if np.any(g_std gyro_std_thresh) or np.any(a_std accel_std_thresh): print(警告数据抖动过大请重新静止采集) print(陀螺标准差:, g_std) print(加速度标准差:, a_std) return None bias_accel accel.mean(axis0) bias_gyro gyro.mean(axis0) return bias_accel, bias_gyro if __name__ __main__: data load_csv(imu_raw.csv) result estimate_bias(data) if result: ba, bg result print(加速度零偏:, ba) print(陀螺零偏:, bg)这里的判据是血泪经验很多人采集时手还没离开桌子就按了开始数据里混进了扰动均值直接偏掉。gyro_std_thresh设 0.5单位取决于你的量程如果是 °/s 就偏松rad/s 就偏紧accel_std_thresh设 0.05g 左右。超了就重采别硬用。3.3 六面法标定位置定义与矩阵求解六面法的核心是建立「每个位置下重力在机体坐标系的理论投影」和「实测加速度均值」之间的方程。设校正模型为 a_corrected M × (a_raw - b)其中 M 是 3×3 矩阵b 是零偏。六个位置提供 18 个方程未知数 12 个M 9 个 b 3 个超定可解。# six_face_calib.py import numpy as np # 六个面的理论重力方向机体坐标系单位 g # 顺序X朝下, -X朝下, Y朝下, -Y朝下, Z朝下, -Z朝下 GRAVITY_DIRS np.array([ [ 1, 0, 0], [-1, 0, 0], [ 0, 1, 0], [ 0,-1, 0], [ 0, 0, 1], [ 0, 0,-1], ], dtypefloat) def solve_calibration(measured_means): # measured_means: 6 x 3每个面实测加速度均值 # 模型: measured S (true) b求 S 和 b # 用最小二乘对每个轴独立拟合 A np.hstack([GRAVITY_DIRS, np.ones((6,1))]) # 6 x 4 S np.zeros((3,3)) b np.zeros(3) for axis in range(3): y measured_means[:, axis] coef, *_ np.linalg.lstsq(A, y, rcondNone) S[axis, :] coef[:3] b[axis] coef[3] return S, b if __name__ __main__: # 示例六个面实测均值实际从采集数据算 measured np.array([ [ 0.98, 0.02, 0.01], [-0.97, 0.03, 0.00], [ 0.01, 0.99, 0.02], [ 0.00,-0.98, 0.01], [ 0.02, 0.01, 0.97], [ 0.01, 0.00,-0.99], ]) S, b solve_calibration(measured) print(校正矩阵 S:\n, S) print(零偏 b:, b) # 理想情况下 S 接近单位阵非对角元反映轴间耦合GRAVITY_DIRS的顺序必须和你的实际翻转顺序严格对应错一个面结果全废。solve_calibration对每个轴做最小二乘得到该轴的标度因数和耦合系数。理想情况下 S 接近单位阵非对角元越小越好。如果非对角元超过 0.05说明夹具垂直度不够或者某个面没放平。3.4 导出参数并在固件里加载标定完要把 S 和 b 导出成固件能读的格式。常见做法是存成头文件或二进制结构体。// imu_calib.h 自动生成勿手改 #ifndef IMU_CALIB_H #define IMU_CALIB_H // 校正矩阵行优先 static const float IMU_CALIB_M[3][3] { { 1.0021f, 0.0032f, -0.0011f}, { 0.0028f, 0.9987f, 0.0015f}, {-0.0009f, 0.0012f, 1.0013f}, }; // 零偏单位与原始输出一致 static const float IMU_CALIB_B[3] { 0.012f, -0.008f, 0.021f }; // 运行时校正函数 static inline void imu_apply_calib(const float raw[3], float out[3]) { for (int i 0; i 3; i) { out[i] 0.0f; for (int j 0; j 3; j) { out[i] IMU_CALIB_M[i][j] * (raw[j] - IMU_CALIB_B[j]); } } } #endifimu_apply_calib先减零偏再乘矩阵顺序不能反。矩阵和零偏的数值由 Python 脚本生成直接替换即可。注意浮点精度嵌入式上如果用的是单精度标定矩阵的非对角元太小可能被截断建议保留四位小数。4. 标定结果验证怎么判断标对了没有4.1 静态验证残差与角度误差标定完不能只看参数好不好看要用独立数据验证。静态验证最简单重新采集一段静止数据用校正矩阵处理后看加速度模长是否接近 1g角速度是否接近零。# verify_static.py import numpy as np from bias_calib import load_csv def verify(data, S, b): accel data[:, :3] gyro data[:, 3:] accel_corr (S (accel - b).T).T gyro_corr gyro - b # 陀螺零偏单独处理 a_norm np.linalg.norm(accel_corr, axis1) print(校正后加速度模长均值: %.4f g (理想 1.0) % a_norm.mean()) print(校正后加速度模长标准差: %.4f % a_norm.std()) print(校正后陀螺均值: , gyro_corr.mean(axis0)) if __name__ __main__: data load_csv(verify_raw.csv) S np.array([[1.0021,0.0032,-0.0011], [0.0028,0.9987,0.0015], [-0.0009,0.0012,1.0013]]) b np.array([0.012,-0.008,0.021]) verify(data, S, b)判据加速度模长均值应在 1.0 ± 0.01g 内标准差小于 0.005g。陀螺均值应接近零残余超过 0.1°/s 说明零偏没标干净。这些数字不是绝对的取决于你的传感器等级但趋势必须对。4.2 动态验证转台或手动旋转的一致性没有转台的话手动绕各轴匀速旋转看积分角度和实际转过的角度是否一致。比如绕 Z 轴转 90 度陀螺积分应该也是 90 度左右。误差在 2% 以内算合格超过 5% 说明标度因数没标好。这一步能暴露静态验证发现不了的问题建议必做。5. 避坑与排查标定翻车的五个典型场景5.1 现象标定后漂移反而变大原因通常是六面法的位置顺序和代码里的GRAVITY_DIRS不匹配或者某个面采集时设备没放稳。解决方法是重新核对顺序并在每个面采集后检查标准差超限就重采。别嫌麻烦这一步省不得。5.2 现象加速度模长始终偏离 1g如果校正后模长还是偏先查量程设置。很多 IMU 默认 ±2g但代码里按 ±8g 换算比例直接错四倍。确认寄存器配置和换算系数一致再查标定矩阵是否被正确加载。5.3 现象温度变化后零偏重新出现常温标定的零偏在温度变化后失效说明需要温度补偿。检查 imu-utils 是否开启了温度记录拟合 k 时数据点是否覆盖了工作温度范围。只标一个温度点是不够的。5.4 现象串口数据丢包导致标定结果抖动高采样率下串口丢包很常见表现为数据行数比预期少或者时间戳不连续。解决方法是降低采样率、提高波特率或者在固件里加缓冲。采集脚本里可以加一个行数校验少于预期 90% 就报警。5.5 现象校正矩阵非对角元异常大非对角元超过 0.1 通常不是传感器问题而是夹具或操作问题。检查夹具是否真的垂直翻面时有没有碰到线缆。如果反复标都是这样可能是传感器本身轴间耦合严重这种模块建议直接换。6. 进阶技巧把标定参数做成运行时自检标定不是一劳永逸的温度、震动、老化都会让参数漂移。我一般会在固件里加一个运行时自检设备静止超过 5 秒时自动计算当前零偏和标定值的偏差超过阈值就置一个标志位提示需要重新标定。// runtime_check.c #include math.h #define BIAS_DRIFT_THRESH 0.05f // 零偏漂移阈值单位同原始输出 #define STATIC_TIME_MS 5000 static float bias_ref[3] {0.012f, -0.008f, 0.021f}; static uint32_t static_start 0; void imu_runtime_check(const float gyro[3], uint32_t now_ms) { float mag sqrtf(gyro[0]*gyro[0] gyro[1]*gyro[1] gyro[2]*gyro[2]); if (mag 0.02f) { // 认为静止 if (static_start 0) static_start now_ms; if (now_ms - static_start STATIC_TIME_MS) { for (int i 0; i 3; i) { if (fabsf(gyro[i] - bias_ref[i]) BIAS_DRIFT_THRESH) { // 置标志位上层决定是否重新标定 imu_need_recalib 1; } } } } else { static_start 0; } }这个自检逻辑很轻不占多少资源但能帮你提前发现参数失效。阈值BIAS_DRIFT_THRESH根据你的传感器噪声水平调太紧会误报太松没意义。我一般设成标定残差标准差的 3 倍左右。另一个技巧是把标定参数和温度绑定做成查找表。如果工作温度范围窄线性补偿够用范围宽的话分段线性或查表更稳。imu-utils 本身不强制你怎么存但导出格式留了扩展字段可以塞温度系数。最后说个习惯每次改完硬件换模块、改安装方式、动线缆重新标一次。别信「上次标过了」安装应力变化会直接改变零偏。这个习惯帮我省了很多次现场调试的后悔药。希望帮到你。本文还有配套的精品资源点击获取