数值计算方法(第2章) 深圳大学 科学与工程计算 数值分析 课件.ppt

迭代法收敛的阶 证明： 满足压缩映像原理 迭代法收敛的阶 敛速是线性的 线性收敛到 。 Steffensen迭代格式 由线性收敛知 当n充分大时有 即 Steffensen迭代格式 展开有： Steffensen迭代格式 已知 ，则 ， 改成 n=0，1，2，… Steffensen迭代格式 也可以改写成 其中迭代函数 Steffensen迭代法收敛的充要条件 定理2.2.3 Steffensen迭代法收敛的充要条件 证明：必要性 Steffensen迭代法收敛的充要条件 充分性 Steffensen算法的收敛速度 Steffensen算法的收敛速度 定理2.2.5 在定理2.2.3假设下，若 产生的序列 至少平方收敛到 。 Steffensen算法的收敛速度 Steffensen算法的收敛速度 Steffensen算法的收敛速度 Steffensen算法的收敛速度 由定理2.2.4知 至少以平方速度收敛到 。 也就是说：简单迭代法是线性收敛；Steffensen迭代至少平方以上收敛（加速收敛）。 例题 例2.2.3试用Steffensen算法求解方程 解法一、取 ，由 n = 0，1，2，… 例题 取初值 ，计算结果如下： N Xn Yn Zn 0 1.5 1.357208808 1.330860959 1 1.324899181 1.324752379 1.324724496 2 1.324717957 1.324717957 1.324717957 例题 解法二、取 ，由 对于该迭代函数在一般迭代法中是发散的，而Steffensen格式却是收敛的。 n=0，1，2，… 例题 取初值 ，计算结果如下： N Xn Yn Zn 0 1.5 2.375 1.239648437 1 1.416292975 1.840921915 5.238872769 2 1.355650442 1.491398279 2.317270699 3 1.328948777 1.347062883 1.444351224 4 1.324804489 1.325173544 1.327117281 5 1.324717944 1.324718152 1.324718980 6 1.324717957 Steffensen迭代格式几何解释 Steffensen迭代算法 Steffensen迭代算法 2.3 Newton迭代法 设x * 是方程f (x ) = 0的根，又x0 为x * 附近的一个值 ，将f (x ) 在x0附近做泰勒展式 令 ，则 Newton迭代法 去掉 的二次项，有： 即 以x1代替x0重复以上的过程，继续下去得： Newton迭代法 以此产生的序列{Xn}得到f(x)=0的近似 解，称为Newton法，又叫切线法。 Newton迭代法几何解释 几何意义 例题 例2.3.1 用Newton法求 的近似解。 解：由零点定理。 等步长扫描法求有根区间 用计算机求有根区间：等步长扫描法。 设h0是给定的步长，取 ， 若 则扫描成功；否则令 ，继续上述方法，直到成 功。如果 则扫描失败。再将h 缩小， 继续以上步骤。 等步长扫描算法 算法：（求方程 的有根区间） (1) 输入 ； (2) ； (3) ，若 输出失败信息，停机。 (4)若 。输出 ，已算出方程的一个根，停机。 等步长扫描算法 (5) 若 。输出 为有根区间，停机 (6) ，转 3） 注：如果对足够小的步长h扫描失败。说明： 在 内无根 在 内有偶重根 二分法 用二分法(将