第2章迭代法选编.ppt

第2章 方程求根; 图 2.3 ;预备定理;§2 迭代法 ;计算结果见下表; 对于一般形式的方程(2.1),首先我们设法将其化为下列等价形式 x=g(x) (2.7) 然后按(2.7)构造迭代公式 (2.8) 从给定的初始近似根x0出发,按迭代公式(2.8)可以得到一个数列 ?x0,x1,x2,…,xk,… 若这个数列｛xk｝有极限,则迭代公式(2.8)是收敛的。此时数列的极限; 就是原方程(2.1)的根。 虽然迭代法的基本思想很简单,但效果并不总是令人满意的。对于上例,若按方程写成另一种等价形式 x=x3-1 (2.9) 建立迭代公式 ? xk+1=x3k-1, k=0,1,2,… 仍取初始值x0=1.5, 则迭代结果为 x1=2.375 x2=12.3976 ; 迭代法的几何意义:把方程(2.1)求根的问题改写成(2.7)变为求数列｛xn｝的极限,实际上是把求根问题转化为求 ; 图 2.4 ; 迭代过程(2.8)就是在x轴取初始近似值x0,过x0作y轴的平行线交曲线y=g(x)于p0,p0的横坐标为x0,纵坐标为g(x0)(g(x0)=x1),也即 p0(x0,x1) 再在x轴上取x1作为新的近似值,过x1作y轴的平行线交曲线y=g(x)于p1,p1的横坐标为x1,纵坐标为 g(x1)(g(x1)=x2),也即 p1(x1,x2) 而这相当于过p0引平行于x轴的直线交y=x于 ? Q1(x1,x2); 再过Q1引平行于y轴的直线交曲线y=g(x)于 ? p1(x1,x2) ?仿此可得到点列 ?p0(x0,x1),p1(x1,x2),p2(x2,x3),… ?若 ; 则迭代法收敛,见图2.4(a);否则迭代法发散,见图2.4(b)。 必须说明两点: ①若g(x)可微,可用充分条件 ;这里q＜1是非常重要的条件,否则不能保证迭代收敛。 ②对于收敛的迭代过程,误差估计式(2.11)说明迭代值的偏差｜xk-xk-1｜相当小,就能保证迭代误差｜x-xk｜足够小。因此在具体计算时常常用条件 ｜xk-x k-1｜＜ε (2.15) 来控制迭代过程结束。 ; 迭代法的突出优点是算法的逻辑结构简单,且在计算时,中间结果若有扰动,仍不会影响计算结果。其计算步骤为: (1)确定方程f(x)=0的等价形式x=g(x),为确保迭代过程的收敛,要求g(x)满足｜g′(x)｜≤q＜1 (2)选取初始值x0,按公式 ? x k+1=g(xk), k=0,1,2,… ?进行迭代; (3)若｜x k+1-xk｜＜ε,则停止计算,x≈x k+1。; 例2 求方程 x= e-x 在x=0.5附近的一个根。按五位小数计算,计算结果 的精度要求为ε=10-3。