(泛函分析)压缩映射
在泛函分析中,压缩映射(Contraction Mapping),也称为压缩算子或压缩函数,是一个非常重要的概念,特别是在不动点理论中。一个映射 T : X → X T:X\to X T:X→X在一个度量空间 ( X , d ) (X,d) (X,d)上被称为压缩映射,如果存在一个常数 k ∈ [ 0 , 1 ) k\in[0,1) k∈[0,1)使得对于所有 x , y ∈ X x,y\in X x,y∈X都有:
d ( T ( x ) , T ( y ) ) ≤ k ⋅ d ( x , y ) . d(T(x),T(y))\leq k\cdot d(x,y). d(T(x),T(y))≤k⋅d(x,y).
这意味着映射 T T T将任何两点之间的距离缩小至少一个比例因子 k < 1 k<1 k<1。
核心性质
- Lipschitz连续性:压缩映像必然是Lipschitz连续的,且Lipschitz常数 k < 1 k < 1 k<1。
- 唯一不动点:在完备度量空间中,压缩映像存在唯一不动点(巴拿赫不动点定理)。
压缩映射原理(Banach不动点定理)
压缩映射原理是关于完备度量空间中的一个重要定理,它指出:如果 ( X , d ) (X,d) (X,d)是一个完备的度量空间,并且 T : X → X T:X\to X T:X→X是一个压缩映射,则 T T T恰好有一个不动点 x ∗ ∈ X x^*\in X x∗∈X,即满足 T ( x ∗ ) = x ∗ T(x^*)=x^* T(x∗)=x∗。此外,从 X X X中任意一点 x 0 x_0 x0出发,通过迭代序列 x n + 1 = T ( x n ) x_{n+1}=T(x_n) xn+1=T(xn)可以得到这个不动点 x ∗ x^* x∗,并且该序列收敛于 x ∗ x^* x∗。
示例问题
考虑求解方程 x 2 − x − 1 = 0 x^2-x-1=0 x2−x−1=0的正根。我们可以使用压缩映射原理来近似找到这个根。
解决方案
首先,注意到方程可以重写为 x = 1 + 1 x x=1+\frac{1}{x} x=1+x1。定义映射 T ( x ) = 1 + 1 x T(x)=1+\frac{1}{x} T(x)=1+x1。我们的目标是证明 T T T在某个区间内是一个压缩映射,并利用Banach不动点定理找到其不动点作为原方程的解。
-
选择合适的区间:假设我们关注 x > 0 x>0 x>0的情况,特别是 x ≥ 1 x\geq1 x≥1。这是因为我们知道方程的正根大于1。
-
验证 T T T是否为压缩映射:我们需要检查是否存在 k < 1 k<1 k<1使得对于所有的 x , y ≥ 1 x,y\geq1 x,y≥1有
∣ T ( x ) − T ( y ) ∣ = ∣ 1 + 1 x − ( 1 + 1 y ) ∣ = ∣ 1 x − 1 y ∣ = ∣ y − x x y ∣ ≤ k ∣ x − y ∣ . |T(x)-T(y)|=\left|1+\frac{1}{x}-(1+\frac{1}{y})\right|=\left|\frac{1}{x}-\frac{1}{y}\right|=\left|\frac{y-x}{xy}\right|\leq k|x-y|. ∣T(x)−T(y)∣= 1+x1−(1+y1) = x1−y1 = xyy−x ≤k∣x−y∣.注意到 1 x y \frac{1}{xy} xy1在 x , y ≥ 1 x,y\geq1 x,y≥1时最大值为1(当 x = y = 1 x=y=1 x=y=1时取等号),但实际上由于 x , y > 1 x,y>1 x,y>1, 1 x y < 1 \frac{1}{xy}<1 xy1<1,因此我们可以选取 k = 1 min ( x , y ) 2 k=\frac{1}{\min(x,y)^2} k=min(x,y)21作为一个合理的估计,但为了简化,通常我们会尝试直接找到一个固定的 k < 1 k<1 k<1适用于整个区间。在这里,考虑到 x , y ≥ 1 x,y\geq1 x,y≥1,我们可以选取 k = 1 2 k=\frac{1}{2} k=21作为一个合理的估计,因为 1 x y \frac{1}{xy} xy1实际上会更小。
-
应用压缩映射原理:既然我们已经表明 T T T在区间 [ 1 , + ∞ ) [1,+\infty) [1,+∞)上是压缩的(尽管上面的分析需要更精确地确定 k k k的具体范围),根据Banach不动点定理,我们知道存在唯一的 x ∗ > 1 x^*>1 x∗>1使得 T ( x ∗ ) = x ∗ T(x^*)=x^* T(x∗)=x∗。我们可以通过迭代法 x n + 1 = T ( x n ) x_{n+1}=T(x_n) xn+1=T(xn)从任一初始值开始逼近 x ∗ x^* x∗。
例如,如果我们从 x 0 = 2 x_0=2 x0=2开始,那么迭代过程将是:
- x 1 = T ( 2 ) = 1 + 1 2 = 1.5 x_1=T(2)=1+\frac{1}{2}=1.5 x1=T(2)=1+21=1.5
- x 2 = T ( 1.5 ) = 1 + 2 3 = 1.666... x_2=T(1.5)=1+\frac{2}{3}=1.666... x2=T(1.5)=1+32=1.666...
- 继续此过程直到数值稳定。
通过这种方式,我们可以逼近方程 x 2 − x − 1 = 0 x^2-x-1=0 x2−x−1=0的正根,这也是黄金分割比 ϕ = 1 + 5 2 \phi=\frac{1+\sqrt{5}}{2} ϕ=21+5。
示例问题:在区间 [ 0 , 1 ] [0, 1] [0,1] 上,求方程 x = cos ( x ) x = \cos(x) x=cos(x) 的解。
解析:
- 定义映射:设 T ( x ) = cos ( x ) T(x) = \cos(x) T(x)=cos(x),空间为 X = [ 0 , 1 ] X = [0, 1] X=[0,1],度量 d ( x , y ) = ∣ x − y ∣ d(x, y) = |x - y| d(x,y)=∣x−y∣。
- 验证压缩性:
∣ T ( x ) − T ( y ) ∣ = ∣ cos ( x ) − cos ( y ) ∣ ≤ ∣ x − y ∣ ( 由中值定理,导数 ∣ sin ( c ) ∣ ≤ 1 ) . |T(x) - T(y)| = |\cos(x) - \cos(y)| \leq |x - y| \quad (\text{由中值定理,导数 } |\sin(c)| \leq 1). ∣T(x)−T(y)∣=∣cos(x)−cos(y)∣≤∣x−y∣(由中值定理,导数 ∣sin(c)∣≤1).
但此处 k = 1 k = 1 k=1,不满足严格压缩。需缩小定义域:- 在 X = [ 0 , π 4 ] X = [0, \frac{\pi}{4}] X=[0,4π] 上, ∣ cos ′ ( x ) ∣ = ∣ sin ( x ) ∣ ≤ sin ( π 4 ) = 2 2 < 1 |\cos'(x)| = |\sin(x)| \leq \sin(\frac{\pi}{4}) = \frac{\sqrt{2}}{2} < 1 ∣cos′(x)∣=∣sin(x)∣≤sin(4π)=22<1,故 T T T 为压缩映像。
- 应用定理:由巴拿赫不动点定理,存在唯一解 x ∗ ∈ [ 0 , π 4 ] x^* \in [0, \frac{\pi}{4}] x∗∈[0,4π],且迭代 x n + 1 = cos ( x n ) x_{n+1} = \cos(x_n) xn+1=cos(xn) 收敛于 x ∗ x^* x∗。
示例问题 微分方程边值问题
问题:求解边值问题:
{ u ′ ′ ( x ) = f ( x , u ( x ) ) , x ∈ [ 0 , 1 ] , u ( 0 ) = u ( 1 ) = 0. \begin{cases} u''(x) = f(x, u(x)), & x \in [0, 1], \\ u(0) = u(1) = 0. \end{cases} {u′′(x)=f(x,u(x)),u(0)=u(1)=0.x∈[0,1],
解析:
- 转化为积分方程:通过格林函数法,问题等价于:
u ( x ) = ∫ 0 1 G ( x , t ) f ( t , u ( t ) ) d t ≜ T ( u ) ( x ) , u(x) = \int_0^1 G(x, t) f(t, u(t)) \, dt \triangleq T(u)(x), u(x)=∫01G(x,t)f(t,u(t))dt≜T(u)(x),
其中 G ( x , t ) G(x, t) G(x,t) 为格林函数。 - 验证压缩性:在适当函数空间(如 C [ 0 , 1 ] C[0, 1] C[0,1])中,若 f f f 满足Lipschitz条件:
∣ f ( t , u ) − f ( t , v ) ∣ ≤ L ∣ u − v ∣ , |f(t, u) - f(t, v)| \leq L |u - v|, ∣f(t,u)−f(t,v)∣≤L∣u−v∣,
则 T T T 为压缩映像(压缩常数 k = L ⋅ max x , t ∣ G ( x , t ) ∣ k = L \cdot \max_{x,t} |G(x, t)| k=L⋅maxx,t∣G(x,t)∣)。 - 存在唯一解:由巴拿赫不动点定理,方程存在唯一解 u ∗ ∈ C [ 0 , 1 ] u^* \in C[0, 1] u∗∈C[0,1]。
示例问题:经济学中的均衡问题
问题:在供需模型中,求市场均衡价格 p ∗ p^* p∗,使得供给量等于需求量: S ( p ∗ ) = D ( p ∗ ) S(p^*) = D(p^*) S(p∗)=D(p∗)。
解析:
- 定义映射:设 T ( p ) = p − α ( S ( p ) − D ( p ) ) T(p) = p - \alpha (S(p) - D(p)) T(p)=p−α(S(p)−D(p)),其中 α > 0 \alpha > 0 α>0 为调整参数。
- 验证压缩性:若供需函数满足:
∣ S ′ ( p ) − D ′ ( p ) ∣ ≤ 1 α , |S'(p) - D'(p)| \leq \frac{1}{\alpha}, ∣S′(p)−D′(p)∣≤α1,
则 T T T 为压缩映像。 - 迭代求解:通过迭代 p n + 1 = T ( p n ) p_{n+1} = T(p_n) pn+1=T(pn),可数值逼近均衡价格 p ∗ p^* p∗。