ch非线性方程的数值解法.pptVIP

* * 第二章 非线性方程的数值解法 2.1 二分法 2.2 一般迭代法 2.3 牛顿迭代法 2.4 弦截法 （1）确定初始含根区间 数值计算方法主要分为两大类。 第一类是区间收缩法。 （2）收缩含根区间 第二类是迭代法。 （1）选定根的初始近似值 （2）按某种原则生成收敛于根的近似点列 2.1 二分法 基本假设 : 2.1.1 二分法的计算步骤 常用终止原则为: 2.1.2 二分法的收敛性与事前误差估计 所以，二分法总是收敛的 例 2.1 试用二分法求 的一个正根，使误差小于10-3 故可取初始区间 解 2.1.3?? 二分法评述 优点：简单可靠，易于编程实现，它对函数要求低，适 用于的奇数重根情形。 缺点：不能直接用于求偶重根，不能用于求复根，也难 以向方程组推广使用，收敛速度慢。 2.2 一般迭代法 迭代法的算法思想为： (1) 把（2-2）等价变换为如下形式 (2) 建立迭代格式 或更一般地建立迭代格式 (3) 适当选取初始值，递推计算出所需的解。 2.2.1 迭代法的算法思想 2.2.2 迭代法的收敛性 则称在 内 李普希兹连续。 定义2.1 设在某个区间 ? 内，函数 ? 满足下述李普希兹条件： 则 在 内李普希兹连续。 命题2.1 若 在闭区间 内连续且 命题得证。 证 定理2.1 设 x*= g(x*), g(x) 在闭区间 ?： 内李普希兹连续，则对任何初值 由迭代格式 xk+1= g(xk) 计算得到的解序列 收敛于 x* ( 这时我们称迭代格式 xk+1= g(xk)在 x* 的邻域 ? 上局部收敛）。 （1）首先用数学归纳法证明： 由假设知 又设 ，则 综上，由归纳法原理知，结论成立。 证 因此， ，定理得证。 反设存在 矛盾。所以结论成立。 2) 迭代函数在 x* 附近李普希兹连续从而收敛的迭代格式统称为皮卡（Picard）迭代 (2) 由（1）的结论和 g(x) 在 ? 内李普希兹连续的假设，可递推得到 注 1) g(x) 在?内李普希兹连续的条件保证了x* 为 f (x)=0 在 ? 内的唯一根。 证 推论 设 x*= g(x*) ， 若 g(x) 在 x* 附近连续可微且 ，则迭代格式 xk+1= g(xk) 在 x* 附近局部收敛。 注 由于 x* 事先未知，故实际应用时，代之以近似判则 。 但需注意，这实际上是假设了x0充分接近 x* ，若 x0 离 x* 较远，迭代格式可能不收敛。 定理2.2 （非局部收敛定理）如果 在 上连续可微且以下条件满足： 注 虽然定理2.1的条件是充分条件，但其条件并不很强，实际上，我们易证如下命题。 命题2.2 若在区间 内 ，则对任何 ，迭代格式 不收敛。 2.2.3 迭代法的误差估计 故对正整数 p，有 （2）事后误差估计 由此，对给定的精度 可进行 (1) 事前误差估计 简单地代之以 或 例 2.2 试建立收敛的迭代格式求解 在区间(2,3)内满足精度要求 的根。 首先可简单的把 等价化为 由此建立迭代格式 所以该迭代格式在内不收敛，不可取。 为建立收敛的迭代格式，我们把 等价化为 从而建立迭代格式 解 易知在 x 0时 g(x) 单调增，故有 2 g(2) g(x) g(3) 3 故由定理2.2得：任取x0 ?[2,3]，该迭代格式收敛。 取 x0=2 计算，结果见表 2-2（书17页）。 2.2.4 迭代法的收敛速度与加速收敛