计算方法第3章非线性方程的数值解法.ppt

上一页 下一页 返回首页 第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 对切线法的要点可描述如下：