ESC
输入关键词搜索文章
目录

第四章 位移算子与指数映射

戴建生 · 高等教育出版社 2014
Rodrigues公式与指数映射:从李代数到李群
学习目标
  • 理解坐标变换与位移算子的区别
  • 掌握 Rodrigues 公式的推导(从指数映射幂级数展开)
  • 掌握 Rodrigues 向量与 Euler-Rodrigues 参数
  • 理解指数映射 $so(3) \to SO(3)$$se(3) \to SE(3)$
  • 理解四元数与对偶四元数作为 SE(3) 表示的优势
章节位置
本节在课程中的位置

第三章给出了旋量的定义和互易积运算,本章将这些运算嵌入刚体运动学框架:如何用旋量轴参数化刚体的旋转和平移。

上一章留下了问题:旋量描述的是刚体运动的"静态"几何特征,但刚体从一个位姿到另一个位姿的变换如何用旋量表示?本章回答这个问题。

本章之后:第五章将有限位移旋量放到李群 SE(3) 的框架下,引入 $6 \times 6$ 伴随矩阵。

4.1
坐标变换与位移算子

坐标变换:将向量从某一坐标系变换至另一坐标系。刚体坐标变换包括旋转、平移和镜像。

关于 $z$ 轴旋转角度 $\alpha$ 的变换(式 4.1):

$$\mathbf{p} = R \mathbf{p}_l = \begin{pmatrix} \cos\alpha & -\sin\alpha & 0 \\ \sin\alpha & \cos\alpha & 0 \\ 0 & 0 & 1 \end{pmatrix} \mathbf{p}_l$$

$R$ 为正交矩阵,满足 $R^T R = I$(式 4.3)。

PDF图 4.1:纯旋转运动引起的坐标变换p.88

pdf/screw-algebra/screw-algebra-lie-groups.pdf · p.88

打开原文

4.1.1
反对称矩阵与 so(3)

定义 4.1:满足 $A = -A^T$ 的方阵为反对称矩阵(斜对称矩阵)。

向量叉积的矩阵形式:给定 $\mathbf{s} = (l, m, n)^T$,其反对称矩阵表示为:

$$[\hat{s}] = \begin{pmatrix} 0 & -n & m \\ n & 0 & -l \\ -m & l & 0 \end{pmatrix}$$

性质:$[\hat{s}] \mathbf{v} = \mathbf{s} \times \mathbf{v}$$[\hat{s}]$ 是李代数 $so(3)$ 的元素。

反对称矩阵的关键性质:

$$[\hat{s}]^3 = -[\hat{s}], \quad [\hat{s}]^2 = \mathbf{s}\mathbf{s}^T - \|\mathbf{s}\|^2 I$$

这是推导 Rodrigues 公式的核心。

4.1.2
齐次变换矩阵

旋转矩阵扩展为 $4 \times 4$ 齐次变换矩阵(式 4.7):

$$H = \begin{pmatrix} R & \mathbf{d} \\ \mathbf{0}^T & 1 \end{pmatrix} \in SE(3)$$

其中 $R \in SO(3)$$\mathbf{d} \in \mathbb{R}^3$$SE(3)$ 维数为 6。

PDF图 4.2:齐次变换矩阵p.89

pdf/screw-algebra/screw-algebra-lie-groups.pdf · p.89

打开原文

前置知识回顾

  • so(3)$3 \times 3$ 反对称矩阵构成的向量空间,李代数。去哪里补:线性代数教材或李群初步。
  • SE(3):特殊欧几里得群,刚体运动群。去哪里补:机器人学基础教材。
4.2
Rodrigues 公式

Rodrigues(1840)给出了绕单位轴 $\mathbf{u}$ 旋转角度 $\theta$ 的旋转矩阵。

指数映射推导:从 $\exp([\hat{u}]\theta)$ 出发,展开幂级数:

$$e^{[\hat{u}]\theta} = I + [\hat{u}]\sin\theta + [\hat{u}]^2(1-\cos\theta)$$
PDF图 4.3:Rodrigues 向量p.96

pdf/screw-algebra/screw-algebra-lie-groups.pdf · p.96

打开原文

Rodrigues 公式

$$R = I + [\hat{u}]\sin\theta + [\hat{u}]^2(1-\cos\theta)$$
推导关键:利用 $[\hat{u}]^3 = -[\hat{u}]$ 将无穷级数化简为闭式。
PDF图 4.4:Rodrigues 公式推导p.98

pdf/screw-algebra/screw-algebra-lie-groups.pdf · p.98

打开原文

4.3-4.4
四元数与对偶四元数

Rodrigues 参数等价于 Hamilton 四元数的 Rodrigues 参数化。单位四元数 $q = (q_0, \mathbf{q})$ 描述旋转。

对偶四元数 $\hat{q} = q_0 + \epsilon q$:实部 $q_0$ 表示旋转,虚部 $\epsilon q$ 表示平移。Clifford(1873)发明。

$$\hat{q} = q_0 + \epsilon \mathbf{q}, \quad \epsilon^2 = 0$$

对偶四元数的乘法规则:$(q_0 + \epsilon\mathbf{q})(r_0 + \epsilon\mathbf{r}) = q_0 r_0 + \epsilon(q_0\mathbf{r} + r_0\mathbf{q})$

对偶四元数完整描述 SE(3) 的元素,比 $4 \times 4$ 齐次矩阵更简洁,在轨迹规划中有广泛应用。

PDF图 4.5:对偶四元数p.104

pdf/screw-algebra/screw-algebra-lie-groups.pdf · p.104

打开原文

4.5-4.7
指数映射:se(3) → SE(3)

se(3) 中的旋量 $\hat{\mathcal{S}}$ 到 SE(3) 的指数映射:

$$\exp(\hat{\mathcal{S}} \theta) = \begin{pmatrix} \exp([\hat{s}]\theta) & (I - \exp([\hat{s}]\theta))(\mathbf{s} \times \mathbf{s}_0)/\|\mathbf{s}\|^2 + \mathbf{s}\mathbf{s}^T \mathbf{s}_0\theta/\|\mathbf{s}\|^2 \\ \mathbf{0}^T & 1 \end{pmatrix}$$

$\mathbf{s} \neq \mathbf{0}$ 时,旋量表示绕轴 $\mathbf{s}$ 旋转角度 $\theta$ 同时平移距离 $h\theta$$h$ 为旋距)。当 $\mathbf{s} = \mathbf{0}$ 时,退化为纯平移 $\exp([\hat{v}]) = I + [\hat{v}]$

物理意义:指数映射将瞬时旋量(se(3))映射为有限位移旋量(SE(3))。这是连接微分运动学与有限位移分析的核心工具。
复习速查
  • 齐次变换 $H = [R\; \mathbf{d};\; \mathbf{0}^T\; 1]$:刚体位姿的 $4 \times 4$ 矩阵表示
  • 反对称矩阵 $[\hat{u}]$:向量叉积的矩阵形式,$[\hat{u}]^3 = -[\hat{u}]$
  • Rodrigues 公式 $R = I + [\hat{u}]\sin\theta + [\hat{u}]^2(1-\cos\theta)$
  • Rodrigues 向量 $\tan(\theta/2)\mathbf{u}$:旋转的最小参数化
  • 对偶四元数 $q_0 + \epsilon \mathbf{q}$:统一表示旋转与平移
  • 指数映射 $\exp : se(3) \to SE(3)$

参考来源