四元数向量旋转：微软实现与自定义实现结果不一致，求排查错误原因

首先咱们直接抓核心问题：你和微软实现的差异根源在于四元数旋转的乘法顺序和共轭的位置，这两种写法对应了完全相反的旋转方向——不是谁对谁错，而是要看你需要的是正向还是反向旋转。

先理清楚四元数旋转的标准数学定义 #

对于代表旋转的单位四元数 q，要旋转向量 v，正确的步骤是：

1. 把向量 v 转换成纯四元数 v_q = (v.x, v.y, v.z, 0)
2. 计算旋转后的四元数：q * v_q * q⁻¹，其中 q⁻¹ 是 q 的逆四元数。因为单位四元数的逆等于它的共轭（q⁻¹ = q_conj），所以公式也可以写成 q * v_q * q_conj。

你的实现与微软实现的核心差异 #

• 微软的优化代码：是标准公式 q * v_q * q_conj的展开计算，完全符合正向旋转的数学定义，这也是游戏引擎、图形库等行业通用的四元数旋转实现方式。
• 你的自定义实现：用的是 q_conj * v_q * q，等价于 q⁻¹ * v_q * q，本质上是反向旋转——相当于用原旋转四元数的逆来做正向旋转。

为什么你的测试场景能正常工作？ #

大概率是这两种情况之一：

1. 你在传入 rotation 参数时，实际传入的是原旋转的逆四元数，这样 q_conj * v_q * q 就刚好等价于标准的正向旋转；
2. 你的测试场景对旋转方向不敏感，或者本身需要的就是反向旋转（比如相机视角的旋转，有时会用反向旋转来匹配观察逻辑）。

简单测试验证差异 #

可以用一个明确的例子快速验证：

• 构造绕Y轴旋转90度的单位四元数：q = (0, sin(π/4), 0, cos(π/4))（即X=0，Y≈0.7071，Z=0，W≈0.7071）
• 待旋转向量 v = (1, 0, 0)
  • 用微软的实现：旋转后结果为 (0, 0, 1)（正向绕Y轴转90度，X轴正方向转到Z轴正方向）
  • 用你的实现：旋转后结果为 (0, 0, -1)（反向旋转，X轴正方向转到Z轴负方向）

总结 #

如果你的需求是标准正向旋转，微软的实现是符合数学定义和行业通用逻辑的；如果你的场景需要反向旋转，或者你传入的参数本身就是逆四元数，那你的实现是适配你的特定场景的。

内容的提问来源于stack exchange，提问作者Fries of Doom