数值分析非线性方程求根.PPT

* * Newton下山法 （续） Remark1：由数值试验可知，用Newton下山法求解方程的根，需要一定的经验。当函数形状复杂时， Newton下山法可能收敛很慢。只有当xk充分接近x*时，才具有二阶收敛速度。 Remark3： Newton法是一种特定的不动点迭代法，因此它可以用于求方程f(x)=0的复根。 Remark2：对形状复杂的函数使用Newton法和Newton下山法时，程序中应该加上 是否为零的判断，以排除 不能进行迭代的可能。 * * 三、弦割法 为了得到超线性收敛且避免导数计算的迭代格式，常用函数增量与自变量增量之比来近似替代Newton迭代格式中的导数项。 取 用上述增量比代替Newton迭代格式中的一阶导数，则有 * * 用弦AB代替 ，用弦AB与x轴交点C的横坐标 作为 的新的近似根。再用弦AC与X轴交点的横坐标 作为 的新的近似根。 弦割法求 要用到 开始必须先给出 。 此迭代过程为弦割法，又称为拟牛顿法。弦割法求 的过程，相当于过A、B两点作直线：　　 用来描述 弦割法的迭代格式（续） * * 局部收敛性定理 设 在 的某邻域内连续， 则存在 的一个邻域　　 ， 当 时，弦割法产生的序列 收敛 于 ，且收敛阶为 。　 定理： * * 设 的等价方程为 对 在 使用弦割法，求出的值记为 ，则 这就是Steffensen迭代法。 弦割法与Steffensen迭代法有密切的关系。 对于弦割法的修正 * * 修正的弦割法（续） 牛顿法和弦割法的另一种修改方式为：对于给定的 令 ，则 与 等价。对如此特定的 应用Steffensen迭代，有 该迭代式也称为Steffensen迭代，它也可看作是对Newton法和弦割法的一种修改。这一迭代公式既不要计算Newton法中的导数 ，也不需弦割法中的两个初值，它是一种单点迭代法。 * * 当 在 附近满足 连续的条件时，则修正的弦割法在 附近是平方收敛的。该迭代法也称为Steffensen迭代法。 代入Steffensen迭代格式　 即修正的Newton迭代法。 即 修正的弦割法（续） 事实上令 ，则有　 * * 局部收敛性及收敛阶（续） 而对于平方收敛的 ，由 可得 ，因此大约需要7次迭代。 结论：平方收敛的序列收敛速度要比线性收敛的序列大的多。 Remark:不动点迭代的方法是以求解非线性方程的实根来讨论的，但类似的结果完全可以推广到求方程的复数根的情形。 * * 四、不动点迭代的加速 设方程为 ， 的某个预测值为 校正两次 由于 ?1在 与 之间 ?2在 与 之间。 在 邻域中，消去导数信息，即 即 因此可以得下述Steffensen迭代方法： * * 给出 ，校正 ，再校正 或 若记 ,k=0,1,2,…..,上式也可以写成。 上式称为Steffensen迭代格式。 改进 k=0,1,2,…. 不动点迭代的加速（续） * * 定理：设函数 有不动点 ，且在 的邻域内具有二阶连续导数，若 且 则Steffensen迭代式为二阶收敛，而且其极限为 。 证明 设Aitken或Steffensen迭代式为 k=0,1,2,… 则 对于Steffensen迭代，有以下定理： 不动点迭代的加速（续） * * (1)验证: