线代 A1 笔记 第 10 部分(矩阵逆/分块 习题课)

Last updated on November 21, 2023 3:49 PM

P1

Statement

已知 AA 可逆,问 A,B,C,DA,B,C,D 满足什么条件的时候 [ABCD]\begin{bmatrix} A & B \\ C & D \\\end{bmatrix} 可逆?求之。

Solution

类似于普通矩阵求逆,一步步初等变换吧。先用第一行左乘 CA1-CA^{-1} 再加到第二行上,得到

[ABI00DCA1BCA1I]\begin{bmatrix} A & B & I & 0 \\ 0 & D-CA^{-1}B & -CA^{-1} & I \\\end{bmatrix}

再把 AA 变成单位阵,即第一行左乘 A1A^{-1}

[IA1BA100DCA1BCA1I]\begin{bmatrix} I & A^{-1}B & A^{-1} & 0 \\ 0 & D-CA^{-1}B & -CA^{-1} & I \\\end{bmatrix}

接下来需要消去 A1BA^{-1}B,还要让 E=DCA1BE=D-CA^{-1}B 变成单位阵,这就需要 E=DCA1BE = D-CA^{-1}B 可逆。第二行乘上 E1E^{-1},得到

[IA1BA100IE1CA1E1]\begin{bmatrix} I & A^{-1}B & A^{-1} & 0 \\ 0 & I & -E^{-1}CA^{-1} & E^{-1} \\\end{bmatrix}

第二行左乘上 A1B-A^{-1}B 加到第一行上:

[I0A1+A1BE1CA1A1BE10IE1CA1E1]\begin{bmatrix} I & 0 & A^{-1}+A^{-1}BE^{-1}CA^{-1} & -A^{-1}BE^{-1} \\ 0 & I & -E^{-1}CA^{-1} & E^{-1} \\\end{bmatrix}

所以

[ABCD]1=[A1+A1BE1CA1A1BE1E1CA1E1]=[A1+A1B(DCA1B)1CA1A1B(DCA1B)1(DCA1B)1CA1(DCA1B)1]\begin{aligned} \begin{bmatrix} A & B \\ C & D \\\end{bmatrix}^{-1} &= \begin{bmatrix} A^{-1}+A^{-1}BE^{-1}CA^{-1} & -A^{-1}BE^{-1} \\ -E^{-1}CA^{-1} & E^{-1} \\\end{bmatrix} \\ &= \begin{bmatrix} A^{-1}+A^{-1}B(D-CA^{-1}B)^{-1}CA^{-1} & -A^{-1}B(D-CA^{-1}B)^{-1} \\ -(D-CA^{-1}B)^{-1}CA^{-1} & (D-CA^{-1}B)^{-1} \\\end{bmatrix} \end{aligned}

P2

Statement

求如下矩阵的逆:

A=[123n012n100120001]A= \begin{bmatrix} 1 & 2 & 3 & \cdots & n \\ 0 & 1 & 2 & \cdots & n-1 \\ 0 & 0 & 1 & \cdots & \vdots \\ \vdots & \vdots & \vdots & & 2 \\ 0 & 0 & \cdots & 0 & 1 \\\end{bmatrix}

Solution

观察发现若设 JJnn 阶若当块,则

A=I+2J+3J2++nJn1A= I+2J+3J^2+ \cdots +nJ^{n-1}

我们用错位相减法:

AJ=J+2J2+3J3++(n1)Jn1AAJ=I+J+J2++Jn1A(IJ)2=(IJ)(I+J+J2++Jn1)=I\begin{aligned} AJ &= J + 2J^2 + 3J^3 + \cdots + (n-1)J^{n-1} \\ A - AJ &= I + J + J^2 + \cdots + J^{n-1}\\ A(I-J)^{2} &= (I-J)(I+J+J^{2}+ \cdots +J^{n-1})\\ &= I \end{aligned}

所以 A1=(IJ)2A^{-1} = (I-J)^2

P3

Statement

A,B,C,DA,B,C,D 都为 nn 阶方阵,且 AC=CAAC=CA,问是否一定有

