上一讲中,我们介绍了vector space的定义,我们知道了vector space必须对线性组合(加,数乘)闭合。在这一讲里,Prof. Strang重点介绍了矩阵的两个subspace:column space and null space.
Subspace的union和intersection
在正式介绍之前,我们先看一个subspace的性质。假设 VVV, WWW 分别是 Rn\mathbb{R}^nRn 的subspace,那么
-
V∪WV\cup WV∪W 并不是一个vector space
-
V∩WV\cap WV∩W 是一个vector space
第一点很容易得出,不再赘述。我们着重看第二点。?x,y∈V∩W\forall~x,y\in V\cap W? x,y∈V∩W, 则有
c1x+c2y∈Vc_1 x+c_2 y \in Vc1?x+c2?y∈V
c1x+c2y∈Wc_1 x+c_2 y\in Wc1?x+c2?y∈W
?c1,c2\forall~c_1,c_2? c1?,c2?. 这很显然因为 x,y∈Vx, y \in Vx,y∈V and x,y∈Wx, y \in Wx,y∈W.
As a result,
c1x+c2y∈V∩Wc_1 x+c_2 y \in V\cap Wc1?x+c2?y∈V∩W
V∩WV\cap WV∩W 是一个subspace。
Column space of A: C(A)C(A)C(A)
上一讲中我们已经touch了矩阵的column space,我们已经知道了矩阵的列的线性组合可以张成Rn\mathbb{R}^nRn的一个子空间(nnn是矩阵的行数)。这一节中我们要探讨的事column space到底有什么用。
好,给定一个矩阵
A=[112213314415]A=\begin{bmatrix} 1&1&2 \\ 2&1&3\\ 3&1&4\\ 4&1&5\\ \end{bmatrix}A=?????1234?1111?2345??????
显然,它只有三列,他的column space最多是 R4\mathbb{R}^4R4 中的一个三维空间,而且显然他的列空间只是一个二维空间,因为第三列等于前两列的和。
我们思考这样一个问题:
Ax=b\bm{Ax=b}Ax=b
[112213314415][x1x2x3]=[b1b2b3b4]\begin{bmatrix} 1&1&2 \\ 2&1&3\\ 3&1&4\\ 4&1&5\\ \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3\\ \end{bmatrix}= \begin{bmatrix} b_1 \\ b_2 \\ b_3 \\ b_4 \\ \end{bmatrix} ?????1234?1111?2345?????????x1?x2?x3?????=?????b1?b2?b3?b4???????
随便取一个vector b=[b1,b2,b3,b4]?\bm{b}=[b_1,b_2,b_3,b_4]^\topb=[b1?,b2?,b3?,b4?]? 这个方程他有解吗?直观上来看,我们有4个方程,三个未知数。方程个数多于未知数的时候,方程之间很容易有冲突。因此,大多数的 b\bm{b}b 这个方程是无解的。那,他什么时候有解尼?
其实答案也是呼之欲出的,因为实际上 x\bm{x}x 就是对 AAA 的各列进行线性组合,因此只要b\bm{b}b 在 AAA 的column space中,这个方程一定有解!
Null space of A: N(A)N(A)N(A)
好,我们来看一下另一个space – null space. 不同于column space (defined for b\bm{b}b), null space 是为了 x\bm{x}x 定义的。回到上一个例子,我们思考以下问题:如果我们把 b\bm{b}b 置成0,所有满足条件的 x\bm{x}x 是什么样子的尼?
Ax=0\bm{Ax=0}Ax=0
[112213314415][x1x2x3]=[0000]\begin{bmatrix} 1&1&2 \\ 2&1&3\\ 3&1&4\\ 4&1&5\\ \end{bmatrix} \begin{bmatrix} x_1\\ x_2\\ x_3\\ \end{bmatrix}= \begin{bmatrix} 0 \\ 0 \\ 0 \\ 0 \\ \end{bmatrix} ?????1234?1111?2345?????????x1?x2?x3?????=?????0000??????
我先说答案,所有把 A\bm{A}A 置零的 x\bm{x}x 也构成一个 vector space,被称为 A\bm{A}A 的null space
好,现在开始分析。我们来看看什么样的 x\bm{x}x 满足这个公式。首先我们能找到一个零解对吧,即 x=0\bm{x}=\bm{0}x=0, 其次我们还能找到一个非零解 [1,1,?1]?[1,1,-1]^\top[1,1,?1]?,那么实际上x\bm{x}x的解就是
N(A)=c[11?1]N(A)=c\begin{bmatrix} 1\\ 1\\ -1\\ \end{bmatrix}N(A)=c???11?1????
好们这里有一个问题,为什么A\bm{A}A的null space是一个vector space?其实很好理解,因为如果Ax1=0\bm{Ax_1=0}Ax1?=0 and Ax2=0\bm{Ax_2=0}Ax2?=0, 那么显然 A(c1x1+c2x2)=0\bm{A}(c_1\bm{x_1}+c_2\bm{x_2})=\bm{0}A(c1?x1?+c2?x2?)=0.
注意,AAA 的 null space 的维度是与 AAA 的列的个数一致的。
齐次线性方程组的基础解系
其实我们可以回顾一下国内线性代数学习中讲的齐次线性方程组的基础解系的内容,对照来看。
- 齐次线性方程组的基础解系
Am×nx=0\bm{A}_{m\times n}\bm{x}=\bm{0}Am×n?x=0
-
若 rank(Am×n)=n\text{rank}(\bm{A}_{m\times n})=nrank(Am×n?)=n (列满秩),方程有唯一解 (零解)。
-
若 rank(Am×n)≠n\text{rank}(\bm{A}_{m\times n})\neq nrank(Am×n?)??=n (列不满秩),方程有无穷解。
这两个结论很好理解了,因为当 A\bm{A}A 列满秩时候,它的列线性无关: column space C(A)C(\bm{A})C(A)是 Rn\mathbb{R}^nRn,null space N(A)N(\bm{A})N(A) 只有一个全零解。另一方面,如果A\bm{A}A 列不满秩,它的null space是一个space,当然是无穷解了。