第7章节.非线性方程迭代法.pptVIP

非线性方程迭代法 邹昌文 * 迭代法 / Fixed-Point Iteration / f (x) = 0 x = g (x) 等价变换 f (x) 的根 g (x) 的不动点 思路 从一个初值 x0 出发，计算 x1 = g(x0), x2 = g(x1), …, xk+1 = g(xk), … 若 收敛，即存在 x* 使得 ，且 g 连续，则由 可知 x* = g(x* )，即x* 是 g 的不动点，也就是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 ? 定理 （充分条件）考虑方程 x = g(x), g(x)?C[a, b], 若 ( I ) 当 x?[a, b] 时， g(x)?[a, b]； ( II ) ? 0 ? L 1 使得 | g’(x) | ? L 1 对 ? x?[a, b] 成立。 则任取 x0?[a, b]，由 xk+1 = g(xk) 得到的序列 收敛于g(x) 在[a, b]上的唯一不动点。并且有误差估计式： ? ? ( k = 1, 2, … ) 且存在极限 ? ? ? ? ? ? 证明：① g(x) 在[a, b]上存在不动点？ 令 有根 ② 不动点唯一？ 反证：若不然，设还有 ，则 在 和 之间。 而 ③ 当k ? ? 时， xk 收敛到 x* ？ ? ? ? ④ ? ⑤ ? ⑥ ? 可用 来控制收敛精度 L 越 收敛越快 小 注：定理条件非必要条件，可将[a, b]缩小，定义局部收敛性：若在 x* 的某? 领域 B? = { x | | x ? x* | ? ? } 有 g?C1[a, b] 且 | g’(x*) | 1，则由?x0?B? 开始的迭代收敛。即调整初值可得到收敛的结果。 设在区间[a, b]上方程 x= φ (x) 有根x*,且对一切x∈[a, b] 都有|φ’(x)|≥ 1，则对于该区间上任意x0(≠x*), 迭代公式xk+1= φ (xk )一定发散。 证明: 不可能收敛于0。 定理 改进、加速收敛 / accelerating convergence / ? 待定参数法： 若 | g’(x) | ? 1，则将 x = g(x) 等价地改造为 求K，使得 例：求 在 (1, 2) 的实根。 如果用 进行迭代，则在(1, 2)中有 现令 希望 ，即 在 (1, 2) 上可取任意 ，例如K = ?0.5，则对应 即产生收敛序列。 ? Aitken 加速： x y y = x y = g(x) x* x0 P(x0, x1) x1 x2 P(x1, x2) 一般地，有： 比 收敛得略快。 ? Steffensen 加速： 详见P273 定理 （牛顿法收敛的充分条件）设 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； 则Newton’s Method产生的序列{ xk } 收敛到f (x) 在 [a, b] 的唯一根。 有根 根唯一 产生的序列单调有界，保证收敛。 定理 （局部收敛性）设 f ?C2[a, b]，若 x* 为 f (x) 在[a, b]上的根，且 f ’(x*) ? 0，则存在 x* 的邻域 使得任取初值 ，Newton’s Method产生的序列{ xk } 收敛到x*，且满足 证明：Newton’s Method 事实上是一种特殊的不动点迭代 其中 ，则 收敛 由 Taylor 展开： 只要 f ’(x*) ? 0，则令 可得结论。 在单根 /simple root / 附近收敛快 ? 牛顿迭代法的改进与推广 ? 重根 / multiple root / 加速收敛法： Q1: 若　　　 ，Newt