2.3.1Newton迭代方法的计算公式.docxVIP

PAGE PAGE # 3 牛顿( Newton)迭代方法 Newton 迭代方法的计算公式 牛顿迭代法计算公式的推导过程 本节所讨论的是： f (x) 0 。 设 x*是 f (x) 0的根， f (x)在 x*的邻域内 具有二阶连续导数，在 x* 的邻域内取一点 x0， 使 f (x0 ) 0，将它在 x0点二阶 Taylor 展开 得： f(x) f(x0) f (x0)(x x0) f 2(! )(x x0)2f (x0) f (x0)(x x0) 又 f (x) 0 ，则有： f (x0 ) f (x0)(x x0) 0 f (x) 0的近似解 f (x) 0的近似解 x x0 f (x0)f (x0) 记 x1 x f(x0 ) x0 0 f (x0) 类似，在点 x1处 Taylor 展开，可得： x x1f (x1)f (x x x1 f (x1)f (x1) x2 x1 f(x1)f (x1) 依次往下做，可得一般的迭代格式： xn 1 x xn 1 xn f (xn ), f (xn), (n 0,1, ) 上述迭代格式称为求 f (x) 0的解的牛顿迭 代法。 几何意义 在点 (x0, f (x0)) 处作 f ( x)的切线，交 x 轴于一点，求该点的横坐标。此切线方程为： y f (x0) f (x0)(x x0)，当 y 0时， f(x0 ) 得： x x0 0 ，正是 x1的值。 0 f (x0) 1 依次类推， 在点 (xn, f (xn))作函数 f (x) 的 切 线 ， 交 x 轴 于 一 点 ， 切 线 方 程 为 ： y f (xn) f (xn)(x xn)，当 y 0时， 得： x xn f(xn) ，正是 xn 1的值 n f (xn) n 1 牛顿迭代法又称为切线求根法。 迭代法收敛的条件与收敛速度(针对单根而 言) 定理：设 f(x*) 0,f (x*) 0，且 f (x)在 x*的 邻域内具有二阶连续导数， 则由牛顿迭代法产 生的迭代序列 xn 1 xn f (xn) , (n 0,1, )f (xn) 局部收敛于 x* ，且为平方收敛。 证明：在牛顿迭代法的迭代格式中，迭代函数 为： (x)f(x)f (x) (x) f (x)在 x* 的邻域内具有二阶连续导 数，即 (x) 1( f (x))2 f (x) f(x)(f (x))2 (x) 1 f (x) f (x) (f (x))2 又 f(x*) 0， (x*) 0 1， 牛顿迭代法局部收敛于 x* (由定理 2) 又 (x)f (x) 2 f (x)f (x) f(x)f (x) 2f(x) f (x) 2 f (x) (x) 但是，充分的程度没有具体的描述， 而且，若 x0 的值没有取好，有可能得不到收敛的结果。 牛顿迭代法收敛的一个保证条件。 补充定理：设 f ( x)在[ a, b]区间上的二阶 导数存在，且满足： ① f(a)f(b) 0 ； ② x [a,b], f ( x)不变号； ③ x [a,b], f ( x)保持符号不变。 ④初始值 x0 [a,b], f (x0)f (x0) 0， 则牛顿迭代法产生的迭代序列 xn 收敛于 f (x) 0在[a,b]区间的唯一根。 证明：由①②知方程 f(x) 0在区间有且只有一个 根，记为 x* 。 不失一般性，设 f (a) 0,f(b) 0,f (x) 0,f (x) 0. 其他情形可类似证明。按④应取 x0 [ a,b], f(x0) 0 由 f (x) 0知 f (x) 为 单 调 增 函 数 ， 从 而 知 x以 x0 x 以 x0 为初值，迭代一次 x0 x1 x0 f (x0) x0 f (x0) 0 另一方面，将 f(x*)在 x0处作泰勒展开，得 f(x* ) f(x0) f (x0 )( x* x0) 1 * 2 2 f ( 0)(x* x0)2 0 其中 0介于 x*和 x0之间。将上式两边除以 f (x0)， 得 f ( 0)(x* x0 f ( 0)(x* x0)2 0 2f (x0) 0 ff((xx00)) (x x0) f (x0) 移项得 因而x0f (x0)f (x0)* 因而 x0 f (x0) f (x0) * x x1 x1 f ( 0)(x* x0)2 0 2f (x0) 0 x0 般的，若 x* xk ，同理可证 x xk 1 xk x。这就说明 xk单调下降有下界 x* ，因此必收敛。 设 lim xk x x。 kk 再对迭代格式 xk 1 再对迭代格式 xk 1 f (xk ) 两边取极限， f (xk ) f(x) f (x) 即 x为 f (x) 即 x为 f (x