绕Z轴复合变换矩阵如何构造？

Question

在三维图形变换中，如何构造一个绕Z轴的复合变换矩阵，使其同时包含旋转、缩放和偏移操作？常见问题在于变换顺序对最终结果的影响：若先旋转再平移，与先平移后旋转，物体位置明显不同。特别是在局部坐标系绕Z轴旋转时，平移分量需考虑旋转后的坐标方向。此外，缩放是否关于原点进行也会影响视觉效果。因此，如何正确组合旋转、缩放和平移矩阵，并按合理顺序左乘或右乘，成为构建绕Z轴复合变换矩阵的关键技术难点。初学者常因忽略变换顺序或矩阵乘法的非交换性而导致模型错位或变形。

Answer 1

三维图形变换中的绕Z轴复合变换矩阵构造技术详解

1. 基础概念：理解基本变换矩阵的数学表达

在三维图形学中，任何复杂的几何变换都可以分解为三种基本操作：平移（Translation）、旋转（Rotation）和缩放（Scaling）。这些操作均可用4×4齐次坐标变换矩阵表示。

• 平移矩阵：
  $$ T = \begin{bmatrix} 1 & 0 & 0 & t_x \\ 0 & 1 & 0 & t_y \\ 0 & 0 & 1 & t_z \\ 0 & 0 & 0 & 1 \end{bmatrix} $$
• 绕Z轴旋转矩阵（角度θ）：
  $$ R_z(\theta) = \begin{bmatrix} \cos\theta & -\sin\theta & 0 & 0 \\ \sin\theta & \cos\theta & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{bmatrix} $$
• 缩放矩阵：
  $$ S = \begin{bmatrix} s_x & 0 & 0 & 0 \\ 0 & s_y & 0 & 0 \\ 0 & 0 & s_z & 0 \\ 0 & 0 & 0 & 1 \end{bmatrix} $$

2. 变换顺序的重要性与非交换性分析

矩阵乘法不具备交换律，即 \( A \times B \neq B \times A \)，因此变换顺序直接影响最终结果。以下表格展示了不同顺序下对单位立方体顶点 (1,1,0) 的变换效果（设旋转45°，平移(2,3,0)，缩放(2,1,1)）：

变换顺序	数学表达	最终位置近似值
S → R → T	T × R × S	(3.41, 4.12, 0)
R → S → T	T × S × R	(4.12, 3.41, 0)
T → R → S	S × R × T	(4.83, 6.24, 0)
R → T → S	S × T × R	(6.24, 4.83, 0)
S → T → R	R × T × S	(−0.71, 5.66, 0)
T → S → R	R × S × T	(−3.54, 5.66, 0)

可见，仅改变顺序即可导致显著不同的空间位置，尤其在涉及局部坐标系旋转时更为敏感。

3. 构造绕Z轴复合变换的标准流程

为了实现“以物体自身中心为原点进行缩放、绕其局部Z轴旋转、再整体移动到目标位置”的需求，推荐使用如下标准顺序：

1. 将物体平移到世界原点（若缩放需关于质心）
2. 执行缩放操作
3. 绕Z轴旋转（局部Z轴此时与世界Z轴重合）
4. 平移回原始位置或指定目标位置

对应的复合变换矩阵为：

$$ M_{final} = T_{world} \times R_z(\theta) \times S \times T_{-center} $$

其中 \( T_{-center} \) 是将物体中心移至原点的逆平移，\( T_{world} \) 是最终定位平移。

4. 实际代码实现示例（OpenGL风格列主序）

#include <glm/glm.hpp>
#include <glm/gtc/matrix_transform.hpp>

glm::mat4 CreateCompositeTransform_Z(
    glm::vec3 scale,
    float rotationZ,
    glm::vec3 origin,     // 物体本地原点
    glm::vec3 worldPos    // 最终世界位置
) {
    glm::mat4 transform = glm::mat4(1.0f);
    
    // 步骤1：反向平移至原点
    transform = glm::translate(transform, -origin);
    
    // 步骤2：缩放
    transform = glm::scale(transform, scale);
    
    // 步骤3：绕Z轴旋转
    transform = glm::rotate(transform, rotationZ, glm::vec3(0,0,1));
    
    // 步骤4：平移到世界位置
    transform = glm::translate(transform, worldPos);
    
    return transform;
}

5. 局部坐标系下的变换陷阱与规避策略

当物体已存在初始旋转时，其局部Z轴可能不再与世界Z轴对齐。此时直接应用 \( R_z(\theta) \) 将导致绕世界Z轴而非局部轴旋转。解决方案包括：

• 使用四元数或欧拉角累积旋转状态，动态计算当前局部轴方向
• 采用层级变换结构（Scene Graph），父节点变换影响子节点基向量
• 在模型空间完成旋转后再叠加世界变换

graph TD
    A[原始顶点坐标] --> B[T⁻¹: 平移至原点]
    B --> C[S: 应用缩放]
    C --> D[Rz: 绕Z轴旋转]
    D --> E[T: 平移至世界位置]
    E --> F[最终屏幕投影]
    
    style A fill:#f9f,stroke:#333
    style F fill:#bbf,stroke:#333

6. 高级应用场景：动画系统中的复合变换优化

在骨骼动画或机器人运动学中，常需频繁重构复合变换。建议采用以下优化策略：

• 缓存中间变换矩阵，避免重复计算
• 使用SRT（Scale-Rotate-Translate）分解存储对象属性，便于插值
• 利用SIMD指令加速矩阵乘法（如Intel SSE或ARM NEON）
• 在GPU端通过Uniform Buffer Object批量上传变换矩阵

现代渲染引擎（如Unreal、Unity）内部均采用类似的SRT组合逻辑，并自动处理变换顺序一致性问题。