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

第八章 旋量系零空间构造理论

戴建生 · 高等教育出版社 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$$
PDF图 8.1:旋量系零空间构造p.203

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

打开原文

其中 $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$$
PDF图 8.2:矩阵增广法p.205

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

打开原文

其中 $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$(多维零空间)时,对旋量矩阵进行分块处理:

  1. 将旋量矩阵分块为若干子矩阵
  2. 对每个子矩阵分别应用代数余子式法,得到低维零空间基
  3. 逐级增广,逐步构造高维零空间

每一步增广都保持零空间维数的正确性,最终得到 $k-n$ 个线性无关的零空间基向量。

复习速查
  • 零空间问题$J_S B = 0$,求互易旋量系
  • 矩阵增广法:构造增广向量(式 8.4),利用代数余子式
  • 代数余子式法$k-n$ 维零空间由代数余子式直接构造
  • 多维零空间:分块 + 逐级增广

参考来源

  • 戴建生 (2014)《旋量代数与李群、李代数》第八章,P.202-P.232。教材:media/screw-algebra/screw-algebra-lie-groups.pdf