第二讲 方程求根.ppt

§3 简单迭代法 因此简单迭代法又称为不动点迭代法。 §4 迭代过程的加速 §5 Newton迭代法及其变形 用一系列切线与x轴交点的横坐标逼近函数f与x轴的交点 习题二 P Q 一、 Newton迭代法及其局部收敛性 对于方程 令： 该格式称为Newton迭代格式。 几何意义： 数学上：借助于Taylor公式（级数）将非线性问题线性化处理； 几何上：利用直线代替曲线； 迭代上：追求更快收敛速度 Newton-Raphson 局部收敛性 Th：若方程 在根 附近具有二阶连续的导数， 且 则Newton迭代法在根的附近至少是平方收敛的。 证明： Newton迭代法算法流程 输入：初始近似值x0；精度要求eps；最大迭代次数N； 输出：迭代次数i与近似解x或失败信息 Step1 Set i=1; Step2 while i≤N do Steps 3—6; Step3 Set x=x0-f(x0)/f’(x0); Step4 if abs(x-x0)<eps then output(i,x);(算法成功) stop Step5 Set i=i+1; Step6 Set x0=x; Step7 output(‘The method failed after N iterations,N=’,N) 例8 用Newton迭代法求方程xex-1=0在0.5附近的根,精度要求?=10-5. 解 Newton迭代格式为 0000-0000.0000000003 0.0000000003 0.5 00000 1 2 3 4 |xk-xk-1| ?(xk) xk k exp8文件夹 称为下山因子的下界。 二、Newton迭代法的改进与变形 1、Newton下山法： Newton迭代法的局部收敛性对于初值的选取方式依赖性很大，在事先无法给出较好的初值情况下，可用如下迭代格式求方程的根的近似值。 称为下山因子，它的作用是保证函数列 从而保证迭代法的收敛性。 单调递减， 通常，下山因子取小于1大于 的常数。 事实上，Newton下山法可以看作是由如下过程得到的： 将Newton迭代法得到的近似值 与上一步的近似值 做如下线性组合 思想：综合前两步的结果产生一个新的近似值（一般来说，可以加快收敛）。 2、求重根（重数大于等于2）的Newton迭代法 （1）重数m已知的情况： （2）重数未知的情况：（Newton-Raphson方法） 令 可以证明若 是方程 的m重根， 则 一定是方程 的单根， 从而可以用Newton迭代法求根的近似值， 且收敛速度至少是2阶的。 （3）重数的估计 3、Newton割线法（快速弦截法） 几何意义：用割线PQ与x轴交点的横坐标逼近根。 （需要两个初值）收敛阶为1.618 P Q Newton割线法算法流程 输入：初值x0，x1，精度要求eps，最大迭代次数N 输出：迭代次数k，近似解x或失败信息 Step1 Set k=2 p0=f(x0);p1=f(x1); Step2 while k≤N Do Steps 3—6 Step3 Set x=x1-p1(x1-x0)/(p1-p0); Step4 if abs(x-x1)<eps then output(k,x);(算法成功) Stop； Step5 Set k=k+1； Step6 Set x0=x1；p0=p1；x1=x；p1=f(x1); Step7 output(‘The method failed after N iterations,N=’,N)(算法失败) 4、抛物线法 设已知方程f(x)=0的三个近似根 则可以过如下三点 即构造一个二次多项式P2(x)近似函数f(x)，从而可以用二次方 构造一条抛物线， 程P2(x)=0的根近似原方程的根。这种方法称为抛物线法或密勒 法。 具体格式如下：首先得到二次多项式P2(x) 其中， 从而可解得一元二次方程P2(x)=0的两个根，并将其看作方程f(x)=0的新的根的近似值 其中， 问题：计算结果中会产生两个根的近似值，如何取舍？ 在方程f(x)=0的三个近似根 中，我们自然