计算方法第二章分解.pptVIP

则 从而 由于f(x)连续，因此，存在 r 0 ，使得对一切 有 且 　位于 x 与 p 之间.由于 ，因此 g (x) 据 §3 定理1 推论知初始值　　　　　　　　，Newton法收敛．由于 §4. Newton-Raphson方法 31 故 Newton 法仅为线性收敛． 当 p 是 f(x) =0 的 q 重根时，为提高迭代法的收敛速度，我们作如下考虑． 将 f(x) 和 在点 p 按 Taylor公式展开，得到 其中， 和 位于 x 与 p 之间，因此 令 因为 ，因此当 x 充分接近于 p 时， 均不为 零．据（4．4）和（4．3）知，p 是方程 f(x) =0 的单根.以 F(x) 代替Newton §4. Newton-Raphson方法 32 法中的 f(x)，便得到迭代公式 如果迭代函数 具有即需的连续性条件，那么无论根重数多少，这个方法至少是二阶收的．理 论上，这个方法的缺点是增加二阶导数 的计算以及迭代过程中的计算 更复杂．但是，实际上重根的出现将产生严重的舍入误差问题． 例2 在例1中，方程 在区间（1，2）中有一个单根 p ．取初始值 ，应用Newton法迭代四次 §4. Newton-Raphson方法 33 得 ．现在改用迭代公式（4．5）来计算它．据（4．5）式 取 ，得 例3 方程 有一个二重根十了 ．我们应用 Newton法和（4．5）来算它 的近似值．现在，Newton迭代公式（4．2）及其修改公式（4．5）分别为 §4. Newton-Raphson方法 34 取初始值 ，用 迭代三次得结果如下： 由此看出，在 中 准确到 ．若用 Newton法（4．2）来计算，要 达到这个精确度则需迭代20次。 　 §4. Newton-Raphson方法 35 　　Newton法有明显的几何意义（见图2.1）。函数方程 f(x) =0 的根 p 是曲线 y=f(x) 与 X 轴的交点的横坐标。过点　　　　　　　作曲线的切线，则切线 方程为 切线　　与 X 轴的交点的横坐标为 a p 图　２.１ 因此，Newton法又叫做切线法． §4. Newton-Raphson方法 36 前面我们指出了，若 f(x) 具有足够阶连续导数，且选取初始近似值 充分接近于方程 f(x)＝0的根p，则Newton迭代序列 收敛于 p ．在实际应用 中，有的实际问题本身可以提供接近于根的初始值，但有的问题却难以 确定接近于根的初始值．关于初始值的选取，我们有下面的定理． 定理2 设函数 f(x) 在有限区间[a，b]上存在二阶导数，且满足条件： 则　Newton法（４.２）对任意的初始值 ，都收敛于方程f(x)=0 的唯一解 p，且收敛阶数为２。 §4. Newton-Raphson方法 37 在定理2中，条件（ ）保证方程 f(x) ＝0 在(a，b)内至少有一个根．条件（ ） 表明函数 f(x) 不是严格单调增大（ ）就是严格单调减小（ ），因 而 f(x) = 0在(a，b)内有唯一根．条件（ ）说明 f 的图形不是凹向上 就是凹向下 ．条件（ ）保证，当 时，Newton序列 在（a，b）中．例如，从图2．2看到，取