数值计算方法第2版 第2章 非线性方程的数值解法.ppt

单根和重根 有根区间 预备定理 迭代法的收敛与发散 迭代的计算步骤 迭代法计算框图的说明 2.2.3 迭代过程的收敛速度 2.3.3 牛顿迭代法的修正 2. 牛顿下山法 2.3 牛顿迭代法 2.3.1 迭代公式的建立 2.3.2 牛顿迭代法的收敛情况 2.3.3 牛顿迭代法的修正法 2.3.1 迭代公式的建立 3.几何意义 过曲线上的点pk(xk , f(xk))作切线，切线方程 y=f(xk)+f ?(xk)(x – xk) 切线方程和横轴的交点(xk+1 , 0) ，即 0= f(xk)+f ?(xk)(xk+1 – xk) 若 f ?(xk )≠0，解出xk+1，则得Newton迭代公式 例 用牛顿迭代法求方程 xex-1=0 在x=0.5附近的根。 解 牛顿迭代法 取x0=0.5，经计算可得 普通迭代法18次才能得到的计算结果。 ,则x2-a=0, 求 等价于求方程 令 例 造平方根表。用牛顿迭代法计算 ( 其中a0) 解 的正实根。因为 f′(x)=2x ，由牛顿迭代公式得 当a=115时，取初值 x0=10，迭代4次可得 10，10.7500，10.723837，10.723805， 10.723805 ≈10.723805 是否还能用牛顿法计算一个正数的立方根？ ,则x3-3=0, 求 等价于求方程 令 例 用牛顿迭代法求 解 的正实根。由牛顿迭代公式得 当a=4111.7910时，取初值 x0=8，迭代4次可得 7.48，7.439977, 7.439760, 7.439760 , 令 例 用牛顿迭代法造倒数表，计算 解 3、牛顿迭代法的计算步骤 (1)给出x0 , ε，N (2)计算 (3)若 则转(4)；否则 ，转(2); (4)输出x1 ,结束。 牛顿迭代法局部收敛: 2.4.2 牛顿迭代法的收敛情况 1.局部收敛性 结论： f(x)=0的单根x* 附近存在着连续的二阶导数，当初值在单根x*附近时，牛顿法具有平方收敛速度。 证 牛顿法迭代函数 当 f(x*)=0，而 f?(x*)≠0 ，则x*是f(x)=0的单根,单根x*附近存在着连续的二阶导数，有 牛顿法至少具有平方收敛速度。 1. 简化牛顿法(平行弦法) 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*|≤