关于姿态解算与融合的代码注释篇（三）

https://blog.csdn.net/sinat_23338865/article/details/53886114

加速度计和陀螺仪都能计算出姿态，但为何要对它们融合呢，是因为加速度计对振动之类的扰动很敏感，但长期数据计算出的姿态可信，而陀螺仪虽然对振动这些不敏感，但长期使用陀螺仪会出现漂移，因此我们要进行互补，短期相信陀螺，长期相信加计。不过，其实加计无法对航向角进行修正，修正航向角需要磁力计，在融合之前先要对传感器原始数据进行一些处理。理想情况下，加速度计水平放置时，XY轴应该是0输出的，仅Z轴输出1个G，因此，我们需要对加速度计进行XY轴的零点校准；同样的，陀螺仪在水平静止放置时各轴输出应为0，因此需对陀螺仪进行三轴的校准。方法就是把机体标准水平静止放置时采集它个一两百次数据求个平均作为校准值保存起来，然后工作状态下各轴输出的数据就是采集来的数据减去校准值。仅此还不够，陀螺仪不进行滤波还可以接受，但加速度计噪声比较大，所以需要来个滑动窗口滤波，窗口深度为8-10即可。

接下来上六轴数据融合代码：（基础上再加些注释）#defineKp 10.0f                        // 这里的KpKi是用于调整加速度计修正陀螺仪的速度#defineKi 0.008f                        #definehalfT 0.001f             // 采样周期的一半，用于求解四元数微分方程时计算角增量floatq0 = 1, q1 = 0, q2 = 0, q3 = 0;    // 初始姿态四元数，由上篇博文提到的变换四元数公式得来floatexInt = 0, eyInt = 0, ezInt = 0;    //当前加计测得的重力加速度在三轴上的分量//与用当前姿态计算得来的重力在三轴上的分量的误差的积分voidIMUupdate(float gx, float gy, float gz, float ax, float ay, float az)   //g表陀螺仪，   a表加速度计{float q0temp,q1temp,q2temp,q3temp;//四元数暂存变量，求解微分方程时要用float norm; //矢量的模或四元数的范数float vx, vy, vz;//当前姿态计算得来的重力在三轴上的分量float ex, ey, ez;//当前加计测得的重力加速度在三轴上的分量//与用当前姿态计算得来的重力在三轴上的分量的误差// 先把这些用得到的值算好float q0q0 = q0*q0;float q0q1 = q0*q1;float q0q2 = q0*q2;float q1q1 = q1*q1;float q1q3 = q1*q3;float q2q2 = q2*q2;float q2q3 = q2*q3;float q3q3 = q3*q3;if(ax*ay*az==0)//加计处于自由落体状态时不进行姿态解算，因为会产生分母无穷大的情况return;norm = sqrt(ax*ax + ay*ay + az*az);//单位化加速度计，ax = ax /norm;// 这样变更了量程也不需要修改KP参数，因为这里归一化了ay = ay / norm;az = az / norm;//用当前姿态计算出重力在三个轴上的分量，//参考坐标n系转化到载体坐标b系的用四元数表示的方向余弦矩阵第三列即是（博文二中有提到）vx = 2*(q1q3 - q0q2);                                                                                                             vy = 2*(q0q1 + q2q3);vz = q0q0 - q1q1 - q2q2 + q3q3 ;//计算测得的重力与计算得重力间的误差，向量外积可以表示这一误差//原因我理解是因为两个向量是单位向量且sin0等于0//不过要是夹角是180度呢~这个还没理解ex = (ay*vz - az*vy) ;                                                                  ey = (az*vx - ax*vz) ;ez = (ax*vy - ay*vx) ;exInt = exInt + ex * Ki;                                           //对误差进行积分eyInt = eyInt + ey * Ki;ezInt = ezInt + ez * Ki;// adjusted gyroscope measurementsgx = gx + Kp*ex + exInt;  //将误差PI后补偿到陀螺仪，即补偿零点漂移gy = gy + Kp*ey + eyInt;gz = gz + Kp*ez + ezInt;    //这里的gz由于没有观测者进行矫正会产生漂移，表现出来的就是积分自增或自减//下面进行姿态的更新，也就是四元数微分方程的求解q0temp=q0;//暂存当前值用于计算q1temp=q1;//网上传的这份算法大多没有注意这个问题，在此补上q2temp=q2;q3temp=q3;//采用一阶毕卡解法q0 = q0temp + (-q1temp*gx - q2temp*gy -q3temp*gz)*halfT;q1 = q1temp + (q0temp*gx + q2temp*gz -q3temp*gy)*halfT;q2 = q2temp + (q0temp*gy - q1temp*gz +q3temp*gx)*halfT;q3 = q3temp + (q0temp*gz + q1temp*gy -q2temp*gx)*halfT;//单位化四元数在空间旋转时不会拉伸，仅有旋转角度，这类似线性代数里的正交变换norm = sqrt(q0*q0 + q1*q1 + q2*q2 + q3*q3);q0 = q0 / norm;q1 = q1 / norm;q2 = q2 / norm;q3 = q3 / norm;//四元数到欧拉角的转换，公式推导见博文二//其中YAW航向角由于加速度计对其没有修正作用，因此此处直接用陀螺仪积分代替Q_ANGLE.Z = GYRO_I.Z; // yawQ_ANGLE.Y = asin(-2 * q1 * q3 + 2 * q0* q2)*57.3; // pitchQ_ANGLE.X = atan2(2 * q2 * q3 + 2 * q0 * q1,-2 * q1 * q1 - 2 * q2* q2 + 1)* 57.3; // roll}

