数值分析对分法和一般迭代法.pptVIP

数值分析对分法和一般迭代法第一页，共四十一页，2022年，8月28日 非线性方程求根1．根的存在性。方程有没有根？如果有根，有几个根？定理1：设函数 f (x) 在区间[a, b]上连续，如果f (a) ? f (b) 0， 则方程 f (x) = 0 在[a, b]内至少有一实根x*。 2．这些根大致在哪里？如何把根隔离开来？3．根的精确化第二页，共四十一页，2022年，8月28日 abx*f(x)1．画出 f(x) 的略图，从而看出曲线与x 轴交点的位置。2．从左端点x = a出发，按某个预先选定的步长h一步一步地向右跨，每跨一步都检验每步起点x0和终点x0 + h的函数值，若那么所求的根x*必在x0与x0+h之间，这里可取x0或x0+h作为根的初始近似。第三页，共四十一页，2022年，8月28日 ??????????????????开 始读入a, ha ? x0f (x0) ? y0x0 + h ?x0f (x0) ? y00打 印结 束否是继续扫描?第四页，共四十一页，2022年，8月28日 例1：考察方程 的含根区间x00.51.01.5f (x)符号－－－＋可见，含根区间为 [1, 1.5]第五页，共四十一页，2022年，8月28日 §4.1 对分区间法 (Bisection Method )原理：若 f(x) ?C[a, b]，且 f (a) · f (b) 0，则f(x) 在 (a, b) 上必有一根。第六页，共四十一页，2022年，8月28日 abx1x2a1b2x*b1a2或不能保证 x 的精度?2xx*第七页，共四十一页，2022年，8月28日 执行步骤1．计算f (x)在有解区间[a, b]端点处的值，f (a)，f (b)。2．计算f (x)在区间中点处的值f (x1)。3．判断若f (x1) = 0，则x1即是根，否则检验：（1）若f (x1)与f (a)异号，则知解位于区间[a, x1]， b1=x1, a1=a;（2）若f (x1)与f (a)同号，则知解位于区间[x1, b]， a1=x1, b1=b。反复执行步骤2、3，便可得到一系列有根区间： (a, b), (a1, b1), …, (ak, bk), …第八页，共四十一页，2022年，8月28日 4、当时，则即为根的近似。①简单; ② 对f (x) 要求不高(只要连续即可) .①无法求复根及偶重根 ② 收敛慢 第九页，共四十一页，2022年，8月28日 误差 分析：第1步产生的有误差第 k 步产生的 xk 有误差对于给定的精度 ? ,可估计二分法所需的步数 k ：第十页，共四十一页，2022年，8月28日 定义f (x)f (a)? f (b)0f (a)? f (b)=0f (a) =0打印b, k打印a, k结束是是是否否否m=(a+b)/2|a-b|?f(a)?f(b)0打印m, ka=mb=m结束k=k+1是是否否输入 k = 0第十一页，共四十一页，2022年，8月28日 例2 用二分法求 在(1,2)内的根，要求绝对误差不超过 解： f(1)=-50 有根区间 中点 f(2)=140 -(1,2)+ f(1.25)0 (1.25,1.5) f(1.375)0 (1.25,1.375) f(1.313)0 (1.313,1.375) f(1.344)0 (1.344,1.375) f(1.360)0 (1.360,1.375) f(1.368)0 (1.360,1.368) f(1.5)0 (1,1.5) 第十二页，共四十一页，2022年，8月28日 12 例3，求方程f(x)= x 3 –e-x =0的一个实根。 因为 f(0)0,f(1)0。 故f(x)在(0,1)内有根用二分法解之，(a,b）=（0,1)’计算结果如表：k a bk xk f(xk)符号 0 0 1 0.5000 －