矩阵的四个基本空间, 不了解下吗?

秩-零化度定理告诉我们: m×nm\times n 阶矩阵 AA 的零空间(Nullspace) N(A)N(A) 和列空间(Column sapce) C(A)C(A) 的关系:

n=dimN(A)+dimC(A) n = \dim N(A) + \dim C(A)

本次依据秩-零化度定理, 介绍四个基本空间, 并证明它们的正交性关系, 最后, 给出一道经典例题.

四个基本空间

将矩阵 AA 行列互换, 称为 AA 的转置, 记为 ATA^{\mathsf T}. 则对应的: C(AT)C(A^{\mathsf T}) 表示 AA 的行空间(Row space), 对应地将 N(AT)N(A^{\mathsf T}) 称为 AA 的左零空间(Left nullspace).

首先, 因为行秩 = 列秩, 所以

n=dimN(A)+dimC(AT) n = \dim N(A) + \dim C(A^{\mathsf T})

然后, 转置并参考秩-零化度定理

m=dimN(AT)+dimC(A) m = \dim N(A^{\mathsf T}) + \dim C(A)

综上, 对于 m×nm \times n 阶矩阵 AA 有:

空间符号矩阵含义空间维数 向量维数
行空间(Row Space)C(AT)C(A^{\mathsf T})ATyA^{\mathsf T}yrrnn
列空间(Column Space)C(A)C(A)AxAxrrmm
零空间(Nullspace) 或核空间(Kernel Space)N(A)N(A)Ax=0Ax=0nrn-rnn
左零空间(Left Nullspace)N(AT)N(A^{\mathsf T})ATy=0 A^{\mathsf T}y=0mrm-rmm

正交关系

实际上, 行空间和零空间在 Rn\mathbb{R}^n 中是直和关系, 相应的, 列空间和左零空间在 Rm\mathbb{R}^m 中是直和关系, 如下图所示

四个基本空间正交关系

下面, 证明这两个直和(正交)关系

