数值分析-求根方程.pptVIP

非线性方程求根/* Solutions of Nonlinear Equations */ §1　 引言 §2 二分法 §3 迭代法 §4 牛顿法 §5 劈因子法 * * 引言 /* Introduction */ 求 f (x) = 0 的根 方程f (x)=0,如果存在某个x*，使得f (x*)=0。则称x*为方程f (x) = 0 的根或者称为函数f (x)的零点。 特别，如果 则称x*为f(x)=0的m重根，或者f(x)的m重零点。当m=1,称为单根。 1. 基本概念 2. 方程求根一般步骤 求出有根区间----方程在此区间有且仅有一根。 将根细化----将已知的近似根用各种方法进一步提高精确， 使其满足给定的精确要求。 3. 隔离区间求法 原理：若 f (x)?C[a, b]，且 f (a) · f (b) 0，则 f (x)在 (a, b) 上必有一根。 作出f (x)函数的草图，由f (x)与x轴交点横轴标大概确定隔 离区间。 将f(x)=0等价改写为f1(x)=f2(x)，作出f1(x)、f2(x)的草 图，由两者交点横轴标大概确定隔离区间。 逐步搜索法，得到满足原理条件的各个区间。 例1：求方程 的有根区间。 解法1：图解法。 由左图所示，得到有根区间 为[1.6, 1.7]。 解法2：分析法。 严格单调。 故函数的有根区间为[0,2] 事实上， 解：当x0 时，f(x)0，故有根区间肯定在 例2：求方程 的有根区间。 分别作出f1(x), f2(x)的图像， 由两图像交点横坐标看出有 根区间为：[0.5,1]。 实际上， 假设区间 内有唯一的实根。 二分法 /* Bisection Method */ 于是有根区间就变为 ，长度只有原区间 的一半。 取 的中点 取 的中点 区间 中点 显然每个小区间都有单根且有如下关系： 可得一系列的小区间和中点 如果取 点列 a b x1 x2 a b When to stop? 或 不能保证 x 的精度 x* ?2 x x* §2 Bisection Method 误差 分析： 第1步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： ①简单; ② 对f (x) 要求不高(只要连续即可) 。 ①无法求复根及偶重根; ② 收敛慢 。 例：用二分法求方程 在区间[0,1]内的一个实根，要求精确到3位有效数字。 解： 号 号 号 1 0(+) 0.5(-) 0.25(+) 0 0(+) 1.0(-) 0.5(-) 2 0.25(+) 0.5(-) 0.375(+) 3 0.375(+) 0.5(-) 0.4375(+) 4 0.4375(+) 0.5(-) 0.46875(+) 迭代法 /* Fixed-Point Iteration */ f (x) = 0 x =φ(x) 等价变换 思路 从一个初值 x0 出发，计算 x1 =φ(x0), x2 =φ(x1), …, xk+1 =φ(xk), … 。若 收敛，即存在 x* 使得 ，且φ连续，则由 可知 x* =φ(x* )，即x* 是f 的根。 称xn+1 =φ(xn)为迭代格式(迭代公式)；称φ(x)为迭代函数；求得序列 为迭代序列；上述求根方法称为迭代法。 迭代格式的构造及其收敛性条件 例：给定方程 f (x)=x2-2x-3=0，显然它在区间[2,4]内有唯一实根x*=3。按不同迭代格式来求此解。 解： 可以看出，当k－＞+∞时，xk≈3.0 可以看出，当k－＞ +∞时，xk －＞ +∞ 收敛 ( I ) 当 x?[a, b] 时，φ(x)?[a, b]； ( II ) ? 0 ? L 1 使得 |φ’(x) | ? L 1 对 ? x?[a, b] 成立。 其中L为正常数。则有如下结论成立： 定理 考虑方程 x =φ(x), φ(x)?C1[a, b], 若φ