在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空间的旋转。直接使用欧拉角计算会导致旋转轴混乱和物体“翻转”的失真现象。

解决方案是利用四元数增量式地更新旋转状态：

1. 根据当前帧与上一帧的触摸点位移，计算出一个屏幕空间的旋转向量。
2. 将该向量映射到一个基于当前摄像机朝向的3D旋转轴上。
3. 根据该轴和位移大小（作为角度因子）构造一个增量旋转四元数。
4. 将增量四元数与物体当前的总旋转四元数相乘（四元数乘法，表示旋转的组合），得到新的旋转状态。
5. 将最终的四元数转换为矩阵并应用。

这种方法避免了欧拉角的顺序依赖性和奇点问题，使得物体旋转始终顺滑且符合直觉。

总结而言，四元数是OpenGL中处理复杂3D旋转不可或缺的数学工具。它通过抽象的复数形式，高效、稳定地解决了旋转的表示、组合和插值问题，是现代实时图形学领域的一项基础性技术。

---

参考来源

• OPENGL ES旋转问题的解决
• OpenGL 四元数旋转
• 四元数在opengl坐标转换的使用
• 旋转矩阵、欧拉角、四元数及四元数插值
• 三维旋转之四元数
• 四元数旋转实现：MATLAB中的向量旋转实践