Last updated on November 15, 2023 12:13 AM
矩阵的逆
设 A∈Mm×n(K),若 ∃B 使得 AB=I,则称 B 为 A 的右逆;若 ∃B 使得 BA=I,则称 B 为 A 的左逆。
何时有右逆?当且仅当 A 的列能表出 Km 的标准基 ε1,⋯,εm,故 A 行满秩,即 rank(A)=m。同理,矩阵有左逆当且仅当 A 列满秩,rank(A)=m。注意:左逆和右逆可能是不唯一的。
若方阵 A 的左/右逆均存在且相等,则将其记作 A−1,称 A 可逆。下面证明逆的唯一性。
假设 A,B∈Mn(K),满足 AB=I,下面证明 BA=I。ABA=A,A(BA−I)=0。由于 A 满秩,故 BA−I=0,即 BA=I。
番外:投影矩阵
现在考虑对于列满秩矩阵 A,构造其的一个左逆。
由于 rank(ATA)=rank(A),故 ATA 满秩,可逆。又由 (ATA)−1ATA=I,故 (ATA)−1AT 为 A 的一个左逆,令其为 A+。
这个东西的用处在于其可以将 β 投到 C(A) 中。比如对于一个无解的方程组 AX=β,我们要求他的一个近似解,自然我们需要找到一个“最近”的 γ 使得 γ∈C(A),然后解 AX=γ。而 AA+ 就是一个很好的投影矩阵。AA+β 即为 β 在 C(A) 上的正交投影。而取 X=A+β 作为解亦是自然的。
下面证明为何其为正交投影:AA+β∈C(A) 是自然的,然后我们用 A 中的列向量与 β−γ 做内积,这相当于用 AT 乘上他,就是 AT(β−AA+β)=AT(I−A(ATA)−1AT)β=AT(I−AA−1(AT)−1AT)β=0,证毕。
类似的,A 行满秩时,A+=AT(AAT)−1 为其最好的一个右逆。AA+=I,A+A 是向 A 行空间做投影的矩阵。
伴随矩阵法
将方阵 A 的每个元素 aij 替换为对应位置的代数余子式 Aij 再转置得到的方阵称作 A 的伴随矩阵,记作 A∗。
我们用矩阵乘法可以轻松看出 AA∗=A∗A=∣A∣I,故可以得到如下定理:
当 ∣A∣=0 时,A 可逆且 A−1=∣A∣1A∗。
比如说对于二元的情况,就有
[acbd]−1=ad−bc1[d−c−ba]
主对角线交换,副对角线取反。
这能帮助我们证明 Cramer 法则:
X=A−1β=∣A∣1⎣⎢⎢⎡A11⋮An1⋯⋯A1n⋮Ann⎦⎥⎥⎤⎣⎢⎢⎡b1⋮bn⎦⎥⎥⎤=∣A∣1⎣⎢⎢⎡A11b1+⋯+A1nbn⋮An1b1+⋯+Annbn⎦⎥⎥⎤=∣A∣1⎣⎢⎢⎡B1⋮Bn⎦⎥⎥⎤
并且还能得到一个有趣的小结论:整系数且行列式为 1 的矩阵的逆阵也是整系数。因为 ∣A∣∣A−1∣=∣AA−1∣=∣I∣=1,
初等变换求逆
在继续之前,先引出可逆方阵的一个乘法性质:
(AB)−1=B−1A−1
证明略过。
一般而言,计算矩阵的逆不使用伴随矩阵法,计算量太大。我们基于任何满秩方阵都能分解为 P1⋯Ps 的事实,我们若能找到一系列初等矩阵满足 P1P2⋯PsA=I,那么是不是就有 A−1=P1P2⋯Ps。
实际操作的时候将 I 和 A 写在一起:[A∣I],然后将 A 消成单位阵,消完后右半边就是 A−1。
同理,解矩阵方程 AX=B 的时候,显然 X=A−1B,将 B 放在 A 右侧一起消元,消完后 B 就变成 A−1B 了。
比如求 XB=C,做一个转置便有 BTXT=CT,按照上述方法解即可。再比如 AXB=C,先把 XB 合在一起,求出 XB,然后再解 X 即可。
最后我们可以得到如下总结:
定理:对于 A∈Mn(K),以下五个条件等价:
- A 可逆;
- ∃B s. t. AB=I(或 BA=I);
- rank(A)=n;
- ∣A∣=0;
- A 可以分解为若干初等矩阵的乘积。
下面是一些例题。
例题:已知 n 阶方阵 A 和 B,且 A+B=AB,证明 A 与 B 可交换。
证明:先移项,然后配凑成两个因式之积:AB−A−B=0,所以 AB−A−B+I=I,即 (A−I)(B−I)=I,所以 (A−I) 与 (B−I) 互为逆矩阵,自然可交换,即 (B−I)(A−I)=I,展开得到 BA=A+B=AB,证毕。