第八章 旋量系零空间构造理论
戴建生 · 高等教育出版社 2014
旋量系零空间构造:从几何到代数的求解法则
学习目标
- 理解旋量系零空间的数学表示(齐次线性方程组 $JB=0$)
- 掌握矩阵增广法(构造一维零空间)
- 掌握代数余子式法(定理 8.1)的原理和步骤
- 理解多维零空间的递归分块与逐级增广
- 理解与 Gauss-Seidel 消元法的比较
章节位置
第七章给出了旋量系关联关系定理,本章回到一个具体的计算问题:给定一组约束旋量(旋量系 $S$),满足这些约束的允许运动旋量(互易旋量系 $S_r$)如何构造?
这等价于求解齐次线性方程组 $JB = 0$。
8.1
旋量系零空间的数学表示
旋量系 $S = \{S_1, \cdots, S_n\}$ 与其互易旋量系 $S_r$ 的关系可以表示为齐次线性方程组:
$$J_S \cdot B = 0$$
其中 $J_S$ 是旋量矩阵($n \times 6$),$B$ 是待求的互易旋量系矩阵($6 \times (6-n)$)。
定义 8.1:$n \times k$ 矩阵 $J$ 的零空间是齐次线性方程组 $JB = 0$ 所有解向量的集合。
核心问题:给定约束旋量($n$ 个线性无关旋量),求满足约束的允许运动旋量——等价于求 $6 \times n$ 旋量矩阵的零空间。
前置知识回顾
- 线性代数基础:矩阵的秩、零空间、行空间、列空间。去哪里补:线性代数教材。
- 引理 8.1:行空间 $R(J^T)$ 与零空间 $N(J)$ 正交互补,$\dim R(J^T) + \dim N(J) = k$。
8.2
矩阵增广法(构造一维零空间)
当 $k - n = 1$(一维零空间)时,可以构造一个与矩阵 $J$ 其他行向量线性无关的增广向量 $v_a$:
$$v_a = ((-1)^{r+2}\det J_{c1}, (-1)^{r+3}\det J_{c2}, \cdots, (-1)^{r+k+1}\det J_{ck})^T$$
其中 $J_{cj}$ 是消去第 $j$ 列后的子矩阵的行列式(代数余子式)。
8.3
代数余子式法(定理 8.1)
定理 8.1:若齐次线性方程组含 $k$ 个未知数和 $n$ 个方程($n < k$),其系数矩阵为 $n \times k$,则 $k-n$ 维零空间的基向量可由代数余子式直接构造。
核心思想:零空间向量与矩阵的每个代数余子式直接相关,不需要逐行消元。
与 Gauss-Seidel 消元法相比:运算效率更高,误差阶数更低,利用了旋量系的几何结构(而非通用数值代数)。
8.5
多维零空间构造:递归分块与逐级增广
当 $k - n > 1$(多维零空间)时,对旋量矩阵进行分块处理:
- 将旋量矩阵分块为若干子矩阵
- 对每个子矩阵分别应用代数余子式法,得到低维零空间基
- 逐级增广,逐步构造高维零空间
每一步增广都保持零空间维数的正确性,最终得到 $k-n$ 个线性无关的零空间基向量。
复习速查
- 零空间问题:$J_S B = 0$,求互易旋量系
- 矩阵增广法:构造增广向量(式 8.4),利用代数余子式
- 代数余子式法:$k-n$ 维零空间由代数余子式直接构造
- 多维零空间:分块 + 逐级增广
参考来源
- 戴建生 (2014)《旋量代数与李群、李代数》第八章,P.202-P.232。教材:
media/screw-algebra/screw-algebra-lie-groups.pdf。