第三节牛顿迭代法.pptVIP

三 简化牛顿法与牛顿下山法 （2） 牛顿下山法. 四 重根情形 五 弦截法与抛物线法 2 抛物线法 7.6 解非线性方程组的牛顿迭代法 拟Newton方程 1. 非线性方程组求根问题是前面介绍的方程（即 ） 求根的直接推广，只要把前面介绍的单变量函数 看 成向量函数 则可将单变量方程求根方法推广到方程 组（6.2）. 若已给出方程（6.2）的一个近似根 ， 将函数 的分量 在 用多元函数泰 勒展开，并取其线性部分，则可表示为 令上式右端为零，得到线性方程组 （6.3） 其中 （6.4） 称为 的雅可比(Jacobi)矩阵. 求解线性方程组（6.3）并记解为 ，则得 （6.5） 这就是解非线性方程组（6.2）的牛顿迭代法. 例12 求解方程组 给定初值 , 用牛顿法求解. 解 先求雅可比矩阵 由牛顿法（6.5）得 即 由 逐次迭代得到 的每一位都是有效数字. * * 一 牛顿法及其收敛性 牛顿法是一种线性化方法，其基本思想是将非线性方 程 逐步归结为某种线性方程来求解. 设已知方程 有近似根 （假定 )， 将函数 在点 展开，有 于是方程 可近似地表示为 （1） 这是个线性方程，记其根为 ，则 的计算公式为 10.4 牛顿迭代法 （2） 这就是牛顿(Newton)法. 牛顿法的几何解释. 方程 的根 可解释为曲线 与 轴 的交点的横坐标（图7-3). 设 是根 的某个近似值， 过曲线 上横坐标为 的点 引切线，并将该切线与 轴的交点的横坐标 作为 的新的近似值. 图7-3 注意到切线方程为 这样求得的值 必满足(1），从而就是牛顿公式（2） 的计算结果. 由于这种几何背景，牛顿法亦称切线法. 牛顿法（2）的收敛性，可直接由上节定理得到，对（2） 其迭代函数为 由于 假定 是 的一个单根，即 ， 则由上式知 ，于是依据可以断定， 牛顿法在根 的邻近至少是平方收敛的. 又因 故 可得 （3.3） 例7.3.1 用牛顿法解方程 （3.4） 解 这里牛顿公式为 取迭代初值 ，迭代结果列于表7-5中. 所给方程（3.4）实际上是方程 的等价形式. 若用不动点迭代到同一精度要迭代28次，可见牛顿法的收敛速度是很快的. 对于给定的正数 ，应用牛顿法解二次方程 可导出求开方值 的计算程序 （3.5） 这种迭代公式对于任意初值 都是收敛的. 事实上，对（3.5）式施行配方手续，易知 二 牛顿法应用举例 以上两式相除得 据此反复递推有 （3.6） 记 整理（3.6）式，得 对任意 ，总有 ，故由上式推知，当 时 ，即迭代过程恒收敛. 解 取初值 ，对 按（3.5）式迭代3次 便得到精度为 的结果 （见表7-6). 由于公式（3.5）对任意 初值 均收敛，并且收 敛的速度很快，因此可取确定 的初值如 编成通用程序. 例7.3.2 求 . 牛顿法的优点 收敛快， 牛顿法的缺点 一 每步迭代要计算 及 ，计算量较大 且有时 计算较困难， 二是初始近似 只在根 附近才能保证收敛， 如 给的不合适可能不收敛. 为克服这两个缺点，通常可用下述方法. （1） 简化牛顿法，也称平行弦法. 其迭代公式为 （3.7） 迭代函数 若在根 附近成立 ，即取 ，则迭代法（3.7）局部收敛. 在（3.7）中取 ，则称为简化牛顿法，这 类方法计算量省，但只有线性收敛，其几何意义是用平行 弦与 轴交点作为 的近似. 如图7-4所示. 图7-