【NA】不动点原理

基本思路.

  • 设映射 F : D ⊂ R n → R n F:D⊂R^n→R^n F:DRnRn,考虑方程 F ( x ) = 0 F(x)=0 F(x)=0,求解思路是将其转化为等价形式 x = G ( x ) x=G(x) x=G(x),其中映射 G : D ⊂ R n → R n G:D⊂R^n→R^n G:DRnRn,称为不动点方程。若存在点 x ∗ ∈ D x^*∈D xD 使得 x ∗ = G ( x ∗ ) x^*=G(x^*) x=G(x) 成立,那么 x ∗ x^* x 称为映射 G G G 的不动点,由等价性可知 x ∗ x^* x 就是原方程 F ( x ) = 0 F(x)=0 F(x)=0 的解
  • 于是求解方程 F ( x ) = 0 F(x)=0 F(x)=0 转化为求解映射 G ( x ) G(x) G(x) 的不动点,关键在于判断 G ( x ) G(x) G(x) 是否存在不动点以及如何求解该不动点。
  • 不动点原理】是用于判断映射 G G G 在定义域 D D D 上是否存在不动点的理论。

压缩映射.

  • 考虑连续映射 G : D ⊂ R n → R n G:D⊂R^n→R^n G:DRnRn,对于 ∀ x , y ∈ D ∀x,y∈D x,yD,恒有 ∣ ∣ G ( x ) − G ( y ) ∣ ∣ ≤ α ∣ ∣ x − y ∣ ∣ ||G(x)-G(y)||≤\alpha||x-y|| G(x)G(y)αxy 成立。
  • α ∈ ( 0 , 1 ) \alpha∈(0,1) α(0,1) 时, G ( x ) G(x) G(x) D D D 上的压缩映射;另外 α = 1 \alpha=1 α=1 时, G ( x ) G(x) G(x) D D D 上的非膨胀映射
  • 定理一】对于压缩映射 G : R n → R n G:R^n→R^n G:RnRn 而言, G G G R n R^n Rn 中有唯一的不动点。
  • 证明】任取 x 0 ∈ R n x_0∈R^n x0Rn,构造点列 x k + 1 = G ( x k ) , k = 0 , 1 , 2... x_{k+1}=G(x_k),k=0,1,2... xk+1=G(xk),k=0,1,2...,由于 G G G 是压缩映射,那么有 ∣ ∣ x k + 1 − x k ∣ ∣ = ∣ ∣ G ( x k ) − G ( x k − 1 ) ∣ ∣ ≤ α ∣ ∣ x k − x k − 1 ∣ ∣ ≤ . . . ≤ α k ∣ ∣ x 1 − x 0 ∣ ∣ ||x_{k+1}-x_k||=||G(x_k)-G(x_{k-1})||≤\alpha||x_k-x_{k-1}||≤...≤\alpha^k||x_1-x_0|| xk+1xk=G(xk)G(xk1)αxkxk1...αkx1x0.
  • 根据三角不等式可知 ∣ ∣ x k + p − x k ∣ ∣ = ∣ ∣ ( x k + p − x ) + ( x − x k ) ∣ ∣ ≤ ∣ ∣ x k + p − x ∣ ∣ + ∣ ∣ x − x k ∣ ∣ ||x_{k+p}-x_k||=||(x_{k+p}-x)+(x-x_k)||≤||x_{k+p}-x||+||x-x_k|| xk+pxk=(xk+px)+(xxk)xk+px+xxk,因此 ∣ ∣ x k + p − x k ∣ ∣ ≤ ∣ ∣ x k + p − x k + p − 1 ∣ ∣ + ∣ ∣ x k + p − 1 − x k + p − 2 ∣ ∣ + . . . + ∣ ∣ x k + 1 − x k ∣ ∣ ≤ ( α k + p − 1 + α k + p − 2 + . . . + α k ) ∣ ∣ x 1 − x 0 ∣ ∣ = α k ( 1 − α p ) 1 − α ∣ ∣ x 1 − x 0 ∣ ∣ ||x_{k+p}-x_k||≤||x_{k+p}-x_{k+p-1}||+||x_{k+p-1}-x_{k+p-2}||+...+||x_{k+1}-x_k||≤(\alpha^{k+p-1}+\alpha^{k+p-2}+...+\alpha^k)||x_1-x_0||=\frac{\alpha^k(1-\alpha^p)}{1-\alpha}||x_1-x_0|| xk+pxkxk+pxk+p1+xk+p1xk+p2+...+xk+1xk(αk+p1+αk+p2+...+αk)x1x0=1ααk(1αp)x1x0,由于 α ∈ ( 0 , 1 ) \alpha∈(0,1) α(0,1) ,所以 k → ∞ k→∞ k ∣ ∣ x k + p − x k ∣ ∣ → 0 ||x_{k+p}-x_k||→0 xk+pxk0,因此 x k {x_k} xk 是柯西序列。
  • 因为 R n R^n Rn 是完备的赋范线性空间,因此 x k x_k xk 的极限 x ∗ ∈ R n x^*∈R^n xRn,又因为 G G G 是连续映射,所以有 x ∗ = l i m k → ∞ ( x k + 1 ) = l i m k → ∞ [ G ( x k ) ] = l i m x → x ∗ G ( x ) = G ( x ∗ ) x^*=lim_{k→∞}(x_{k+1})=lim_{k→∞}[G(x_k)]=lim_{x→x^*}G(x)=G(x^*) x=limk(xk+1)=limk[G(xk)]=limxxG(x)=G(x),于是 G G G 的不动点存在性证明完成。
  • 关于唯一性,假设还有另外一不动点满足 y ∗ = G ( y ∗ ) y^*=G(y^*) y=G(y),所以有 ∣ ∣ x ∗ − y ∗ ∣ ∣ = ∣ ∣ G ( x ∗ ) − G ( y ∗ ) ∣ ∣ ≤ α ∣ ∣ x ∗ − y ∗ ∣ ∣ , α ∈ ( 0 , 1 ) ||x^*-y^*||=||G(x^*)-G(y^*)||≤\alpha||x^*-y^*||,\alpha∈(0,1) xy=G(x)G(y)αxy,α(0,1),显然只有 ∣ ∣ x ∗ − y ∗ ∣ ∣ ≡ 0 ||x^*-y^*||≡0 xy0 满足上式,根据范数正定性可知 x ∗ ≡ y ∗ x^*≡y^* xy,是同一个不动点。
  • 上述定理的条件过于苛刻,要求 G G G 在整个 R n R^n Rn 上为压缩映射,考虑将条件放松,简单地判断不动点是否存在。

