求方程根的近似方法.pptVIP

§2.2 迭代法 2.收敛定理 * f(x)=0的根(或f(x)的零点)，当f(x)复杂时，很难求（需要找到有效简单的近似方法去求）。 第二章 求方程根的近似方法 §2.1 二分法 理 论: f(x) ∈ C[a,b], 单调, f(a)f(b)0 f(x)=0在(a,b)中有惟一根。 a b x1 x2 a b 何时停下来? 或 不能保证 x 的精度 x* ?2 x x* 解： f(1)=-50 有根区间 中点xn f(2)=140 -(1,2)+ x1=1.5 f(1.5)0 (1,1.5) x2=1.25 f(1.25)0 (1.25,1.5) x3=1.375 f(1.375)0 (1.25,1.375) x4=1.313 f(1.313)0 (1.313,1.375) x5=1.344 f(1.344)0 (1.344,1.375) x6=1.360 f(1.360)0 (1.360,1.375) x7=1.368 f(1.368)0 (1.360,1.368) x8=1.364 例2.1.1 用二分法求 在(1,2)内 的根，要求绝对误差不超过 ，则 （事后估计） 误差 分析： 第1步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? , 可估计二分法所需的步数 k ： 缺点：收敛速度慢， 不易求偶数重根. 如图 注：用二分法求根，最好先给出 f (x) 草图以确定根的大概位置。或用搜索程序，将[a, b]分为若干小区间，对每一个满足 f (ak)·f (bk) 0 的区间调用二分法程序，可找出区间[a, b]内的多个根，且不必要求 f (a)·f (b) 0 。 y x 优点：条件和方法简单(只要求f(x)连续即可)，方法收敛； 一. 迭代法的建立与收敛性 所以, ?为f的根的充要条件是?为?的不动点。 前者收敛：1.5; 1.35721; 1.33086; 1.32588; 1.32494; 1.32476; 1.32473; 1.32472; 1.32472;… 后者发散: 1.5; 2.375; 12.39; … 问题：何时收敛？ x y y = x x y y = x x y y = x x y y = x ? ? ? ? y=?(x) y= ?(x) y= ?(x) y= ?(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 定理2.2.1 注1：L越小，收敛越快。 由定理结论(3)或(2.2.2)，只要前后两次迭代值的差值足 够小，就可使近似值 达到任意的精度。在实际计算 中，一般用 来控制迭代过程结束。 注2：定理条件非必要条件，可将[a, b]缩小，定义局部收敛性： 定义2.2.1 若存在 ? 的某? 邻域 B? = { x | | x ? ? | ? ? } , 使由?x0?B? 开始的迭代都收敛, 则称迭代法具有局部收敛性。 定理2.2.2 设?(x)在?的某?邻域内具有连续的一阶导数， 且 | ?(?) | 1， 则迭代法xn+1 = ?(xn)具有局部收敛性。 证明 省略。 3．编程停机判断 时，由（2.2.2）式知 比较小，此时停机， （取定初值x0）计算，当 由 二.迭代法的收敛阶(收敛速度) 则称xn p阶收敛,相应的迭代法称为p阶方法. 特别, p=1时叫线性收敛,此时要求0C1; p=2时叫平方收敛; p越大越好.(why?) 若有实数 p0,使 定义2.2.2: 设 定理2.2.3 设?(x)在?的某邻域 内有充分多阶连续导数，则迭代法xn+1 = ?(xn)为p阶收敛的充要条件是 ??(?