《四元数解算姿态完全解析及资料汇总.docx》由会员分享，可在线阅读，更多相关《四元数解算姿态完全解析及资料汇总.docx（19页珍藏版）》请在金锄头文库上搜索。
1、四元数完全解析及资料汇总本文原帖出自匿名四轴论坛，附件里的资源请到匿名论坛下载：http:/ 好，废话到此为止，开始说主体。四元数和姿态角怎么说呢？先得给和我一样的小菜鸟们理一理思路，小鸟我在此画了一个“思维导图”（我承认我画的丑），四元数解算姿态首先分为两部分理解：第一部分先理解什么是四元数，四元数与姿态角间的关系；第二部分要理解怎么由惯性单元测出的加速度和角速度求出四元数，再由四元数求出欧拉角。图1 渣渣思维导图在讲解什么是四元数时，小弟的思维是顺着说的，先由四元数的定义说起，说到四元数与姿态角间的关系。但在讲解姿态解算时，小弟的思维是逆向的，就是反推回来的，从欧拉角一步步反推回到惯性器件
2、的测量数据，这样逆向说是因为便于理解，因为实际在工程应用时和理论推导有很大差别。实际应用时正确的求解顺序应该为图1中序号顺序，即1-2-3-.但在笔者讲解姿态求解时思路是如图2的。图2 逆向讲解思路大家在看四元数时最好结合着代码一块看，小弟看的是匿名四轴的代码，感觉写的非常好也非常清晰，粘出来大家一块观摩。红色部分是核心代码，总共分为八个步骤，和图1中的八个步骤是一一对应的。讲解介绍时也是和代码对比起来讲解的。代码可以去匿名官网上下载，都是开源的，不是小弟的，所以小弟不方便加在附件中。/四元数更新姿态#define Kp 2.0f /加速度权重，越大则向加速度测量值收敛越快#define Ki
3、 0.001f/误差积分增益void ANO_IMU:Quaternion_CF(Vector3f gyro,Vector3f acc, float deltaT)Vector3f V_gravity, V_error, V_error_I;/1.重力加速度归一化acc.normalize();/2.提取四元数的等效余弦矩阵中的重力分量Q.vector_gravity(V_gravity);/3.向量叉积得出姿态误差V_error = acc % V_gravity;/4.对误差进行积分V_error_I += V_error * Ki;/5.互补滤波，姿态误差补偿到角速度上，修正角速度积分漂
4、移Gyro += V_error * Kp + V_error_I;/6.一阶龙格库塔法更新四元数Q.Runge_Kutta_1st(Gyro, deltaT);/7.四元数归一化Q.normalize();/8.四元数转欧拉角Q.to_euler(&angle.x, &angle.y, &angle.z);好的，下面搬砖开始！。嘿咻嘿咻！一 什么是四元数？关于四元数的定义摘自秦永元的惯性导航，里面有非常好的讲解，大家可以直接看绪论和第九章就可以。下面我粘贴了部分原文，粘贴的比较多比较详细，应为本人比较笨还爱较真，所以按本人的风格就要详尽一点，大牛们都可以自动忽略。四元数定义、表达方式及运算方
5、法摘自惯性导航-秦永元P289-292好，关于四元数定义就搬这么多，其他的大家去附件下载惯性导航的pdf自己看吧。下面开始搬四元数与姿态解算关系的。嘿咻嘿咻二、四元数与姿态阵间的关系从上面我们知道了四元数的定义，可这四个数和我们要求的三个姿态角有什么关系呢？下面是详细的推导，同样摘自惯性导航-秦永元P292-297。四元数与姿态阵间的关系摘自惯性导航-秦永元P292-297 呃，粘了这么多其实就是为了想知道推导过程小伙伴好理解，真正有用的就是（9.2.34）这个公式。这个公式把四元数转换成了方向余弦矩阵中的几个元素，再用这几个元素转换为欧拉角。就求解除了姿态！先从四元数q0q3转成方向余弦矩阵
6、：xbybzb=q02+q12-q22-q322(q1q2-q0q3)2(q1q3+q0q2)2(q1q2+q0q3)q02-q12+q22-q322(q2q3-q0q1)2(q1q3-q0q2)2(q2q3+q0q1)q02-q12-q22+q32xnynzn再从方向余弦矩阵转换为欧拉角：=sin-1（T32）俯仰角pitch=tan-1（-T31T33）横滚角roll=tan-1（T12T22）航向角yaw好的，公式原理都讲清楚了，下面来看一下匿名的代码：/四元数转欧拉角，这里四元数是q1q4 和公式里q0q3相对应。void Quaternion:to_euler(float *roll
7、, float *pitch, float *yaw) if (roll) *roll = degrees(atan2f(2.0f*(q1*q2 + q3*q4),1 - 2.0f*(q2*q2 + q3*q3);/*roll = degrees(atan2f(-2.0f*(q2*q4 - q1*q3),1 - 2.0f*(q2*q2 + q3*q3); if (pitch) / 使用safe_asin()来处理pitch接近90/-90时的奇点 *pitch = degrees(safe_asin(2.0f*(q1*q3 - q2*q4);/*pitch = degrees(safe_asi
8、n(2.0f*(q3*q4 + q1*q2); if (yaw) *yaw = degrees(atan2f(2.0f*(q2*q3 - q1*q4), 2.0f*(q1*q1 + q2*q2) - 1);/*yaw = degrees(atan2f(2.0f*(q2*q3 - q1*q4), 2.0f*(q1*q1 + q3*q3) - 1); 对比余弦矩阵转换为欧拉角的公式很容易理解了吧，注意一下，红色是匿名原版的代码，和公式还是有出入的，自己细心观察一下吧。被注释了的黑色代码是我根据上面的公式写的。但黑色的求解出来的欧拉角反映出来的姿态是不对的，具体表现为俯仰（pitch）和横滚（rol
9、l）是相反的，为啥根据公式写的是不对的？群里的小伙伴给与了我热心的解答。这个错误主要是由于方向余弦矩阵的旋转顺序不一样，也就是公式（9.2.39）不一样，这是由于旋转的先后顺序不同引起的，具体大家参考惯性导航绪论来看就能明白，因为这一点小弟还有点混乱，就点到这为止。以上就是四元数求解欧拉角的方法。通过公式可以看到，要求欧拉角需要求得四元数的方向余弦矩阵；要求得四元数方向余弦矩阵，需要求得四元数q0q3，那么如何求得q0q3？接下来详细介绍。三、四元数微分及龙格库塔求Q0Q3 首先我们先来看一下在程序里如何求解的q0q3：/一阶龙格库塔法更新四元数void Quaternion:Runge_Ku
10、tta_1st(Vector3f &g, float deltaT) q1 += 0.5 * (-q2*g.x - q3*g.y - q4*g.z)* deltaT; q2 += 0.5 * (q1*g.x + q3*g.z - q4*g.y)* deltaT; q3 += 0.5 * (q1*g.y - q2*g.z + q4*g.x)* deltaT; q4 += 0.5 * (q1*g.z + q2*g.y - q3*g.x)* deltaT;这就是一阶龙格库塔法求解q的微分方程，传入参数只需要这个周期的角速度g.x、g.y、g.z和周期时间deltaT。下面一张是从某位大神的贴吧上盗的
11、图，描绘的是一阶龙格库塔的计算式。相信很多人和我一样，单看上图很难理解其中的意思和其由来，于是我又找了很多帖子，感谢前人做出的贡献，小弟在这里再次整理大神的四元数微分方程推导公式，便于大家理解。摘自附件中推導_四元數.pdf虽然在下也不是很懂，不过粘出来还是能起到理解的作用，这样大家就不会觉得这是凭空变出来的，本人数学功底薄弱，没有对推导进行过验证，如果有不对的地方欢迎指正。接着使用一阶龙格库塔（Runge-Kutta）发求出q0q3，这一点很多人不知道一阶龙格库塔怎么推导的，下面也是这位网友的推导，大家参考着理解吧。这里的角速度wx、wy、wz是由捷联陀螺的输出（对机械转子陀螺必须经过误差补
12、偿，将在下面介绍）。对比着匿名四轴的代码看一看（g.x、g.y、g.z是捷联陀螺的输出），代码的意思就比较清楚了。在往上一步步推，我们就要求陀螺输出了，并且还要对数据进行互补滤波处理。四、惯性单元测量值融合这部分看似很简单，但是也有让笔者难以理解的地方，希望后人能补充修正进行更好的讲解。有了上一步的龙格库塔方程，我们现在需要的就是角速度的测量值。在四轴上安装陀螺仪，可以测量四轴倾斜的角速度，由于陀螺仪输出的是四轴的角速度，不会受到四轴振动影响。因此该信号中噪声很小。四轴的角度又是通过对角速度积分而得，这可进一步平滑信号，从而使得角度信号更加稳定。因此四轴控制所需要的角度和角速度可以使用陀螺仪所
13、得到的信号。由于从陀螺仪的角速度获得角度信息，需要经过积分运算。如果角速度信号存在微小的偏差，经过积分运算之后，变化形成积累误差。这个误差会随着时间延长逐步增加，最终导致电路饱和，无法形成正确的角度信号。如何消除这个累积误差呢？可以通过上面的加速度传感器获得的角度信息对此进行校正。利用加速度计所获得的角度信息 g 与陀螺仪积分后的角度 进行比较，将比较的误差信号经过比例Tg 放大之后与陀螺仪输出的角速度信号叠加之后再进行积分。对于加速度计给定的角度g ，经过比例、积分环节之后产生的角度必然最终等于g 。由于加速度计获得的角度信息不会存在积累误差，所以最终将输出角度中的积累误差消除了。加速度计所
14、产生的角度信息g 中会叠加很强的有四轴运动加速度噪声信号。为了避免该信号对于角度 的影响，因此比例系数 Tg 应该非常小。这样，加速度的噪声信号经过比例、积分后，在输出角度信息中就会非常小了。由于存在积分环节，所以无论比例Tg多么小，最终输出角度必然与加速度计测量的角度g相等，只是这个调节过程会随着Tg 的减小而延长。先把这个过程的代码粘出来，看着代码一步步理解：#define Kp 2.0f /加速度权重，越大则向加速度测量值收敛越快#define Ki 0.001f /误差积分增益/1.重力加速度归一化acc.normalize();/2.提取四元数的等效余弦矩阵中的重力分量Q.vector_gravity(V_gravity);/3.向量叉积得出姿态误差V_error
