15-矩阵

2023 字
10 分钟
15-矩阵

公式#

操作顺序可交换#

(AT)1=(A1)T(A)1=(A1)(A)T=(AT)\begin{aligned} & (A^T)^{-1} = (A^{-1})^{T} \\ & (A^*)^{-1} = (A^{-1})^{*} \\ & (A^*)^{T} = (A^{T})^{*} \end{aligned}

T,*,-1的运算顺序可任意交换

操作倒置#

(AB)1=B1A1(AB)=BA(AB)T=BTAT\begin{aligned} & (AB)^{-1} = B^{-1}A^{-1} \\ & (AB)^{*} = B^{*}A^{*} \\ & (AB)^{T} = B^{T}A^{T} \end{aligned}

本身为自己的逆运算#

(AT)T=A(A1)1=A\begin{aligned} & (A^T)^T = A \\ & (A^{-1})^{-1} = A \end{aligned}

转置的优良性#

AT{AT=A(kA)T=kAT(A+B)T=AT+BTA^T \Rightarrow \left\{ \begin{aligned} & |A^T| = |A| \\ & (kA)^T = kA^T \\ & (A + B)^T = A^T + B^T \end{aligned} \right.

逆矩阵#

A1A=E{A1=1A(kA)1=1kA1A^{-1}A = E \Rightarrow \left\{ \begin{aligned} & |A^{-1}| = \frac{1}{|A|} \\ & (kA)^{-1} = \frac{1}{k}A^{-1} \end{aligned} \right.

伴随矩阵#

A=AA1{(kA)=kn1AA=An1(A)=An2AA^{*} = |A|A^{-1} \Rightarrow \left\{ \begin{aligned} & (kA)^{*} = k^{n-1}A^{*} \\ & |A^*| = |A|^{n-1} \\ & (A^*)^* = |A|^{n-2}A \end{aligned} \right.

分块矩阵#

把矩阵按若干行、若干列切开,每一块当成一个整体参与运算。 例如

A=(A11A12A21A22)A= \begin{pmatrix} A_{11}&A_{12}\\ A_{21}&A_{22} \end{pmatrix}

其中每个 AijA_{ij} 本身可以是矩阵。

Note

分块以后要先检查每一块的阶数。块与块之间能不能加、能不能乘,本质上还是看普通矩阵的维数是否匹配。

加法与数乘#

两个矩阵按相同方式分块时,才能按块相加:

(A11A12A21A22)+(B11B12B21B22)=(A11+B11A12+B12A21+B21A22+B22)\begin{pmatrix} A_{11}&A_{12}\\ A_{21}&A_{22} \end{pmatrix} + \begin{pmatrix} B_{11}&B_{12}\\ B_{21}&B_{22} \end{pmatrix} = \begin{pmatrix} A_{11}+B_{11}&A_{12}+B_{12}\\ A_{21}+B_{21}&A_{22}+B_{22} \end{pmatrix}

数乘直接乘到每一块:

kA=(kA11kA12kA21kA22)kA= \begin{pmatrix} kA_{11}&kA_{12}\\ kA_{21}&kA_{22} \end{pmatrix}

乘法#

分块乘法和普通矩阵乘法形式一样,但每个元素换成了矩阵块:

(A11A12A21A22)(B11B12B21B22)=(A11B11+A12B21A11B12+A12B22A21B11+A22B21A21B12+A22B22)\begin{pmatrix} A_{11}&A_{12}\\ A_{21}&A_{22} \end{pmatrix} \begin{pmatrix} B_{11}&B_{12}\\ B_{21}&B_{22} \end{pmatrix} = \begin{pmatrix} A_{11}B_{11}+A_{12}B_{21}&A_{11}B_{12}+A_{12}B_{22}\\ A_{21}B_{11}+A_{22}B_{21}&A_{21}B_{12}+A_{22}B_{22} \end{pmatrix}

前提是中间维数匹配,例如 A11B11A_{11}B_{11}A12B21A_{12}B_{21} 都能相乘,并且结果阶数相同才能相加。

转置#

分块矩阵转置时,块的位置转置,每一块也要转置:

(ABCD)T=(ATCTBTDT)\begin{pmatrix} A&B\\ C&D \end{pmatrix}^{T} = \begin{pmatrix} A^T&C^T\\ B^T&D^T \end{pmatrix}

分块对角矩阵#

若每个 AiA_i 都是方阵,且

A=(A1000A2000As)A= \begin{pmatrix} A_1&0&\cdots&0\\ 0&A_2&\cdots&0\\ \vdots&\vdots&\ddots&\vdots\\ 0&0&\cdots&A_s \end{pmatrix}

A=A1A2As|A|=|A_1||A_2|\cdots |A_s|

若每个 AiA_i 都可逆,则

A1=(A11000A21000As1)A^{-1}= \begin{pmatrix} A_1^{-1}&0&\cdots&0\\ 0&A_2^{-1}&\cdots&0\\ \vdots&\vdots&\ddots&\vdots\\ 0&0&\cdots&A_s^{-1} \end{pmatrix}

分块三角矩阵#

若主对角线上的块都是方阵,则

AB0D=AD,A0CD=AD\begin{vmatrix} A&B\\ 0&D \end{vmatrix} =|A||D|,\qquad \begin{vmatrix} A&0\\ C&D \end{vmatrix} =|A||D|

对应的逆矩阵:

(AB0D)1=(A1A1BD10D1)\begin{pmatrix} A&B\\ 0&D \end{pmatrix}^{-1} = \begin{pmatrix} A^{-1}&-A^{-1}BD^{-1}\\ 0&D^{-1} \end{pmatrix}(A0CD)1=(A10D1CA1D1)\begin{pmatrix} A&0\\ C&D \end{pmatrix}^{-1} = \begin{pmatrix} A^{-1}&0\\ -D^{-1}CA^{-1}&D^{-1} \end{pmatrix}

其中 A,DA,D 均需可逆。

Schur 补#

对二阶分块矩阵

M=(ABCD)M= \begin{pmatrix} A&B\\ C&D \end{pmatrix}

AA 可逆,则

M=ADCA1B|M|=|A|\cdot |D-CA^{-1}B|

其中

DCA1BD-CA^{-1}B

叫做关于 AA 的 Schur 补。

DD 可逆,也可以写成

M=DABD1C|M|=|D|\cdot |A-BD^{-1}C|
Tip

不要把 M|M| 直接写成 ADBC|A||D|-|B||C|。一般矩阵块不能这样拆,只有低阶数值行列式才有类似形式。

分块初等变换#

当块大小匹配时,可以像普通行列变换一样做分块变换:

(EOCA1E)(ABCD)=(ABODCA1B)\begin{pmatrix} E&O\\ -CA^{-1}&E \end{pmatrix} \begin{pmatrix} A&B\\ C&D \end{pmatrix} = \begin{pmatrix} A&B\\ O&D-CA^{-1}B \end{pmatrix}

这就是 Schur 补公式的来源:左边乘的分块矩阵行列式为 11,所以行列式不变。

常见用途:

  • 化成分块三角矩阵
  • 计算含零块的行列式
  • 求分块矩阵的逆
  • 处理线性方程组中的变量消元

分块矩阵方程组#

当系数矩阵是分块矩阵时,可以把未知向量和常数向量也按相同方式分块。

例如

(ABCD)(x1x2)=(b1b2)\begin{pmatrix} A&B\\ C&D \end{pmatrix} \begin{pmatrix} x_1\\ x_2 \end{pmatrix} = \begin{pmatrix} b_1\\ b_2 \end{pmatrix}

等价于方程组

{Ax1+Bx2=b1Cx1+Dx2=b2\left\{ \begin{aligned} Ax_1+Bx_2&=b_1\\ Cx_1+Dx_2&=b_2 \end{aligned} \right.
Note

分块以后每个乘积都要满足维数匹配,例如 Ax1Ax_1Bx2Bx_2 必须是同型列向量,才能相加等于 b1b_1

分块三角方程组#

(ABOD)(x1x2)=(b1b2)\begin{pmatrix} A&B\\ O&D \end{pmatrix} \begin{pmatrix} x_1\\ x_2 \end{pmatrix} = \begin{pmatrix} b_1\\ b_2 \end{pmatrix}

A,DA,D 都可逆,则先解第二个方程:

Dx2=b2Dx_2=b_2

得到

x2=D1b2x_2=D^{-1}b_2

再代回第一个方程:

Ax1+Bx2=b1Ax_1+Bx_2=b_1

得到

x1=A1(b1BD1b2)x_1=A^{-1}(b_1-BD^{-1}b_2)

所以

{x2=D1b2x1=A1(b1BD1b2)\left\{ \begin{aligned} x_2&=D^{-1}b_2\\ x_1&=A^{-1}(b_1-BD^{-1}b_2) \end{aligned} \right.

同理,若是下三角分块矩阵

(AOCD)(x1x2)=(b1b2)\begin{pmatrix} A&O\\ C&D \end{pmatrix} \begin{pmatrix} x_1\\ x_2 \end{pmatrix} = \begin{pmatrix} b_1\\ b_2 \end{pmatrix}

则先解

Ax1=b1Ax_1=b_1

再解

Cx1+Dx2=b2Cx_1+Dx_2=b_2

得到

{x1=A1b1x2=D1(b2CA1b1)\left\{ \begin{aligned} x_1&=A^{-1}b_1\\ x_2&=D^{-1}(b_2-CA^{-1}b_1) \end{aligned} \right.

Schur 补消元#

对一般二阶分块方程组

{Ax1+Bx2=b1Cx1+Dx2=b2\left\{ \begin{aligned} Ax_1+Bx_2&=b_1\\ Cx_1+Dx_2&=b_2 \end{aligned} \right.

AA 可逆,可以先由第一个方程表示

x1=A1(b1Bx2)x_1=A^{-1}(b_1-Bx_2)

代入第二个方程:

CA1(b1Bx2)+Dx2=b2CA^{-1}(b_1-Bx_2)+Dx_2=b_2

整理得

(DCA1B)x2=b2CA1b1(D-CA^{-1}B)x_2=b_2-CA^{-1}b_1

其中

DCA1BD-CA^{-1}B

就是关于 AA 的 Schur 补。

如果 DCA1BD-CA^{-1}B 也可逆,则

{x2=(DCA1B)1(b2CA1b1)x1=A1(b1Bx2)\left\{ \begin{aligned} x_2&=(D-CA^{-1}B)^{-1}(b_2-CA^{-1}b_1)\\ x_1&=A^{-1}(b_1-Bx_2) \end{aligned} \right.

DD 可逆,也可以先消去 x2x_2。由第二个方程

x2=D1(b2Cx1)x_2=D^{-1}(b_2-Cx_1)

代入第一个方程,得到

(ABD1C)x1=b1BD1b2(A-BD^{-1}C)x_1=b_1-BD^{-1}b_2

其中

ABD1CA-BD^{-1}C

是关于 DD 的 Schur 补。

如果 ABD1CA-BD^{-1}C 也可逆,则

{x1=(ABD1C)1(b1BD1b2)x2=D1(b2Cx1)\left\{ \begin{aligned} x_1&=(A-BD^{-1}C)^{-1}(b_1-BD^{-1}b_2)\\ x_2&=D^{-1}(b_2-Cx_1) \end{aligned} \right.
Tip

解分块方程组的核心不是把每一块当成数,而是用可逆块做消元。选哪个块消元,通常看 AADD 哪个更容易求逆。

#

证明一般由

  1. 矩阵子式不等于0的最高阶
  2. 线性表示
  3. 极大线性无关组
  4. 转化为分块矩阵使用二级结论

进行证明,

子式#

在矩阵里任选k行和k列组成的行列式为k阶子式,存在等于矩阵秩的非0子式,大于矩阵秩的子式全为0

不等式证明#

下面默认矩阵阶数都满足运算要求。记 r(A)r(A) 为矩阵 AA 的秩。

基本不等式#

AAm×nm\times n 矩阵,则

0r(A)min{m,n}0\le r(A)\le \min\{m,n\}

并且

r(A)=r(AT)r(A)=r(A^T)

因为行秩等于列秩,转置只是在交换行空间和列空间。

加法不等式#

A,BA,B 同型,则

r(A+B)r(A)+r(B)r(A+B)\le r(A)+r(B)

证明

A,B,A+BA,B,A+B 按列写成

A=(α1,α2,,αn),B=(β1,β2,,βn)A=(\alpha_1,\alpha_2,\cdots,\alpha_n),\qquad B=(\beta_1,\beta_2,\cdots,\beta_n)

A+B=(α1+β1,α2+β2,,αn+βn)A+B=(\alpha_1+\beta_1,\alpha_2+\beta_2,\cdots,\alpha_n+\beta_n)

每一列 αi+βi\alpha_i+\beta_i 都能由向量组

α1,,αn,β1,,βn\alpha_1,\cdots,\alpha_n,\beta_1,\cdots,\beta_n

线性表示。

所以 A+BA+B 的列向量组能由 A,BA,B 的全部列向量组成的向量组线性表示,因此

r(A+B)r(AB)r(A)+r(B)r(A+B)\le r\begin{pmatrix}A&B\end{pmatrix}\le r(A)+r(B)

乘法上界#

AAm×nm\times n 矩阵,BBn×pn\times p 矩阵,则

r(AB)min{r(A),r(B)}r(AB)\le \min\{r(A),r(B)\}

证明一:证明 r(AB)r(A)r(AB)\le r(A)

BB 按列写成

B=(β1,β2,,βp)B=(\beta_1,\beta_2,\cdots,\beta_p)

AB=(Aβ1,Aβ2,,Aβp)AB=(A\beta_1,A\beta_2,\cdots,A\beta_p)

AβiA\beta_iAA 的列向量的线性组合,所以 ABAB 的每一列都能由 AA 的列向量组线性表示。

因此

r(AB)r(A)r(AB)\le r(A)

证明二:证明 r(AB)r(B)r(AB)\le r(B)

转置后

r(AB)=r((AB)T)=r(BTAT)r(AB)=r((AB)^T)=r(B^TA^T)

由证明一可知

r(BTAT)r(BT)=r(B)r(B^TA^T)\le r(B^T)=r(B)

所以

r(AB)r(B)r(AB)\le r(B)

分块上三角矩阵秩不等式#

r(ACOB)r(A)+r(B)r\begin{pmatrix} A&C\\ O&B \end{pmatrix}\ge r(A)+r(B)

证明

r(A)=ra,r(B)=rbr(A)=r_a,\qquad r(B)=r_b

A,BA,B 分别存在一个 ra,rbr_a,r_b 级满秩子矩阵 A1,B1A_1,B_1,即

A10,B10|A_1|\ne 0,\qquad |B_1|\ne 0

D=(ACOB)D= \begin{pmatrix} A&C\\ O&B \end{pmatrix}

中,取出 A1A_1 所在的行列和 B1B_1 所在的行列,同时在右上角取对应的子块 C1C_1,得到一个 ra+rbr_a+r_b 阶子矩阵

(A1C1OB1)\begin{pmatrix} A_1&C_1\\ O&B_1 \end{pmatrix}

它是分块上三角矩阵,所以

A1C1OB1=A1B10\begin{vmatrix} A_1&C_1\\ O&B_1 \end{vmatrix} =|A_1||B_1|\ne 0

因此 DD 存在一个 ra+rbr_a+r_b 阶非零子式,所以

r(D)ra+rb=r(A)+r(B)r(D)\ge r_a+r_b=r(A)+r(B)

Sylvester 不等式#

AAm×nm\times n 矩阵,BBn×pn\times p 矩阵,则

r(AB)r(A)+r(B)nr(AB)\ge r(A)+r(B)-n

这是乘法秩的下界,和上界合起来是

r(A)+r(B)nr(AB)min{r(A),r(B)}r(A)+r(B)-n\le r(AB)\le \min\{r(A),r(B)\}

证明

构造分块矩阵

D=(AOEnB)D= \begin{pmatrix} A&O\\ E_n&B \end{pmatrix}

其中 EnE_nnn 阶单位矩阵。因为

r(D)=r(DT)=r(ATEnOBT)r(AT)+r(BT)=r(A)+r(B)r(D)=r(D^T) =r\begin{pmatrix} A^T&E_n\\ O&B^T \end{pmatrix} \ge r(A^T)+r(B^T)=r(A)+r(B)

所以

r(D)r(A)+r(B)r(D)\ge r(A)+r(B)

另一方面,对 DD 做分块初等变换:

(EmAOEn)(AOEnB)=(OABEnB)\begin{pmatrix} E_m&-A\\ O&E_n \end{pmatrix} \begin{pmatrix} A&O\\ E_n&B \end{pmatrix} = \begin{pmatrix} O&-AB\\ E_n&B \end{pmatrix}

再对列做分块初等变换:

(OABEnB)(EnBOEp)=(OABEnO)\begin{pmatrix} O&-AB\\ E_n&B \end{pmatrix} \begin{pmatrix} E_n&-B\\ O&E_p \end{pmatrix} = \begin{pmatrix} O&-AB\\ E_n&O \end{pmatrix}

这些变换都不改变秩。把最后一个矩阵交换分块行、分块列,可化成

(EnOOAB)\begin{pmatrix} E_n&O\\ O&-AB \end{pmatrix}

所以

r(D)=n+r(AB)r(D)=n+r(AB)

于是

r(A)+r(B)r(D)=n+r(AB)r(A)+r(B)\le r(D)=n+r(AB)

r(AB)r(A)+r(B)nr(AB)\ge r(A)+r(B)-n

Frobenius 不等式#

ABCABC 有意义,则

r(ABC)r(AB)+r(BC)r(B)r(ABC)\ge r(AB)+r(BC)-r(B)

证明

考虑分块对角矩阵:

r(ABC)+r(B)=r(ABCOOB)=r(OABBCB)=r(BCBOAB)=r(BCBOAB)r(BC)+r(AB)\begin{aligned} r(ABC)+r(B) &= r\begin{pmatrix} ABC&O\\ O&B \end{pmatrix}\\ &= r\begin{pmatrix} O&-AB\\ BC&B \end{pmatrix}\\ &= r\begin{pmatrix} BC&B\\ O&-AB \end{pmatrix}\\ &= r\begin{pmatrix} BC&B\\ O&AB \end{pmatrix}\\ &\ge r(BC)+r(AB) \end{aligned}

其中等号来自分块初等行列变换和交换分块行,秩不变;最后一步使用分块上三角矩阵秩不等式。

于是

r(ABC)+r(B)r(BC)+r(AB)r(ABC)+r(B)\ge r(BC)+r(AB)

r(ABC)r(AB)+r(BC)r(B)r(ABC)\ge r(AB)+r(BC)-r(B)

分块矩阵中的秩#

横向拼接时:

max{r(A),r(B)}r(AB)r(A)+r(B)\max\{r(A),r(B)\}\le r\begin{pmatrix}A&B\end{pmatrix}\le r(A)+r(B)

证明(AB)\begin{pmatrix}A&B\end{pmatrix} 的列向量组包含 AA 的全部列向量,也包含 BB 的全部列向量,所以它的秩至少不小于 r(A),r(B)r(A),r(B) 中较大的一个。

另一方面,取 AA 的一个极大线性无关组和 BB 的一个极大线性无关组,合在一起最多有 r(A)+r(B)r(A)+r(B) 个向量。A,BA,B 的全部列向量都能由这组向量线性表示,所以拼接矩阵的列向量组也能由它线性表示。因此

r(AB)r(A)+r(B)r\begin{pmatrix}A&B\end{pmatrix}\le r(A)+r(B)

纵向拼接时:

max{r(A),r(B)}r(AB)r(A)+r(B)\max\{r(A),r(B)\}\le r\begin{pmatrix}A\\B\end{pmatrix}\le r(A)+r(B)

证明:转置后变成横向拼接:

r(AB)=r(ATBT)r\begin{pmatrix}A\\B\end{pmatrix} =r\begin{pmatrix}A^T&B^T\end{pmatrix}

再用横向拼接结论。

分块对角矩阵:

r(A00B)=r(A)+r(B)r\begin{pmatrix} A&0\\ 0&B \end{pmatrix} =r(A)+r(B)

因为两个块分别作用在互不相交的行列位置上,互不产生线性组合关系。

常用推论#

P,QP,Q 可逆,则

r(PAQ)=r(A)r(PAQ)=r(A)

因为

r(PAQ)r(A)r(PAQ)\le r(A)

A=P1(PAQ)Q1A=P^{-1}(PAQ)Q^{-1}

所以

r(A)r(PAQ)r(A)\le r(PAQ)

两边合起来得到等号。

AA 可逆,则

r(AB)=r(B),r(BA)=r(B)r(AB)=r(B),\qquad r(BA)=r(B)

本质上是左乘或右乘可逆矩阵不改变秩。

文章分享

如果这篇文章对你有帮助,欢迎分享给更多人!

15-矩阵
https://skaco2.com/posts/02-math/15-矩阵/
作者
SKACO2
发布于
2026-05-19
许可协议
CC BY-NC-SA 4.0

评论区

Profile Image of the Author
SKACO2
笼中鸟,何时飞!
公告
欢迎来到我的博客!
音乐
封面

音乐

暂未播放

0:00 0:00
暂无歌词
分类
标签
站点统计
文章
61
分类
10
标签
59
总字数
64,400
运行时长
0
最后活动
0 天前

目录