简要证明

  1. Rn=C(AT)N(A)\mathbb{R}^n = C(A^{\mathsf T}) \oplus N(A)
  • x1C(AT),x2N(A)\forall\, x_1 \in C(A^{\mathsf T}), \forall \, x_2 \in N(A), 成立 x=x1+x2Rnx=x_1+x_2 \in \mathbb{R}^n
  • 假设 xC(AT)N(A)x \in C(A^{\mathsf T}) \cap N(A), 即
{Ax=0xN(A)ATy=xyC(AT) \begin{cases} Ax = 0 & \forall \, x \in N(A) \\\\[3pt] A^{\mathsf T}y = x & \exists\, y \in C(A^{\mathsf T}) \end{cases}

联立即得

AATy=0 AA^{\mathsf T}y = 0

左右同乘 yTy^{\mathsf T}

yTAATy=(ATy)TATy=0 y^{\mathsf T}AA^{\mathsf T}y = (A^{\mathsf T}y)^{\mathsf T}A^{\mathsf T}y=0

得到

x=ATy=0 x = A^{\mathsf T}y = 0

C(AT)N(A)={0} C(A^{\mathsf T}) \cap N(A) = \{0\}

结论得证.


  1. Rm=C(A)N(AT)\mathbb{R}^m = C(A) \oplus N(A^{\mathsf T})

B=ATB=A^{\mathsf T}, 利用证明 1

Rm=C(BT)N(B) \mathbb{R}^m = C(B^{\mathsf T}) \oplus N(B)

BT=(AT)T=AB^{\mathsf T} = (A^{\mathsf T})^{\mathsf T} = A 带入上式, 即证.

经典例题

假设存在 m×nm\times n 阶矩阵 AA, 给定一个向量 bRmb\in\,\mathbb{R}^m, 且已知 Ax=bAx = b 有解. 试证

  1. 存在唯一 y \in\, C(A^{\mathsf T}), 使得 Ay=bAy = b 成立.
  2. 所有解中, 行空间 C(AT)C(A^{\mathsf T}) 中的解 yy 的长度最小.
  3. 假设 AA 行满秩(rank(A)=m{\rm rank}(A) = m), 求 yC(AT)y \in \, C(A^{\mathsf T}) (用 A,bA, b 表示)

证:

1 . 存在唯一 y \in\, C(A^{\mathsf T}), 使得 Ay=bAy = b 成立.

  • 存在性
    假设 xRnx \in \mathbb{R}^n 满足 Ax=bAx=b, 因为 C(AT)C(A^{\mathsf T})N(A)N(A) 是直和关系, 所以 xx唯一地分解为

    x=y+z x = y + z

    其中, yC(AT),zN(A)y\in\, C(A^{\mathsf T}), z \in N(A), 那么

    Ax=A(y+z)=Ay+Az=Ay=b\begin{align*} Ax &= A(y+z) \\\\[3pt]&= Ay + Az \\\\[3pt]&= Ay = b \end{align*}

    成立.

  • 唯一性
    假设 yC(AT)y^{\prime} \in\, C(A^{\mathsf T}), 也满足 Ay=bAy^{\prime} = b. 则

    AyAy=A(yy)=0 Ay - Ay^{\prime} = A(y-y^{\prime}) =0

    yyN(A)y-y^{\prime} \in\, N(A), 又因为 yyC(AT)y-y^{\prime} \in \, C(A^{\mathsf T}), 所以

    yy=0 y-y^{\prime} = 0

    y=yy = y^{\prime}.


2 . 所有解中, 行空间 C(AT)C(A^{\mathsf T}) 中的解 yy 的长度最小.

联系问 1, 解 x=y+zx = y + z, 其中 yC(AT),zN(A)y\in\, C(A^{\mathsf T}), z \in N(A). 因为二者正交, 所以

x2=(y+z)T(y+z)=y2+z2y2\begin{align*} \lVert x \rVert^2 &= (y+z)^{\mathsf T}(y+z) \\\\[3pt]&= \lVert y \rVert^2 + \lVert z \rVert^2 \geq \lVert y \rVert^2\end{align*}

z=0z = 0 时, 等号成立.


3 . 假设 AA 行满秩(rank(A)=m{\rm rank}(A) = m), 求 yC(AT)y \in \, C(A^{\mathsf T}) (用 A,bA, b 表示)

因为 yC(AT)y \in \, C(A^{\mathsf T}), 则 yy 可表示为 y=ATxy=A^{\mathsf T}x.则

Ay=AATx=b Ay = AA^{\mathsf T} x = b

AATAA^{\mathsf T} 满秩, 则

x=(AAT)1b x = (AA^{\mathsf T})^{-1}b

可得到

y=AT(AAT)1by=A^{\mathsf T}(AA^{\mathsf T})^{-1}b

下面证明 rankA=rankAAT{\rm rank} A ={\rm rank } AA^{\mathsf T}, 只需证明

{Ax=0ATAx=0 \begin{cases}Ax =0 \\\\[3pt] A^{\mathsf T}Ax=0 \end{cases}

同解. 显然, Ax=0ATAx=0Ax = 0 \Longrightarrow A^{\mathsf T}Ax=0.

若满足 ATAx=0A^{\mathsf T}Ax=0, 左右同乘 xTx^{\mathsf T}

xTATAx=(Ax)TAx=0 x^{\mathsf T}A^{\mathsf T}Ax=(Ax)^{\mathsf T}Ax=0

即得 Ax=0Ax = 0. 所以 rankA=rankATA=rankAAT{\rm rank} A = {\rm rank } A^{\mathsf T}A={\rm rank } AA^{\mathsf T}.

总结

从空间的角度理解线代才是掌握线代的不二法门.

矩阵的四个基本空间是很重要的概念, 可以帮助我们从空间的角度理解线性方程组解的结构, 甚至是最小二乘法.

由于篇幅原因, 下次介绍最小二乘法(ATAx=ATbA^{\mathsf T}Ax = A^{\mathsf T}b).