第五章非线性方程求根.pptVIP

2. 根的搜索 (1) 图解法（利用作图软件如 Matlab) 定理4 如果 在 附近的某个领域内有 ( )阶 连续导数，且 §3 牛顿法 一、牛顿法的迭代公式 考虑非线性方程 五、求m重根的牛顿法 1、迭代格式 可取 ，若记 则（5.3.4）式可改写为 称为松弛因子，这种方法称为松弛法。为使迭代速度加快， 需要边计算边调整松弛因子。由于计算松弛因子需要用到微商， 在实际应用中不便使用，具有一定局限性。若迭代法是线性收 敛的，当计算 不方便时，可以采用埃特金加速公式。 (2) 埃特金加速公式 设迭代法是线性收敛，由定义知 成立，故当 时有 由此可得 的近似值 （5.3.5） 由此获得比 和 更好的近似值 ，利用（5.3.5） 序列 的方法称为 (3) Steffensen 加速法 将Aitken加速公式与不动点迭代相结合，可得 （5.3.6） 式构造 埃特金（Aitken）加速方法。 利用（5.3.6）式构造序列 的方法称为Steffensen加速方法。 即每进行两次不动点迭代，就执行一次Aitken加速。 例2 试用简单迭代法和Steffensen加速法求方程 在 附近的根，精确至四位有效数。 解：记 ，简单迭代法公式为： 计算得 0.56719 13 0.56486 6 0.56707 12 0.57117 5 0.56728 11 0.56006 4 0.56691 10 0.57970 3 0.56714 16 0.56756 9 0.54524 2 0.56716 15 0.56641 8 0.60653 1 0.56712 14 0.55844 7 0.5 0 xk k xk k xk k Aitken加速公式 计算得 所以, 。 4．迭代过程的局部收敛 定义1： 若存在 的某一邻域 ，迭代过程 对任意初值 均收敛，则 称迭代过程 在根 邻近具有局 部收敛性。 定理3 设 为方程 的根， 在 的邻近 连续，且 ，则迭代过程 在 的邻近具有局部收敛性。 5．迭代过程的收敛速度 设由某方法确定的序列{xk}收敛于方程的根x*， 如果存在正实数p，使得 （C为非零常数） 定义2 则称序列{xk}收敛于x*的收敛速度是p阶的，或称该方法 具有p 阶收敛速度。当p = 1时，称该方法为线性（一次）收敛； 当p = 2时，称方法为平方（二次）收敛；当1 p 2或 C=0,p=1时，称方法为超线性收敛。 则迭代格式 在 附近是 阶局部收敛的，且有 原理：将非线性方程线性化— Taylor 展开 取 x0 ? x*，将 f (x)在 x0 做一阶Taylor展开: ，? 在 x0 和 x 之间。 将 (x* ? x0)2 看成高阶小量，则有： 只要 f ?C1，且每步迭代都有 ， 而且 则 x*就是 f (x)的根。公式（9.4.1）称为牛顿迭代公式。 (9.4.1) 构造迭代公式 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 则由（9.4.1）确定的牛顿迭代序列{xk}二阶收敛于f (x)在[a, b]上的唯一单根x*。 注：Newton法的收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 四. 牛顿迭代法的局部收敛性与收敛速度 ， ,且 设 在包含 的一个区间二阶连续可 导，则Newton迭代法至少二阶收敛，即 值得注意的是，当 充分光滑且 是 的重根时，牛顿 法在 的附近是线性收敛的。且Newton迭代法在 上的 收敛性依赖于初值 的选取。即初值 的选取充分靠近 时， 一般可保证Newton迭代法收敛。 并得