2 非线性方程求根.pptVIP

? /* Descent Method */ ——Newton’s Method 局部微调： 下山法 原理：若由 xk 得到的 xk+1 不能使 | f | 减小，则在 xk 和 xk+1 之间找一个更好的点 ，使得 。 xk xk+1 注：? = 1 时就是Newton’s Method 公式。 当 ? = 1 代入效果不好时，将 ? 减半计算。 求复根 ? /* Finding Complex Roots */ —— Newton 公式中的自变量可以是复数 记 z = x + i y, z0 为初值，同样有 设 p.28 代入公式，令实、虚部对应相等，可得 定义 　　　设迭代 xk+1 = g(xk) 收敛到g(x) 的不动点 x*。 设 ek = xk ? x*，若　　　　　　　，则称该迭代为p 阶收敛，其中 C 称为渐进误差常数。 ? 一般 Fixed-Point Iteration 有　　　　　　　 ，称为线 性收敛 ，这时 p = 1，0 C 1。 ? Aitken 加速有　 　　　　 。 称为超线性收敛 §5 迭代法的收敛阶 ? Newton’s Method 有 ，只要 就有 p ? 2。重根是线性收敛的。 Q: 如何实际确定收敛阶和渐进误差常数？ 定理 设 x* 为x = g(x) 的不动点，若 ，p ? 2； ，且 ，则 xk+1 = g(xk) 在 内 p 阶收敛。 证明： x* k ? ? C 道铁09-12班第2周1 信息工程09-123班第2周1 第二周5 * 第2章 非线性方程求根 求 f (x) = 0 的根 §1 多项式及代数方程根的界 （自习） 依据：若 f ?C[a, b]，且 f (a) · f (b) 0，则 f 在 (a, b) 上必有一根。此时[a, b]称为有根区间。 §2 二分法 原理 设[a, b]是方程f(x)=0的一个有根区间，取 是[a1, b1]的中点。 若f(x1) = 0，则x1是f(x) = 0的一个根， 若f(x1) f(a1) 0，则取a2=x1，b2=b1，否则取a2=a1，b2=x1 得到[a2, b2]，满足： 以[a2, b2]取代[a1, b1]，继续以上过程，得到[a3, b3]……, 直到某一xk时， ，使f(xk)=0，则xk是f(x)=0 的根，若不存在这样的xk，能得到一系列闭区间： 表明存在x*， f(x*)=0 ， x* ?[ak , bk]，因此{ak }单调上升，有上界x* ， {bk }单调下降，有下界x* ，且这二个序列均存在极限： － － 定 理2-4 若 f ?C[a, b]，且 f (a) · f (b) 0，则由二分法产生的序列{xk} 收敛于f (x) =0的一个根x*，且 a b x1 x2 a b When to stop? 或 不能保证 x 的精度 x* ?2 x x* 第1步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： ①简单; ② 对f (x) 要求不高(只要连续即可) . ①无法求复根及偶重根 ② 收敛慢 注：用二分法求根，最好先给出 f (x) 草图以确定根的大概位置。或用搜索程序，将[a, b]分为若干小区间，对每一个满足 f (ak)·f (bk) 0 的区间调用二分法程序，可找出区间[a, b]内的多个根，且不必要求 f (a)·f (b) 0 。 误差 分析： 二分法求根算法 例 例 a b (a+b)/2 x* (a, f (a)) (b, f (b)) Is it really better than Bisection Method? 注：试位法每次迭代比二分法多算一次乘法，而且不保证收敛。 ? 试位法（比例求根法） f (x) = 0 x = g (x) 等价变换 f (x) 的根 g (x) 的不动点 思路 从一个初值 x0 出发，计算 x1 = g(x0), x2 = g(x1), …,