16-矩阵

1656 字
8 分钟
16-矩阵

公式#

操作顺序可交换#

(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,所以行列式不变。

常见用途:

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

#

子式#

在矩阵里任选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)

证明

BBm×nm\times n 矩阵,且

r(B)=sr(B)=s

由初等行变换、初等列变换可知,存在可逆矩阵 P,QP,Q,使得

PBQ=(EsOOO)PBQ= \begin{pmatrix} E_s&O\\ O&O \end{pmatrix}

记这个矩阵为 JJ

A=AP1,C=Q1CA'=AP^{-1},\qquad C'=Q^{-1}C

可逆矩阵相乘不改变秩,所以

r(ABC)=r(AJC),r(AB)=r(AJ),r(BC)=r(JC)r(ABC)=r(A'JC'),\qquad r(AB)=r(A'J),\qquad r(BC)=r(JC')

AA' 按列分块,把 CC' 按行分块:

A=(A1A2),C=(C1C2)A'=\begin{pmatrix}A_1&A_2\end{pmatrix}, \qquad C'=\begin{pmatrix}C_1\\C_2\end{pmatrix}

其中 A1A_1 取前 ss 列,C1C_1 取前 ss 行。于是

AJ=(A1O),JC=(C1O),AJC=A1C1A'J=\begin{pmatrix}A_1&O\end{pmatrix},\qquad JC'=\begin{pmatrix}C_1\\O\end{pmatrix},\qquad A'JC'=A_1C_1

所以

r(AB)=r(A1),r(BC)=r(C1),r(ABC)=r(A1C1)r(AB)=r(A_1),\qquad r(BC)=r(C_1),\qquad r(ABC)=r(A_1C_1)

A1C1A_1C_1 使用 Sylvester 不等式,其中中间矩阵的列数是 ss,得到

r(A1C1)r(A1)+r(C1)sr(A_1C_1)\ge r(A_1)+r(C_1)-s

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)

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

文章分享

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

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

评论区

Profile Image of the Author
SKACO2
Hello……
公告
欢迎来到我的博客!
音乐
封面

音乐

暂未播放

0:00 0:00
暂无歌词
分类
标签
站点统计
文章
53
分类
8
标签
54
总字数
58,255
运行时长
0
最后活动
0 天前

目录