计算机图形学模型变换应用正确性及变换矩阵相关问题咨询

咱们一步步把这个问题掰明白，搞清楚变换矩阵顺序、逆矩阵的来龙去脉，还有旋转角度的计算方法～

你提到的逆矩阵和变换顺序混乱，本质是没分清这两种场景：

1. 物体变换（模型变换）：直接按操作顺序构建矩阵 #

如果是给物体做变换：先缩放S，再绕某个轴旋转R₁，接着绕另一个轴旋转R₂，最后平移T——这是把物体从局部空间转到世界空间的标准模型变换流程。

这里要注意图形学里的两种向量表示习惯：

• 列向量右乘（OpenGL、Unity等主流API常用）：变换矩阵的顺序是先做的变换放在乘法右侧，整体矩阵为：

  M = T × R₂ × R₁ × S
  

  原因很简单：列向量v' = M × v，计算时会先执行S×v（缩放），再R₁×(S×v)（旋转R₁），接着R₂×(R₁×S×v)（旋转R₂），最后T×(R₂×R₁×S×v)（平移），完全符合你想要的操作顺序。
• 行向量左乘：这种情况矩阵顺序反过来，整体矩阵是M = S × R₁ × R₂ × T，因为行向量的变换是v' = v × M，计算顺序是从左到右依次执行变换。

2. 为什么会用到逆矩阵？——坐标系反向变换的场景 #

你提到的逆矩阵，大概率是在做坐标系之间的反向转换，比如：

• 把世界空间的点转换到物体的局部空间（比如碰撞检测时要判断点是否在物体内部）
• 相机的观察变换：把世界空间转到相机视角空间（相当于把整个世界向相机移动的反方向变换）

矩阵乘法有个规则：(A×B×C)⁻¹ = C⁻¹ × B⁻¹ × A⁻¹，也就是逆矩阵的顺序是原矩阵顺序的反转。

比如刚才的物体变换矩阵是M_obj = T × R₂ × R₁ × S，那么反过来，把世界点转到物体局部空间的逆变换矩阵就是：

M_obj⁻¹ = S⁻¹ × R₁⁻¹ × R₂⁻¹ × T⁻¹

这里旋转矩阵的逆有个特殊点：旋转矩阵是正交矩阵，它的逆等于自身的转置（R⁻¹ = R^T），因为旋转矩阵的行/列都是互相垂直的单位向量，转置后就满足逆矩阵的条件，计算起来特别方便。

3. 你的两种矩阵选项：到底哪个对？ #

你纠结的M = T × R₁⁻¹ × R₂⁻¹ × S和M = T × R₂⁻¹ × R₁⁻¹ × S，都不符合标准的变换逻辑，大概率是把变换方向和矩阵乘法顺序搞混了：

• 如果是物体变换，根本不需要逆矩阵，顺序应该是T × R₂ × R₁ × S（列向量场景）
• 如果是逆变换（世界→局部），顺序应该是S⁻¹ × R₁⁻¹ × R₂⁻¹ × T⁻¹

有人告诉你第一种正确，可能是混淆了行/列向量的乘法方向，或者把坐标系变换的顺序搞反了。

4. 从旋转矩阵反推角度的具体公式 #

不同轴的旋转矩阵，提取角度的方式不同：

绕X轴旋转θ的矩阵 #

[1,  0,        0       ]
[0, cosθ,  -sinθ]
[0, sinθ,   cosθ]

• 取cosθ = R[1][1]，sinθ = R[2][1]
• 用θ = atan2(sinθ, cosθ)计算角度（范围是[-π, π]，自动处理象限问题）

绕Y轴旋转θ的矩阵 #

[cosθ,  0,  sinθ]
[0,     1,  0   ]
[-sinθ, 0, cosθ]

• 取cosθ = R[0][0]，sinθ = R[0][2]
• 同样用θ = atan2(sinθ, cosθ)计算

绕Z轴旋转θ的矩阵 #

[cosθ, -sinθ, 0]
[sinθ,  cosθ, 0]
[0,     0,    1]

• 取cosθ = R[0][0]，sinθ = R[1][0]
• θ = atan2(sinθ, cosθ)

任意轴的旋转矩阵（Rodrigues公式生成） #

可以通过矩阵的迹（对角线元素之和）计算：

• trace(R) = 1 + 2cosθ → cosθ = (trace(R) - 1)/2
• θ = arccos((trace(R)-1)/2)
• 旋转轴n的分量：n_x = (R[2][1]-R[1][2])/(2sinθ)，n_y = (R[0][2]-R[2][0])/(2sinθ)，n_z = (R[1][0]-R[0][1])/(2sinθ)
• 注意：当θ=0或π时，sinθ=0，需要特殊处理（比如θ=π时，旋转轴可以从矩阵的特征值为1的特征向量提取）

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