跳到主要内容

行列式

1. 线性方程组的通解​

对于二元一次方程组:

{a11x+a12y=b1a21x+a22y=b2\begin{cases} a_{11}x + a_{12}y = b_1 \\ a_{21}x + a_{22}y = b_2 \end{cases}

如果它有唯一解,那么利用高斯消元法,可以将它的解为:

{x=b1a22−a12b2a11a22−a12a21y=b2a11−a21b1a11a22−a12a21\begin{cases} x = \frac{b_1 a_{22} - a_{12}b_2}{a_{11}a_{22} - a_{12}a_{21}} \\ y = \frac{b_2 a_{11} - a_{21}b_1}{a_{11}a_{22} - a_{12}a_{21}} \end{cases}

这个解的分子和分母的形式非常相似,写成方阵形式时都是交叉相乘再相减,以分母为例:

a11a22−a12a21=∣a11a12a21a22∣a_{11}a_{22} - a_{12}a_{21} =\begin{vmatrix} a_{11} & a_{12} \\ a_{21} & a_{22} \end{vmatrix}

这样,上面方程组的解就可以表示为:

{x=∣b1a12b2a22∣∣a11a12a21a22∣y=∣a11b1a21b2∣∣a11a12a21a22∣\begin{cases} x = \frac{\begin{vmatrix} b_1 & a_{12} \\ b_2 & a_{22} \end{vmatrix}}{\begin{vmatrix} a_{11} & a_{12} \\ a_{21} & a_{22} \end{vmatrix}} \\ y = \frac{\begin{vmatrix} a_{11} & b_1 \\ a_{21} & b_2 \end{vmatrix}}{\begin{vmatrix} a_{11} & a_{12} \\ a_{21} & a_{22} \end{vmatrix}} \end{cases}

同样的,用高斯消元法可以求解三元一次方程组的解:

{x=b1a22a33−b1a23a32−a12b2a33+a12a23b3+a13b2a32−a13a22b3a11a22a33−a11a23a32−a12a21a33+a12a23a31+a13a21a32−a13a22a31y=a11b2a33−a11a23b3−b1a21a33+b1a23a31+a13a21b3−a13b2a31a11a22a33−a11a23a32−a12a21a33+a12a23a31+a13a21a32−a13a22a31z=a11a22b3−a11b2a32−a12a21b3+a12b2a31+b1a21a32−b1a22a31a11a22a33−a11a23a32−a12a21a33+a12a23a31+a13a21a32−a13a22a31\begin{cases} x = \frac{b_1 a_{22} a_{33} - b_1 a_{23} a_{32} - a_{12} b_2 a_{33} + a_{12} a_{23} b_3 + a_{13} b_2 a_{32} - a_{13} a_{22} b_3}{a_{11} a_{22} a_{33} - a_{11} a_{23} a_{32} - a_{12} a_{21} a_{33} + a_{12} a_{23} a_{31} + a_{13} a_{21} a_{32} - a_{13} a_{22} a_{31}} \\ y = \frac{a_{11} b_2 a_{33} - a_{11} a_{23} b_3 - b_1 a_{21} a_{33} + b_1 a_{23} a_{31} + a_{13} a_{21} b_3 - a_{13} b_2 a_{31}}{a_{11} a_{22} a_{33} - a_{11} a_{23} a_{32} - a_{12} a_{21} a_{33} + a_{12} a_{23} a_{31} + a_{13} a_{21} a_{32} - a_{13} a_{22} a_{31}} \\ z = \frac{a_{11} a_{22} b_3 - a_{11} b_2 a_{32} - a_{12} a_{21} b_3 + a_{12} b_2 a_{31} + b_1 a_{21} a_{32} - b_1 a_{22} a_{31}}{a_{11} a_{22} a_{33} - a_{11} a_{23} a_{32} - a_{12} a_{21} a_{33} + a_{12} a_{23} a_{31} + a_{13} a_{21} a_{32} - a_{13} a_{22} a_{31}} \end{cases}

显然解的各项的分子与分母的形式都相似,如果将这种形式的参数用方阵表示,就是:

{x=∣b1a12a13b2a22a23b3a32a33∣∣a11a12a13a21a22a23a31a32a33∣y=∣a11b1a13a21b2a23a31b3a33∣∣a11a12a13a21a22a23a31a32a33∣z=∣a11a12b1a21a22b2a31a32b3∣∣a11a12a13a21a22a23a31a32a33∣\begin{cases} x = \frac{\begin{vmatrix} b_1 & a_{12} & a_{13} \\ b_2 & a_{22} & a_{23} \\ b_3 & a_{32} & a_{33} \end{vmatrix}}{\begin{vmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{vmatrix}} \\ y = \frac{\begin{vmatrix} a_{11} & b_1 & a_{13} \\ a_{21} & b_2 & a_{23} \\ a_{31} & b_3 & a_{33} \end{vmatrix}}{\begin{vmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{vmatrix}} \\ z = \frac{\begin{vmatrix} a_{11} & a_{12} & b_1 \\ a_{21} & a_{22} & b_2 \\ a_{31} & a_{32} & b_3 \end{vmatrix}}{\begin{vmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{vmatrix}} \end{cases}

同样的,扩展到nn元时,解可以表示为:

xi=∣a11⋯b1⋯a1na21⋯b2⋯a2n⋮⋱⋮an1⋯bn⋯ann∣ith∣a11a12⋯a1na21a22⋯a2n⋮⋮⋱⋮an1an2⋯ann∣x_{i} = \frac{\stackrel{ith}{\begin{vmatrix} a_{11} & \cdots &b_1 & \cdots & a_{1n} \\ a_{21} & \cdots &b_2 & \cdots & a_{2n} \\ \vdots & & \ddots & & \vdots \\ a_{n1} & \cdots & b_{n} & \cdots & a_{nn} \end{vmatrix}}}{\begin{vmatrix} a_{11} & a_{12} & \cdots & a_{1n} \\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{vmatrix}}

这种方阵表示为解析式时,有n!n!项,每一项由前面的符号和nn项系数的乘积组成,即:

(−1)ta1p1a2p2⋯anpn(-1)^{t} a_{1p_1} a_{2p_2} \cdots a_{np_n}

其中p1p2⋯pnp_1 p_2 \cdots p_n为{ 1,2,⋯ ,n }\Set{1,2,\cdots,n}的全排列,aija_{ij}为参数方阵的元素,t由排列的逆序数来确定。

逆序数

在一个排列中,如果一对数的前后位置与大小顺序相反,即前面的数大于后面的数,则称为一个逆序对(英语:reverse pair)。一个排列中逆序对的个数称为该排列的逆序数(英语:reverse count)。

例如,对于数列A={ 3,2,5,4,1 }A = \Set{3,2,5,4,1},用ti{t_i}分别表示数列里第ii个元素左侧比它大的元素的个数,那么:

{t1=0t2=1t3=0t4=1t5=4\begin{cases} t_1 &= 0 \\ t_2 &= 1 \\ t_3 &= 0 \\ t_4 &= 1 \\ t_5 &= 4 \end{cases}

该数列的逆序数为tA=∑i=15ti=6t_A=\sum_{i=1}^5 t_i = 6。

上面的数列,交换一下5和4的顺序,即B={ 3,2,4,5,1 }B = \Set{3,2,4,5,1},用ti{t_i}分别表示数列里第ii个元素左侧比它大的元素的个数,那么:

{t1=0t2=1t3=0t4=0t5=4\begin{cases} t_1 &= 0 \\ t_2 &= 1 \\ t_3 &= 0 \\ t_4 &= 0 \\ t_5 &= 4 \end{cases}

此时,该数列的逆序数为tB=∑i=15ti=5t_B=\sum_{i=1}^5 t_i = 5。

提示

交换排列的任意两个元素,排列的逆序数会改变奇偶性。

求解线性方程组时,总结出来的方阵即为行列式(Determinant)。

2. 行列式的定义​

n阶行列式

方阵AA的行列式(英语:determinant),是一种数学运算,结果是一个数,记作∣A∣|A|或det⁡(A)\det(A)。 n×nn \times n的方阵AA,其行列式可以表示为:

det⁡(A)=∣A∣=∣a11a12⋯a1na21a22⋯a2n⋮⋮⋱⋮an1an2⋯ann∣\det(A) = |A| = \begin{vmatrix} a_{11} & a_{12} & \cdots & a_{1n} \\ a_{21} & a_{22} & \cdots & a_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{vmatrix}

且

det⁡(A)=∑(−1)ta1p1a2p2⋯anpn\det(A) = \sum(-1)^{t} a_{1p_1} a_{2p_2} \cdots a_{np_n}

或

det⁡(A)=∑(−1)tap11ap22⋯apnn\det(A) = \sum(-1)^{t} a_{p_1 1} a_{p_2 2} \cdots a_{p_n n}

其中tt为排列p1p2⋯pnp_1 p_2 \cdots p_n的逆序数,aija_{ij}为AA的元素。

2.1. 二阶行列式​

det⁡[abcd]=∣abcd∣=ad−bc\det\begin{bmatrix} a & b \\ c & d \end{bmatrix} = \begin{vmatrix} a & b \\ c & d \end{vmatrix} = ad - bc

等于列向量 [ac]\begin{bmatrix} a \\ c \end{bmatrix} 和 [bd]\begin{bmatrix} b \\ d \end{bmatrix} 张成的平行四边形面积。

2.2. 三阶行列式​

∣a11a12a13a21a22a23a31a32a33∣=a11∣a22a23a32a33∣−a12∣a21a23a31a33∣+a13∣a21a22a31a32∣\begin{aligned} \begin{vmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{vmatrix} &= a_{11} \begin{vmatrix} a_{22} & a_{23} \\ a_{32} & a_{33} \end{vmatrix} - a_{12} \begin{vmatrix} a_{21} & a_{23} \\ a_{31} & a_{33} \end{vmatrix} + a_{13} \begin{vmatrix} a_{21} & a_{22} \\ a_{31} & a_{32} \end{vmatrix} \end{aligned}

等于三个列向量张成的平行六面体体积,这与向量的混合积一致: det⁡[abc]=a⋅(b×c)\det\begin{bmatrix} \bm{a} & \bm{b} & \bm{c} \end{bmatrix} = \bm{a} \cdot (\bm{b} \times \bm{c}) 。

2.3. 行列式的计算复杂度​

直接根据其定义进行行列式计算,以3阶行列式为例,它共有3!3!种排列方式,每种排列需要22次乘法,所以总共需要2×3!=122\times 3! = 12次乘法。同时,它还需要3!−1=53!-1 = 5次加/减法,所以总的计算量为12+5=1712+5 = 17次。

对于nn阶,它共有n!n!种排列方式,每种排列需要n−1n-1次乘法,所以总共需要(n−1)n!(n-1)n!次乘法。 同时,它还需要n!−1n!-1次加/减法,所以总的计算量为(n−1)n!+n!−1=n×n!−1(n-1)n! + n! - 1 = n\times n! - 1次。

注意

直接计算nn阶行列式的计算复杂度为O(n!)O(n!),对于复杂的方阵,计算行列式的代价是很大的。

3. 拉普拉斯展开​

对于一个3×33 \times 3的方阵A=[a11a12a13a21a22a23a31a32a33]A = \begin{bmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{bmatrix},将其行列式按第一行的项合并,有:

det⁡(A)=a11a22a33−a11a23a32−a12a21a33+a12a23a31+a13a21a32−a13a22a31=a11(a22a33−a23a32)−a12(a21a33−a23a31)+a13(a21a32−a22a31)=a11×∣a22a23a32a33∣−a12×∣a21a23a31a33∣+a13×∣a21a22a31a32∣=a11∣A11∣−a12∣A12∣+a13∣A13∣=∑j=13(−1)1+ja1jdet⁡(A1j)\begin{aligned} \det(A) &= a_{11} a_{22} a_{33} - a_{11} a_{23} a_{32} - a_{12} a_{21} a_{33} + a_{12} a_{23} a_{31} + a_{13} a_{21} a_{32} - a_{13} a_{22} a_{31} \\ &= a_{11} (a_{22} a_{33} - a_{23} a_{32}) - a_{12}( a_{21} a_{33} - a_{23} a_{31}) + a_{13}( a_{21} a_{32} - a_{22} a_{31}) \\ & = a_{11} \times \begin{vmatrix} a_{22} & a_{23} \\ a_{32} & a_{33} \end{vmatrix} - a_{12} \times \begin{vmatrix} a_{21} & a_{23} \\ a_{31} & a_{33} \end{vmatrix} + a_{13} \times \begin{vmatrix} a_{21} & a_{22} \\ a_{31} & a_{32} \end{vmatrix}\\ &= a_{11}|A_{11}| - a_{12}|A_{12}| + a_{13}|A_{13}|\\ &=\sum_{j=1}^3 (-1)^{1+j} a_{1j} \det(A_{1j}) \end{aligned}

其中AijA_{ij}表示AA的去掉第ii行第jj列后的子阵(英语:sub matrix),其对应的行列式det⁡(Aij)\det(A_{ij}),称作余子式(英语:minor)。

将上述形式推广到n阶方阵,有:

det⁡(A)=∣A∣=∑j=1n(−1)i+jaijdet⁡(Aij)=∑j=1naijCij\begin{aligned} \det(A) &= |A| \\ &= \sum_{j=1}^n (-1)^{i+j} a_{ij}\det(A_{ij}) \\ &= \sum_{j=1}^n a_{ij} C_{ij} \end{aligned}

其中,CijC_{ij}为代数余子式(英语:cofactor),表示的是余子式和它前面的符号结合起来,即:

Cij=(−1)i+jdet⁡(Aij)C_{ij} = (-1)^{i+j} \det(A_{ij})

除了可以按第一行展开,还可以按第一列展开,甚至可以按任意一行或列展开,这就定义了一种递归的计算行列式的方法,即拉普拉斯展开(英语:Laplace expansion)。

拉普拉斯展开

对于任意n×nn \times n的方阵AA,n≥2n \geq 2,都有:

  1. 以第ii行展开:
det⁡(A)=∣A∣=∑j=1n(−1)i+jaijdet⁡(Aij)=∑j=1naijCij\det(A) = |A| = \sum_{j=1}^n (-1)^{i+j} a_{ij}\det(A_{ij}) = \sum_{j=1}^n a_{ij} C_{ij}
  1. 以第jj列展开:
det⁡(A)=∣A∣=∑i=1n(−1)i+jaijdet⁡(Aij)=∑i=1naijCij\det(A) = |A| = \sum_{i=1}^n (-1)^{i+j} a_{ij}\det(A_{ij}) = \sum_{i=1}^n a_{ij} C_{ij}

其中,Cij=(−1)i+jdet⁡(Aij)C_{ij} = (-1)^{i+j} \det(A_{ij}),为代数余子式。

例如,∣A∣=∣123456789∣ |A| = \begin{vmatrix} 1 & 2 & 3 \\ 4 & 5 & 6 \\ 7 & 8 & 9 \end{vmatrix}的计算过程如下:

∣A∣=1×(−1)1+1×∣5689∣+2×(−1)1+2×∣4679∣+3×(−1)1+3×∣4578∣=1×∣5689∣−2×∣4679∣+3×∣4578∣=1×(5×9−6×8)−2×(4×9−6×7)+3×(4×8−5×7)=1×(−3)−2×(−6)+3×(−3)=0\begin{aligned} |A| &= 1 \times (-1)^{1+1} \times \begin{vmatrix} 5 & 6 \\ 8 & 9 \end{vmatrix} + 2 \times (-1)^{1+2} \times \begin{vmatrix} 4 & 6 \\ 7 & 9 \end{vmatrix} + 3 \times (-1)^{1+3} \times \begin{vmatrix} 4 & 5 \\ 7 & 8 \end{vmatrix}\\ &= 1 \times \begin{vmatrix} 5 & 6 \\ 8 & 9 \end{vmatrix} - 2 \times \begin{vmatrix} 4 & 6 \\ 7 & 9 \end{vmatrix} + 3 \times \begin{vmatrix} 4 & 5 \\ 7 & 8 \end{vmatrix}\\ &= 1 \times (5 \times 9 - 6 \times 8) - 2 \times (4 \times 9 - 6 \times 7) + 3 \times (4 \times 8 - 5 \times 7)\\ &= 1 \times (-3) - 2 \times (-6) + 3 \times (-3)\\ &= 0 \end{aligned}

3.1 计算复杂度​

令M(n)M(n)为nn阶行列式需要的乘法次数,A(n)A(n)为nn阶行列式需要的加法/减法次数。显然,M(2)=2M(2)=2,A(2)=1A(2)=1。

对于n≥3n \geq 3,有:

A(n)=n−1+nA(n),M(n)=n+nM(n−1) A(n) = n - 1 + nA(n), \quad M(n) = n + nM(n-1)

先化简A(n)A(n):

A(n)+1=n(A(n−1)+1)=n(n−1)(A(n−2)+1)=n(n−1)⋯3(A(2)+1)=n⋅(n−1)⋯3⋅2⋅1=n!\begin{aligned} A(n) + 1 &= n(A(n-1) + 1) \\ & = n(n-1)(A(n-2) + 1) \\ & = n(n-1)\cdots 3 (A(2) + 1) \\ & = n\cdot(n-1) \cdots 3 \cdot 2 \cdot 1 \\ & = n! \end{aligned}

那么,A(n)=n!−1A(n) = n! - 1。

再化简M(n)M(n):

M(n)=n+nM(n−1)=n+n(⋯(5+5⋅(4+4⋅(3+3⋅M(2)))))=n+n(⋯(5+5⋅(4+4⋅(3+3⋅2))))=n+n(⋯(5+5⋅4+5⋅4⋅3+5⋅4⋅3⋅2))=n+n⋅(n−1)+n⋅(n−1)⋅(n−2)+⋯+n!=n!(1+12!+⋯+1(n−1)!)\begin{aligned} M(n) &= n + nM(n-1) \\ & = n + n(\cdots (5 + 5\cdot(4 + 4\cdot(3 + 3\cdot M(2))))) \\ & = n + n(\cdots (5 + 5\cdot(4 + 4\cdot(3 + 3\cdot 2)))) \\ & = n + n(\cdots (5 + 5\cdot4 + 5\cdot4\cdot3 + 5\cdot4\cdot3\cdot2)) \\ & = n + n \cdot (n-1) + n \cdot (n-1) \cdot (n-2) + \cdots + n! \\ & = n!(1 + \frac{1}{2!} + \cdots + \frac{1}{(n-1)!}) \end{aligned}

那么,有:n!<M(n)<(1−e)n!n! < M(n) < (1-e)n!

因此,通过拉普拉斯定理展开行列式计算复杂度为A(n)+M(n)A(n)+M(n),即O(n!)O(n!)。

注意

使用拉普拉斯展开计算行列式的计算复杂度为O(n!)O(n!),对于复杂的矩阵,计算行列式的代价是很大的。

4. 行列式的性质​

4.1 基本性质​

定理

对于一个n×nn \times n的方阵AA,其行列式有如下性质:

  1. 如果A某行或列的元素全为0,则行列式为0。
  2. 如果A的某两行或两列相等,则行列式为0。
  3. 交换AA的两行或两列得到BB,则∣B∣=−∣A∣|B| = -|A|,即交换两行或两列会改变行列式的符号。
  4. AA的某行或某列乘以kk得到BB,则∣B∣=k∣A∣|B| = k|A|。
  5. ∣kA∣=kn∣A∣|kA| = k^n|A|
  6. 将某一行/列的k倍,加减到另一行/列,得到的行列式不变。

4.2 三角阵的行列式​

三角阵的行列式

对于一个n×nn \times n的方阵AA,如果它是上三角阵或下三角阵,则其行列式可以用对角线元素的乘积得到,即:

∣A∣=∏i=1naii|A| = \prod_{i=1}^n a_{ii}

通过不断进行行/列变换得到三角阵的复杂度为O(n3)O(n^3),三角阵的行列式计算复杂度是O(n)O(n)。先化为三角阵再计算行列式,可大幅降低行列式的计算复杂度。

例如计算A=∣3−169−52047∣A = \begin{vmatrix} 3 & -1 & 6 \\ 9 & -5 & 2 \\ 0 & 4 & 7 \end{vmatrix}的行列式:

∣A∣=∣3−169−52047∣(第二行减去第一行的3倍)=∣3−160−2−16047∣(第三行加第二行的2倍)=∣3−160−2−1600−25∣=3×(−2)×(−25)=150\begin{aligned} |A| &= \begin{vmatrix} 3 & -1 & 6 \\ 9 & -5 & 2 \\ 0 & 4 & 7 \end{vmatrix}\text{(第二行减去第一行的3倍)}\\ &=\begin{vmatrix} 3 & -1 & 6 \\ 0 & -2 & -16 \\ 0 & 4 & 7 \end{vmatrix} \text{(第三行加第二行的2倍)}\\ &=\begin{vmatrix} 3 & -1 & 6 \\ 0 & -2 & -16 \\ 0 & 0 & -25 \end{vmatrix}\\ &= 3 \times (-2) \times (-25) \\ &= 150 \end{aligned}

4.3 矩阵与行列式​

定理
  1. 对于方阵AA,当且仅当det⁡(A)≠0\det(A) \neq 0时,AA可逆,且有det⁡(A−1)=1det⁡(A)\det(A^{-1}) = \frac{1}{\det(A)}
  2. det⁡(AB)=det⁡(A)det(B)\det(AB) = \det(A)det(B)
  3. det⁡(A+B)≠det⁡(A)+det⁡(B)\det(A+B) \neq \det(A) + \det(B)
  4. det⁡(AT)=det⁡(A)\det(A^T) = \det(A)

5. 行列式的几何意义​

将方阵 AA 视为一个线性变换,det⁡(A)\det(A) 的绝对值表示变换后单位体积被缩放了多少倍,其符号表示变换是否改变了空间的定向(英语:orientation):

  • det⁡(A)>0\det(A) > 0:变换保持空间定向不变(右手系仍为右手系)。
  • det⁡(A)<0\det(A) < 0:变换反转了空间定向(右手系变为左手系)。
  • det⁡(A)=0\det(A) = 0:变换将空间压缩到了更低维度(降维),即 AA 为奇异矩阵,不可逆。

6. 克拉默法则​

以三元一次方程组的解为例:

x=∣b1a12a13b2a22a23b3a32a33∣∣a11a12a13a21a22a23a31a32a33∣,y=∣a11b1a13a21b2a23a31b3a33∣∣a11a12a13a21a22a23a31a32a33∣,z=∣a11a12b1a21a22b2a31a32b3∣∣a11a12a13a21a22a23a31a32a33∣ x = \frac{\begin{vmatrix} b_1 & a_{12} & a_{13} \\ b_2 & a_{22} & a_{23} \\ b_3 & a_{32} & a_{33} \end{vmatrix}}{\begin{vmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{vmatrix}} ,\quad y = \frac{\begin{vmatrix} a_{11} & b_1 & a_{13} \\ a_{21} & b_2 & a_{23} \\ a_{31} & b_3 & a_{33} \end{vmatrix}}{\begin{vmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{vmatrix}} ,\quad z = \frac{\begin{vmatrix} a_{11} & a_{12} & b_1 \\ a_{21} & a_{22} & b_2 \\ a_{31} & a_{32} & b_3 \end{vmatrix}}{\begin{vmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{vmatrix}}

显然,分母是参数方阵的行列式,分子是b\mathbf{b}替换参数方阵的第ii列得到的矩阵的行列式。

推广到n阶时,AA是一个n×nn\times n矩阵,b\mathbf{b}是一个n×1n\times 1列向量。Ai(b)A_i(\mathbf{b})表示用b\mathbf{b}替代AA在第ii列得到的新矩阵,即:

Ai(b)=[a1⋯b⋯an]ithA_i(\mathbf{b}) = \stackrel{ith}{\begin{bmatrix} \mathbf{a_1} \cdots \mathbf{b} \cdots \mathbf{a_n} \end{bmatrix}}

单位矩阵InI_n的每一列都是一个单位向量,即In=[e1e2⋯en]I_n = \begin{bmatrix} e_1 & e_2 & \cdots & e_n \end{bmatrix},其中eie_i是ii号单位向量。

如果Ax=bA\mathbf{x}=\mathbf{b},那么有:

AIi(x)=A[e1⋯x⋯en]=[Ae1⋯Ax⋯Aen]=[a1⋯b⋯an]=Ai(b)\begin{aligned} AI_i(\mathbf{x}) &= A\begin{bmatrix} \mathbf{e_1} & \cdots & \mathbf{x} & \cdots & \mathbf{e_n} \end{bmatrix} \\ &= \begin{bmatrix} A\mathbf{e_1} & \cdots & A\mathbf{x} & \cdots & A\mathbf{e_n} \end{bmatrix} \\ &= \begin{bmatrix} \mathbf{a_1} & \cdots & \mathbf{b} & \cdots & \mathbf{a_n} \end{bmatrix} \\ &= A_i(\mathbf{b}) \end{aligned}

因此,

det⁡(AIi(x))=det⁡(A)det⁡(Ii(x))=det⁡(Ai(b))\det(AI_i(\mathbf{x})) = \det(A)\det(I_i(\mathbf{x})) = \det(A_i(\mathbf{b}))

而

det⁡(Ii(x))=[10⋯x1⋯0001⋯x2⋯00⋮⋮⋱⋮⋮⋮00⋯xi⋯00⋮⋮⋮⋱⋮⋮00⋯xn−1⋯0000⋯xn⋯01]=xi\det(I_i(\mathbf{x})) = \begin{bmatrix} 1 & 0 & \cdots & x_1 & \cdots & 0 & 0 \\ 0 & 1 & \cdots & x_2 & \cdots & 0 & 0 \\ \vdots & \vdots & \ddots & \vdots & & \vdots & \vdots \\ 0 & 0 & \cdots & x_i & \cdots & 0 & 0 \\ \vdots & \vdots & & \vdots & \ddots & \vdots & \vdots \\ 0 & 0 & \cdots & x_{n-1} & \cdots & 0 & 0 \\ 0 & 0 & \cdots & x_n & \cdots & 0 & 1 \end{bmatrix} = x_i

于是,

xi=det⁡(Ai(b))det⁡(A)x_i = \frac{\det(A_i(\mathbf{b}))}{\det(A)}

这就是克拉默法则(英语:Cramer's rule)。

克拉默法则

设AA是一个n×nn\times n矩阵,b\mathbf{b}是一个n×1n\times 1列向量,那么Ax=bA\mathbf{x}=\mathbf{b}的解x\mathbf{x}满足:

x=1det⁡(A)[det⁡(A1(b))⋮det⁡(An(b))]\mathbf{x} = \frac{1}{\det(A)}\begin{bmatrix} \det(A_1(\mathbf{b})) \\ \vdots \\ \det(A_n(\mathbf{b})) \end{bmatrix}

即:

xi=det⁡(Ai(b))det⁡(A)x_i = \frac{\det(A_i(\mathbf{b}))}{\det(A)}

其中i∈{ 1,2,⋯ ,n }i \in \Set{1, 2, \cdots, n}。

6.1 应用​

设AA是一个n×nn\times n矩阵,如果AA可逆,那么存在矩阵XX,使得AX=InAX=I_n,令xj\mathbf{x_j}为XX的第jj列。每次求解一列,即Axj=ejA\mathbf{x_j}=\mathbf{e_j},应用克拉默法则,有:

xij=det⁡(Ai(ej))det⁡(A)x_{ij} = \frac{\det(A_i(\mathbf{e_j}))}{\det(A)}

而,

det⁡(Ai(ej))=[a11⋯0⋯a1n⋮⋱⋮⋮aj1⋯1⋯ajn⋮⋮⋱⋮an1⋯0⋯ann]=(−1)j+idet⁡(Aji)=Cji\det(A_i(\mathbf{e_j})) = \begin{bmatrix} a_{11} & \cdots & 0 & \cdots & a_{1n} \\ \vdots & \ddots & \vdots & & \vdots \\ a_{j1} & \cdots & 1 & \cdots & a_{jn} \\ \vdots & & \vdots & \ddots & \vdots \\ a_{n1} & \cdots & 0 & \cdots & a{nn} \end{bmatrix} = (-1)^{j+i}\det(A_{ji}) = C_{ji}

于是,xij=Cjidet⁡(A)x_{ij} = \frac{C_{ji}}{\det(A)}。因此X=1det⁡(A)[Cji]=1det⁡(A)[Cij]TX = \frac{1}{\det(A)}[ C_{ji} ] = \frac{1}{\det(A)}[ C_{ij} ]^T。

其中

[Cji]=[Cij]T=[C11⋯Cn1⋮⋱⋮C1n⋯Cnn][C_{ji}] = [C_{ij}]^T = \begin{bmatrix} C_{11} & \cdots & C_{n1} \\ \vdots & \ddots & \vdots \\ C_{1n} & \cdots & C_{nn}\\ \end{bmatrix}

称为AA的伴随矩阵(英语:adjoint matrix),记作adj(A)\mathrm{adj}(A)。

矩阵的逆

设AA是一个n×nn\times n矩阵,如果AA可逆,那么有:

A−1=1det⁡(A)adj(A)A^{-1} = \frac{1}{\det(A)}\mathrm{adj}(A)

其中,adj(A)\mathrm{adj}(A)是AA的伴随矩阵。

例如,求解A=[2121−1312−2]A=\begin{bmatrix} 2 & 1 & 2 \\ 1 & -1 & 3 \\ 1 &2 & -2 \end{bmatrix}的逆。

首先,AA的各代数余子式为:

C11=+∣−132−2∣=−4,C12=−∣131−2∣=5,C13=+∣1−112∣=3,C21=−∣122−2∣=6,C22=+∣221−2∣=−6,C23=−∣2112∣=−3,C31=+∣12−13∣=5,C32=−∣2213∣=−4,C33=+∣211−1∣=−3\begin{aligned} C_{11} &= +\begin{vmatrix} -1 & 3 \\ 2 & -2 \end{vmatrix} = -4, \quad C_{12} = -\begin{vmatrix} 1 & 3 \\ 1 & -2 \end{vmatrix} = 5, \quad C_{13} = +\begin{vmatrix} 1 & -1 \\ 1 & 2 \end{vmatrix} = 3, \quad \\ C_{21} &= -\begin{vmatrix} 1 & 2 \\ 2 & -2 \end{vmatrix} = 6, \quad C_{22} = +\begin{vmatrix} 2 & 2 \\ 1 & -2 \end{vmatrix} = -6, \quad C_{23} = -\begin{vmatrix} 2 & 1 \\ 1 & 2 \end{vmatrix} = -3, \quad \\ C_{31} &= +\begin{vmatrix} 1 & 2 \\ -1 & 3 \end{vmatrix} = 5, \quad C_{32} = -\begin{vmatrix} 2 & 2 \\ 1 & 3 \end{vmatrix} = -4, \quad C_{33} = +\begin{vmatrix} 2 & 1 \\ 1 & -1 \end{vmatrix} = -3 \end{aligned}

因此,

det⁡(A)=2×(−4)+1×5+2×3=3\det(A) = 2 \times (-4) + 1 \times 5 + 2 \times 3 = 3

同时,AA的伴随矩阵为:

adj(A)=[−4536−6−35−4−3]T=[−4655−6−43−3−3]\mathrm{adj}(A) = \begin{bmatrix} -4 & 5 & 3 \\ 6 & -6 & -3 \\ 5 & -4 & -3 \end{bmatrix}^T = \begin{bmatrix} -4 & 6 & 5 \\ 5 & -6 & -4 \\ 3 & -3 & -3 \end{bmatrix}

因此,AA的逆为:

A−1=13[−4655−6−43−3−3]=[−4325353−2−431−1−1]A^{-1} = \frac{1}{3}\begin{bmatrix} -4 & 6 & 5 \\ 5 & -6 & -4 \\ 3 & -3 & -3 \end{bmatrix} = \begin{bmatrix} \frac{-4}{3} & 2 & \frac{5}{3} \\ \frac{5}{3} & -2 & -\frac{4}{3} \\ 1 & -1 & -1 \end{bmatrix}