chapter2非线性方程的数值解法.ppt

0 ( ) ( ) f x g x x ? ? ? ? ? ? § 4 牛顿法 /* Newton - Raphson Method */ 一、牛顿迭代公式的推导 1 、 待定参数法 不动点迭代的 关键 是构造满足 收敛条件 的 迭代函数 ( ) g x 一种 自然的选择 是令 0 ( ) ( ) ( ) g x x cf x c ? ? ? 为了加速不动点迭代的收敛过程，应尽可能使迭代函数 在 处有 更多阶导数等于零 （定理 2.3 ）。 ( ) g x x x ? ? 1 ( ) c f x ? ? ? ? ? 1 ( ) ( ) g x cf x ? ? ? ? ? ? 令 0 ? 1 1 0 ( ) ( ) ( ) ( ) ( ) ( ) ( ) g x h x f x h x f x h x f x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 0 ( ) ( ( ) ) ( ) h x f x f x ? ? ? ? ? ? ? ? ? 现设 ( ) ( ) ( ) g x x h x f x ? ? 1 0 ( ) ( ( ) ) ( ) h x f x f x ? ? ? ? ? 取 0 ( ) g x ? ? ? 满足 ? 因此，选取迭代函数 ( ) ( ) ( ) f x g x x f x ? ? ? Newton – Raphson 迭代格式 1 0 1 2 ( ) , , , ( ) k k k k f x x x k f x ? ? ? ? ? 称之为 牛顿 — 拉夫森 方法，简称 牛顿法 原理： 将非线性方程 线性化 取 x 0 ? x* ，将 f ( x ) 在 x 0 做一阶 Taylor 展开 : 2 0 0 0 0 ) ( ! 2 ) ( ) )( ( ) ( ) ( x x f x x x f x f x f ? ? ? ? ? ? ? ? ? ， ? 在 x 0 和 x 之间 2 、 Taylor 展开法 /* Taylors expansion Method */ 将 ( x* ? x 0 ) 2 看成 高阶小量 ，则有： ) * )( ( ) ( *) ( 0 0 0 0 x x x f x f x f ? ? ? ? ? ) ( ) ( * 0 0 0 x f x f x x ? ? ? ? x y x* x 0 1 0 1 2 ( ) , , , ( ) k k k k f x x x k f x ? ? ? ? ? 只要 f ? C 1 ，每一步迭代都有 而且 ，则 x * 就是 f 的根。 lim k k x x ? ?? ? 0 ( ) k f x ? ? 1 x ? 0 0 0 ( ) ( )( ) y f x f x x x ? ? ? ? 与 x 轴交点的横坐标 1 x 2 x 无开方运算，又无除法运算。 例 1 ： ? 写出求 的 Newton 迭代格式； ? 写出求 的 Newton 迭代格式 , 要求公式中既 0 ( ) a a ? 1 0 ( ) a a ? 解： ? 等价于求方程 的正根 2 0 0 ( ) ( ) f x x a a ? ? ? ? 2 1 1 0 1 2 2 2 ( ) ( ) , , , ( ) k k k k k k k k k f x x a a x x x x k f x x x ? ? ? ? ? ? ? ? ? ? 2 ( ) f x x ? ? ? 解法一： 等 价于求方程 的正根 2 1 0 0 ( ) ( ) ( ) f x x a a ? ? ? ? 1 2 ( ) ( ) f x x a ? ? ? 2 1 1 1 2 ( ) ( ) ( ) ( ) k k k k k k k x f x a x x x f x x a ? ? ? ? ? ? ? ? 1 1 0 1 2 2 ( ) , , , k x k a ? ? ? ? 3 2 ( ) f x x ? ? 解法二： 等价于求方程