ABCD=ADCB\begin{vmatrix} A & B \\ C & D \\\end{vmatrix} = |AD-CB|

Solution

这题没说 AA 一定可逆,所以会比较棘手。

先考虑 AA 可逆的情形,此时我们对其施加初等变换可以知道

[ICA1I][ABCD]=[ABDCA1B]ABCD=ABDCA1B=ADCA1B=ADACA1B=ADCB\begin{aligned} \begin{bmatrix} I & \\ -CA^{-1} & I \\\end{bmatrix} \begin{bmatrix} A & B \\ C & D \\\end{bmatrix} &= \begin{bmatrix} A & B \\ & D-CA^{-1}B \\\end{bmatrix} \\ \begin{vmatrix} A & B \\ C & D \\\end{vmatrix} &= \begin{vmatrix} A & B \\ & D-CA^{-1}B \\\end{vmatrix}\\ &= |A||D-CA^{-1}B|\\ &= |AD - ACA^{-1}B|\\ &= |AD - CB| \end{aligned}

最后一步用到了 AC=CAAC=CA 的性质。

AA 不可逆的情形就比较麻烦了,考虑摄动法。考虑 A+tIA+tI(A+tI)C=AC+tC=CA+tC=C(A+tI)(A+tI)C = AC + tC = CA + tC = C(A+tI),可交换性仍满足。A+tI|A+tI| 为关于 ttnn 阶多项式,根是有限个的。

那么令 δ=minti>0{ti}\delta = \displaystyle \min_{t_i> 0}\left\{ t_i \right\},其中 tit_iA+tI=0A+tI=0 的根。

则在 t(0,δ)t\in(0,\delta) 这个区间内,对于任意小的 tt 我们都有 A+tIBCD=(A+tI)DCB\begin{vmatrix} A+tI & B \\ C & D \\\end{vmatrix} = |(A+tI)D - CB| 成立。注意到当 t0t\to 0 时,两边分别会收敛到 ABCD\begin{vmatrix} A & B \\ C & D \\\end{vmatrix}ADCB|AD-CB|,得证。

P4

Statement

EE 的逆,其中 i=1nai0\displaystyle \prod_{i=1}^{n} a_i\ne 0

E=[1+a11111+a21111+an]E = \begin{bmatrix} 1+a_1 & 1 & \cdots & 1 \\ 1 & 1+a_2 & \cdots & 1 \\ \vdots & \vdots & & \vdots \\ 1 & 1 & \cdots & 1+a_n \\\end{bmatrix}

Solution

注意到 EE 可以拆成:

E=diag{a1,,an}(In+[1/a11/a21/an][11])E = \operatorname{diag}\left\{ a_1, \cdots ,a_n \right\}\left( I_n + \begin{bmatrix} 1 / a_1 \\ 1 / a_2 \\ \vdots \\ 1 / a_n \\\end{bmatrix}\begin{bmatrix} 1 & \cdots & 1 \\\end{bmatrix} \right)

右边的式子符合 IABI - AB 的形式,所以改号拉开填逆,我们注意到 I+[11][1/a11/a21/an]I + \begin{bmatrix} 1 & \cdots & 1 \\\end{bmatrix}\begin{bmatrix} 1 / a_1 \\ 1 / a_2 \\ \vdots \\ 1 / a_n \\\end{bmatrix} 的逆就是 (1+1a1++1an)1\displaystyle \left( 1 + \frac{1}{a_1}+ \cdots +\frac{1}{a_n} \right) ^{-1},于是

E1=(In(1+1a1++1an)1[1/a11/a21/an][11])diag{1a1,,1an}E^{-1} = \left( I_n - \left( 1 + \frac{1}{a_1}+ \cdots +\frac{1}{a_n} \right) ^{-1} \begin{bmatrix} 1 / a_1 \\ 1 / a_2 \\ \vdots \\ 1 / a_n \\\end{bmatrix}\begin{bmatrix} 1 & \cdots & 1 \\\end{bmatrix}\right) \operatorname{diag}\left\{ \frac{1}{a_1}, \cdots ,\frac{1}{a_n} \right\}

