方程求根的迭代法.ppt

第四章 方程求根的迭代法 高 云 方程求根需要考虑的问题 迭代法的基本思想 转换例子 几何含义 几何含义 压缩映像定理 压缩映像定理证明 压缩映像定理证明 全局收敛与局部收敛 迭代过程的收敛速度 P阶收敛 必要性的证明 迭代过程的加速 埃特金算法 一点注记 Newton迭代 Newton迭代 Newton迭代 收敛性定理 全局收敛性定理 举例（一） 重根情形 提高收敛阶 降低初始点的要求 割线法 割线法公式 误差估计 局部收敛性定理 本章结束 谢谢！ * * 求 f (x) = 0 的根 代数方程： f (x) = a0 + a1x + . . . + anxn 超越方程： f (x) 含超越函数，如 sin(x), ex, lnx 等 实根与复根 根的重数 f (x) = ( x – x*)m · g(x) 且 g(x*) ? 0, 则 x* 为 f (x) 的 m 重根 有根区间：[a, b] 上存在 f (x) = 0 的一个实根 在有根的前提下求出方程的近似根。 研究 内容： 基本思路 同解 迭代公式 给定初值 存在 等价于 迭代函数 ？转换是否唯一 几何 意义 (1) x = ?1(x) = x3-6x2+10x-2 ; (2) ; (3) ; (4) ; 例： 已知方程 x3-6x2+9x-2=0 在 [3,4] 内有一根，考虑迭代 ？哪种转换方法好 x y y = x x y y = x x* x* y=g(x) y=g(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x y y = x x y y = x x* x* y= ?(x) y=?(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? 定理 设 ?(x)?C[a, b] 且可导，若 (2) ? 0 ? L 1，使得 | ?’(x) | ? L 对 ? x?[a, b] 成立 (1) a ? ?(x) ? b 对一切 x?[a, b] 都成立 则有 (a) 对任意 x0?[a, b]，由 xk+1 = ?(xk) 产生的迭代序列 均收敛到 ? (x) 在 [a, b] 中的唯一不动点 x*。 (b) 有如下的误差估计 可用| x k+1-xk | 来控制收敛精度 L 越小收敛越快 (a) 由压缩映像定理可知，不动点 x* 存在且唯一。 (b) 又 定理的条件保证了不动点迭代的全局收敛性。 即迭代的收敛性与初始点的选取无关。 这种在 x* 的邻域内具有的收敛性称为局部收敛性。 定理中的条件 | ?’(x) | ? L 1 可以适当放宽 (2’) ?’(x) 在 x* 的某个邻域内连续，且 | ?’(x*) | 1 由 ?’(x) 的连续性及 | ?’(x*) | 1 即可推出： 存在 x* 的某个 ? 邻域 N(x*) =[x*-? , x* +? ], 使得对 ? x? N(x*)都有| ?’(x) | ? L 1, 则由 ?x0? N(x*) 开始的迭代都收敛。 定义 则称该迭代为 r 阶收敛。 (1) 当 r =1 时称为线性收敛，此时 C 1； (2) 当 r =2 时称为二次收敛，或平方收敛； (3) 当 r 1 时称为超线性收敛。 二分法线性收敛 不动点迭代中，若 ?’(x*) ? 0，则 取极限得 (C为常数) 线性收敛 设迭代 xk+1 = ?(xk) ，若 ?(p)(x) 在 x* 的某邻域内连续，则该迭代法具有 p 阶收敛的充要条件是 定理 并且有 证明： 充分性. 根据泰勒展开有 必要性. 设迭代 xk+1 = ?(xk) 是 p 阶收敛。 迭代两边取极限 ，由 ?(x) 的连续性可知 x* = ?(x*) 。 设 p0 是满足 的最小正整数。 由充分性的证明过程可知迭代 p0 阶收敛。 又 若 p0 p， 与迭代 p 阶收敛矛盾 p0 = p 设有不动点迭代： 设： 缺点：每次迭代需计算 设： Aitken 加速 当 x 收敛到 x* 时，修正项分子趋于零。 基本思想： 将非线性方程线性化 设 xk 是 f (x)=0 的近似根， 将 f (x) 在 xk 一阶 Taylor 展开: ，? 在 xk 和 x 之间。 x y x* xk xk+1 条件： f’(x) ? 0 Newton 法可以看作下面的不动点迭代： 其中