数值分析7-4,5(牛顿法,弦截法).pptVIP

7.4 -7.5 牛顿法及其推广 /* Newton Method */ 一、牛顿迭代法的公式 二、牛顿迭代法的改进与推广 原理：将非线性方程线性化——泰勒展开 /* Taylor’s expansion */ 取 x0 ? x* ，将 f(x*) 在 x0 做一阶泰勒展开: ? 在 x0 和 x*之间。 将 (x* ? x0)2看成高阶小量，则有： 一、牛顿迭代法的公式 线性 /* linear */ x y x* x0 只要 每一步迭代都有f ’( xk ) ? 0， 而且 ，则 x*就是 f 的根。 牛顿迭代法的基本思想 将非线性方程 f(x)=0 的求根问题归结为计算一系列线性方程的求根问题。 牛顿迭代法的计算步骤 (1)给出初始近似根 x0 及精度ε； (3)若 ，转向(4)， 否则 ，转向(2); (4)输出满足精度的根 x1 ，结束。 (2)计算 例 用牛顿迭代法求方程 在 x=0.5 附近的根。取 解 其牛顿迭代公式为 取初值 x0=0.5 ，迭代结果见下表 易见 故 0.5 0.57102 0.56716 0.56714 xk 0 1 2 3 k k 0 1 2 3 xk 0.880000 0.884688 0.884675 0.884675 例 2 计算 的近似值, ?=10-6 x0=0.88 解：令 x= 问题转化为求f(x)=x2-0.78265=0 的正根 由牛顿迭代公式 xk+1= xk-?(xk)/?(xk)= xk/2+0.78265/2xk 迭代结果 满足了精度要求,故 ≈0.884675 设 f ?C2[a, b]，若 x* 为 f (x) 在[a, b]上的根，且 f ’(x*) ? 0，则存在 x* 的邻域 定理 （局部收敛性） Newton’s Method产生的序列{ xk } 收敛到x*，且满足 使得任取初值 ? Newton’s Method 有 ，只要 就有 p ? 2。重根是线性收敛的。 证明： Newton’s Method 事实上是一种特殊的不动点迭代，其中 收敛 则 ? 由泰勒展开： 在单根 /*simple root */ 附近收敛快 只要 f ’(x*) ? 0，则令 可得结论。 注： Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 注 (1) 牛顿法要求初值充分接近根以保证局部收敛性。 (2)牛顿迭代法的主要优点是收敛较快，是平方收敛的缺点是公式中需要求 f(x) 的导数。若 f(x)比较复杂，则使用牛顿公式就大为不便。。 ? 重根 /* multiple root */ 加速收敛法： 问题1: 若　　　 ，Newton’s Method 是否仍收敛？ 设 x*是 f 的 n 重根，则： 因为 Newton’s Method 事实上是一种特殊的不动点迭， 其中 二、牛顿迭代法的改进与推广 /* improvement and generalization */ 且 K1: 有局部收敛性，但重数 n 越高，收敛越慢。 则 问题2: 如何加速重根的收敛？ K2: 将求 f 的重根转化为求另一函数的单根。 ? 令　　　　　 则 f 的重根 = ? 的单根。 ? 下山法 /* Descent Method */ ——Newton’s Method 局部微调： 原理：若由 xk 得到的 xk+1 不能使 | f | 减小，则在 xk 和 xk+1 之间找一个更好的点 ，使得 xk xk+1 注： ? = 1 时就是Newton’s Method 公式。 当 ? = 1 代入效果不好时，将 ? 减半计算。 。 ? 弦截法 /* Secant Method */ Newton’s Method 一步要计算 f 和 f ’，相当于2个函数值，比较费时。现用 f 的值近似 f ’，可少算一个函数值。 x0 x1 切线 /* tangent line */ 割线 /*