发现 (InBA)1=In+B(ImAB)1A(I_n - BA)^{-1} = I_n + B(I_m - AB)^{-1}A 还是很好用的。

事已至此先睡觉吧。

P5

Statement

AAn (n2)n~(n\ge 2) 级方阵,AA^{*}AA 的伴随矩阵,证明:

rank(A)={n,rank(A)=n1,rank(A)=n10,rank(A)<n1\operatorname{rank}(A^{*})=\begin{cases} n, & \operatorname{rank}(A)=n \\ 1, & \operatorname{rank}(A)=n-1\\0,&\operatorname{rank}(A)<n-1 \end{cases}

Solution

rank(A)=n\operatorname{rank}(A)=n 时,由 AA=AIA A^{*} = |A|I,知 AA=An|A||A^{*}| = |A|^n,即 A=An10|A^{*}|=|A|^{n-1}\ne 0,故 rank(A)=n\operatorname{rank}(A^{*}) = n

rank(A)<n1\operatorname{rank}(A)<n-1 时,AA 的任意 n1n-1 阶子式均为 00,所以 A=0A^{*} = 0rank(A)=0\operatorname{rank}(A^{*}) = 0

比较棘手的是 rank(A)=n1\operatorname{rank}(A) = n -1 的情况。显然存在至少一个 n1n-1 阶子式不为 00,故 rank(A)1\operatorname{rank}(A^{*})\ge 1,接下来证明 rank(A)1\operatorname{rank}(A^{*})\le 1 即可。

由于 AA=AI=0A A^{*}=|A|I = 0,故 rank(A)+rank(A)n\operatorname{rank}(A) + \operatorname{rank}(A^{*})\le n,由于 rank(A)=n1\operatorname{rank}(A) = n-1,故 rank(A)1\operatorname{rank}(A^{*})\le 1,证毕。

P6

Statement

证明 Frobenius 不等式:

A,B,CA,B,C 分别为 m×n,n×k,k×sm\times n,n\times k,k\times s 矩阵,则 rank(ABC)+rank(B)rank(AB)+rank(BC)\operatorname{rank}(ABC)+\operatorname{rank}(B)\ge \operatorname{rank}(AB) + \operatorname{rank}(BC)

注:当 B=InB=I_n 时,上式退化为 rank(A)+rank(C)rank(AC)+n\operatorname{rank}(A)+\operatorname{rank}(C)\le \operatorname{rank}(AC) + n

A=CA=C 为方阵且 B=Ak1B =A^{k-1} 时,上式变为 rank(Ak1)rank(Ak)rank(Ak)rank(Ak+1)\operatorname{rank}(A^{k-1})-\operatorname{rank}(A^k)\ge \operatorname{rank}(A^k)-\operatorname{rank}(A^{k+1})。即,次数越高,相邻两项的秩的差越小

Solution

直接用分块初等变换即可。我们有

rank(ABC)+rank(B)=rank[ABC00B]\operatorname{rank}(ABC)+\operatorname{rank}(B)=\operatorname{rank} \begin{bmatrix} ABC & 0 \\ 0 & B \\\end{bmatrix}

第二行左乘 A-A 加到第一行,然后第二列右乘 CC 加到第一列:

[ImAIn][ABCB][IsCIk]=[ABBCB]\begin{bmatrix} I_m & A \\ & I_n \\\end{bmatrix} \begin{bmatrix} ABC & \\ & B \\\end{bmatrix} \begin{bmatrix} I_s & \\ C & I_k \\\end{bmatrix} = \begin{bmatrix} & -AB \\ BC & B \\\end{bmatrix}

注意到 rank[ABBCB]rank(BC)+rank(AB)\operatorname{rank} \begin{bmatrix} & -AB \\ BC & B \\\end{bmatrix}\ge \operatorname{rank}(BC) + \operatorname{rank}(AB),证毕。


线代 A1 笔记 第 10 部分(矩阵逆/分块 习题课)
https://blog.imyangty.com/note-linear-algebra-a1/10/
Author
YangTY
Posted on
November 20, 2023
Licensed under