预备概念.

  • 仅出于理解概念简述,非严格定义。

有界集.

  • 指那些可以被有界区间包含的实数集,通俗来说就是可以被长度有限区间包含的集合,例如区间 [ 2 , 5 ] [2,5] [2,5].
  • ∃ k ∃k k 使得对 ∀ s ∈ S ∀s∈S sS 都有 s ≤ k s≤k sk,那么 k k k 称为集合 S S S 的上界,类似地可以定义下界。
  • 具有上下界的集合 S S S 称为有界集。

闭集.

  • 如果一个集合不含边界点或没有边界点,那么它称为开集,例如区间 ( 2 , 5 ) (2,5) (2,5).
  • 补集为开集的集合是闭集,例如区间 [ 2 , 5 ] [2,5] [2,5],其实数域补集为 ( − ∞ , 2 ) ∪ ( 5 , + ∞ ) (-∞,2)∪(5,+∞) (2)(5,+).

凸集.

  • 以二维平面直观来看,凸集是一个具有这样性质的集合:集合中任意两点之间的连线都在集合内。
  • 使用数学语言表述为 ∀ x , y ∈ S , t ∈ ( 0 , 1 ) ∀x,y∈S,t∈(0,1) x,yS,t(0,1),都有 t x + ( 1 − t ) y ∈ S tx+(1-t)y∈S tx+(1t)yS 成立。

