线代 A1 笔记 第 8 部分(矩阵的逆)

Last updated on November 15, 2023 12:13 AM

矩阵的逆

AMm×n(K)A\in M_{m\times n}(K),若 B\exists B 使得 AB=IAB=I,则称 BBAA 的右逆;若 B\exists B 使得 BA=IBA=I,则称 BBAA 的左逆。

何时有右逆?当且仅当 AA 的列能表出 KmK^m 的标准基 ε1,,εm\varepsilon_1, \cdots ,\varepsilon_m,故 AA 满秩,即 rank(A)=m\operatorname{rank}(A) = m。同理,矩阵有左逆当且仅当 AA 满秩,rank(A)=m\operatorname{rank}(A)=m。注意:左逆和右逆可能是不唯一的。

若方阵 AA 的左/右逆均存在且相等,则将其记作 A1A^{-1},称 AA 可逆。下面证明逆的唯一性。

假设 A,BMn(K)A,B\in M_n(K),满足 AB=IAB=I,下面证明 BA=IBA=IABA=AABA=AA(BAI)=0A(BA-I)=0。由于 AA 满秩,故 BAI=0BA-I=0,即 BA=IBA=I

番外:投影矩阵

现在考虑对于列满秩矩阵 AA,构造其的一个左逆。

由于 rank(ATA)=rank(A)\operatorname{rank}(A^TA)=\operatorname{rank}(A),故 ATAA^TA 满秩,可逆。又由 (ATA)1ATA=I(A^TA)^{-1}A^TA = I,故 (ATA)1AT(A^TA)^{-1}A^TAA 的一个左逆,令其为 A+A^+

这个东西的用处在于其可以将 β\beta 投到 C(A)\mathcal{C}(A) 中。比如对于一个无解的方程组 AX=βAX=\beta,我们要求他的一个近似解,自然我们需要找到一个“最近”的 γ\gamma 使得 γC(A)\gamma\in \mathcal{C}(A),然后解 AX=γAX=\gamma。而 AA+A A^+ 就是一个很好的投影矩阵。AA+βA A^+\beta 即为 β\betaC(A)\mathcal{C}(A) 上的正交投影。而取 X=A+βX=A^+\beta 作为解亦是自然的。

下面证明为何其为正交投影:AA+βC(A)A A^+\beta\in \mathcal{C}(A) 是自然的,然后我们用 AA 中的列向量与 βγ\beta - \gamma 做内积,这相当于用 ATA^T 乘上他,就是 AT(βAA+β)=AT(IA(ATA)1AT)β=AT(IAA1(AT)1AT)β=0A^T(\beta - A A^+\beta)=A^T(I - A(A^TA)^{-1}A^T)\beta=A^T(I - A A^{-1}(A^T)^{-1}A^T)\beta = 0,证毕。

类似的,AA 行满秩时,A+=AT(AAT)1A^+ = A^T(AA^T)^{-1} 为其最好的一个右逆。AA+=IA A^+=IA+AA^+A 是向 AA 行空间做投影的矩阵。

伴随矩阵法

将方阵 AA 的每个元素 aija_{ij} 替换为对应位置的代数余子式 AijA_{ij} 再转置得到的方阵称作 AA伴随矩阵,记作 AA^*

我们用矩阵乘法可以轻松看出 AA=AA=AIA A^* = A^*A= |A|I,故可以得到如下定理:

A0|A|\ne 0 时,AA 可逆且 A1=1AA\displaystyle A^{-1}=\frac{1}{|A|}A^*

比如说对于二元的情况,就有

[abcd]1=1adbc[dbca]\begin{bmatrix} a & b \\ c & d \\\end{bmatrix}^{-1} = \frac{1}{ad-bc} \begin{bmatrix} d & -b \\ -c & a \\\end{bmatrix}

主对角线交换,副对角线取反

这能帮助我们证明 Cramer 法则:

X=A1β=1A[A11A1nAn1Ann][b1bn]=1A[A11b1++A1nbnAn1b1++Annbn]=1A[B1Bn]\begin{aligned} X &= A^{-1}\beta = \frac{1}{|A|}\begin{bmatrix} A_{11} & \cdots & A_{1n} \\ \vdots & & \vdots \\ A_{n1} & \cdots & A_{nn} \\\end{bmatrix}\begin{bmatrix} b_1 \\ \vdots \\ b_n \\\end{bmatrix}\\ &= \frac{1}{|A|}\begin{bmatrix} A_{11}b_1 + \cdots + A_{1n}b_n \\ \vdots \\ A_{n 1}b_1 + \cdots + A_{nn}b_n \\\end{bmatrix} = \frac{1}{|A|}\begin{bmatrix} B_1 \\ \vdots \\ B_n \\\end{bmatrix} \end{aligned}

并且还能得到一个有趣的小结论:整系数且行列式为 11 的矩阵的逆阵也是整系数。因为 AA1=AA1=I=1|A||A^{-1}| = |A A^{-1}|=|I|=1

初等变换求逆

在继续之前,先引出可逆方阵的一个乘法性质:

(AB)1=B1A1(AB)^{-1} = B^{-1}A^{-1}

证明略过。

一般而言,计算矩阵的逆不使用伴随矩阵法,计算量太大。我们基于任何满秩方阵都能分解为 P1PsP_1\cdots P_s 的事实,我们若能找到一系列初等矩阵满足 P1P2PsA=IP_1P_2\cdots P_sA = I,那么是不是就有 A1=P1P2PsA^{-1} = P_1P_2\cdots P_s

实际操作的时候将 IIAA 写在一起:[AI]\begin{bmatrix} A \mid I \\\end{bmatrix},然后将 AA 消成单位阵,消完后右半边就是 A1A^{-1}

同理,解矩阵方程 AX=BAX=B 的时候,显然 X=A1BX = A^{-1}B,将 BB 放在 AA 右侧一起消元,消完后 BB 就变成 A1BA^{-1}B 了。

比如求 XB=CXB=C,做一个转置便有 BTXT=CTB^TX^T=C^T,按照上述方法解即可。再比如 AXB=CAXB=C,先把 XBXB 合在一起,求出 XBXB,然后再解 XX 即可。

最后我们可以得到如下总结:

定理:对于 AMn(K)A\in M_n(K),以下五个条件等价:

  1. AA 可逆;
  2. B\exists B s. t. AB=IAB=I(或 BA=IBA=I);
  3. rank(A)=n\operatorname{rank}(A) = n
  4. A0|A|\ne 0
  5. AA 可以分解为若干初等矩阵的乘积。

下面是一些例题。

例题:已知 nn 阶方阵 AABB,且 A+B=ABA+B=AB,证明 AABB 可交换。
证明:先移项,然后配凑成两个因式之积:ABAB=0AB-A-B=0,所以 ABAB+I=IAB-A-B+I=I,即 (AI)(BI)=I(A-I)(B-I)=I,所以 (AI)(A-I)(BI)(B-I) 互为逆矩阵,自然可交换,即 (BI)(AI)=I(B-I)(A-I)=I,展开得到 BA=A+B=ABBA=A+B=AB,证毕。


线代 A1 笔记 第 8 部分(矩阵的逆)
https://blog.imyangty.com/note-linear-algebra-a1/8/
Author
YangTY
Posted on
November 15, 2023
Licensed under