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

* * 第二章 非线性方程的数值解法 2.1 二分法 2.2 一般迭代法 2.3 牛顿迭代法 2.4 弦截法 （1）确定初始含根区间 数值计算方法主要分为两大类。 第一类是区间收缩法。 （2）收缩含根区间 第二类是迭代法。 （1）选定根的初始近似值 （2）按某种原则生成收敛于根的近似点列 2.1 二分法 一、根的隔离 将含根区间一个个隔开，找到根的范围，使每个区间只有一个根。 定理：对f(x)=0,f(x)在[a,b]上连续, f(a)f(b)0且f(x) 严格单调上升（或严格单调下降），则f(x)在[a,b]内仅有一根。 1。利用零点存在定理。 2.作图 x y y=tgx y=1/x 找出交点 3.隔离法： 二、对分法 设f(x)在(a,b)上连续且f(x)=0在(a,b)只有一个根 1.算法 2.收敛性 3.误差控制 例2.1：试用二分法求 的非零实根，使其误差小于 0 (+)1.5 (-)2 1.75 + 1.75 2 1.875 + 1.875 2 1.9375 - 1.875 1.9375 1.90265 + 1.90265 1.9375 1.921875 + 5 1.92188 1.9375 1.92688 + 2.1.3?? 二分法评述 优点：简单可靠，易于编程实现，它对函数要求低，适用于的奇数重根情形。 缺点：不能直接用于求偶重根，不能用于求复根，也难以向方程组推广使用，收敛速度慢。 2.2 一般迭代法 迭代法的算法思想为： (1) 把（2-2）等价变换为如下形式 (2) 建立迭代格式 或更一般地建立迭代格式 (3) 适当选取初始值，递推计算出所需的解。 一 迭代法的算法思想 二 迭代法的收敛性 则称在 内 李普希兹连续。 定义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 若在区间 内 ，则对任何 ，迭代格式 不收敛。 三、 迭代法的误差估计 故对正整数 p，有 （2）事后误差估计 由此，对给定的精度