研究生数值分析第2章B.doc

§2.4牛顿(Newton)迭代法 一、牛顿法的构造 在求解方程时,设可导,且在附近,则在附近一次泰勒展式 如果是的解,则代入上式解得: 将左边作的一次近似,,反复运用以上公式,得到迭代公式: 得数列,如收敛,则。 二、几何意义 是曲线,其一次泰勒展开式,是过曲线上的切线，而满足,即是该切线与x轴的交点。而是曲线与x轴的交点，因而每次近似均是过的切线与x轴的交点来近似逼近。 三、牛顿法计算步骤 选取初值 对于 计算和 判断 如果是转出循环 如果不是继续循环，输出 四、数值例子 例1:用牛顿法解方程,精确到0.01． 解:,,,, 选取初值 ，此时迭代结束，． 用两分法求解方程达到0.00065的精度需要10次达到要求。 用一般迭代法构造 七次迭代可达到以上精度,而牛顿迭代只需三次。 五、收敛性问题 牛顿法是一种迭代法，令，则牛顿迭代就是而此时，因此如果f(x)在x*的某个邻域内f’(x)≠0时，，所以由定理2.3知，．此即定理4.1的结果。 六、牛顿下山法 牛顿法是一种局部收敛方法，初值要取在解附近才能保证迭代序列收敛。为扩大收敛范围，使对任意初值迭代序列都收敛，我们可引入一个参数，并将牛顿法改为： 其中，称为下山因子，此迭代法称为牛顿下山法。 计算时可取，直到满足为止。由此得到的序列由于满足下山条件，从而单调减，故收敛。但只是线性收敛。 七、多重根的情况 如果x*是f(x)=0的重根，将迭代公式修改为：，而实际情况往往是不知道的，而如x*是f(x)的重根，则x*是f’(x)的重根，则x*是的单根了，因此对k(x)可以使用牛顿迭代的方法来求解，即． 八、离散牛顿法（也叫弦截法或割线法） 1、离散牛顿法的思想 设方程f(x)=0在区间[x0, x1]内有唯一实根，其解就是曲线与x轴的交点，而弦截法是选取曲线上两点，用过两点的直线（弦）与x轴的交点作为x*的一次近似值。 不妨假定f(x0)0, f(x1)0,，连接P0(x0 , f(x0)), P1(x1 , f(x1))的直线P0 P1(弦)与x轴的交点为x2，用为x*的一次近似，过P0 P1的直线为： ， 该直线与x轴交点（令y=0）得： 由此得到，再计算，点。 (1)、如，就是方程的解； (2)、如，再作弦继续求与x轴的交点； (3)、如，再与相邻作弦与x轴交点为；如此反复计算。 一般地，离散牛顿法的计算公式为： 每一步均加以判断。 2、离散牛顿法的计算流程 假定f(x)在 [a,b]上连续, (异号)，不妨设，置 对于n=1,2,... 并计算如或转出循环； 如，； 否则置继续循环 输出 由于离散牛顿法与第2节的单点迭代不同，其收敛性要复杂一些，但可以证明离散牛顿法是超线性收敛的，它的收敛阶为．即有下列定理： 定理4.2:设在方程的根的某个邻域内二阶连续可微,且,则存在,当时,由弦截法产生的数列收敛于，且该方法具有超线性收敛速度，它的收敛阶为．（证明参见关治，陆金甫《数值分析基础》第348页） 若为弦的固定一端时，则上述离散牛顿法变为 此迭代称为单点弦截法。可以证明此方法具有局部收敛性，且只有线性收敛速度。 x y X x0 x* x1 x2 Y y=f(x) X x0 x* x1 x3 x2 Y f(x)0 P0 P2 P1 X x0 x* x1 x2 x3 Y f(x)0 P0 P2 P1