第五章迭代法52Newton法弦截法讲解.ppt

第五章 方程求根的迭代法 Newton迭代 Newton迭代 举例（一） 重根情形 提高收敛阶 降低初始点的要求 第五章 方程求根的迭代法 快速弦截法 弦截法 局部收敛性定理 * * Newton迭代法 基本思想： 将非线性方程线性化 设 xk 是 f (x)=0 的近似根， 将 f (x) 在 xk Taylor 展开: , ? 在 xk 和 x 之间. x y x* xk xk+1 x* 是单根时： f’(x* ) ? 0 Newton 法可以看作下面的不动点迭代： 其中 ?’(x*) = 0 Newton 法至少 二阶 局部收敛 定理 设 f(x) 在其零点 x* 的某个邻域内二阶连续可导且 x*是单根 ，则存在 x* 的某个 ? 邻域 N(x*) =[x*-? , x* +? ], 使得对 ? x0? N(x*)，Newton 法产生的序列以不低于二阶的收敛速度收敛到 x* 。 (局部收敛定理) 例：只用加减乘除运算，设计一个二阶收敛算法计算 (a 0)。 解：转化为求 x2-a = 0 的正根 Newton 迭代： 二阶收敛 设 x* 是 f(x) 的 m(m?2) 重根，Newton法是否收敛？ Taylor 展式 Newton 迭代： 线性收敛。 且重数 m 越高，收敛越慢。 Taylor 展式 Taylor 展式 提高收敛速度 但 m 通常无法预先知道！ 法一：取 二阶收敛 法二：将求 f(x) 的重根转化为求 另一个函数 的单根。 构造针对 ?(x) 的具有二阶收敛的 Newton 迭代： 令 ，则 x* 是 ?(x) 的单重根。 例7：求 x3-x-1=0 在 x = 1.5 附近的正根。(Ch4exp7.m) Newton 下山法： ?k ?k 为数列 中满足 的最大数。 算法 (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法稍慢，但对初始点要求同样高。 两步弦截法 (单步)弦截法 快速弦截法的特点: 不用计算导数; 计算xk+1 时, 要用到 xk, xk-1 . 也称为两步弦截法. 线性收敛 超线性收敛 定理 设 x* 是 f(x) 的单重零点，f”(x) 在 x* 的某个邻域内连续，则存在 x* 的一个邻域 N(x*) =[x*-? , x* +? ]，使得当 x0 , x1 ? N(x*) 时，由两步弦截法产生的序列收敛到 x*，且收敛阶为 超线性收敛 * * *