在OpenGL图形编程中，处理3D物体旋转是一个核心且复杂的课题。为了解决传统欧拉角方法中存在的万向节死锁（Gimbal Lock）问题，并使用更平滑、更稳定的方式进行旋转插值，四元数（Quaternion）已成为现代图形引擎中表示和计算旋转的标准工具。其数学本质是一个四维复数，形式为q = [w, (x, y, z)]，其中w是标量部分，(x, y, z)是矢量部分，代表了旋转轴的方向。
一、核心原理：为什么使用四元数？
相较于欧拉角和旋转矩阵，四元数在表示3D旋转时具有显著优势，如下表所示：
| 表示方法 | 优点 | 缺点 | 主要应用场景 |
|---|---|---|---|
| 欧拉角 | 直观（俯仰、偏航、滚转），易于人类理解。 | 存在万向节死锁；插值不平滑。 | 简单的摄像机控制、飞行器姿态显示。 |
| 旋转矩阵 | 表示明确，可直接用于顶点变换。 | 9个参数，冗余度高；插值困难。 | 基本的坐标变换。 |
| 四元数 | 仅有4个参数，无冗余；无万向节死锁；可实现平滑的球面线性插值(SLERP)。 | 数学概念抽象，不易直观理解。 | 复杂的3D动画、骨骼蒙皮、摄像机环绕、物理仿真。 |
万向节死锁是欧拉角的一个致命缺陷，当俯仰角（Pitch）为±90度时，偏航（Yaw）和滚转（Roll）会失去一个自由度，导致旋转行为异常。四元数从数学结构上避免了这一问题，它通过绕一个单一轴旋转任意角度来定义，从而保证了旋转空间的完备性。
二、关键操作与实现
在OpenGL中，我们需要将四元数转换为4x4旋转矩阵，然后与模型视图矩阵（Model-View Matrix）相乘，才能最终作用于顶点着色器中的顶点位置。
1. 从轴-角到四元数
给定一个旋转轴v = (x, y, z)（需为单位向量）和一个旋转角度θ（弧度制），可以构造相应的四元数：
w = cos(θ/2) x = sin(θ/2) * axis.x y = sin(θ/2) * axis.y z = sin(θ/2) * axis.z此公式的几何意义是将旋转信息编码到四维空间的一个点上。
2. 四元数转换为旋转矩阵
OpenGL的着色器最终处理的是矩阵。一个归一化四元数q = [w, (x, y, z)]可以转换为以下3x3旋转矩阵R（可以轻松扩展为4x4齐次坐标矩阵）：
R = [ 1 - 2y² - 2z², 2xy - 2wz, 2xz + 2wy, 2xy + 2wz, 1 - 2x² - 2z², 2yz - 2wx, 2xz - 2wy, 2yz + 2wx, 1 - 2x² - 2y² ]3. 四元数插值 (SLERP)
这是四元数最强大的特性之一，用于在两个旋转之间生成平滑的过渡动画。球面线性插值（Spherical Linear Interpolation, SLERP）公式如下：
SLERP(q0, q1, t) = ( sin((1-t)Ω) * q0 + sin(tΩ) * q1 ) / sin(Ω)其中，t是插值参数（0到1），Ω是q0与q1之间的夹角（通过点积求得）。这保证了旋转在四维球面上沿最短路径匀速移动，动画效果自然。
三、代码实现示例（使用C++和GLM库）
GLM（OpenGL Mathematics）是一个被广泛使用的数学库，它提供了完整的四元数支持。
#include <glm/glm.hpp> #include <glm/gtc/matrix_transform.hpp> #include <glm/gtc/quaternion.hpp> #include <glm/gtx/quaternion.hpp> int main() { // 1. 创建四元数：绕Y轴旋转90度 glm::vec3 rotationAxis(0.0f, 1.0f, 0.0f); // 旋转轴：Y轴 float rotationAngle = glm::radians(90.0f); // 旋转角度：转换为弧度 glm::quat myQuaternion = glm::angleAxis(rotationAngle, rotationAxis); // 2. 将四元数转换为旋转矩阵 glm::mat4 rotationMatrix = glm::mat4_cast(myQuaternion); // 关键转换函数 // 3. 组合模型矩阵：将旋转应用到模型上 glm::mat4 model = glm::mat4(1.0f); // 初始化为单位矩阵 model = model * rotationMatrix; // 应用旋转 // 4. 四元数插值示例 (SLERP) glm::quat startQuat = glm::angleAxis(glm::radians(0.0f), glm::vec3(0, 1, 0)); glm::quat endQuat = glm::angleAxis(glm::radians(180.0f), glm::vec3(0, 1, 0)); float t = 0.5f; // 中间点（50%进度） glm::quat interpolatedQuat = glm::slerp(startQuat, endQuat, t); // 执行SLERP插值 // 5. 在渲染循环中，将最终的model矩阵传递给着色器 // shader.setMat4("model", model); return 0; }以上代码演示了使用GLM库进行四元数创建、转换和插值的基本流程。glm::mat4_cast函数是执行四元数到旋转矩阵转换的核心。
四、应用场景与问题解决
四元数在实际开发中解决了诸多难题，一个典型应用是触屏控制3D物体旋转。当用户用手指在屏幕上滑动时，需要将2D触摸位移转换为3D空间的旋转。直接使用欧拉角计算会导致旋转轴混乱和物体“翻转”的失真现象。
解决方案是利用四元数增量式地更新旋转状态：
- 根据当前帧与上一帧的触摸点位移，计算出一个屏幕空间的旋转向量。
- 将该向量映射到一个基于当前摄像机朝向的3D旋转轴上。
- 根据该轴和位移大小（作为角度因子）构造一个增量旋转四元数。
- 将增量四元数与物体当前的总旋转四元数相乘（四元数乘法，表示旋转的组合），得到新的旋转状态。
- 将最终的四元数转换为矩阵并应用。
这种方法避免了欧拉角的顺序依赖性和奇点问题，使得物体旋转始终顺滑且符合直觉。
总结而言，四元数是OpenGL中处理复杂3D旋转不可或缺的数学工具。它通过抽象的复数形式，高效、稳定地解决了旋转的表示、组合和插值问题，是现代实时图形学领域的一项基础性技术。
参考来源
- OPENGL ES旋转问题的解决
- OpenGL 四元数旋转
- 四元数在opengl坐标转换的使用
- 旋转矩阵、欧拉角、四元数及四元数插值
- 三维旋转之四元数
- 四元数旋转实现：MATLAB中的向量旋转实践
