这部分内容探讨零空间、列空间和行空间及其与线性变换的关系。零空间帮助识别线性方程的解,列空间反映了矩阵的映射范围,而行空间揭示了方程间的关系。
在线性代数的应用中, R n ~\mathbb{R^n}~ R n 的子空间产生的方式通常有两种:
( 1 ) (1)~ ( 1 ) 方程解的集合 :作为齐次线性方程组的所有解的集合,这种子空间通常称为 零空间 。
( 2 ) (2)~ ( 2 ) 线性组合的集合 :作为一组特定向量的所有线性组合的集合,这种子空间称为 列空间 或 行空间 。
接下来我们将介绍并使用这些新的术语重新讨论之前学过的内容,以更清晰地理解子空间的性质。
1. 零空间 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A
1.1 零空间的定义
在线性代数中,零空间是我们重新理解齐次线性方程组解的结构和性质的一个核心概念。我们来看下面这个齐次方程组:
x 1 − 3 x 2 − 2 x 3 = 0 − 5 x 1 + 9 x 2 + x 3 = 0 \begin{align*}
x_1 - 3x_2 - 2x_3 &= 0 \\
-5x_1 + 9x_2 + x_3 &= 0
\end{align*} x 1 − 3 x 2 − 2 x 3 − 5 x 1 + 9 x 2 + x 3 = 0 = 0
它可以用矩阵方程 A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 来表示,其中
A = [ 1 − 3 − 2 − 5 9 1 ] \mathbf{A} = \begin{bmatrix} 1 & -3 & -2 \\ -5 & 9 & 1 \end{bmatrix} A = [ 1 − 5 − 3 9 − 2 1 ]
我们称满足方程 A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 的所有 x ~\mathbf{x}~ x 的集合(解集)为矩阵 A ~\mathbf{A}~ A 的零空间。严格的定义如下:
定义
零空间
一个 m × n ~m\times n~ m × n 矩阵 A ~\mathbf{A}~ A 的零空间 ( n u l l s p a c e ) ~(\mathbf{null ~ space})~ ( null space ) ,记作 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A ,是齐次方程 A x = 0 ~\mathbf{Ax}=\mathbf{0}~ Ax = 0 的所有解的集合。用集合表示法,可以写作:
Nul A = { x ∣ x ∈ R n 且 A x = 0 } \text{Nul} \, \mathbf{A} = \{ \mathbf{x} \mid \mathbf{x} \in \mathbb{R}^n \text{ 且 } \mathbf{A} \mathbf{x} = \mathbf{0} \} Nul A = { x ∣ x ∈ R n 且 Ax = 0 } 1.2 判定向量是否属于零空间
我们判断一个向量是否属于矩阵 A ~\mathbf{A}~ A 的零空间,只需要判断该向量是否为齐次方程 A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 的解即可。例如,判断向量 u ~\textcolor{#ff8200}{\mathbf{u}}~ u 是否存在于矩阵 A ~\textcolor{#2196f3}{\mathbf{A}}~ A 的零空间,只需要判断 A u = 0 ~\mathbf{A}\mathbf{u}=\mathbf{0}~ Au = 0 是否成立即可。经计算:
A u = [ 1 − 3 − 2 − 5 9 1 ] [ 5 3 − 2 ] = [ 5 − 9 + 4 − 25 + 27 − 2 ] = [ 0 0 ] \mathbf{A} \mathbf{u} = \textcolor{#2196f3}{\begin{bmatrix} 1 & -3 & -2 \\ -5 & 9 & 1 \end{bmatrix}} \textcolor{#ff8200}{\begin{bmatrix} 5 \\ 3 \\ -2 \end{bmatrix}} = \begin{bmatrix} 5 - 9 + 4 \\ -25 + 27 - 2 \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \end{bmatrix} Au = [ 1 − 5 − 3 9 − 2 1 ] 5 3 − 2 = [ 5 − 9 + 4 − 25 + 27 − 2 ] = [ 0 0 ]
所以, u ~\mathbf{u}~ u 存在于矩阵 A ~\mathbf{A}~ A 的零空间内,即: u ∈ Null A ~\mathbf{u}\in \text{Null} \, \mathbf{A}~ u ∈ Null A 。我们可以对矩阵 A ~\mathbf{A}~ A 进行行化简并得到齐次方程组的解集 x = t v , 其中, v = [ 1 2 − 3 2 1 ] T \mathbf{x}=t\mathbf{v},~\text{其中,}~\mathbf{v}=\begin{bmatrix} \frac{1}{2} & -\frac{3}{2} & 1 \end{bmatrix}^T x = t v , 其中 , v = [ 2 1 − 2 3 1 ] T ,所以零空间是向量 v ~\mathbf{v}~ v 的所有线性组合,可以用 Span { v } ~\text{Span}\{\mathbf{v}\}~ Span { v } 表示,几何上可以表示为 R 3 ~\mathbb{R^3}~ R 3 中一条经过原点的直线。
1.3 零空间性质
零空间是一个向量的集合,集合中的向量必要符合矩阵乘法 A x ~\mathbf{A}\mathbf{x}~ Ax 的定义,即: A m × n x n × 1 = 0 ~\mathbf{A}_{m\times n}\mathbf{x}_{n\times 1}=\mathbf{0}~ A m × n x n × 1 = 0 。有如下定理:
定理 2
零空间的性质
一个 m × n ~m\times n~ m × n 矩阵 A ~\mathbf{A}~ A 的零空间是 R n ~\mathbb{R^n}~ R n 的一个子空间。等价地, m ~m~ m 个方程、 n ~n~ n 个未知数的齐次线性方程组 A x = 0 ~\mathbf{A} \mathbf{x} = \mathbf{0}~ Ax = 0 的所有解的集合是 R n ~\mathbb{R^n}~ R n 的一个子空间。
这个定理在验证特定条件下的向量集合是否为子空间时非常有用。在实际问题中,我们通常会得到来自系统(工程、物理和计算机科学等领域)中不同变量之间的约束或平衡关系构成的条件方程,例如下面的条件方程: a − 2 b + 5 c = d ~a-2b+5c=d~ a − 2 b + 5 c = d 和 c − a = b ~c-a=b~ c − a = b 。满足此方程的所有向量 ( a , b , c , d ) ~(a,b,c,d)~ ( a , b , c , d ) 构成的集合 H ~H~ H 是 R 4 ~\mathbb{R^4}~ R 4 的一个子空间。这是因为上面的条件方程可以写成如下的齐次方程组:
a − 2 b + 5 c − d = 0 − a − b + c = 0 \begin{align*}
a - 2b + 5c - d &= 0 \\
-a - b + c \quad \quad &= 0
\end{align*} a − 2 b + 5 c − d − a − b + c = 0 = 0
这是一个包含 2 ~2~ 2 个方程、 4 ~4~ 4 个未知数的齐次方程组,根据定理 2 ~2~ 2 可知它的解集 H ~H~ H 是 R 4 ~\mathbb{R^4}~ R 4 的子空间。
1.4 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 的显示表示
上面对零空间的定义并不是显而易见的,在 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 中的向量与矩阵 A ~\mathbf{A}~ A 中的元素没有明显的直接关系,我们称 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 是被 隐式定义 的。隐式定义虽然简洁,但缺乏对集合结构的直观描述,难以直接进行计算和分析。在实际应用中通常会将隐式定义转换为显式表示。为此,我们可以通过求解矩阵方程 A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 的解集来获得显式表示。请观察下面求解 A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 的过程:
将通解中的自由变量提取,并表示为一组向量的线性组合,这些向量 { u , v , w } ~\{\mathbf{u},\mathbf{v},\mathbf{w}\}~ { u , v , w } 就构成了 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 的一个生成集合。上面的示例展示了如何构建 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 的生成集合,并可以得出以下两个关键性质:
( 1 ) ~(1)~ ( 1 ) 线性无关 :生成集合中的向量是线性无关的,因为自由变量作为生成向量的权重,只有当所有权重(即自由变量)均为零时,整个组合结果才为零。
( 2 ) ~(2)~ ( 2 ) 生成集合的维数 :当 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 包含非零向量时,生成集合的向量个数等于自由变量的数量,即方程 A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 中的自由变量的数量。
2. 列空间 Col A ~\text{Col}\,\mathbf{A}~ Col A
列空间 ( C o l u m n S p a c e ) ~(\mathbf{Column ~ Space})~ ( Column Space ) 提供了矩阵 A ~\mathbf{A}~ A 的列向量所能覆盖的所有向量集合的范围,常用于研究矩阵的映射范围、线性独立性以及其他与线性变换相关的性质。与零空间不同,列空间是通过线性组合显式定义:
定义
列空间
列空间 Col A ~\text{Col}\,\mathbf{A}~ Col A 是一个 m × n ~m\times n~ m × n 矩阵 A ~\mathbf{A}~ A 的列向量的线性组合的集合。具体来说,如果矩阵 A ~\mathbf{A}~ A 的列向量是 a 1 , a 2 , … , a n ~\mathbf{a}_1,\mathbf{a}_2,\dots,\mathbf{a}_n~ a 1 , a 2 , … , a n ,那么列空间定义为:
Col A = Span { a 1 , … , a n } \text{Col}\,\mathbf{A}=\text{Span}\,\{\mathbf{a}_1,\dots,\mathbf{a}_n\} Col A = Span { a 1 , … , a n }
也就是说,列空间是由矩阵 A ~\mathbf{A}~ A 的列向量所张成的所有线性组合的集合。
因为 Span { a 1 , … , a n } \text{Span}\,\{\mathbf{a}_1,\dots,\mathbf{a}_n\} Span { a 1 , … , a n } 是一个子空间,根据上一节的 定理 1 可以得出列空间是定义在 R m ~\mathbb{R^m}~ R m 中的子空间,因为 A ~\mathbf{A}~ A 的列向量位于 R m ~\mathbb{R^m}~ R m 中。
定理 3
矩阵的列空间
一个 m × n ~m\times n~ m × n 矩阵 A ~\mathbf{A}~ A 的列空间是 R m ~\mathbb{R^m}~ R m 的一个子空间。
由于 A x ~\mathbf{A}\mathbf{x}~ Ax 表示矩阵 A ~\mathbf{A}~ A 的列向量的线性组合,所以在列空间 Col A ~\text{Col}\,\mathbf{A}~ Col A 中的向量可以表示为 A x ~\mathbf{A}\mathbf{x}~ Ax 的形式,其中 x ~\mathbf{x}~ x 是某个来自 R n ~\mathbb{R^n}~ R n 的向量。
那么,列空间有如下等价定义:
Col A = { b : b = A x , x ∈ R n } \text{Col}\mathbf{A} = \{\mathbf{b}: \mathbf{b} = \mathbf{A}\mathbf{x} ~,~ \mathbf{x} \in \mathbb{R^n}\} Col A = { b : b = Ax , x ∈ R n }
在实际应用中,构造一个矩阵使其列空间等于给定的向量空间的情况是很常见的需求。例如我们前面介绍过的 透视投影矩阵 P ~\mathbf{P}~ P ,它的列空间被定义为所有能够通过矩阵 P ~\mathbf{P}~ P 作用后生成的向量集合,也就是在透视投影下所有可以投影到二维平面上的点。从这个角度来看,矩阵 A ~\mathbf{A}~ A 的列空间 Col A ~\text{Col}\,\mathbf{A}~ Col A 也是线性变换 x ↦ A x ~\mathbf{x}\mapsto \mathbf{A}\mathbf{x}~ x ↦ Ax 的 值域 ( R a n g e ) ~(\mathbf{Range})~ ( Range ) 。
3. 行空间 Row A ~\text{Row}\,\mathbf{A}~ Row A
对于一个 m × n ~m\times n~ m × n 的矩阵 A ~\mathbf{A}~ A ,每一行都可以视为 R n ~\mathbb{R^n}~ R n 中的一个向量。矩阵 A ~\mathbf{A}~ A 的所有行向量的线性组合的集合被称为 A ~\mathbf{A}~ A 的 行空间 ( R o w S p a c e ) ~(\mathbf{Row~Space})~ ( Row Space ) ,记作 Row A ~\text{Row}\,\mathbf{A}~ Row A 。因为每一行都有 n ~n~ n 个元素,行空间 Row A ~\text{Row}\,\mathbf{A}~ Row A 是 R n ~\mathbb{R^n}~ R n 的一个子空间。此外,矩阵 A ~\mathbf{A}~ A 的行向量与转置矩阵 A T ~\mathbf{A}^T~ A T 的列向量相同,因此行空间 Row A ~\text{Row}\,\mathbf{A}~ Row A 可以表示为 Col A T ~\text{Col}\,\mathbf{A}^T~ Col A T 。
4. 对比 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 和 Col A ~\text{Col}\,\mathbf{A}~ Col A
零空间和列空间代表了矩阵的不同几何属性。列空间 Col A ~\text{Col}\,\mathbf{A}~ Col A 包含矩阵 A ~\mathbf{A}~ A 的所有列向量的线性组合,描述了 A ~\mathbf{A}~ A 能够映射到的输出空间;而零空间 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 包含了被 A ~\mathbf{A}~ A 映射到零向量的所有输入向量,即 A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 的所有解。例如下面的矩阵 A ~\mathbf{A}~ A 它的列空间 Col A ~\text{Col}\,\mathbf{A}~ Col A 是 R 3 ~\mathbb{R^3}~ R 3 的一个子空间,而零空间是 R 4 ~\mathbb{R^4}~ R 4 的子空间:
接下来,我们可以从零空间和列空间的结构入手,进一步探讨它们在实际应用中的不同处理方法。在线性代数中, 寻找 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 或 Col A ~\text{Col}\,\mathbf{A}~ Col A 中的非零向量和 判断 一个向量是否属于 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 或 Col A ~\text{Col}\,\mathbf{A}~ Col A 是线性代数中非常常见的两个问题。下面我们分别针对零空间和列空间来讨论这两个问题。
4.1 寻找 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 和 Col A ~\text{Col}\,\mathbf{A}~ Col A 中的非零向量
寻找 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 中的向量,就是求解齐次方程 A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 中的 x ~\mathbf{x}~ x ,这需要对 [ A 0 ] \begin{bmatrix}\mathbf{A} & \mathbf{0}\end{bmatrix} [ A 0 ] 做行变换,以上面的矩阵 A ~\mathbf{A}~ A 为例,行变换过程如下(补充说明:把 A ~\mathbf{A}~ A 化简为 简化阶梯形 总共需要 7 ~7~ 7 个步骤,下面的动画过程对此进行了简化):
上面的示例,零空间 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 中有非零向量,这意味着矩阵 A ~\mathbf{A}~ A 存在线性依赖关系,即方程组有无穷多解。寻找列空间 Col A ~\text{Col}\,\mathbf{A}~ Col A 中的向量比较容易,它就是由矩阵 A ~\mathbf{A}~ A 的列向量线性组合得到的,即 Col A ~\text{Col}\,\mathbf{A}~ Col A 中的向量可以由 b = x 1 a 1 + x 2 a 2 + ⋯ + x n a n \mathbf{b}=x_1\mathbf{a}_1+ x_2\mathbf{a}_2+\dots+x_n\mathbf{a}_n b = x 1 a 1 + x 2 a 2 + ⋯ + x n a n 表示。
4.2 判断 v ~\mathbf{v}~ v 是否存在于 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 和 Col A ~\text{Col}\,\mathbf{A}~ Col A 中
判断 v ∈ Nul A ~\mathbf{v}\in \text{Nul}\,\mathbf{A}~ v ∈ Nul A 很简单。根据对 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 的定义,若 A u = 0 ~\mathbf{A}\mathbf{u}=\mathbf{0}~ Au = 0 成立,那么 v ∈ Nul A ~\mathbf{v}\in \text{Nul}\,\mathbf{A}~ v ∈ Nul A 。判断 v ∈ Col A ~\mathbf{v}\in \text{Col}\,\mathbf{A}~ v ∈ Col A 是否成立,就要尝试能否找到 A ~\mathbf{A}~ A 中列向量的一种组合方式,使得 v = x 1 a 1 + x 2 a 2 + ⋯ + x n a n ~\mathbf{v}=x_1\mathbf{a}_1+x_2\mathbf{a}_2+\dots+x_n\mathbf{a}_n~ v = x 1 a 1 + x 2 a 2 + ⋯ + x n a n 成立,这其实是在判断矩阵方程 A x = v ~\mathbf{A}\mathbf{x}=\mathbf{v}~ Ax = v 中是否有解,这个过程也可以通过对 [ A v ] \begin{bmatrix}\mathbf{A} & \mathbf{v} \end{bmatrix} [ A v ] 进行行化简确定。
上面的增广矩阵的最右列不是主元列,所以根据 线性方程组的存在性定理 可知,方程组 A x = v ~\mathbf{A}\mathbf{x}=\mathbf{v}~ Ax = v 有解,那么 v ∈ Col A ~\mathbf{v}\in \text{Col}\,\mathbf{A}~ v ∈ Col A 成立。
4.3 整体对比 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A 和 Col A ~\text{Col}\,\mathbf{A}~ Col A 表 4-1 m × n ~\quad m\times n~ m × n 矩阵 A ~\mathbf{A}~ A 的Nul A \text{Nul}\,\mathbf{A}~ Nul A 和Col A \text{Col}\,\mathbf{A}~ Col A 的对比
Nul A \text{Nul}\,\mathbf{A}~ Nul A
Col A \text{Col}\,\mathbf{A}~ Col A
1 ~1~ 1
Nul A \text{Nul}\,\mathbf{A}~ Nul A 是 R n ~\mathbb{R^n}~ R n 的子空间。
Col A \text{Col}\,\mathbf{A}~ Col A 是 R m ~\mathbb{R^m}~ R m 的子空间。
2 ~2~ 2
Nul A \text{Nul}\,\mathbf{A}~ Nul A 是
隐式定义 的。
Col A \text{Col}\,\mathbf{A}~ Col A 是
显式定义 的。
3 ~3~ 3
求Nul A \text{Nul}\,\mathbf{A}~ Nul A 中的向量,需要对[ A 0 ] \begin{bmatrix}\mathbf{A} & \mathbf{0}\end{bmatrix} [ A 0 ] 做行变换。
Col A \text{Col}\,\mathbf{A}~ Col A 中的向量可由 A ~\mathbf{A}~ A 的列组合出来。
4 ~4~ 4
Nul A \text{Nul}\,\mathbf{A}~ Nul A 与 A ~\mathbf{A}~ A 中元素之间没有明显关系。
A ~\mathbf{A}~ A 中的列是Col A \text{Col}\,\mathbf{A}~ Col A 中的向量。
5 ~5~ 5
Nul A \text{Nul}\,\mathbf{A}~ Nul A 中的任意向量 v ~\mathbf{v}~ v ,满足 A v = 0 ~\mathbf{A}\mathbf{v}=\mathbf{0}~ Av = 0 。
Col A \text{Col}\,\mathbf{A}~ Col A 中的向量
v ~\mathbf{v}~ v 满足方程
A x = v ~\mathbf{A}\mathbf{x}=\mathbf{v}~ Ax = v 相容 。
6 ~6~ 6
判断 v ∈ Nul A ~\mathbf{v}\in \text{Nul}\,\mathbf{A}~ v ∈ Nul A 是否成立,只需计算 A v ~\mathbf{A}\mathbf{v}~ Av 。
判断 v ∈ Col A ~\mathbf{v}\in \text{Col}\,\mathbf{A}~ v ∈ Col A 是否成立,需对[ A v ] \begin{bmatrix}\mathbf{A} & \mathbf{v}\end{bmatrix} [ A v ] 做行变换。
7 ~7~ 7
Nul A = { 0 } \text{Nul}\,\mathbf{A}=\{\mathbf{0}\} Nul A = { 0 } 当且仅当
A x = 0 ~\mathbf{A}\mathbf{x}=\mathbf{0}~ Ax = 0 仅有一个
平凡解 。
Col A = R m \text{Col}\mathbf{A}=\mathbb{R^m} Col A = R m 当且仅当A x = b \mathbf{A}\mathbf{x}=\mathbf{b} Ax = b 对任意b ∈ R m \mathbf{b}\in \mathbb{R^m} b ∈ R m 有解。
8 ~8~ 8
Nul A = R n \text{Nul}\mathbf{A}=\mathbb{R^n} Nul A = R n 当且仅当x ↦ A x \mathbf{x}\mapsto \mathbf{Ax} x ↦ Ax 是一对一的。
Col A = R m \text{Col}\mathbf{A}=\mathbb{R^m} Col A = R m 当且仅当x ↦ A x \mathbf{x}\mapsto \mathbf{Ax} x ↦ Ax 将 R n ~\mathbb{R^n}~ R n 映射到 R m ~\mathbb{R^m}~ R m 上。
5. 线性变换
前面我们学习了零空间和列空间,它们都是通过矩阵定义的子空间。现在我们将视角从矩阵提升到更一般的 线性变换 ,这是一个能够统一描述各种"保持线性结构"的映射的核心概念。
5.1 线性变换的定义 定义
线性变换的定义
从向量空间 V ~V~ V 到向量空间 W ~W~ W 的线性变换 T ~T~ T 是这样一种规则:它将 V ~V~ V 中每个向量 x ~\mathbf{x}~ x 映射为 W ~W~ W 中唯一的向量 T ( x ) ~T(\mathbf{x})~ T ( x ) ,并满足以下条件:
i. 对于 V ~V~ V 中的任意向量 u ~\mathbf{u}~ u 和 v ~\mathbf{v}~ v ,有 T ( u + v ) = T ( u ) + T ( v ) ~T(\mathbf{u}+\mathbf{v})=T(\mathbf{u}) + T(\mathbf{v})~ T ( u + v ) = T ( u ) + T ( v ) 。
ii. 对于 V ~V~ V 中的任意向量 u ~\mathbf{u}~ u 和任意标量 c ~c~ c ,有 T ( c u ) = c T ( u ) ~T(c\mathbf{u})=cT(\mathbf{u})~ T ( c u ) = c T ( u ) 。
线性变换 T ~T~ T 有两个重要的关联子空间: 核 ( kernel ) ~(~\text{kernel}~)~ ( kernel ) ,也称为零空间,是 V ~V~ V 中所有满足 T ( u ) = 0 ~T(\mathbf{u})=\mathbf{0}~ T ( u ) = 0 的向量集合; 值域 ( range ) ~(~\text{range}~)~ ( range ) ,是 W ~W~ W 中所有形如 T ( x ) ~T(\mathbf{x})~ T ( x ) 的向量集合。当线性变换是矩阵变换 T ( x ) = A x ~T(\mathbf{x})=\mathbf{A}\mathbf{x}~ T ( x ) = Ax 时,核就是矩阵 A ~\mathbf{A}~ A 的零空间 Nul A ~\text{Nul}\,\mathbf{A}~ Nul A ,值域就是列空间 Col A ~\text{Col}\,\mathbf{A}~ Col A 。
5.2 线性变换与求导运算
线性变换的概念不仅适用于矩阵,还能推广到函数空间。在微积分中, 求导运算 本身就是一个线性变换。这个视角将微分方程的求解问题纳入了线性代数的统一框架。
例 9 (需要微积分基础):设 V ~V~ V 是区间 [ a , b ] ~[a,b]~ [ a , b ] 上所有可微函数的集合, W ~W~ W 是 [ a , b ] ~[a,b]~ [ a , b ] 上所有连续函数的集合。定义变换 D : V → W ~D:V\to W~ D : V → W 为求导运算,即 D ( f ) = f ′ ~D(f)=f'~ D ( f ) = f ′ 。证明 D ~D~ D 是线性变换,并确定其核与值域。
明确向量空间
设 V ~V~ V 是区间 [ a , b ] ~[a,b]~ [ a , b ] 上所有可微函数的集合(且导数连续), W ~W~ W 是 [ a , b ] ~[a,b]~ [ a , b ] 上所有连续函数的集合 C [ a , b ] ~C[a,b]~ C [ a , b ] 。定义变换 D : V → W ~D:V\to W~ D : V → W 为求导运算: D ( f ) = f ′ ~D(f)=f'~ D ( f ) = f ′ 。
验证加法保持性
对于 V ~V~ V 中任意两个函数 f ~f~ f 和 g ~g~ g ,根据求导的和法则:
D ( f + g ) = ( f + g ) ′ = f ′ + g ′ = D ( f ) + D ( g ) D(f+g) = (f+g)' = f' + g' = D(f) + D(g) D ( f + g ) = ( f + g ) ′ = f ′ + g ′ = D ( f ) + D ( g )
这说明 D ~D~ D 保持向量加法运算。
验证标量乘法保持性
对于 V ~V~ V 中任意函数 f ~f~ f 和任意标量 c ~c~ c ,根据求导的常数倍法则:
D ( c f ) = ( c f ) ′ = c ⋅ f ′ = c D ( f ) D(cf) = (cf)' = c \cdot f' = cD(f) D ( c f ) = ( c f ) ′ = c ⋅ f ′ = cD ( f )
这说明 D ~D~ D 保持标量乘法运算。
确定核(零空间)
D ~D~ D 的核是所有满足 D ( f ) = 0 ~D(f)=0~ D ( f ) = 0 的函数,即 f ′ = 0 ~f'=0~ f ′ = 0 。根据微积分基本定理,导数恒为零的函数必为常数函数。因此:
ker ( D ) = { f ∈ V ∣ f ( x ) = c , c ∈ R } \ker(D) = \{ f \in V \mid f(x) = c, \, c \in \mathbb{R} \} ker ( D ) = { f ∈ V ∣ f ( x ) = c , c ∈ R }
核是 [ a , b ] ~[a,b]~ [ a , b ] 上所有常数函数的集合。
确定值域
D ~D~ D 的值域是所有形如 D ( f ) = f ′ ~D(f)=f'~ D ( f ) = f ′ 的函数集合。由于 V ~V~ V 中的函数可微且导数连续,其导数 f ′ ~f'~ f ′ 是 [ a , b ] ~[a,b]~ [ a , b ] 上的连续函数。反过来,对于 W ~W~ W 中任意连续函数 g ~g~ g ,根据微积分基本定理, F ( x ) = ∫ a x g ( t ) d t F(x)=\int_a^x g(t)\,dt F ( x ) = ∫ a x g ( t ) d t 是 g ~g~ g 的一个原函数,且 F ′ = g F'=g F ′ = g 。因此:
range ( D ) = W = C [ a , b ] \text{range}(D) = W = C[a,b] range ( D ) = W = C [ a , b ]
值域是 [ a , b ] ~[a,b]~ [ a , b ] 上所有连续函数的集合。
例 10 (需要微积分基础):微分方程 y ′ ′ + ω 2 y = 0 y'' + \omega^2 y = 0 y ′′ + ω 2 y = 0 描述了弹簧振动、钟摆运动等物理系统。证明该微分方程的解集是某个线性变换的核。
理解微分方程的物理背景
微分方程 y ′ ′ + ω 2 y = 0 ~y''+\omega^2 y=0~ y ′′ + ω 2 y = 0 描述了多种物理系统的振动行为:弹簧的简谐振动、钟摆的小角度摆动、LC电路中的电压振荡等。其中 ω ~\omega~ ω 是角频率常数。
定义线性变换
设 V ~V~ V 是所有二阶可微函数的集合, W ~W~ W 是所有连续函数的集合。定义变换 T : V → W ~T:V\to W~ T : V → W 为:
T ( f ) = f ′ ′ + ω 2 f T(f) = f'' + \omega^2 f T ( f ) = f ′′ + ω 2 f
这个变换将函数 f ~f~ f 映射到 f ′ ′ + ω 2 f ~f''+\omega^2 f~ f ′′ + ω 2 f 。
验证 T 是线性变换
对于任意函数 f , g ∈ V ~f,g\in V~ f , g ∈ V 和标量 c ~c~ c :
T ( f + g ) = ( f + g ) ′ ′ + ω 2 ( f + g ) = f ′ ′ + g ′ ′ + ω 2 f + ω 2 g = ( f ′ ′ + ω 2 f ) + ( g ′ ′ + ω 2 g ) = T ( f ) + T ( g ) \begin{align*} T(f+g) &= (f+g)'' + \omega^2(f+g) \\[1ex] &= f'' + g'' + \omega^2 f + \omega^2 g \\[1ex] &= (f'' + \omega^2 f) + (g'' + \omega^2 g) \\[1ex] &= T(f) + T(g) \end{align*} T ( f + g ) = ( f + g ) ′′ + ω 2 ( f + g ) = f ′′ + g ′′ + ω 2 f + ω 2 g = ( f ′′ + ω 2 f ) + ( g ′′ + ω 2 g ) = T ( f ) + T ( g )
T ( c f ) = ( c f ) ′ ′ + ω 2 ( c f ) = c f ′ ′ + c ω 2 f = c ( f ′ ′ + ω 2 f ) = c T ( f ) \begin{align*} T(cf) &= (cf)'' + \omega^2(cf) \\[1ex] &= cf'' + c\omega^2 f \\[1ex] &= c(f'' + \omega^2 f) = cT(f) \end{align*} T ( c f ) = ( c f ) ′′ + ω 2 ( c f ) = c f ′′ + c ω 2 f = c ( f ′′ + ω 2 f ) = c T ( f )
因此 T ~T~ T 满足线性变换的两个条件。
微分方程的解是 T 的核
微分方程 y ′ ′ + ω 2 y = 0 ~y''+\omega^2 y=0~ y ′′ + ω 2 y = 0 的解集正是 T ~T~ T 的核:
ker ( T ) = { f ∈ V ∣ T ( f ) = 0 } = { f ∣ f ′ ′ + ω 2 f = 0 } \ker(T) = \{ f \in V \mid T(f) = 0 \} = \{ f \mid f'' + \omega^2 f = 0 \} ker ( T ) = { f ∈ V ∣ T ( f ) = 0 } = { f ∣ f ′′ + ω 2 f = 0 }
这意味着 求解微分方程等价于找出线性变换的核 。
求解核(微分方程的通解)
该微分方程的特征方程为 r 2 + ω 2 = 0 ~r^2+\omega^2=0~ r 2 + ω 2 = 0 ,解得 r = ± i ω ~r=\pm i\omega~ r = ± iω 。因此通解为:
y = c 1 cos ( ω t ) + c 2 sin ( ω t ) y = c_1 \cos(\omega t) + c_2 \sin(\omega t) y = c 1 cos ( ω t ) + c 2 sin ( ω t )
其中 c 1 , c 2 ~c_1,c_2~ c 1 , c 2 是任意常数。这说明 T ~T~ T 的核是由 cos ( ω t ) \cos(\omega t) cos ( ω t ) 和 sin ( ω t ) \sin(\omega t) sin ( ω t ) 张成的二维向量空间。
5.3 线性变换与移动平均
线性变换在金融领域也有重要应用。在股票技术分析中, 移动平均 ( moving average ) ~(~\text{moving average}~)~ ( moving average ) 是一种常用的指标,它通过平滑价格波动来揭示整体趋势。移动平均本质上也是一个线性变换。
例 11 :设 { p k } ~\{p_k\}~ { p k } 表示某只股票在第 k ~k~ k 天的价格序列。定义两日移动平均变换 M 2 ~M_2~ M 2 为:
M 2 ( { p k } ) = { p k + p k − 1 2 } M_2(\{p_k\}) = \left\{ \frac{p_k + p_{k-1}}{2} \right\} M 2 ({ p k }) = { 2 p k + p k − 1 }
证明 M 2 ~M_2~ M 2 是线性变换,并求其核。
理解信号空间
设 S ~\mathbb{S}~ S 是所有离散信号 { p k } ~\{p_k\}~ { p k } 的集合,其中 p k ~p_k~ p k 表示第 k ~k~ k 天的股票价格。假设在研究时间段之外 p k = 0 ~p_k=0~ p k = 0 。两日移动平均变换 M 2 : S → S ~M_2:\mathbb{S}\to\mathbb{S}~ M 2 : S → S 定义为:
M 2 ( { p k } ) = { p k + p k − 1 2 } M_2(\{p_k\}) = \left\{ \frac{p_k + p_{k-1}}{2} \right\} M 2 ({ p k }) = { 2 p k + p k − 1 }
验证加法保持性
对于 S ~\mathbb{S}~ S 中任意两个信号 { p k } ~\{p_k\}~ { p k } 和 { q k } ~\{q_k\}~ { q k } :
M 2 ( { p k } + { q k } ) = M 2 ( { p k + q k } ) = { ( p k + q k ) + ( p k − 1 + q k − 1 ) 2 } = { p k + p k − 1 2 } + { q k + q k − 1 2 } = M 2 ( { p k } ) + M 2 ( { q k } ) \begin{align*} M_2(\{p_k\} + \{q_k\}) &= M_2(\{p_k + q_k\}) \\[1ex] &= \left\{ \frac{(p_k + q_k) + (p_{k-1} + q_{k-1})}{2} \right\} \\[1ex] &= \left\{ \frac{p_k + p_{k-1}}{2} \right\} + \left\{ \frac{q_k + q_{k-1}}{2} \right\} \\[1ex] &= M_2(\{p_k\}) + M_2(\{q_k\}) \end{align*} M 2 ({ p k } + { q k }) = M 2 ({ p k + q k }) = { 2 ( p k + q k ) + ( p k − 1 + q k − 1 ) } = { 2 p k + p k − 1 } + { 2 q k + q k − 1 } = M 2 ({ p k }) + M 2 ({ q k })
验证标量乘法保持性
对于任意信号 { p k } ~\{p_k\}~ { p k } 和标量 c ~c~ c :
M 2 ( c { p k } ) = M 2 ( { c p k } ) = { c p k + c p k − 1 2 } = c { p k + p k − 1 2 } = c M 2 ( { p k } ) \begin{align*} M_2(c\{p_k\}) &= M_2(\{cp_k\}) \\[1ex] &= \left\{ \frac{cp_k + cp_{k-1}}{2} \right\} \\[1ex] &= c \left\{ \frac{p_k + p_{k-1}}{2} \right\} \\[1ex] &= cM_2(\{p_k\}) \end{align*} M 2 ( c { p k }) = M 2 ({ c p k }) = { 2 c p k + c p k − 1 } = c { 2 p k + p k − 1 } = c M 2 ({ p k })
因此 M 2 ~M_2~ M 2 是线性变换。
求核的条件
信号 { p k } ~\{p_k\}~ { p k } 在 M 2 ~M_2~ M 2 的核中,当且仅当对所有 k ~k~ k :
p k + p k − 1 2 = 0 ⇒ p k = − p k − 1 \frac{p_k + p_{k-1}}{2} = 0 \quad \Rightarrow \quad p_k = -p_{k-1} 2 p k + p k − 1 = 0 ⇒ p k = − p k − 1
这意味着相邻两项互为相反数。
推导核的形式
由 p k = − p k − 1 ~p_k=-p_{k-1}~ p k = − p k − 1 递推可得:
p k = − p k − 1 = ( − 1 ) 2 p k − 2 = ( − 1 ) 3 p k − 3 = ⋯ = p 0 ⋅ ( − 1 ) k p_k = -p_{k-1} = (-1)^2 p_{k-2} = (-1)^3 p_{k-3} = \cdots = p_0 \cdot (-1)^k p k = − p k − 1 = ( − 1 ) 2 p k − 2 = ( − 1 ) 3 p k − 3 = ⋯ = p 0 ⋅ ( − 1 ) k
因此核中的任意信号都可以写成:
{ p k } = p 0 ⋅ { ( − 1 ) k } \{p_k\} = p_0 \cdot \{(-1)^k\} { p k } = p 0 ⋅ {( − 1 ) k }
即 交替信号 { ( − 1 ) k } = { … , 1 , − 1 , 1 , − 1 , … } ~\{(-1)^k\}=\{\dots,1,-1,1,-1,\dots\}~ {( − 1 ) k } = { … , 1 , − 1 , 1 , − 1 , … } 的标量倍。
核的物理意义
核由所有交替信号的倍数组成。这说明两日移动平均能够 完全消除 日间的交替波动(如每日涨跌交替的噪声),同时保留整体趋势。这正是移动平均在技术分析中的核心作用:平滑短期波动,揭示长期趋势。