2019-牛顿迭代法-文档资料.ppt

10.4 牛顿迭代法 一 牛顿法及其收敛性 牛顿法是一种线性化方法，其基本思想是将非线性方 f ( x ) 逐步归结为某种线性方程来求解 ? 0 程 . f ( x ) ? 0 f ? x k ) ， 设已知方程 有近似根 （假定 ( x k ) ? 0 x k f ( 在点 x ) 展开，有 将函数 ? f ( x ) ? f ( x ) ? f ( x )( x ? x ), k k k f ( x ) 可近似地表示为 ? 0 于是方程 ? f ( x ) ? f ( x )( x ? x ) ? 0 . k k k （ 1 ） x k ? 1 x k ? 这是个线性方程，记其根为 ，则 的计算公式为 1 1 f ( x ) k x ? x ( k ? 0 , 1 , ? ), k ? 1 k ? ? f ( x ) k （ 2 ） 这就是 牛顿 (Newton) 法 . 牛顿法的几何解释 . f ( x ) ? 0 y ? ) x * 方程 的根 可解释为曲线 与 轴 f ( x 的交点的横坐标（图 7-3). x x * x k 的某个近似值， 设 是根 y ? f ( x ) 过曲线 上横坐标为 x k 的点 x P 引切线，并将该切线与 k x k 作为 ? 1 轴的交点的横坐标 x * 的新的近似值 . 图 7-3 2 注意到切线方程为 ? y ? f ( x ) ? f ( x )( x ? x ). k k k x k ? 这样求得的值 必满足 (1 ），从而就是牛顿公式（ 2 ） 1 的计算结果 . 由于这种几何背景，牛顿法亦称 切线法 . 牛顿法（ 2 ）的收敛性，可直接由上节定理得到，对（ 2 ） 其迭代函数为 f ( x ) g ( x ) ? x ? , f ? ( x ) 由于 ? f ( x ) f ? ( x ) ? g ( x ) ? . 2 [ f ? ( x )] ? f ( x *) ? 0 , f ( x *) ? 0 f 的一个单根，即 ( x ) 假定 ， x 是 * ? g ( x *) ? 0 则由上式知 ，于是依据可以断定， x * 牛顿法在根 的邻近至少是平方收敛的 . 3 又因 e k ? 1 g ( p ) ( x *) ? . 故 可得 p e k p ! ? ( x f ? *) ? ? ( x g *) ? , f ? ( x *) e k ? 1 g ( p ) ( x *) ? . p e k p ! ? x x * f ? ( x *) k ? 1 ? l i m ? . 2 k ? ? ( x x *) 2 f ? ( x *) k ? （ 3.3 ） 例 7.3.1 用牛顿法解方程 x x e ? 1 ? 0 . （ 3.4 ） ? x k 解 这里牛顿公式为 x e k ? x k ? 1 ? x k ? 1 ? x k , x 0 ? ，迭代结果列于表 0 . 5 取迭代初值 7-5 中 . 4 x ? e 所给方程（ 3.4 ）实际上是方程 的等价形式 . 若用不 动点迭代到同一精度要迭代 28 次，可见牛顿法的收敛速度 是很快的 . 表 7 k 0 1 2 3 ? 5 计算结果 x k 0 . 5 0 . 57102 0 . 56716 0 . 56714 ? x 5 二 牛顿法应用举例 对于给定的正数 ，应用牛顿法解二次方程 C x ? C ? 0 , 2 可导出求开方值 的计算程序 C 1 C x ( x ). k ? 1 ? k ? 2 x k x 0 ? 0 这种迭代公式对于任意初值 都是收敛的 . （ 3.5 ） 事实上，对（ 3.5 ）式施行配方手续，易知 x k ? 1 ? x k ? 1 ? 1 C ? ( x k ? 2 x k 1 C ? (