数值计算方法 第2章节 第二节.pptVIP

预备定理 迭代法的收敛与发散 迭代的计算步骤 2.2.3 迭代过程的收敛速度 * 2.2 迭代法 2.2.1 迭代原理 2.2.2 迭代的收敛性 2.2.3 迭代的收敛速度 2.2.4 迭代的加速 2.2.1 迭代原理 计算结果见下表 方程f(x)=0化为等价形式的方程 x=φ(x) ， 构造迭代公式 xk+1= φ(xk ) , k=0,1,2,…… 取初始近似根x0 ,进行迭代计算 x1= φ (x0), x2= φ(x1),…….. 则有x1 , x2 ,, ……. , xk , …….,得到迭代序列{xk }. 如果这个序列有极限，则迭代公式是收敛的。这时 则 ， x* 为不动点，等价地有 f(x*)=0, x* 即为方程的根。连续函数φ (x)称为迭代函数。 实际计算到 |xk – xk-1 |ε（ε是预定的精度），取x*≈xk 。 迭代公式收敛指迭代序列{xk }收敛，迭代公式发散指迭代序列{xk }不收敛，即发散。迭代公式不一定总是收敛。 例如求方程 f(x)=x3-x-1=0的一个根。 对应的迭代公式为 取初值 迭代序列{xk }发散. x1= φ (x0 ) x2= φ (x1 ) 迭代法收敛与发散的图示 收敛的情形 发散的情形 2.2.2 迭代的收敛性 迭代法的收敛条件及误差估计式 定理（充分性条件） 设函数 φ (x) 在[a, b]上连续，且 (1)对 x∈[a, b],有 φ (x) ∈[a, b] (2) 存在0 L 1，使对任意 x ∈[a, b]有 | φ ′(x)| ≤ L 1 则方程x= φ (x)在 [a, b]上的根x*存在且唯一；对初值 x0 ∈[a, b] ,迭代过程 xk+1= φ (xk)均收敛于方程的根x*。 定理中的 (1)对 x∈[a, b],有 φ (x) ∈[a, b]，称为适定性（映内性）。 证明 先证根的存在性。 作连续函数ψ(x)=x-φ(x),由条件（1）x∈[a, b], φ (x) ∈[a, b]，即a≤ φ (x)、x≤b ,于是 ψ(a)=a- φ (a) ≤0 ψ(b)=b- φ(b) ≥0 由于ψ (x)是连续函数，故必存在 x* ∈[a, b] 使ψ (x*)=0.即 ψ (x*)= x*-φ (x*)=0. 于是 x*= φ (x*)即x*为方程 x= φ (x)的根。 其次，证根的唯一性 。 设y*也是方程的根，则x*= φ (x*), y*= φ (y*), x*- y*= φ (x*) –φ (y*)= φ ′(ξ)(x*- y*) x*- y*–φ ′(ξ)(x*- y*)=0， （x*- y*）[1- φ ′(ξ)]=0 由条件（2）| φ ′(x)| ≤ L 1，故有x*-y*=0，即x*=y* 所以方程在[a,b]的根唯一。 再证迭代的收敛性 。 由xk= φ (xk-1 ) , x*= φ (x*), 有 |xk - x*| = | φ ′(ξ)(xk-1-x*)|≤L |x k-1- x*| ≤L2|xk-2-x*| ≤L3| xk-3-x*|≤…… ≤Lk |x0-x*| →0 （k→∞） 所以，对[a,b]上任取的x0 ,迭代公式xk+1= φ (xk )都收敛于x*。 L越小收敛得越快。 定理是充分性条件 xk - x* = φ (xk-1) – φ(x*)= φ ′(ξ)(xk-1-x*) 推论：在定理的条件下，有误差估计式 验后误差估计式 验前误差估计式 证明： |xk - x*| ≤L |x k-1- x*|= L |x k-1- x k + x k -x* | ≤L （|x k- x*|+ |x k-1- x k |） （1-L) |x k- x*| ≤L