关于姿态解算与融合的代码注释篇（三）相关推荐

imu姿态解算+卡尔曼滤波融合JAVA版（此版本卡拉曼滤波奇点有错误）
原版地址:IMU9轴卡尔曼滤波增加mpu6050 陀螺仪零飘矫正,imu算法优化 KalmAndAndIMU 类: import java.util.Vector;public class Kalm ...
基于STM32F407四旋翼无人机 --- 姿态解算讲解（四元数）（叉积法融合陀螺仪数据和加速度数据）（五）
基于STM32F407四旋翼无人机 --- 姿态解算讲解(五) 姿态解算姿态解算定义欧拉角四元数四元数性质方向余弦矩阵四元数方向余弦矩阵叉积法融合陀螺仪数据和加速度数据叉积运算一阶龙 ...
四周无人机的姿态解算（2）
接上一篇 4.尽量直观地分析一下程序直接把小飞机的程序考过来了... 说明: 函数名:姿态更新输入:mpu6050结构体,也就是测得数据+初始偏置输出:三个姿态角数据函数变量的结构体定义: typ ...
基于STM32的四旋翼无人机项目（二）：MPU6050姿态解算（含上位机3D姿态显示教学）
前言:本文为手把手教学飞控核心知识点之一的姿态解算--MPU6050 姿态解算(飞控专栏第2篇).项目中飞行器使用 MPU6050 传感器对飞行器的姿态进行解算(四元数方法),搭配设计的卡尔曼滤波器与 ...
姿态解算(用于飞行器或ROV等)的基本思路和数学原理
文章目录 (一)前言 (二)方向余弦矩阵DCM (三)四元数 (四)求四元数更新四元数误差消除 (五)总结 (一)前言数学是算法的门槛我们就用数学的角度去探寻姿态解算,不谈代码主要是 ...
matlab求逆矩阵_MPU6050姿态解算2-欧拉角amp;旋转矩阵
1 IMU姿态解算 IMU,即惯性测量单元,一般包含三轴陀螺仪与三轴加速度计.之前的文章码农爱学习:MPU6050姿态解算方式1-DMPzhuanlan.zhihu.com 已将对MPU6050这 ...
四旋翼无人机飞控系统设计(姿态解算)
姿态解算姿态传感器读出加速度和角速度,而对一个系统的自动控制往往需要更加上层和贴近应用的的一个属性:角度.所以需要通过加速度和角速度进行数据融合转化得到姿态角度. 以MPU6050为例,姿态 ...
姿态解算四元数、方向余弦、欧拉角、Mahony滤波、四轴
姿态解算四元数.方向余弦.欧拉角.Mahony滤波说明:本文只是做了一些总结,需要一些对这方面的基础概念的了解. 一般人千万不要试图去深入探讨四元数 1. 方向余弦矩阵方向余弦矩阵是使用欧拉角( ...
四元数姿态解算及多传感器融合详细解析
代码路径ardupolit/modules/PX4Firmware/src/modules/attitude_estimator_so3/attitude_estimator_so3_main.cpp ...

关于姿态解算与融合的代码注释篇（三）

关于姿态解算与融合的代码注释篇（三）相关推荐

最新文章

热门文章