奇异矩阵及广义逆矩阵

再介绍奇异矩阵以及矩阵的广义逆之前，先复习几个概念，然后我一步步给出奇异矩阵和广义逆的介绍。

伴随矩阵

1. 余子式

定义1 余子式：在 n n 阶行列式 |A| | A | 中，划去元素 aij a i j 所在的行和列，剩下的元素（显然，剩下了 (n−1)2 ( n − 1 ) 2 个元素），按照原先的排列生成的新行列式的值称为元素 aij a i j 的余子式，记为 Mij M i j 。如下所示：

|A|=∣∣∣∣∣∣∣∣∣a11⋮ai1⋮an1………a1j⋮aij⋮anj………a1n⋮ain⋮ann∣∣∣∣∣∣∣∣∣,Mij=∣∣∣∣∣∣∣∣∣∣a11⋮ai−1 1ai−1 1⋮an1…………a1 j−1⋮ai−1 j−1ai+1 j−1⋮an j−1a1 j+1ai−1 j+1ai+1 j+1an j+1…⋮……⋮…a1nai−1 nai−1 nann∣∣∣∣∣∣∣∣∣∣(1) (1) | A | = | a 11 … a 1 j … a 1 n ⋮ ⋮ ⋮ a i 1 … a i j … a i n ⋮ ⋮ ⋮ a n 1 … a n j … a n n | , M i j = | a 11 … a 1   j − 1 a 1   j + 1 … a 1 n ⋮ ⋮ ⋮ a i − 1   1 … a i − 1   j − 1 a i − 1   j + 1 … a i − 1   n a i − 1   1 … a i + 1   j − 1 a i + 1   j + 1 … a i − 1   n ⋮ ⋮ ⋮ a n 1 … a n   j − 1 a n   j + 1 … a n n |

2. 代数余子式

定义2 代数余子式：由上面余子式的定义，可以得到 n n 阶行列式 |A| | A | 中元素 aij a i j 的代数余子式 Aij A i j 。其计算公式如下：

Aij=(−1)i+jMij(2) (2) A i j = ( − 1 ) i + j M i j

代数余子式是一个很有用的概念，因为用它可以直接计算出行列式的值。我们有下面这样一个定理存在：

定理1：行列式等于它某一行的元素分别与它们对应的代数余子式的乘积之和。可以用如下公式表达：

|A|=∑j=1naij⋅Aij,i∈{1,2,…,n}(3) (3) | A | = ∑ j = 1 n a i j ⋅ A i j , i ∈ { 1 , 2 , … , n }

定理1其实就是行列式计算方法的另一种表达。由这个定理，我们可以接着得到下面的定理2：

定理2：行列式中一行的元素与另一行相应元素的代数余子式的乘积之和为0。可以用如下公式表达：

ak1Ai1+⋯+aknAin=0(4) (4) a k 1 A i 1 + ⋯ + a k n A i n = 0

其中， aij=akj,j=1,2,…,n;i≠k a i j = a k j , j = 1 , 2 , … , n ; i ≠ k .这个定理的证明思路可以用行列式性质：“一行的数乘加到另一行，行列式的值不变”。具体的证明我这里略过了。

3. 伴随矩阵

定义3 伴随矩阵：由矩阵 A A 中所有元素的代数余子式按照 A A 中原有的排列的转置之后，构造而成的矩阵称为 A A 的伴随矩阵。也就是下面公式中的 A∗ A ∗ .

A∗=⎛⎝⎜⎜A11⋮A1n…⋱…An1⋮Ann⎞⎠⎟⎟(5) (5) A ∗ = ( A 11 … A n 1 ⋮ ⋱ ⋮ A 1 n … A n n )

注意，伴随矩阵中，元素的排列恰好原矩阵的转置序列是一致的。这样，根据上面说的定理1和定理2，我们可以得到这样的结论： AA∗=A∗A=dE A A ∗ = A ∗ A = d E .

解释一下，上面的公式中， d d 是 A A 的行列式的值， E E 是单位矩阵。

矩阵的逆

根据 AA∗=A∗A=dE A A ∗ = A ∗ A = d E 的性质，不难得到矩阵的逆的计算方式如下：

A−1=1dA∗(6) (6) A − 1 = 1 d A ∗

根据这个公式，也就可以直接得到矩阵可逆的充要条件，即下面的定理3：

