【3D数学基础：图形与游戏开发】笔记第9章矩阵的更多知识

参考资料&原文链接

参考书籍：【3D数学基础：图形与游戏开发】 ISBN7-302-10946XTP.7262

(美) etcher Dun著、(美) an Arberry 清华大学出版社

GAMES101-现代计算机图形学入门-闫令琪

猎豹网校：游戏开发之3D数学基础

矩阵的行列式

在任意方阵中都存在一个标量，称作该方阵的行列式。方阵M的行列式记作|M|或“detM”。非方阵矩阵的行列式是末定义的。

注意，在书写行列式时，两边用竖线将数字块围起来，省略方括号下面的示意图能帮助起记忆。将主对角线和反对角线上的元素各自相乘，然后用主对角线元素的积减去反对角线元素的积。

行列式的重要性质

矩阵积的行列式等于矩阵行列式的积：|AB|=|A||B|。
这可以扩展到多个矩阵的情况：|M₁M₂····M_n-1M_n| = |M₁||M₂|···|M_n-1||M_n|。
矩阵转置的行列式等于原矩阵的行列式：|M|^T = |M|。
如果矩阵的任意行或列全为零，那么它的行列式等于零。
交换矩阵的任意两行或两列，行列式变负。
任意行或列的非零积加到另一行或列上不会改变行列式的值。

行列式的几何意义

2D行列式的几何意义

矩阵的行列式有着非常有趣的几何解释。2D中，行列式等于以基向量为两边的平行四边形的有符号面积(如图9.1所示)。我们曾介绍过如何利用平行四边形形象化显示坐标空间的变换。

有符号面积是指如果平行四边形相对于原来的方位“翻转”那么面积变负。

3D行列式的几何意义

3D中，行列式等于以变换后的基向量为三边的平行六面体的有符号体积。3D中，如果变换使得平行六面体“由里向外”翻转，则行列式变负。

行列式和矩阵变换导致的尺寸改变相关。其中行列式的绝对值和面积(2D)、体积(3D)的改变相关。行列式的符号说明了变换矩阵是否包含镜象或投影。

矩阵的行列式还能对矩阵所代表的变换进行分类：

如果矩阵行列式为零，那么该矩阵包含投影。
如果阵行列式为负，那么该矩阵包含镜象。

代数余子式

假设矩阵M有r行，c列。记法M^|ij|表示从M中除去第i行和第j列后剩下的矩阵。显然，该矩阵有r-1行，c-1列。矩阵M^|ij|称作M的余子式。考虑3×3阶矩阵M：

余子式M^|12|是一个2×2阶矩阵，是从M中除去第1行和第2列的结果：

对方阵M，给定行、列元素的代数余子式等于相应余子式的有符号行列式，见公式9.3。

这个公式可能看不清楚，它是：C_ij = (-1)^i+j |M^|ij||。

代数余子式矩阵

如上，用记法C_ij表示M的第i行，第j列元素的代数余子式。注意余子式是一个矩阵，而代数余子式是一个标量。代数余子式计算式中的项(-1)^(i+j)有以棋盘形式使矩阵的代数余子式每隔一个为负的效果：

首先，从矩阵中任意选择一行或一列，对该行或列中的每个元素，都乘以对应的代数余子式。这些乘积的和就是矩阵的行列式。例如，任意选择一行，如行i行列式的计算过程如公式9.4所示。

矩阵的逆

注意，只有方阵才能够求逆。

方阵M的逆，记作M^-1，也是个矩阵，当M与M^-1乘时，结果是单位矩阵。用公式表示为公式9.6的形式。

并非所有矩阵都有逆。一个明显的例子是若矩阵的某一行或列上的元素都为零，用任何矩阵乘以该矩阵，结果都是一个零矩阵。如果一个矩阵有逆矩阵，那么称它为可逆的或非奇异的。如果一个矩阵没有逆矩阵，则称它为不可逆的或奇异矩阵。奇异矩阵的行列式为零，非奇矩阵的行列式不为零，所以检测行列式的值是判断矩阵是否可逆的有效方法。此外，对于任意可逆矩阵M，当且仅当v=0时，vM=0。

