第七章非线性方程的数值解法Chapter7NumericalMethods.ppt

第七章 非线性方程的数值解法 /*Chapter7 Numerical Methods for Nonlinear Equations*/ * 求 f (x) = 0 的根 原理：若 f ?C[a, b]，且 f (a) · f (b) 0， 则 f 在 (a, b) 上必有一根。 1. 二分法 a b x1 x2 a b When to stop? 或 不能保证 x 的精度 x* ?2 x x* §6.1.1 二分法 误差 分析： 二分初始时刻 有误差 二分k 次产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： ①简单; ② 对f (x) 要求不高(只要连续即可) . ①无法求复根及偶重根 ② 收敛慢 注：用二分法求根，最好先给出 f (x) 草图以确定根的大概位置。或用搜索程序，将[a, b]分为若干小区间，对每一个满足 f (ak)·f (bk) 0 的区间调用二分法程序，可找出区间[a, b]内的多个根，且不必要求 f (a)·f (b) 0 。 §6.1.1 二分法 2. 不动点迭代法 /* 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 的根。 So basically we are done! I can’t believe it’s so simple! What’s the problem? Oh yeah? Who tells you that the method is convergent? §6.1.2 一般迭代法 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 ? §6.1.2 不动点迭代 ( 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]上的唯一不动点。并且有误差估计式： 定理 考虑方程 x = g(x), g(x)?C[a, b], 若 ? ? ( k = 1, 2, … ) 且存在极限 ? ? ? ? ? ? k §6.1.2 不动点迭代 证明：① g(x) 在[a, b]上存在不动点？ 令 有根 ② 不动点唯一？ 反证：若不然，设还有 ，则 在 和 之间。 而 ③ 当k ? ? 时， xk 收敛到 x* ？ ? ? ? §6.1.2 不动点迭代 ④ ? ⑤ ? ⑥ ? 可用 来控制收敛精度 L 越 收敛越快 小 注：定理条件非必要条件，可将[a, b]缩小，定义局部收敛性：若在 x* 的某? 领域 B? = { x | | x ? x* | ? ? } 有 g?C1[a, b] 且 | g’(x*) | 1，则由?x0?B? 开始的迭代收敛。即调整初值可得到收敛的结果。 §6.1.2 不动点迭代 局部收敛性与迭代法的收敛阶 §6.1.2 不动点迭代 定义 　　　设迭代 xk+1 = g(xk) 收敛到g(x) 的不动点 x*。 设 ek = xk ? x*，若　　　　　　　，则称该迭代为p 阶收敛。 ? 一般 Fixed-Point Iteration 有　　　　　　　 ，称为线性收敛 /* linear convergence */，这时 p = 1，0 C 1。 ? p 1，称为超线性收敛， ? p ＝ 2，称为平方收敛。 §6.1.2 不动点迭代 定理 设 x* 为x = g(x) 的不动点，若 ，p ? 2； ，且