（数值分析）第七章 非线性方程求根.ppt

华长生制作 重复上述步骤，第k步就得到根 的近似 值序列 及包含根 的区间套，且有 总结上述讨论有下列定理 定理 给定方程f(x)=0，设f(x)在[a,b]上连 续，且f(a)f(b)<0.则由二分法产生的序列 收敛于方程的根 ，且具有误差估计 局部收敛性和收敛阶 定理 设不动点迭代的迭代函数 在其 不动点 的某邻域内具有二阶连续导 数， 则斯蒂芬森的迭代 技术是二阶收敛的，且极限仍为 定理2 给定方程f(x)=0且 ，如果满 足条件： （1） （2） （3） 则由牛顿迭代法产生的序列 收敛于方程的惟 一实根 ，且有 定理3 设 在[a,b]上连续，且 则对 牛顿迭代序列 收敛 于方程f(x)=0在[a,b]内的唯一实根 这种迭代对原迭代的各近似值在根的两侧摆动的情 形较为有效。 使用三个迭代值的组合方法 中值定理 差商近似代替导数 即 于是可以得到迭代格式: 其中 上组公式称为Altken公式或Altken加速 将上式综合后可得一个解析式表示的迭代法: 或写成Steffensen迭代法 Altken公式与Steffensen公式是等价的 加速效果也是很明显的 例. 对迭代格式 进行加速解方程组 解: x0 = 0.5 x1 = 0.375 x2 = 0.3509115 x3 = 0.3477369 x4 = 0.3473496 x5 = 0.3473028 x6 = 0.3472971 x7 = 0.3472964 (1)直接使用迭代格式 迭代7次,得到满足精度的解 (2)对迭代格式进行Altken加速 x0 = 0.5 x1 = 0.3451613 x2 = 0.3472961 x3 = 0.3472964 迭代3次,得到满足精度的解 从以上3种结果可见,迭代法加速技术效果比较明显 迭代格式 显然不收敛 x0 = 1.5 x1 = 1.5350706 x2 = 1.5321124 x3 = 1.5320889 x4 = 1.5320889 迭代4次,得到满足精度的解 对迭代格式进行Altken加速 可见加速技术可能将不收敛的迭代法加速为收敛 7.4 牛顿迭代法 标准Newton迭代法及其收敛阶 如果将非线性方程 令 化为等价方程 如果 令 即 则 于是取 --------(1) --------(2) (2)式称为Newton迭代法 由前面分析可知有 Newton迭代法至少平方收敛 局部收敛性 对牛顿法可作如下的几何解释： 为函数f(x)在点 处的切线与横坐标轴的交点,见图.因此Newton迭代法也称为切线法. Y 0 y=f(x) X 例1. 用Newton迭代法求方程的根: 解: 由Newton迭代法 x0 =0.5; x1 =0.3333333333 x2 =0.3472222222 x3 =0.3472963532 x4 =0.3472963553 迭代四次 精度达10-8 Newtonddf.m 　　以上讨论的是Newton法的局部收敛性。对于某些非线性方程，Newton法具有全局收敛性。 初始值的选取可更一般化。 　　例 设a>0,对方程 -a=0试证:取任何初值 　>0，Newton迭代法都收敛到算术根 。 由此可知 证 对f(x)= -a, Newton迭代法为 可见,对于任何 >0,都有 ,并且{ }非增.因此{ }是有下界的非增序列,从而有惟一极限x*.在 内 故对任何 有迭代序列 都平方收敛于 即x*= 。 故有 且 对于Newton迭代法 趋于零 Newton迭代法也可能只是线性收敛 此时Newton迭代法可能不收敛 从而， 时只要 ，这时的Newton迭代法线性收敛。 　　为了改善重根时Newton法的收敛性，有如下两种方法。 若改为取 容易验证 ，故此时迭代法至少二阶收敛. 另一方案是令