大连理工大学矩阵与数值分析第4章-2非线性方程迭代解法2.ppt

Newton 迭代法的几何意义 ---切线法 Newton迭代法 由于Newton迭代法带有 的导数 , 使用 起来不太方便。 为了不求导数， 可用导数的近似式替代 ： 将它代入（4-26）代替 ,得 则 （4-27） （4-27）就是弦截法。 由于弦截法采用了导数的近 似值， 故在Newton法和弦截法都收敛的情况下， 弦截法 的收敛阶为 , 低于Newton法， 为超线 性收敛。 弦截法的几何意义 例4 用Newton法和弦截法分别计算方程 在 附近的根 。 解 （1） 使用Newton法， 并取 （4-28） 迭代3次就得到具有6位有效数字的结果。 * * 4.2 非线性方程的迭代解法 工程实际与科学计算中都遇到大量求解非线性 若数 (4-17) ， ， 为方程（4-17） 的零点。 的根， 使 则称 或称函数 常见的非线性方程有代数方程（二次、三次等）、 超越方程（三角方程，指数、对数方程等）。 方程的问题。 对于非线性方程 困难 即使是最基本的代数方程，当次数超过4时，在一般情况下就不能用公式表示方程的根，即难于用解析法求出方程的根，对于超越方程就更难了。 因此，研究用数值方法计算非线性方程的根就显 得非常必要。 在求根时通常假设非线性方程 中的函数 是关于 的连续函数。 若令 则它在平面直角坐标系 下的图像为连续曲线， 可见,求 的根， 如果 在区间 上仅有一个根，则称 为方程的单根区间； 为方程的多根区间。 上有多个根，则称 若方程在 与 轴的交点 就是求 y x o 方程的单根区间和多根区间统称为方程的有根区间。 为了研究方便，我们主要研究方程在单根区间上的求解方法。 （4-18） 4.2.1 简单迭代法 首先将方程 化为一个与它同解的方程 为 其中 的连续函数。 即如果数 使 反之， ， 。 则也有 ， 则也有 ， 若 称迭代法或迭代过程或迭代格式， 通常称（4-19）为求解非线性方程的简单迭代法， 代入（4-18）的右端，得到 再将 代入（4-18）右端得 ， 继为之，得到一个数列， （4-19） 称为迭代函数， 称第 步的迭代值或简称迭代值。 其一般表示形式为 也 ， 任取一个初始值 则称迭代法收敛， 如果由迭代格式产生的数列收敛，即 否则称迭代法发散。 若迭代法收敛于 则 所以 。 即 几何直观： 在曲线 上得到点列 ， 其横坐标分别 为由公式 所确定的迭代值 ， 若迭代法收敛 ， 则点列 将越来越逼近所求的交点 。 p1 p2 p* y x o p0 y=x （1）化方程为等价方程 则迭代值为 取初始值 用迭代法求 的根。 ， ， ， ，… 显然, 当 时， ,即迭代法收敛于1. 就是方程 的根。 例1 解 则迭代格式为 ,所以 （2） 为等价方程 化 ， 同样取初始值 ， 其迭代格式为 ， ， ， 显然, 当 ,故迭代法发散。 上述例子表明，迭代法的收敛与发散，依赖于迭代函数的构造，迭代函数构造的方法很多。 对于同一个方程，由于构造出来的迭代函数不同，有的迭代函数所构成的迭代法收敛，有的迭代函数所构成的迭代法却发散。 迭代函数满足什么条件时，迭代法收敛？ 从而,迭代函数满足条件： 时,迭代法收敛。 从而,当 或 时,迭代法发散。 定理4.5 设迭代函数 满足 （1） （2） ， 存在正数 对任意 均有 则 在 内存在唯一根 , 且对任意初始值 , 迭代法 收敛于 , 且 1． 2． （4-20） （4-21） 当 时， 满足条件（1）、（2）时， 易证方程 在[a,b]内存在唯一根 。 因为 ， 且 ， 根据微分中值定理可得 其中 。由条件（2）得 （4-22） 证 又因为 将上式移项整理后，得 从而 即（4-20）成立。 再反复使用（4-22）的第2式，得 将上式代入（4-20）即得（4-21）成立。 又因为L＜1， 所以根据（4-21）得 故迭代法收敛。 当迭代函数满足定理4.5的条件且 较小时， 根据 （4-20）式可知， 只要相邻两次计算值的偏差 达到事先给定的精度要求 （即 ）时， 过程就可以终止， 迭代 就可作为 的近似值。 因此， （4-20）式也是判断迭代是否可终止的依据。 如果对 的大小可作出估计时， 由（4-21）式就可以大概估计 出迭代过程所需要的迭代次数， 即 时， 的 大小范围。 由于定理4.5的条件一般难于验证， 而且在大区间 上， 这些条件也不一定都成立， 所以在使用迭代法 时往往在根 的附