M的“标准伴随矩阵”记作“adj M”，定义为M的代数余子式矩阵的转置矩阵。下面是一个例子：

当然还有其他方法可以用来计算矩阵的逆，比如高斯消元法。很多线性代数书都断定该方法更适合在计算机上实现，因为它所使用的代数运算较少，这种说法其实是不正确的。对于大矩阵或某特殊矩阵来说，这也许是对的。然而，对于低阶矩阵，比如几何应用中常见的那些低阶矩阵，标准伴随矩阵可能更快些。因为可以为标准伴随矩阵提供无分支( branchless)实现，这种实现方法在当今的超标量体系结构和专用向量处理器上会更快一些。

矩阵的逆的重要性质

如果M是非奇异矩阵，则该矩阵的逆的逆等于原矩阵：(M^-1)^-1=M。
单位矩阵的逆是它本身：I^-1 = I。
矩阵转置的逆等于它的逆的转置：(M^T)^-1=(M^-1)^T。
矩阵乘积的逆等于矩阵的逆的相反顺序的乘积：(AB)^-1=B^-1A^-1。这可扩展到多个矩阵的情况：
(M₁M₂…M_n-1M_n)^-1=M_n^-1M_n-1^-1·······M₂^-1M₁^-1。

矩阵的逆的几何意义

矩阵的逆在几何上非常有用，因为它使得我们可以计算变换的“反向”或“相反”变换。——能“撤销”原变换的变换。所以，如果向量v用矩阵M来进行变换，接着用M的逆M进行变换，将会得到原向量，除了“投影”以外。这很容易通过代数方法验证：

(vM)M^-1=v(MM^-1) = vI=v

正交矩阵

运算法则

若方阵M是正交的，则当且仅当M与它转置M^T的乘积等于单位矩阵。检测矩阵的正交性：M正交 <=> MM^T = I。
回忆9.2节中的定义，矩阵乘以它的逆等于单位矩阵：MM^-1 = I。
所以，如果一个矩阵是正交的，那么它的转置等于它的逆。即：
M正交 <=> M^T=M^-1
这是一条非常有用的性质，因为在实际应用中经常需要计算矩阵的逆，而3D图形计算中正交矩阵出现得又是如此频繁。如果知道矩阵是正交的，就可以全避免计算逆矩阵了，这也将大大减少计算量。

正交矩阵的几何意义

换句话说，对于任意矩阵M，怎样检测它是否正交?当MM^T=I时M是正交的。它的确切含义如下：

现在做一些解释：

当且仅当一个向量是单位向量时，它与它自身的点积结果是1。因此，仅当r₁、r₂、r₃是单位向量时，第1、4、9式才能成立。
当且仅当两个向量互相垂直时，它们的点积为零。因此，仅当r₁、r₂、r₃互柑垂直时其他等式才成立。

所以一个矩阵是正交的，必须满足以下条件：

矩阵的每一行都是单位向量。
矩阵的所有行相垂直。
对矩阵的列也能得到类似的条件。这使得以下结论非常清楚：如果M是正交的，则M也是止交的。

计算逆矩阵时，仅在先知道矩阵是正交的情况下才能利用正交性的优点。如果预先不知道，那么检查正交性经常是浪费时间。即使在最好的情况下，先检查正交性以确定矩阵是否正交再进行转置，和开始就进行求逆运算也将耗费同样多的时间。而如果矩阵不是正交的，那么这种检查完全是浪费时间。

注意，有个术语上的差别可能会导致轻微的混淆。线性代数中，如果一组向量互相垂直，这组向量就被认为是正交基，它只要求所有向量互相垂直，并不要求所有向量都是单位向量。如它们都是单位向量，则称它们为标准正交基。本书所讲的正交矩阵的行或列向量都是指标准正交基向量所以，本书中由组正交基向量构造的矩阵并不一定是正交矩阵(除非基向量是标准正交的)。