@



一、旋转运动学

1.1 线速度与角速度

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_加速度计

1.2 旋转坐标系下的运动学

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_离散化_02

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_03补充:右扰动模型:

\(R[w]_{X}=[R w]_{X} \cdot R\)

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_离散化_04VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_加速度计_05在旋转坐标系下观察,运动的物体(运动方向和旋转轴不为同一个轴时)会受到科氏力的作用。

二、IMU 测量模型及运动模型

2.1 MEMS(Micro-electromechanical Systems微电子机械系统) 加速度计工作原理

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_离散化_06VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_离散化_07

加速度变化,引起质量块位置的变化,引起电阻电容的变化。

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_08

2.2 MEMS(Micro-electromechanical Systems微电子机械系统) 陀螺仪测量原理

• 陀螺仪主要用来测量物体的旋转角速度,按测量原理分有振动陀螺,光纤陀螺等。

• 低端 MEMS 陀螺上一般采用振动陀螺原理,通过测量 Coriolis force(科氏力) 来间接得到角速度。

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_09VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_加速度计_10VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_公式推导_11VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_12加速度计当然也会受到科氏力影响,但是因为加速度计的质量块并不含有主动驱动的旋转运动。因此,一般情况下ω->0, 在精度要求不高的情况下,科氏力对加速度的影响可以忽略不记。

三、IMU 误差模型

3.1 误差分类

• 加速度计和陀螺仪的误差可以分为:确定性误差,随机误差。

• 确定性误差可以事先标定确定,包括: bias, scale ...

• 随机误差通常假设噪声服从高斯分布,确定其方差大小,包括:高斯白噪声, bias随机游走...

3.2 确定性误差

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_13VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_14

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_加速度计_15VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_离散化_16

3.2.1 六面法标定加速度 bias 和 scale factor

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_公式推导_17VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_陀螺仪_183.2.2 六面法标定陀螺仪 bias 和 scale factor

和加速度计六面法类似,只是陀螺仪的真实值由高精度转台提供,这

里的 6 面是指各个轴顺时针和逆时针旋转。

3.2.3 温度相关的参数标定

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_陀螺仪_19

3.3 随机误差

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_加速度计_20

3.3.1 随机误差的离散化

实际上, IMU 传感器获取的数据为离散采样,因此需要进行离散化

1. 高斯白噪声的离散化

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_陀螺仪_21VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_公式推导_22VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_陀螺仪_23

2. bias 随机游走的离散化

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_公式推导_24VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_25更详细的连续到离散的推导参见:

John L Crassidis. “Sigma-point Kalman fltering for integrated GPS and inertial navigation”. In: IEEE Transactions on Aerospace and Electronic Systems 42.2 (2006), pp. 750–756.

3.3.2 IMU 随机误差的标定

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_陀螺仪_26VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_陀螺仪_27详细可参考:

Allan Variance. “Noise Analysis for Gyroscopes”. In: Freescale Semiconductor Document Number: AN5087 Application Note Rev. 0 2 (2015)

3.4 加速度计的误差模型总结

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_28

3.5 陀螺仪的误差模型总结

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_29详细内容参考:MA Shelley. “Monocular visual inertial odometry on a mobile device”. In: Master’s thesis, Institut für Informatik, TUMünchen, Germany (2014)

四、运动模型(预积分模型)离散时间处理

4.1 IMU 模型

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_离散化_30

4.2 连续时间下 IMU 运动模型

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_31

4.3 运动模型的离散积分——欧拉法

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_离散化_32\({\mathbf{q}}_{w b_{t+1}}={\mathbf{q}}_{w b_{t}}+\dot{\mathbf{q}}_{w b_{t}}\times\Delta t=\mathbf{q}_{w b_{t}} \otimes\left[\begin{array}{c}1 \\ 0\end{array}\right]+\mathbf{q}_{w b_{t}} \otimes\left[\begin{array}{c}0 \\ \frac{1}{2} \omega^{b_{t}}\end{array}\right]=\mathbf{q}_{w b_{k}} \otimes\left[\begin{array}{c}1 \\ \frac{1}{2} \omega \delta t \end{array}\right]\)

4.3 运动模型的离散积分——中值法

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_33

4.4 李代数表示的运动模型

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_34VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_工作原理_35

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_加速度计_36VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_离散化_37

VIO(2)—— IMU传感器测量模型、误差模型、运动模型(预积分模型)_公式推导_38

太多了就不一一放上来了,重点看文档《预积分总结与公式推导》