数值分析10-方程求根的迭代法.ppt

一点注记 Newton迭代 基本思想： 将非线性方程线性化 设 xk 是 f (x)=0 的近似根， 将 f (x) 在 xk 一阶 Taylor 展开: ，? 在 xk 和 x 之间。 x y x* xk xk+1 条件： f’(x) ? 0 Newton迭代 Newton 法可以看作下面的不动点迭代： 其中 ?’(x*) = 0 Newton 法至少 二阶 局部收敛 定理 设 f(x) 在其零点 x* 的某个邻域内二阶连续可导且 f’(x) ? 0，则存在 x* 的某个 ? 邻域 N(x*) =[x*-? , x* +? ], 使得对 ? x0? N(x*)，Newton 法产生的序列以不低于二阶的收敛速度收敛到 x* 。 Newton迭代 Newton 法也可以看作一类特殊的加速迭代 取 ?(x) = x+f(x) 收敛性定理 定理 设 f ?C2[a, b]，且 f 满足 (1) f (a) f (b) 0 (2) 对? x?[a, b]， f’(x) ? 0 且 f”(x) ? 0 ； (3) 初始点 x0 ? [a, b] 满足 f (x0) f”(x0) 0； 则 Newton 法产生的序列收敛到 f 在 [a, b] 的唯一零点 x*。 全局收敛性定理 定理 设 f ?C2[a, b]，且 f 满足 (1) f (a) f (b) 0 (2) 对? x?[a, b]， f’(x) ? 0 且 f”(x) ? 0 ； 则对任意初始点 x0 ? [a, b]，Newton 法产生的序列收敛到 f 在 [a, b] 的唯一零点 x*。 (3) 举例（一） 例：设计一个二阶收敛算法计算 (a 0)。 解：转化为求 x2-a = 0 的正根 Newton 迭代： 二阶收敛 重根情形 设 x* 是 f(x) 的 m(m?2) 重根，Newton法是否收敛？ Taylor 展式 Newton 迭代： 线性收敛。 且重数 m 越高，收敛越慢。 提高收敛阶 提高收敛速度 但 m 通常无法预先知道！ 法一：取 二阶收敛 法二：将求 f(x) 的重根转化为求 另一个函数 的单根。 构造针对 ?(x) 的具有二阶收敛的 Newton 迭代： 令 ，则 x* 是 ?(x) 的单重根。 降低初始点的要求 例：求 sin(x)-x/6=0 的正根。 Newton 下山法： ?k ?k 为数列 中满足 的最大数。 算法 7.2 (Newton下山法) 给定初始点 x0，精度要求 ? 1. 如果 | f (xk)| ? , 停机，输出 xk 2. 计算 ，? =1 3. 如果| f (xk+?dk)| | f(xk)|，令 xk+1= xk+?dk ，返回第1步； 否则? 折半，重新计算第3步 Newton 法的收敛依赖于初始点的选取。 割线法 Newton法的缺点：每步迭代都要计算导数值 只需计算函数值，避免计算导数； 切线斜率 ? 割线斜率 xk-1 xk xk+1 xk+1 切线 割线 需要两个初始点； 收敛比Newton法稍慢，但对初始点要求同样高。 割线法公式 两点割线法 单点割线法 误差估计 引理 设 f(x) 在其零点 x* 的某个邻域 N(x*) =[x*-? , x* +? ] 内存在连续的二阶导数，且 f’(x) ? 0，又设 xk-1 , xk ? N(x*) \{x*}，且互不相等，则由两点割线法的误差 ek=xk-x* 满足 其中 局部收敛性定理 定理 设 x* 是 f(x) 的单重零点，f”(x) 在 x* 的某个邻域内连续，则存在 x* 的一个邻域 N(x*) =[x*-? , x* +? ]，使得当 x0 , x1 ? N(x*) 时，由两点割线法产生的序列收敛到 x*，且收敛阶为 本章结束 谢谢！ * 第四章 方程求根的迭代法 远在公元前1700年的古巴比伦人就已有关于一、二次方程的解法。《九章算术》(公元前50~100年)其中“方程术”有联立一次方程组的一般解法。 1535年意大利数学家坦特格里亚(TorTaglia)发现了三次方程的解法，卡当(H·Cardano)从他那里得到了这种解法，于1545年在其名著《大法》中公布了三次方程的公式解，称为卡当算法。 后来卡当的学生弗瑞里(Ferrari)又提出了四次方程的解法。此成果更激发了数学家们的情绪，但在以后的二个世纪中，求索