Brouwer不动点定理.

  • 设映射 G : D ⊂ R n → R n G:D⊂R^n→R^n G:DRnRn 在有界闭凸集 D D D 上连续,且 G ( D ) ⊂ D G(D)⊂D G(D)D,那么 G G G D D D 中必然存在至少一个不动点。
  • 从一维情况直观来看,取 D = [ a , b ] , G ( D ) = [ a , b ] ⊂ D . D=[a,b],G(D)=[a,b]⊂D. D=[a,b],G(D)=[a,b]D.
  • 假设不存在 x ∗ x^* x 满足 G ( x ∗ ) = x ∗ G(x^*)=x^* G(x)=x,即在区间 [ a , b ] [a,b] [a,b] 上,函数 h ( x ) = G ( x ) − x h(x)=G(x)-x h(x)=G(x)x 没有零点,所以 ∀ x ∈ [ a , b ] , h ( x ) > 0 ⋁ h ( x ) < 0 \forall x\in[a,b],h(x)>0\bigvee h(x)<0 x[a,b],h(x)>0h(x)<0,否则根据函数连续性必然存在零点。以 h ( x ) > 0 h(x)>0 h(x)>0 为例,当 x = b x=b x=b 时,必然有 G ( b ) > b G(b)>b G(b)>b,但 G ( D ) ⊂ [ a , b ] G(D)\subset[a,b] G(D)[a,b],因此 m a x G ( x ) ≤ b maxG(x)≤b maxG(x)b,矛盾出现; h ( x ) < 0 h(x)<0 h(x)<0 情况类似,因而至少存在一个不动点 x ∗ x^* x
    在这里插入图片描述

一叨.

  • 很多问题转化为不动点方程后,求解算法都是迭代法,例如比较熟悉的Hopcroft算法求解最简自动机,DV路由算法以及线性、非线性方程数值解的计算。
  • 关于压缩映射的定义式: ∀ x , y ∈ D , 0 < α < 1 , ∣ ∣ G ( x ) − G ( y ) ∣ ∣ ≤ α ∣ ∣ x − y ∣ ∣ ∀x,y∈D,0<\alpha<1,||G(x)-G(y)||≤\alpha||x-y|| x,yD,0<α<1,G(x)G(y)αxy由此想到Lipschitz条件,函数 f ( x ) f(x) f(x) x 0 x_0 x0 δ \delta δ 邻域内的所有点都满足 ∃   L > 0 , ∀ x ∈ [ x 0 − δ , x 0 + δ ] , ∣ f ( x ) − f ( x 0 ) ∣ ≤ L ⋅ ∣ x − x 0 ∣ \exist~L>0,\forall x\in[x_0-\delta,x_0+\delta],|f(x)-f(x_0)|≤L·|x-x_0|  L>0,x[x0δ,x0+δ],f(x)f(x0)Lxx0则称函数 f ( x ) f(x) f(x) x 0 x_0 x0 满足李普希茨条件。
  • 由此引申出的概念有李普希茨连续映射:设函数 f : D ⊂ X → Y f:D\subset X\rightarrow Y f:DXY,若有正常数 L L L 使得 ∀ x , y ∈ D \forall x,y\in D x,yD ∣ ∣ f ( x ) − f ( y ) ∣ ∣ ≤ L ⋅ ∣ ∣ x − y ∣ ∣ ||f(x)-f(y)||≤L·||x-y|| f(x)f(y)Lxy则称上述映射 f f f 为李普希茨连续映射,显然压缩映射必然是Lipschitz连续映射,是 L < 1 L<1 L<1 时的特殊情况。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值