计算方法非线性方程的数值解法.pptVIP

第二章 非线性方程的数值解法 单根和重根 有根区间 预备定理 迭代法的收敛与发散 迭代的计算步骤 迭代法计算框图的说明 2.2.3 迭代过程的收敛速度 收敛的情形 发散的情形 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 |x k-1- x k| 迭代法的终点判断：只要相邻两次迭代值的偏差充分小，就能保证迭代值足够准确，因而用 |x k- x k-1|控制迭代过程的结束。 定理 设在区间[a, b]上方程 x= φ (x) 有根x*,且对一切x∈[a, b] 都有| φ ′(x)| ≥ 1，则对于该区间上任意x0(≠x*), 迭代公式xk+1= φ (xk )一定发散。 证明 不可能收敛于0。 计算结果见下表 取方程的根 2.0946。 由于 ，故取 迭代法的局部收敛性 由于在实际应用中根 x* 事先不知道，故条件 | φ ′(x* )| 1 无法验证。但已知根的初值x0在根 x*邻域，又根据φ ′(x)的连续性，则可采用 | φ ′(x0 )| 1 来代替| φ ′(x* )| 1，判断迭代的收敛性。