[数值分析]第八章 非线性方程求根课件.pptVIP

知当 时， 所以迭代格式（2）是收敛的。 迭代格式（3）的迭代函数为 当 时 由 时， 知当 所以迭代格式（3）也是收敛的。 结论: 通过以上算例可以看出对迭代函数 所得到的 若小于1，则收敛；且上界越小收敛速度越快。 求导， 的上界若是大于1，则迭代格式发散； 4．迭代过程的收敛速度 设由某方法确定的序列{xk}收敛于方程的根x*， 如果存在正实数p，使得 （C为非零常数） 定义： 则称序列{xk}收敛于x*的收敛速度是p阶的，或称该方法 具有p 阶收敛速度。当p = 1时，称该方法为线性（一次）收敛； 当p = 2时，称方法为平方（二次）收敛；当1 p 2或 C=0,p=1时，称方法为超线性收敛。 5.迭代法收敛阶的判别 定理3 如果 在 附近的某个领域内有p( )阶 连续导数，且 则迭代格式 在 附近是p阶局部收敛的。 §3 牛顿法 一、牛顿法的迭代公式 考虑非线性方程 原理：将非线性方程线性化— Taylor 展开 取 x0 ? x*，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间。 将 (x* ? x0)2 看成高阶小量，则有： 只要 f ?C1，且每步迭代都有f ’( xk ) ? 0， 而且 则 x*就是 f (x)的根。公式（2.3）称为牛顿迭代公式。 (2.3) 构造迭代公式 x* x0 x1 x2 x y f(x) 二、牛顿法的几何意义 三、牛顿法的收敛性 定理4： 设f (x)在[a, b]上存在二阶连续导数且满足下列条件： （1）f (a) ? f (b) 0； （2）f’ (x) ? 0； （3）f? (x) 在区间[a,b]上不变号； （4）取x0? [a, b]，使得f? (x)?f (x0) 0 则由（2.3）确定的牛顿迭代序列{xk}二阶收敛于f (x)在[a, b]上 的唯一单根x*。 注：Newton法的收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 四、求m重根的牛顿法 1、迭代格式 (2.4) 2、重数m的确定 3、迭代格式（2.4)的收敛阶（至少2阶收敛） 五、牛顿法的变形--牛顿下山法 1、牛顿下山法的计算步骤： （1）选取初始近似值x0； （2）取下山因子? = 1； （3）计算 （4）计算f (xk+1)，并比较 与 的大小，分以下二种情况 1）若 ，则当 时，取x*? xk+1，计算过程结束；当 时，则把 xk+1 作为新的 近似值，并返回到（3）。 2）若 ，则当?≤??且|f(xk+1)| ，取x*? xk，计算过程结束； 否则若?≤??，而 时，则把xk+1加上一个适当选定的小正数， 即取xk+1+?作为新的xk值，并转向（3）重复计算；当?＞??；且 时 ，则将下山因子缩小一半，取?/2代入，并转向（3）重复计算。 例5：求方程f (x) = x3 – x – 1 = 0 的根。 k ? xk 0 1 0.6 1 1/25 1.14063 2 1 1.36681 3 1 1.32628 4 1 1.32472 牛顿下山法的计算结果： 牛顿法一步要计算 f 和 f ’，相当于2个函数值。现用 f 的值近似 f ’ : x0 x1 切线 割线 切线斜率　?　割线斜率 x2 2、割线法： 3、 双点弦截法 ： 切线斜率　?　割线斜率 初值 x0 和 x1。 x0 x1 x2 加速收敛技术 由定理9.3.2知L越小迭代法的收敛速度越快，因此， 可以从寻找较小的L来改进迭代格式以加快收敛速度。 1. 松弛法 引入待定参数 ，将 作等价变形为 （9.3.4） 将方程右端记为 ，则得到新的迭代格式 由压缩映像定理知 ，为了使新的迭代格式比原来迭 代格式收敛得更快，只要满足 且 越小，所获取的L就越小， 迭代法收敛的就越快，因此我们希望 。 可取 ，若记 ，则（9.3.4）式可改写为 称为松弛因子，这种方法称为松弛法。为使