定理3：矩阵可逆的充要条件是矩阵行列式的值不为0（ |A|≠0 | A | ≠ 0 ）

我们把这种行列式不为0的矩阵称为是“非退化”的，如果是方阵，那也称之为是“非奇异”的。与“非奇异”相对应的“奇异”矩阵则是数据挖掘，机器学习领域中非常重要的一个概念。下面将具体介绍。

奇异矩阵

定义3 奇异矩阵：奇异矩阵就是不满秩的方阵。

通过这个定义，就不难推出以下3点奇异矩阵的性质（记奇异矩阵为 A A ）：

1. 既然 A A 不满秩，那么 |A|=0 | A | = 0 ；
2. 根据上面说的定理3，如果 |A|=0 | A | = 0 ，那么 A A 一定是不可逆的。同理，非奇异矩阵则一定是可逆的；
3. 根据我之前介绍的线性方程组解的情况（线性方程组解的分析：唯一解，无穷多解以及无解）可知，齐次线性方程组 AX=0 A X = 0 有无穷多的解，而线性方程组 AX=b A X = b 有无穷多的解或者无解。同理对于非奇异矩阵 A′ A ′ ，齐次线性方程组 A′X=0 A ′ X = 0 有且只有唯一零解，线性方程组 A′X=b A ′ X = b 有唯一解。

广义逆矩阵

定义与概念

广义逆矩阵在数理统计、系统理论、优化计算和控制论等多领域中有重要应用。它的概念与上面的奇异矩阵密切相关，被如下定义：

定义4 广义逆矩阵：若 A A 是奇异阵或长方阵，根据上面说的奇异矩阵的性质3可知，线性方程组 Ax=b A x = b 有无穷多的解或者无解。若有解，则解为 x=Xb+(E−XA)y x = X b + ( E − X A ) y ，其中у是维数与A的列数相同的任意向量， X X 是满足 AXA=A A X A = A 的任意矩阵，那么我们称 X X 为 A A 的广义逆矩阵，用 Ag,A−,A1 A g , A − , A 1 等符号表示。当 A A 为非奇异矩阵时， A A 的广义逆矩阵就是 A A 的逆。

定义5：满足以下3个性质的广义逆矩阵 M M 称为 A A 的M-P广义逆矩阵矩阵，记为 A+ A + ：

• AMA = A
• MAM = M
• AM与MA均为对称矩阵

实际上，若 A A 是非奇异矩阵， A−1 A − 1 也满足上面这3个性质。也就是说，M-P逆就是通常逆矩阵的推广。

计算方法

广义逆矩阵的计算方法大致可分为三类：以满秩分解和奇异值分解为基础的直接法，迭代法和其他一些常用于低阶矩阵的非平凡方法。我只介绍基于满秩分解和奇异值分解的推导过程。

1. 满秩分解

对于 m×n m × n 的矩阵 A A ，假设其秩为 r r ，则 A A 有满秩分解 A=Fm×rGr×n A = F m × r G r × n 存在。那么 A A 的M-P逆可以通过满秩分解计算。下面我给出推导过程：

A+=(FG)+=G+F+=GT(GGT)−1⋅(FTF)−1FT=GT(FTAGT)−1FT(7) (7) A + = ( F G ) + = G + F + = G T ( G G T ) − 1 ⋅ ( F T F ) − 1 F T = G T ( F T A G T ) − 1 F T

通过上面的公式，可以算得关于 A A 的广义逆。

2. 奇异值分解

对于 m×n m × n 的矩阵 A A ，其奇异值分解如下：

A=UΣVT(8) (8) A = U Σ V T

其中 U U 是 m×m m × m 的正交阵， V V 是 n×n n × n 的正交阵， Σ Σ 是 m×n m × n 的对角阵。关于奇异值分解的详细介绍可以查阅我之前的博客（矩阵的分解：满秩分解和奇异值分解）。

根据奇异值分解，可以通过下面的计算得到广义逆：

A+=(UΣVT)+=(ΣVT)+⋅U−1=V(Σ)+UT(9) (9) A + = ( U Σ V T ) + = ( Σ V T ) + ⋅ U − 1 = V ( Σ ) + U T

其中， Σ Σ 因为是对角阵，所以广义逆就是他所有元素的倒数。