Preparations for the Proof of the Existence and Uniqueness Theorem
#Math285
General Summary
- 问题重述与算子不动点: 将一阶常微分方程的初值问题 (Initial Value Problem, IVP) 转化为等价的积分方程形式,并认识到求解 IVP 等价于寻找某个特定算子 (Operator) 的不动点 (Fixed Point)。
-
高阶方程降阶: 将
阶常微分方程转化为一个等价的一阶常微分方程组,从而可以将存在唯一性的讨论统一到一阶系统上。 - (选学) 牛顿迭代法 (Newton’s Method): 作为求解方程根 (可以看作一种不动点问题) 的经典方法,牛顿迭代法的思想和收敛性分析为我们理解后续的不动点定理提供了重要的类比和启发。
- 度量空间 (Metric Spaces): 引入度量空间的概念,这是讨论函数之间“距离”和收敛性的数学框架。
- 巴拿赫不动点定理 (Banach’s Fixed-Point Theorem): 这是证明存在唯一性定理的核心数学工具,它保证了在特定条件下,“压缩映射”拥有唯一的不动点。
- 矩阵范数 (Matrix Norms): 推广了向量的“长度”概念到矩阵,为后续在线性化分析和估计中度量矩阵“大小”提供工具。
问题重述与算子不动点 (Problem Restatement & Fixed Points of Operators)
微分积分转化
-
核心思想: 将微分问题转化为积分问题。
考虑一阶 IVP:$$ y’ = f(t, y), \quad y(t_0) = y_0
$$
其中
$$
1
\phi(t) = y_0 + \int_{t_0}^t f(s, \phi(s)) ds \quad \forall t \in I \quad \text{(IE)}
$$
反之,满足此积分方程且图像在
引入算子
1. 含义 (Meaning): 什么是算子 T?
-
从积分方程出发: 我们首先看到,求解初值问题 (IVP)
被证明等价于求解积分方程 (IE) 。 -
抽象化右侧操作: 算子
的定义就是将积分方程的右侧抽象出来,看作一个对函数进行的操作。具体来说:$$ (T\psi)(t) = y_0 + \int_{t_0}^t f(s, \psi(s)) ds
$$
1
2
3
这里的 $T$ 是一个**映射 (Map)** 或者说**变换 (Transformation)**。它的特殊之处在于:
* 它的**输入 (Input)** 不是一个数字或向量,而是一个**函数** $\psi$ (这个函数需要是定义在某个包含 $t_0$ 的区间 $I$ 上的连续函数,并且其图像要落在 $f$ 的定义域 $D$ 内)。
* 它的**输出 (Output)** 也是一个**函数** $(T\psi)$ (这个新函数也是定义在 $t$ 上的)。
-
算子的作用: 算子
接收一个函数 ,然后通过以下步骤生成一个新的函数 :- 将输入的函数
代入到 的第二个变量位置,得到 。 - 将这个结果从
积分到 。 - 最后加上初始值
。
- 将输入的函数
2. 目的 (Purpose): 为什么要引入算子 T?
引入算子
-
统一和简化表达: 积分方程
的核心在于,解函数 同时出现在等号两边。使用算子 ,这个方程可以极其简洁地写成:
$$ \phi = T\phi
$$
1
这立刻就揭示了问题的本质:**我们寻找的解 $\phi$ 正是算子 $T$ 的一个不动点 (Fixed Point)**。不动点的意思就是,经过 $T$ 这个变换之后,它保持不变。
-
利用强大的不动点理论: 数学中,关于不动点 (Fixed Point) 的存在性、唯一性以及如何找到它们,已经有非常成熟和强大的理论,其中最核心的就是巴拿赫不动点定理 (Banach’s Fixed-Point Theorem),也叫压缩映射原理 (Contraction Mapping Principle)。这个定理告诉我们,如果在一个合适的函数空间 (Function Space) (具体来说是一个完备度量空间 Complete Metric Space) 中,算子
满足压缩 (Contraction) 条件(即它能把任意两个函数之间的“距离”缩小一个固定的比例),那么 就一定存在且仅存在一个不动点。 -
证明存在唯一性的关键步骤: 将求解 IVP 问题转化为寻找算子
的不动点问题,是证明 ODE 解的存在唯一性定理 (如 Picard-Lindelöf 定理) 的核心策略。证明过程的关键就在于:- 定义一个合适的函数空间 (通常是某个区间上的连续函数空间,配备
度量)。 - 证明这个空间是完备的。
- 证明在一定条件下 (主要是
关于 满足利普希茨条件 Lipschitz Condition),算子 (或者它的某个迭代) 在这个空间的一个子集上是压缩映射。 - 应用巴拿赫不动点定理,直接得出存在唯一不动点
,这个不动点 就是原 IVP 的唯一解。
- 定义一个合适的函数空间 (通常是某个区间上的连续函数空间,配备
-
启发数值解法: 寻找不动点的迭代思想
也直接启发了一种构造解的近似序列的方法,称为皮卡迭代法 (Picard Iteration)。从一个初始猜测函数 (通常取 ) 开始,反复应用算子 :
如果 是压缩映射,这个函数序列 就会收敛到不动点 ,也就是 IVP 的解。
总结:
引入算子
-
简化了问题的表述 (
) - 使得我们可以应用强大的不动点理论 (特别是巴拿赫不动点定理) 来证明解的存在性和唯一性
- 为构造解的近似序列 (如皮卡迭代) 提供了理论基础和方法
高阶方程降阶 (Order Reduction)
-
核心思想: 将高阶问题转化为一阶向量问题。
考虑 阶常微分方程:$$ y^{(n)} = f(t, y, y’, \dots, y^{(n-1)})
$$
初始条件为
-
状态向量 (State Vector) 的引入:
我们定义一个
维的向量 ,它的分量是 及其直到 阶的导数:
$$
1
\mathbf{y}(t) = \begin{pmatrix} y(t) \\\ y'(t) \\\ \vdots \\\ y^{(n-1)}(t) \end{pmatrix} = \begin{pmatrix} y_0(t) \\\ y_1(t) \\\ \vdots \\\ y_{n-1}(t) \end{pmatrix}
$$
1
(这里用 $y_0, y_1, ...$ 表示向量的分量函数,不要与初始值混淆)。
-
一阶系统: 对
求导,我们得到:$$ \mathbf{y}’(t) =
=
$$
1
令
$$
\mathbf{f}(t, \mathbf{y}) =
$$
1
2
3
4
(注意 $\mathbf{f}$ 的定义,它的前 $n-1$ 个分量直接取自输入向量 $\mathbf{y}$ 的第 2 到第 $n$ 个分量,最后一个分量才是由原方程的 $f$ 给出)。
这样,原 $n$ 阶方程就等价于以下的一阶**向量形式 (vectorial form)** 的 ODE 系统:
\$\$
\mathbf{y}' = \mathbf{f}(t, \mathbf{y})
$$
初始条件为
-
结论: 通过这种降阶方法,任何
阶 ODE 的 IVP 都可以转化为一个一阶 ODE 系统的 IVP。因此,我们只需要研究一阶系统 (包括向量形式) 的解的存在唯一性就足够了。同样地,向量形式的 IVP 也可以写成积分方程形式 ,其解也是对应向量算子 的不动点。 -
例子: 对于
,令 。则 , 。 写成向量形式 ,则
$$
1
\mathbf{y}' = \begin{pmatrix} y_0' \\\ y_1' \end{pmatrix} = \begin{pmatrix} y_1 \\\ -y_0 \end{pmatrix} = \mathbf{f}(t, \mathbf{y})
$$
1
这里 $\mathbf{f}(t, \mathbf{y}) = \begin{pmatrix} y_1 \\\ -y_0 \end{pmatrix}$。
牛顿迭代法 (Newton’s Method)
一维情况分析
-
回顾一维情况:
求解 的根 $x^ f(x) x_n$ 的线性近似 (Linear Approximation) 来代替 。令 解出 ,得到新的、更好的近似点 。
迭代公式: (要求 )。
这可以看作是寻找函数 的不动点的迭代过程: 。
如果序列 收敛到 $x^ T x^* = T(x^) f(x^) = 0$。
收敛性分析:-
二次收敛 (Quadratic Convergence): 如果
足够光滑 ( ) 且 $f’(x^) \neq 0 T’(x^) = 0 x^ T(x_n) - T(x^) \approx T’(x^) (x_n - x^) + \frac{T’‘(x^)}{2}(x_n - x^)^2 x_{n+1} - x^* \approx \frac{T’‘(x^)}{2}(x_n - x^)^2$。这意味着误差大约每步平方,收敛非常快 (有效数字位数大约翻倍)。 -
收缩条件 (Contraction Condition): 如果我们不知道根是否存在,可以考虑差分 。如果能保证在某个区间内 $T’(x) \le C < 1 T$ 是一个压缩映射 (Contraction Mapping)),并且迭代值 始终落在这个区间内,那么可以证明 $x_{n+1} - x_n \le C^n x_1 - x_0 (x_n)$ 是一个柯西序列 (Cauchy Sequence)。-
柯西序列: 一个序列 是柯西序列,如果对于任意小的 ,当 足够大时,都有 $a_m - a_n < \epsilon$。直观地说,序列的尾部项都挤在一起了。 -
实数的完备性 (Completeness of
): 中每一个柯西序列都收敛到 中的一个极限。这是 的一个基本性质。 (有理数集) 则不具备这个性质 (例如,逼近 的有理数序列是柯西序列,但在 中不收敛)。 - 技术比喻 (柯西序列 vs 收敛): 收敛序列是“奔向”一个明确的目的地 (极限)。柯西序列是“队员们”自己靠得越来越近,即使我们不知道最终的目的地在哪,但由于空间的“完备性”(没有“洞”),他们必然会聚集到一个点上。
-
-
利用柯西序列和 的完备性,可以证明只要满足 $T’(x) \le C < 1$ 且迭代序列有界,该序列必定收敛到一个根。
-
二次收敛 (Quadratic Convergence): 如果
多维情况
-
多维情况:
求解方程组 ,其中 。
回顾多元函数的可微性 (Differentiability),线性近似 ,其中 是雅可比矩阵 (Jacobi Matrix) (其元素是偏导数 )。
多维牛顿迭代公式: (要求雅可比矩阵可逆)。
这仍然是一个不动点迭代 ,其中 。
收敛性分析更复杂,但基本思想仍然是证明 在不动点 附近是一个压缩映射。这需要衡量向量的“长度”或“距离”(例如欧氏距离) 以及矩阵 (线性变换 ) 的“大小”或“缩放因子”,即矩阵范数 (Matrix Norm) (后续介绍)。
度量空间 (Metric Spaces)
定义
-
定义: 一个度量空间 (Metric Space)
由一个集合 和一个距离函数 (Distance Function) (或称为度量 Metric) 组成,满足以下性质 ( ):- (M1) 非负性 (Non-negativity):
,且 。 - (M2) 对称性 (Symmetry):
。 - (M3) 三角不等式 (Triangle Inequality):
。-
技术比喻: 度量空间就是一个配备了“尺子”
的集合 。这把“尺子”必须满足我们对距离的基本直觉:距离非负、从 A 到 B 和从 B 到 A 距离一样、从 A 绕道 Z 再到 B 不会比直接从 A 到 B 更近。
-
技术比喻: 度量空间就是一个配备了“尺子”
- (M1) 非负性 (Non-negativity):
Example
-
例子:
-
上的欧氏距离 (Euclidean Distance) 。包括 上的 $x-y \mathbb{C}z-w $。 -
上的 (曼哈顿距离 Manhattan Distance) $d_1(\mathbf{x}, \mathbf{y}) = \sum_{i=1}^nx_i - y_i d_\infty$ (最大距离 Max Distance / Chebyshev Distance) $d_\infty(\mathbf{x}, \mathbf{y}) = \max_{1 \le i \le n}x_i - y_i $。 上的法国铁路度量 (French Railway Metric) (有点奇怪,但符合定义)。-
空间 (Hilbert Cube): 无穷复数序列 且 $\suma_n ^2 < \infty d((a_n), (b_n)) = \sqrt{\suma_n - b_n ^2}$。 - 任意集合
上的离散度量 (Discrete Metric): if , if 。 - 图上的最短路径距离 (Shortest Path Distance)。
-
汉明距离 (Hamming Distance): 长度为
的比特串之间不同位的个数。 - 度量空间的子集
继承原度量 构成子空间 (Subspace)。 -
连续函数空间 (Space of Continuous Functions) : 定义在闭区间 上的所有连续实函数,配备一致收敛度量 (Metric of Uniform Convergence) $d_\infty(f, g) = \max_{x \in [a,b]}f(x) - g(x) d_\infty(f_n, f) \to 0 f_n$ 一致收敛 (Uniformly Converges) 到 。 -
广义度量空间 (Generalized Metric Spaces): 允许距离为
。
-
Analysis on Metric Spaces
-
度量空间中的分析概念:
可以在度量空间中定义开球 (Open Ball) ,进而定义开集、闭集、序列收敛、函数连续性等拓扑概念。
注意: 中一些熟悉的性质在一般度量空间中不一定成立,例如波尔查诺 - 魏尔斯特拉斯定理 (Bolzano-Weierstrass Theorem) (有界序列必有收敛子序列) 在离散度量空间或 空间中就不成立 (Slide 38)。
Complete Metric Spaces
-
完备度量空间 (Complete Metric Spaces) (Slides 40-48):
柯西序列 (Cauchy Sequence) 在度量空间中的定义: s.t. 。
完备性 (Completeness): 一个度量空间 称为完备的 (Complete),如果其中每一个柯西序列都收敛到 中的一个点。
例子: (及 ) 在欧氏度量下是完备的。 在 度量下是完备的 (这是分析学中的重要定理:一致收敛的连续函数序列的极限函数也是连续的)。离散度量空间总是完备的。
反例: (有理数集) 在通常度量下不完备。开区间 在通常度量下不完备 (序列 是柯西的但极限 不在区间内)。 -
重要性质: 度量空间
的子空间 是完备的当且仅当 是 中的闭集 (Closed Set) (如果 本身是完备的)。
完备性检验
正式定义 (Formal Definition):
一个度量空间
要理解这个定义,我们需要拆解两个关键概念:
-
柯西序列 (Cauchy Sequence):
-
直观理解: 一个序列
是柯西序列,意味着当 变得越来越大时,序列中的点彼此之间靠得越来越近。它们“挤”在了一起。 -
形式化定义: 对 任意
,都 存在 一个正整数 ,使得对 所有 ,都有 。 - 注意: 柯西序列本身并不保证一定收敛到 空间内 的某个点。它只保证序列的项最终彼此无限接近。
-
直观理解: 一个序列
-
收敛到
中的一个点 (Converge to a Point in ):-
直观理解: 序列
收敛到点 ,意味着当 变得越来越大时, 会无限接近于点 。 -
形式化定义: 存在一个点
,使得对 任意 ,都 存在 一个正整数 ,使得对 所有 ,都有 。 -
关键点: 对于完备性,这个极限点
必须 属于 空间 本身。
-
直观理解: 序列
如何判断一个度量空间
判断方法主要有两种思路:证明它是完备的,或者证明它不是完备的。
A. 证明
-
直接使用定义:
- 取
中一个 任意的 (arbitrary) 柯西序列 。 -
目标: 证明这个序列
收敛,并且其极限点 必须属于 。 -
步骤:
a. 利用
是柯西序列的性质,构造出一个候选的极限点 。这通常是最难的一步,可能需要利用空间的具体结构(比如实数的序结构、向量空间的结构等)。 b. 证明 确实收敛到这个候选点 ,即证明 。 c. 证明这个极限点 确实在空间 中,即 。
- 取
-
利用已知完备空间:
- 如果
可以被看作是一个已知的完备度量空间 的一个闭子集 (Closed Subset),那么 也是完备的。 -
回顾: 一个子集
是闭子集,意味着 包含了它所有的极限点 (Limit Points) / 聚点 (Accumulation Points)。换句话说,如果 中的一个序列 全部落在 中 ( for all ),并且这个序列收敛到 中的某个点 ( in ),那么这个极限点 也必须在 中 ( )。 -
逻辑: 如果
是 中的一个柯西序列,因为它也是 中的序列,并且 是完备的,所以 在 中收敛到一个极限点 。因为 是 的闭子集,且所有 ,所以这个极限点 必须也属于 。因此, 中的任意柯西序列都收敛到 中的点,故 是完备的。 -
例子: 实数空间 配备标准度量 $d(x,y)=x-y [0, 1] \mathbb{R} [0, 1]$ 也是完备的。
- 如果
B. 证明
-
找反例 (Find a Counterexample): 只需要找到 一个 在
中的柯西序列 ,它 不收敛 到 中的任何点。 -
常见的“不收敛到
中”的情况: 这个柯西序列实际上收敛到一个点 ,但是这个点 不在 空间 中。 -
步骤:
a. 构造一个具体的序列
,其中所有 。 b. 证明 是一个柯西序列。 c. 证明这个序列没有极限点在 中。通常是证明它收敛到一个 。
经典例子:
-
完备的例子:
-
实数空间 和复数空间 ,使用标准度量 $d(x,y)=x-y $。这是分析学的基础。 - 欧几里得空间
,使用欧几里得度量 。 - 任何闭区间 (Closed Interval)
。 - 任何有限 (Finite) 度量空间。(因为柯西序列最终各项必然相等)
-
巴拿赫空间 (Banach Space): 定义为完备的赋范向量空间 (Normed Vector Space)。例如,闭区间 上的连续函数空间 配备上确界范数 (Supremum Norm) $|f|\infty = \sup{x \in [a,b]}f(x) d(f, g) = |f-g|_\infty$ 是完备的。
-
-
不完备的例子:
-
有理数空间 ,使用标准度量 $d(x,y)=x-y $。 -
反例: 构造一个有理数序列趋近于
(例如 )。这个序列在 中是柯西序列,但它的极限 不在 中。 有“洞”。
-
反例: 构造一个有理数序列趋近于
- 开区间
,使用标准度量。-
反例: 序列
(对于 )。所有 。这是一个柯西序列(因为它在 中收敛到 0)。但它的极限 0 不在 这个空间里。
-
反例: 序列
-
配备 L1 范数 (L1 Norm) $|f|_1 = \int_a^bf(x) dx d_1(f, g) = \int_a^bf(x)-g(x) dx$。 -
反例: 可以构造一个连续函数序列,它在
度量下收敛到一个阶跃函数 (Step Function) (它在某个点不连续),而阶跃函数不属于 。
-
反例: 可以构造一个连续函数序列,它在
-
完备性例子
好的,我们来详细梳理一下 Slides 中提到的关于度量空间 (Metric Spaces) 的一些具体的正例 (Examples) 和反例 (Counterexamples),并解释它们各自说明了什么重要概念或特性。
A. 标准且重要的度量空间实例
-
欧氏空间 (Euclidean Space) (
) (Slide 27, Ex 1):-
描述: 集合是 维实向量空间 ,距离是标准的欧氏距离 。 (配 $x-y \mathbb{C}z-w n=1 n=2$ 的特例。 - 说明: 这是我们最熟悉的度量空间,许多分析概念的直观来源。
-
关键性质 (Slide 41): 欧氏空间是完备的 (Complete)。这意味着任何在
中看起来“应该”收敛的序列 (即柯西序列) 确实收敛到 中的一个点。这是许多分析定理的基础。
-
-
上的其他度量 ( ):-
描述: 仍然是
空间,但使用不同的距离定义:-
$d_1(\mathbf{x}, \mathbf{y}) = \sum_{i=1}^n x_i - y_i $ (曼哈顿距离) -
$d_\infty(\mathbf{x}, \mathbf{y}) = \max_{1 \le i \le n} x_i - y_i $ (最大距离/切比雪夫距离)
-
-
说明:
- 不同的几何形状: 这些度量下的“单位球”形状不同 (菱形、正方形,见 Slide 36)。
-
拓扑等价性: 尽管几何形状不同,这三种度量 (
) 在 上是强等价的 (Strongly Equivalent) (Slide 45 exercise b)。这意味着它们定义了相同的开集、闭集、收敛序列和柯西序列。因此,在讨论收敛性、连续性等拓扑性质时,它们是等价的。 -
完备性: 由于它们与完备的
强等价,因此 和 也都是完备的。
-
描述: 仍然是
-
连续函数空间 (
):-
描述: 集合是定义在闭区间 上的所有连续实函数,距离是 $d_\infty(f, g) = \max_{x \in [a,b]}f(x) - g(x) $ (一致收敛度量)。 -
说明: 这是泛函分析中的核心例子。
- 距离的含义: 两个函数之间的距离是它们函数值差异的最大值。
-
收敛性: 在此空间中,
意味着函数序列 一致收敛 (Uniformly Converges) 到函数 。 -
球的形状: 以
为中心的球 包含所有图像完全落在 图像上下宽度 形成的“带子”内的连续函数 。
- 关键性质: 这个空间是完备的。这是微积分中一个重要定理的体现:一致收敛的连续函数序列的极限函数仍然是连续的。这个完备性对于证明 ODE 解的存在性 (Picard 迭代法发生在此类空间中) 至关重要。
-
-
希尔伯特立方 (
):-
描述: 集合是所有复数 (或实数) 无穷序列 $(a_n){n=0}^\infty \sum{n=0}^\inftya_n ^2 < \infty d((a_n), (b_n)) = \sqrt{\sum_{n=0}^\inftya_n - b_n ^2}$。 - 说明: 这是无限维空间的一个典型例子。
- 关键性质: 这个空间是完备的。
-
反例性质: 波尔查诺 - 魏尔斯特拉斯定理在此失效。例如,考虑序列 (第 个位置为 1,其余为 0)。这个序列是有界的 (所有点到零向量的距离都是 1),但它没有任何收敛的子序列,因为任意两个不同的项 之间的距离都是 (如果是实数序列) 或 (如果是复数序列, 距离是 $\sqrt{\suma_n-b_n ^2}$),距离不趋于 0。这表明有界性在无限维空间中不一定能保证序列紧致性。
-
-
离散度量空间 (Discrete Metric Space):
-
描述: 任何集合
都可以配备离散度量: 如果 , 如果 。 -
说明:
-
简单的球: 开球
要么只包含 (如果 ),要么包含整个空间 (如果 )。
-
简单的球: 开球
- 关键性质 : 离散度量空间总是完备的。因为在离散度量下,一个序列是柯西序列当且仅当这个序列最终恒定 (即从某项开始,后面所有项都相同),而恒定序列显然收敛。
-
反例性质: 波尔查诺 - 魏尔斯特拉斯定理在此失效 (如果
是无限集)。例如,在 配备离散度量,序列 是有界的 (所有点都在 内),但它没有任何收敛子序列 (因为任何子序列都不会最终恒定)。
-
描述: 任何集合
B. 说明特定概念的反例 (Counterexamples) / 特殊情况
-
完备性的反例 (Incompleteness):
-
有理数空间 (
): 这是最经典的例子。序列 (逼近 ) 是 中的柯西序列,但它的极限 不在 中。说明 不完备。 -
穿孔实线 (
): 这是一个非闭子集,因此不完备。具体地,序列 属于此空间。它在 中是柯西序列 (极限为 0)。因此它在这个子空间中也是柯西序列。但它的极限 0 不属于 。这表明子空间中存在一个“洞”。 -
开区间 (
): 类似地,序列 在 中是柯西序列,但极限 0 不在 中。 -
非闭子空间: 一般地,完备度量空间的子空间是完备的当且仅当这个子空间是闭集。上面的
和 都是 中非闭的子集。 -
有界但不完备的例子: 考虑在 上定义度量 $d(x, y) =\arctan(x) - \arctan(y) d(x,y) = d_E(x/(1+x ), y/(1+ y )) \pi \mathbb{R} a_n = n \arctan(n) \to \pi/2 \mathbb{R}$ 中不收敛 (没有极限点)。这说明完备性不是一个纯粹的拓扑性质,它依赖于度量本身。
-
有理数空间 (
-
波尔查诺 - 魏尔斯特拉斯性质失效的反例 (Failure of B-W):
- 如前所述,离散度量空间 (当集合无限时) 和无限维空间
都是有界序列不一定有收敛子序列的例子。这说明有限维欧氏空间的这个良好性质不能随意推广。
- 如前所述,离散度量空间 (当集合无限时) 和无限维空间
-
拓扑等价但度量不等价 (Topologically Equivalent but Metrically Different):
, , (Slide 36): 它们的球形状不同,度量值也不同,但在收敛性等拓扑层面表现一致。- Slide 45/46 中的练习构造的
上的度量 (如 ): 它与 诱导相同的开集和收敛序列 (拓扑等价),但它不是强等价的 (因为 是有界的而 不是)。这说明拓扑等价比强等价更弱。
巴拿赫不动点定理 (Banach’s Fixed-Point Theorem)
压缩映射
-
压缩映射 (Contraction Mapping) (Slide 49):
设 是度量空间,映射 称为一个压缩映射,如果存在一个常数 满足 ,使得对于所有 ,都有:$$ d(T(x), T(y)) \le C \cdot d(x, y)
$$
1
2
常数 $C$ 称为**压缩常数 (Contraction Constant)**。
压缩映射会使得任意两点在映射后的距离**严格**缩小 (按比例 $C$)。
检验关键点
- 压缩映射
为自同态映射:即映射从 映射到自身 - 压缩常数严格小于1
对于定义在单一区间上的一维连续函数,我们只需要验证:
- 自映射
- 验证其定义域上任意一点的微分都严格小于1
定理
-
定理条件
- 完备的度量空间(柯西序列收敛到度量空间中的点/完备度量空间的闭子空间)
- 压缩映射 (自映射+压缩)
-
定理叙述:
巴拿赫不动点定理 (Banach’s Fixed-Point Theorem / Contraction Mapping Theorem):
设 是一个完备 (Complete) 的度量空间, 是一个压缩映射 (Contraction Mapping)。则: 在 中存在唯一 (Unique) 的不动点 (Fixed Point) $x^ T(x^) = x^*$。- 对于任意初始点
,由迭代 ( ) 生成的序列 收敛 (Converges) 到不动点 。 - 误差估计 (Error Estimates):
$$ d(x_n, x^*) \le \frac{C^n}{1-C} d(x_1, x_0) \quad (\text{先验估计, a priori estimate})
Matrix Norm 矩阵范数
矩阵范数是向量范数 (Vector Norm) 概念到矩阵上的推广。我们不仅想衡量向量的大小,也想衡量矩阵的“大小”或其作为线性变换的“强度”。
-
为什么需要矩阵范数?
- 衡量矩阵元素的大小。
- 衡量线性变换
对向量的“拉伸”程度。 - 分析矩阵运算 (如乘积、求逆) 的敏感性/稳定性 (例如,条件数 Condition Number)。
- 分析迭代算法 (如求解线性方程组
的迭代法) 的收敛性。
-
定义:
一个函数
(其中 是 或 ,表示矩阵空间) 被称为一个矩阵范数,如果它满足以下性质 (对于任意 和标量 ):-
非负性 (Non-negativity):
-
正定性 (Positive Definiteness):
当且仅当 (零矩阵) -
齐次性 (Homogeneity / Absolutely Scalable): $|\alpha A| = \alpha |A|$ -
三角不等式 (Triangle Inequality):
这些与向量范数的定义完全相同。但对于方阵 (Square Matrices) (
),通常还需要一个额外的、非常重要的性质:-
次可乘性 / 相容性 (Submultiplicativity / Consistency):
-
意义: 这个性质将矩阵乘法与范数联系起来,保证了矩阵乘积的范数不会超过各自范数的乘积。这对于分析迭代过程 (如
) 或误差传播至关重要。
-
意义: 这个性质将矩阵乘法与范数联系起来,保证了矩阵乘积的范数不会超过各自范数的乘积。这对于分析迭代过程 (如
-
非负性 (Non-negativity):
常见矩阵范数
-
诱导范数 / 算子范数 (Induced Norm / Operator Norm):
这是最常用的一类矩阵范数,它们是由向量范数“诱导”产生的。给定定义在
和 上的向量范数 (分别记作 ),矩阵 的诱导范数定义为:
$$
1
\|A\| = \sup_{\mathbf{x} \ne \mathbf{0}} \frac{\|A\mathbf{x}\|_v}{\|\mathbf{x}\|_v} = \sup_{\|\mathbf{x}\|_v = 1} \|A\mathbf{x}\|_v
$$
1
2
* **几何意义:** $\|A\|$ 表示线性变换 $A$ 能将单位向量 (在 $\|\cdot\|_v$ 意义下) 拉伸的最大倍数。
* **重要性质:** 由向量范数诱导出的矩阵范数**自动满足**次可乘性 (如果作用空间和目标空间的向量范数相同,对于方阵来说)。并且,它们还满足与向量范数的**相容性 (Compatibility):** $\|A\mathbf{x}\|_v \le \|A\| \|\mathbf{x}\|_v$。
-
常见的诱导范数 (对于方阵
):-
列和范数 / 1- 范数 (
): 由向量 - 范数诱导。$$ |A|_1 = \max_{1 \le j \le n} \sum_{i=1}^n |a_{ij}| \quad (\text{最大列绝对值和 Max column sum})
-
$$
1
2. **谱范数 / 2- 范数 ($p=2$):** 由向量 $l_2$- 范数 (欧氏范数) 诱导。
$$ |A|2 = \sqrt{\lambda{\max}(A^*A)} = \sigma_{\max}(A) \quad (\text{A 的最大奇异值 Max singular value})
$$
1
2
3
4
其中 $A^*$ 是 $A$ 的共轭转置 (对实矩阵即为转置 $A^T$),$\lambda_{\max}$ 表示最大特征值。
3. **行和范数 / $\infty$-范数 ($p=\infty$):** 由向量 $l_\infty$-范数 (最大范数) 诱导。
\$\$
\|A\|\_\infty = \max\_{1 \le i \le n} \sum\_{j=1}^n |a\_{ij}| \quad \(\text{最大行绝对值和 Max row sum})
$$
-
非诱导范数:
也存在不直接由向量范数诱导但满足基本范数公理和次可乘性的矩阵范数。
- 弗罗贝尼乌斯范数 (Frobenius Norm):
$$
1
\|A\|\_F = \sqrt{\sum\_{i=1}^m \sum\_{j=1}^n |a\_{ij}|^2} = \sqrt{\text{trace}\(A^\*A)}
$$
1
* **特点:** 它就像把 $m \times n$ 矩阵看成 $mn$ 维向量后的欧氏范数。计算简单。它满足次可乘性,但它**不是**由任何向量 $l_p$-范数诱导产生的 (除了 $n=1$ 或 $m=1$ 的平凡情况)。
-
范数的等价性 (Equivalence of Norms):
在有限维空间 (如 ) 中,任意两种范数都是等价的 (Equivalent)。这意味着如果一个矩阵序列在一范数下收敛到零,那么它在任何其他范数下也收敛到零。这在理论分析中很有用。但需要注意,不同范数的值可能相差很大,等价性常数可能影响实际计算中的收敛速度或稳定性界限。 -
与巴拿赫不动点定理的联系 (潜在的):
考虑线性迭代 。这个迭代收敛到不动点 的一个充分条件是 ,其中 是某个满足次可乘性的矩阵范数。-
推导: 设 $x^
x^ = Ax^* + c e_k = x_k - x^ e_{k+1} = x_{k+1} - x^ = (Ax_k + c) - (Ax^* + c) = A(x_k - x^*) = Ae_k$。 - 因此,
(需要范数与向量范数相容)。 - 重复这个过程,
。如果 ,那么当 时, ,所以 ,即 。 - 这本质上是巴拿赫不动点定理在线性变换和特定范数下的体现。映射
在范数小于 1 的条件下是压缩映射。
-
推导: 设 $x^
谱半径与谱范数
1. 谱半径 (Spectral Radius) -
-
定义 (Definition):
对于一个方阵 (Square Matrix) (或 ),它的谱半径 (Spectral Radius) 定义为其所有特征值 (Eigenvalues) 的绝对值 (Absolute Value / Modulus) 中的最大值。$$ \rho(A) = \max { |\lambda_1|, |\lambda_2|, \dots, |\lambda_n| } = \max_{\lambda \in \sigma(A)} |\lambda|
$$
1
2
其中 $\sigma(A)$ 是 $A$ 的特征值集合 (谱)。
* ***技术比喻:*** 想象所有特征值在复平面上的点,谱半径就是这些点到原点距离的最大值,也就是包含所有特征值的最小圆盘的半径。
-
计算 (Calculation):
需要先计算出矩阵 的所有特征值 (解特征方程 ),然后取这些特征值绝对值的最大者。 -
性质 (Properties):
。-
$\rho(\alpha A) = \alpha \rho(A) \alpha$。 对任意正整数 。-
关键: 谱半径
不是一个矩阵范数 (Matrix Norm)。因为它不满足范数的所有公理,特别是三角不等式 和正定性 ( ) 可能不成立。-
反例 (Counterexample): 考虑
。它的唯一特征值是 0,所以 。但 不是零矩阵,违反了正定性。 -
反例 (Triangle Inequality): 考虑
和 。 , ,所以 。但是 ,特征值为 ,所以 。这里 。
-
反例 (Counterexample): 考虑
-
主要用途 (Main Use):
谱半径主要用于判断矩阵幂序列 的收敛性以及相关迭代过程的收敛性。-
定理 (Gelfand’s Formula related):
(零矩阵) 当且仅当 (if and only if) 。 -
应用: 在分析迭代法
的收敛性时,该迭代收敛到唯一解的充要条件是 。
-
定理 (Gelfand’s Formula related):
2. 谱范数 (Spectral Norm) -
-
定义 (Definition):
矩阵 (可以是长方阵) 的谱范数 (Spectral Norm),也称为2- 范数 (2-Norm),是由向量的欧几里得范数 (Euclidean Norm / -norm) 诱导 (Induce) 产生的矩阵范数。
$$ |A|2 = \sup{\mathbf{x} \ne \mathbf{0}} \frac{|A\mathbf{x}|2}{|\mathbf{x}|_2} = \sup{|\mathbf{x}|_2 = 1} |A\mathbf{x}|_2
$$
1
* ***技术比喻:*** 它衡量了线性变换 $A$ 对向量的欧氏长度的最大“拉伸”因子。
-
计算 (Calculation):
谱范数等于矩阵 的最大奇异值 (Largest Singular Value) 。奇异值 是矩阵 $A^A AA^$) 特征值的非负平方根 (其中 是 的共轭转置,对实矩阵即为转置 )。$$ |A|2 = \sigma{\max}(A) = \sqrt{\lambda_{\max}(A^*A)}
$$
1
其中 $\lambda_{\max}(A^*A)$ 是半正定矩阵 $A^*A$ 的最大特征值。
-
性质 (Properties):
是一个真正的矩阵范数。它满足所有范数公理:非负性、正定性、齐次性、三角不等式。- 它也满足次可乘性 (Submultiplicativity):
。 - 与向量范数相容:
。 - 如果
是正规矩阵 (Normal Matrix) (即 $A^A = AA^$,包括 Hermitian/对称矩阵, Skew-Hermitian/斜对称矩阵, Unitary/正交矩阵),那么它的谱范数等于它的谱半径: 。这是因为对于正规矩阵,奇异值等于其特征值的绝对值。 - 对于非正规矩阵 (Non-normal Matrix),一般有
。等号不一定成立。-
例子: 还是
。$A^A = = A^A |A|_2 = \sqrt{1} = 1 \rho(A) = 0 \rho(A) < |A|_2$。
-
例子: 还是
-
主要用途 (Main Use):
- 衡量线性变换
对欧氏距离的影响程度,即最大拉伸率。 - 在线性方程组
的条件数 (Condition Number) 计算中使用: 。条件数衡量了问题对输入的微小扰动的敏感程度。 - 在奇异值分解 (Singular Value Decomposition, SVD) 和相关应用 (如低秩逼近、主成分分析 PCA) 中自然出现。
- 提供
的充分条件:如果 ,那么必然有 ,因此 。但这只是充分条件,不是必要的 (可能 但 )。
- 衡量线性变换
总结对比 (Summary Comparison):
特性 (Feature) | 谱半径 (Spectral Radius) |
谱范数 (Spectral Norm) |
||
---|---|---|---|---|
定义基于 (Based on) | 特征值 (Eigenvalues) | 奇异值 (Singular Values) / |
||
定义域 (Applies to) | 方阵 (Square matrices) |
任意矩阵 (Any matrix) |
||
是否范数 (Is a Norm?) | 否 (No) | 是 (Yes) | ||
关键关系 (Relation) | ||||
主要用途 (Main Use) | 迭代收敛性 ( |
最大拉伸因子, 条件数, SVD, 提供收敛的充分条件 ( |