如何沿双轴旋转矩阵且完全避免第三轴产生旋转？

我有一个对象，需要沿两个轴进行旋转（为简化说明，暂称其为X轴和Y轴，但实际可为任意轴）。
示例代码如下：

Matrix aMat;
aMat.RotateAroundAxis(Vector(1,0,0),45.0f);
aMat.RotateAroundAxis(Vector(0,1,0),25.8f);

执行上述操作时，会在我所旋转的两个轴的叉积方向出现不必要的旋转（例如上述示例中，会产生少量Z轴（0,0,1）方向的旋转）。该旋转量虽小，但已足以被察觉，尤其在多次旋转后，误差会不断累积。
是否存在方法可以抑制这种情况？
（补充说明：若有更优方案，我的实际场景是一个带有配重的平台，希望平台能根据配重位置进行旋转，如同支点平衡一般。我当前的实现方式是根据配重相对平台中心的位置，累积绕平台向上方向的两个正交轴的旋转量。）

问题根源 #

这种额外旋转是矩阵旋转的顺序相关性导致的：每次绕局部轴旋转后，后续旋转的轴已经被之前的旋转改变了方向，多次累积后会引入欧拉角漂移（类似万向锁的副作用），也就是你看到的叉积方向的额外旋转。

可行解决方案 #

1. 用四元数替代矩阵累积旋转 #

四元数对旋转的组合更稳定，不会产生累积的轴漂移，是解决这类问题的首选方案。每次计算旋转增量时用四元数表示，再和当前旋转四元数相乘（注意乘法顺序：局部旋转的增量要后乘到当前四元数上），最后转成矩阵用于渲染。

示例伪代码：

Quaternion currentRot; // 初始化为单位四元数（无旋转）

// 计算绕局部X轴的旋转增量
Quaternion rotX = Quaternion::FromAxisAngle(Vector(1,0,0), 45.0f);
// 计算绕局部Y轴的旋转增量
Quaternion rotY = Quaternion::FromAxisAngle(Vector(0,1,0), 25.8f);

// 组合旋转：先应用X轴旋转，再应用Y轴旋转（局部空间）
currentRot = currentRot * rotX;
currentRot = currentRot * rotY;

// 将四元数转换为矩阵用于后续操作
Matrix aMat = currentRot.ToMatrix();

2. 直接计算目标姿态矩阵（推荐用于平衡平台场景） #

你的场景是平台根据配重实现支点平衡，本质是需要平台达到一个明确的目标姿态，而非逐步累积旋转。直接计算目标姿态的矩阵，完全避免增量旋转的累积误差。

具体步骤：

• 根据配重位置和重力方向，计算平台的目标"向上"向量（平衡时平台的姿态应该让配重产生的力矩为0，可通过力矩方向推导目标Up向量）
• 用目标Up向量构建一组正交的基向量（Right、Up、Forward），直接生成旋转矩阵

示例伪代码：

// 假设配重相对平台中心的位置为weightPos，重力方向为gravity=(0,-1,0)
// 计算目标Up向量（需根据你的平衡逻辑调整，这里仅为示例）
Vector targetUp = Vector::Cross(Vector::Cross(weightPos, gravity), weightPos);
targetUp.Normalize();

// 构建正交的Right向量，避免与Up共线的情况
Vector right = Vector::Cross(targetUp, Vector(0,0,1));
if (right.LengthSquared() < 0.0001f) {
    right = Vector::Cross(targetUp, Vector(1,0,0));
}
right.Normalize();

// 构建Forward向量
Vector forward = Vector::Cross(right, targetUp);
forward.Normalize();

// 用正交基向量构建旋转矩阵（列向量依次为Right、Up、Forward）
Matrix aMat;
aMat.SetColumn(0, right);
aMat.SetColumn(1, targetUp);
aMat.SetColumn(2, forward);

3. 矩阵正交化修正（应急方案） #

如果必须保留矩阵累积的方式，可以在每次旋转后对矩阵的轴向量进行正交化，消除因误差累积导致的轴非正交问题，缓解额外旋转。但这种方法只能治标，无法从根源解决问题。

示例伪代码：

// 每次旋转后执行正交化
Vector right = aMat.GetColumn(0);
Vector up = aMat.GetColumn(1);
Vector forward = aMat.GetColumn(2);

// 重新正交化轴向量
right.Normalize();
up = up - Vector::Dot(up, right) * right;
up.Normalize();
forward = Vector::Cross(right, up);
forward.Normalize();

// 更新矩阵
aMat.SetColumn(0, right);
aMat.SetColumn(1, up);
aMat.SetColumn(2, forward);

内容的提问来源于stack exchange，提问作者KiraHoneybee