计算方法非线性方程求根.ppt

? /* Secant Method */ ： 正割法 Newton’s Method 一步要计算 f 和 f ’，相当于2个函数值，比较费时。现用 f 的值近似求 f ’，可少算一个函数值。 x0 x1 切线 /* tangent line */ 割线 /* secant line */ 切线斜率　?　割线斜率 需要2个初值 x0 和 x1。 收敛比Newton’s Method 慢，且对初值要求同样高。 定义 　　　设迭代 xk+1 = g(xk) 收敛到g(x) 的不动点 x*。 设 ek = xk ? x*，若　　　　　　　，则称该迭代为p 阶收敛，其中 C 称为渐进误差常数。 ? 一般 Fixed-Point Iteration 有　　　　　　　 ，称为线 性收敛 ，这时 p = 1，0 C 1。 ? Aitken 加速有　 　　　　 。 称为超线性收敛 §2.6 迭代法的收敛阶 ? Newton’s Method 有 ，只要 就有 p ? 2。重根是线性收敛的。 Q: 如何实际确定收敛阶和渐进误差常数？ 定理 设 x* 为x = g(x) 的不动点，若 ，p ? 2； ，且 ，则 xk+1 = g(xk) 在 内 p 阶收敛。 证明： x* k ? ? C * 软件+交设10-123（第4周五34） * * 第2章 非线性方程求根 求 f (x) = 0 的根 非线性方程求解的基本思想 确定非线性方程实根范围的方法 图解法： 计算机比较难实现，对人使用方便。 逐步扫描法： 便于计算机实现。 2. 对方程根进一步精确化的方法 二分法，迭代法，Newton迭代法 依据：若 f ?C[a, b]，且 f (a) · f (b) 0，则 f 在 (a, b) 上必有一根。此时[a, b]称为有根区间。 §2.1 二分法 原理 设[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 。 误差 分析： 例 例 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 的根。 我不相信如此简单， 问题究竟是什么？ 如何保证它的收敛性？ §2.2 迭代法 迭代法的几何意义 方程x = g (x