数值分析--第7章 非线性方程及方程组及数值解法.ppt

数值分析 第7章 非线性方程与方程组的数值解法 7.1 方程求根与二分法 7.1.2 二分法 7.2 不动点迭代法及其收敛性 7.2.2 不动点的存在性与迭代法的收敛性 7.2.3 局部收敛性与收敛阶 7.3 迭代收敛的加速方法 7.3.2 斯蒂芬森迭代法 7.4 牛顿法 7.4.2 牛顿法应用举例 7.4.3 简化牛顿法与牛顿下山法 （2） 牛顿下山法. 7.4.4 重根情形 7.5 弦截法与抛物线法 7.5.1 弦截法 7.5.2 抛物线法 迭代3次得到的结果 有6位有效数字. 但如果改用 作为迭代初值，则依牛顿法公式 （4.9）迭代一次得 这个结果反而比 更偏离了所求的根 . 为了防止迭代发散，对迭代过程再附加一项要求，即 具有单调性： （4.10） 满足这项要求的算法称下山法. （4.9） 将牛顿法与下山法结合起来使用，即在下山法保证函 数值稳定下降的前提下，用牛顿法加快收敛速度. 将牛顿法的计算结果 与前一步的近似值 适当加权平均作为新的改进值 （4.11） 其中 称为下山因子， （4.12） (4.12)称为牛顿下山法. （4.11）即为 选择下山因子时从 开始，逐次将 减半进行试算， 若用此法解方程（4.8），当 时由（4.9）求得 直到能使下降条件（4.10）成立为止. ，它不满足条件（4.10）. 通过 逐次取半进行试算，当 时可求得 此时有 ， 显然 . 而 由 计算 时 , 均能使条件（4.10） 成立. 计算结果如下 ： （4.10） （4.9） （4.8） 即为 的近似.一般情况只要能使条件（4.10）成立， 则可得到 ，从而使 收敛. （4.10） 设 ，整数 ， 则 为方程 的 重根，此时有 只要 仍可用牛顿法（4.2）计算，此时迭代函数 的导数为 且 ，所以牛顿法求重根只是线性收敛. （4.2） 则 . （4.13） 求 重根，则具有2阶收敛，但要知道 的重数 . 构造求重根的迭代法，还可令 ， 若 是 的 重根，则 若取 用迭代法 从而可构造迭代法 （4.14） 它是二阶收敛的. 故 是 的单根. 对它用牛顿法，其迭代函数为 例9 方程 的根 是二重根， 用上述三种方法求根. 解 （1） 牛顿法 先求出三种方法的迭代公式： （2） 用（4.13）式 （3） 用（4.14）式 取初值 ，计算结果如表7-9. （4.13） 从结果看出，经过三步计算，方法（2）及（3）均达到 10位有效数字，而由于牛顿法只有线性收敛，所以要达到同 样精度需迭代30次. 当函数 比较复杂时，计算 往往较困难， 值 的计算. 为此可以利用已求函数值 来回避导数 还要算 . 用牛顿法求方程 的根，每步除计算 外， 设 是 的近似根， 利用 构造一次插值多项式 ，并用 的根作为新的 近似根 . （5.1） 由于 因此有 （5.2） （5.2）可以看做将牛顿公式 中的导数 用差商 取代的结果. 接着讨论几何意义. 曲线 上横坐标为 的点分别记为 ， 则弦线 的斜率等于差商值 , （5.2） 按（5.2)式求得的 实际上是弦线 与 轴交点的横坐标. 表7-6 这种算法因此而称为弦截法. 其方程为 （5.2） 而弦截法（5.2），在