计算方法 教学课件 作者 江爱民 主编 张有正 副主编 nn03.pptVIP

上一页 下一页 返回首页 第3章 非线性方程的数值解法 3.1 二分法 3.2 简单迭代法 3.3 牛顿迭代法 3.4 插值法 3.1 二分法 3.1.1 根所在范围的确定 ?根存在定理 若 f ?C[a, b]，且 f a · f b 0，则 f 在 a, b 上必有一根。 为了确定根所在区间，一般可以用逐步探索法，即从某点 出发，以选取适当的 作为步长，计算 ，如 ，以 代替 ，重复以上步骤，继续往前，直到 。如 说明 或 即为方程 的解，如 ，则在区间 内方程有一根。 a b 步长过大时扫描可能会出现的情况 h 3.1.2 二分法 a b x1 x2 a b 何时结束? 或 不能保证 x 的精度 x* ?2 x x* 第一步 取 第二步 若 ，则 即为方程的根，求解结束。 第三步 若 ，说明根在区间 内，将 ；否则说明根在区间 内，将 第四步 若 ，则过程结束，以 作为根的近似解，否则重复以上过程。 ?二分法的要点可以描述如下： 例1 用二分法求方程 在区间 内的近似根，要求误差不超过 。 继续重复以上的过程，我们将各次计算的结果列于下表： 3.2 简单迭代法 3.2.1 简单迭代法的基本过程 迭代法是一种非常重要的逐次逼近的方法，所谓迭代法就是用一个公式反复校正根的近似值，使之逐步精确化，直到达到精度要求为止。 ?迭代法 f x 0 x ? x 等价变换 f x 的根 ? x 的不动点 思路 从一个初值 x0 出发，计算 x1 ? x0 , x2 ? x1 , …, xk+1 ? xk , … 若 收敛，即存在 x* 使得 ，且 ? 连续，则由 可知 x* ? x* ，即x* 是 ? 的不动点，也就是f 的根。 x y y x x y y x x y y x x y y x x* x* x* x* y g x y g x y g x y g x x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 3.2.1 简单迭代法的几何解释与收敛性 对于简单迭代法的收敛性,简单说明如下: 解: 等价变换 继续上面的过程，具体计算如下表所示 3.3.1 公式的导出 3.3 牛顿迭代法 原理：将非线性方程线性化 —— Taylor 展开 取 x0 ? x*，将 f x 在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间。 将 x* ? x0 2 看成高阶小量，则有： 线性 x y x* x0 定理 （收敛的充分条件）设 f ?C2[a, b]，若 1 f a f b 0； 2 在整个[a, b]上 f ”不变号且 f ’ x ? 0； 3 选取 x0 ? [a, b] 使得 f x0 f ” x0 0； 则牛顿法产生的序列 xk 收敛到f x 在 [a, b] 的唯一根。 有根 根唯一 产生的序列单调有界，保证收敛。 3.3.2 牛顿迭代法的几何意义及收敛讨论 x y x* x0 注：牛顿迭代法的 收敛性依赖于x0 的选取,定理中的 3 给出了选取x0 一般原则。 x* x0 ? x0 ? x0 对切线法的要点可描述如下：