数值计算方法 (第2章 非线性方程与方程组的数值解法).ppt

* 代数方程的Newton迭代法 * * 代数方程的Newton迭代法 同理 * 代数方程的Newton迭代法 比较x的同次幂系数得： 故代数方程的Newton迭代公式 * 代数方程的Newton迭代法算法 * * 2.4弦截法 Newton迭代法有一个较强的要求是 且存在。因此，用弦的斜率近似的替代 。 * 弦截法 令y=0，解得弦与x轴的交点是坐标x2 * 弦截法 * 弦截法的几何解释 (a) 定端点弦截法 (b) 变端点弦截法 * 例题 例2.4.1 用快速弦截法求方程 在区间（1，2）内的实根。 解：取x0=1,x1=2,代入公式2.4.2计算结果,如表2.4.1所示。 * k xk f(xk) 0 1 -1 1 2 5 2 1.166666667 -03 1.253112023 -04 1.337206444 0.053880579 5 1.323850096 -0.0036981168 6 1.324707936 -4.273521*10E-5 7 1.324717965 3.79*10E-8 * 弦截法收敛定理 * 弦截法收敛定理 * 求解方程f(x)=0的快速弦截法 * 复习 本章的主要内容 二分法 等步长扫描法、误差估计等 一般迭代法、压缩映像原理 牛顿切线法、代数方程的牛顿切线法 公式推到、几何意义、收敛性 弦截法 定端点弦截法、动端点弦截法 * 例1：用二分法求方程x3-x-1=0在[1,2]的近似解，误差若要求小于10-3，至少要需要迭代多少次？ 解：设连续函数f(x)=x3-x-1，则f ’(x)=3x2-1， 在区间[1,2]上f ’(x)0，f(x)单调联系。 且f(1)=-10，f(2)=50，所以f(x)在[1,2]上 有唯一根。 * 由二分法误差估算式 * 例2：证明对任何初始值x0?R，由迭代公式 xk+1=cosxk，k=0,1,2,…所产生的序列 都收敛于方程x=cosx的根。 * * * * 作业(P40) 1、4 * * * * Steffensen算法的收敛速度 * Steffensen算法的收敛速度 * Steffensen算法的收敛速度 * Steffensen算法的收敛速度 由定理2.2.4知 至少以平方速度收敛到 。 也就是说：简单迭代法是线性收敛；Steffensen迭代至少平方以上收敛（加速收敛）。 * 例题 例2.2.3试用Steffensen算法求解方程 解法一、取 ，由 n = 0，1，2，… * 例题 取初值 ，计算结果如下： N Xn Yn Zn 0 1.5 1.357208808 1.330860959 1 1.324899181 1.324752379 1.324724496 2 1.324717957 1.324717957 1.324717957 * 例题 解法二、取 ，由 对于该迭代函数在一般迭代法中是发散的，而Steffensen格式却是收敛的。 n=0，1，2，… * 例题 取初值 ，计算结果如下： N Xn Yn Zn 0 1.5 2.375 1.239648437 1 1.416292975 1.840921915 5.238872769 2 1.355650442 1.491398279 2.317270699 3 1.328948777 1.347062883 1.444351224 4 1.324804489 1.325173544 1.327117281 5 1.324717944 1.324718152 1.324718980 6 1.324717957 * Steffensen迭代格式几何解释 * Steffensen迭代算法 * Steffensen迭代算法 * 2.3 Newton迭代法 设x*是方程f (x ) = 0的根， 又x0 为x* 附近的一个值， 将f (x ) 在x0 附近做泰勒展式： 令 ，则 * Newton迭代法 去掉 的二次项，有： 即 以x1代替x0重复以上的过程，继续下去得： * Newton迭代法 以此产生的序列{Xn}得到f(x)=