无人机滤波算法之——互补滤波原理

兰博

无人机滤波算法之——互补滤波原理 [复制链接]

互补滤波原理：  

     在四轴入门理论知识那节我们说，加速度计和磁传感器都是极易受外部干扰的传感器，都只能得到2维的角度关系，但是测量值随时间的变化相对较小，结合加速度计和磁传感器可以得到3维的角度关系。陀螺仪可以积分得到三维的角度关系，动态性能好，受外部干扰小，但测量值随时间变化比较大。可以看出，它们优缺点互补，结合起来才能有好的效果。那么三者的数据如何融合呢，接下来介绍互补滤波算法。

     互补滤波就是在短时间内采用陀螺仪得到的角度做为最优，定时对加速度采样来的角度进行取平均值来校正陀螺仪的得到的角度。即短时间内用陀螺仪比较准确，以它为主；长时间用加速度计比较准确，这时候加大它的比重，这就是互补了，不过滤波在哪里？加速度计要滤掉高频信号，陀螺仪要滤掉低频信号，互补滤波器就是根据传感器特性不同，通过不同的滤波器（高通或低通，互补的），然后再相加得到整个频带的信号，例如，加速度计测倾角，其动态响应较慢，在高频时信号不可用，所以可通过低通抑制高频；陀螺响应快，积分后可测倾角，不过由于零漂等，在低频段信号不好。通过高通滤波可抑制低频噪声。将两者结合，就将陀螺和加表的优点融合起来，得到在高频和低频都较好的信号，互补滤波需要选择切换的频率点，即高通和低通的频率。

互补滤波原理框图：

加速度计补偿：

    假设n系为地理坐标系，b系为机体坐标系，在地理坐标系中，加速度的输出为： ，经过 矩阵转换后的值为： 。在b系中，加速度的测量值为： ，现在 和 都表示在b系中数值向下的向量，由此，我们对这两个向量做向量积（叉积），得到误差： ，利用这个误差来修正 矩阵，于是乎，我们的四元数就在这样一个过程中被修正了。

但是，由于加速度计无法感知Z轴上的旋转运动，所以还需要用地磁计来进一步补偿。现在我们假设旋转矩阵 是经过加速度计校正后的矩阵，当某个确定的向量（b系中）经过这个矩阵旋转之后（到n系），这两个坐标系在XOY平面上重合，只是在Z轴旋转上会存在一个偏航角的误差。下图表示的是经过 旋转之后的b系和n系的相对关系。可以明显发现加速度计可以把b系通过四元数法从任意角度拉到与n系水平的位置上，这时，只剩下一个偏航角误差。这也是为什么加速度计误差修正偏航的原因。

地磁计补偿：

现在我们反过来从b系到n系，假设地磁计在b系中的输出为： ，经过 的转换后得到: 。由于 是经过加速度计修正过的旋转矩阵，因此该旋转矩阵只在Z轴上存在一个偏航的误差，这就导致 中的hy不为零。如果 不存在误差，这里的hy应该为0。在n系中，地磁方向与x轴呈一个角度，与z轴呈一个角度，这里我们让x轴对准北边，那么地磁向量为： 。在n系的XOY平面上（水平面）， 的投影为：Sqrt（bx * bx）， 的投影为：Sqrt（(hx * hx)+ (hy * hy)）。由于 存在的偏航误差，导致 的hy不为零，这就是说现在得到的 是真实的地磁向量绕Z轴旋转一定的角度后得到的。但由于是绕Z轴旋转，所以该地磁向量在XOY平面上（n系）投影的大小必定相同，所以有bx^2= hx^2+hy^2，我们求出了地磁向量在X轴方向的真实值。而 得到hz就是地磁向量在Z轴方向上的真实值，我们不做改变，令bz=hz即可。经过这样处理之后我们得到 ，这个地磁向量就是地磁的真实值，类似于重力加速度的 。接着把 经过变换后到b系中得到： 。 与 向量积求误差，再次修正 这样就完成了一次地磁计的补偿。

将加速度计没能做到的z轴上的旋转修正，通过地磁计在XOY平面上的地磁力相同原理，得到了修正。于是乎，Pitch和Roll通过加速度计修正，然后在这个基础之上（该地磁计补偿方法必须依靠加速度计修正提供一致的XOY平面，才会有bx^2= hx^2+hy^2等式成立），Yaw通过地磁计来补偿，最终得到了没有偏差的实时姿态（也就是由四元数组成的旋转矩阵）。

月下良缘

谢谢分享！