数值分析课件 第7章4-6节.ppt

7.4 牛顿法 7.4.2 牛顿法应用举例 7.4.3 简化牛顿法与牛顿下山法 （2） 牛顿下山法. 7.4.4 重根情形 7.5 弦截法与抛物线法 7.5.1 弦截法 7.5.2 抛物线法 7.6 解非线性方程组的牛顿迭代法 插值多项式 有两个零点： （5.3） 式中 问题是该如何确定 . 假定在 三个近似根中， 更接近所求的根 . 为了保证精度，选（5.3）中较接近 的一个值作为新的近似根 . 为此，只要取根式前的符号与 的符号相同. 例11 用抛物线法求解方程 解 设用表7-8的前三个值 作为开始值，计算得 （5.3） 故 代入（5.3）式求得 以上计算表明，抛物线法比弦截法收敛得更快. 在一定条件下可以证明，对于抛物线法，迭代误差有 下列渐近关系式 （5.3） 可见抛物线法也是超线性收敛的，其收敛的阶 , 从（5.3）看到，即使 均为实数， 也 可以是复数，所以抛物线法适用于求多项式的实根和复根. 收敛速度比弦截法更接近于牛顿法. （5.3） 考虑方程组 （6.1） 其中 均为 的多元函数. 用向量记号记 , （6.2） (6.1）就可写成 7.4.1 牛顿法及其收敛性 设已知方程 有近似根 （假定 )， 将函数 在点 展开，有 于是方程 可近似地表示为 （4.1） 牛顿法是一种线性化方法， 逐步归结为某种线性方程来求解. 其基本思想是将非线性方程 这是个线性方程，记其根为 ， 则 的计算公式为 （4.2） 这就是牛顿(Newton)法. 牛顿法的几何解释. 方程 的根 可解释为 曲线 与 轴的交点的横坐标 图7-3 （图7-3). 设 是根 的某个近似值，过曲线 上横坐标为 的点 引切线，并将该切线与 轴的交点的横坐标 作为 的新的近似值. 注意到切线方程为 这样求得的值 必满足（4.1），从而就是牛顿公式 （4.2）的计算结果. 由于这种几何背景，牛顿法也称切线法. 由定理4，可以直接得到牛顿法的收敛性. （4.2） （4.1） 由于 假定 是 的一个单根， 即 ， 则由上式知 (4.2)的迭代函数为 在根 的邻近是平方收敛的. 于是依据定理4可以断定，牛顿法 又因 （4.2） 故由（2.9）可得 （4.3） 例7 （4.4） 解 取迭代初值 ，迭代结果列于表7-5中. 用牛顿法解方程 这里牛顿公式为 所给方程（4.4）实际上是方程 的等价形式. （2.9） 牛顿法的计算步骤： 步骤1 准备 选定初始近似值 ， 步骤2 迭代 迭代一次， 得新的近似值 ， 若用不动点迭代到同一精度 可见牛顿法的收敛速度是很快的. 要迭代17次. 计算 按公式 计算 则终止迭代，以 作为所求的根； 此处 是允许误差，而 其中 是取绝对误差或相对误差的控制常数， 步骤4 修改 如果迭代次数达到预先指定的次数 ， 步骤3 控制 如果 满足 或 ， 否则转步骤4. 一般可取 或者 ，则方法失败； 否则以 代替 转步骤2继续迭代. 对于给定的正数 ，应用牛顿法解二次方程 可导出求开方值 的计算程序 （4.5） 这种迭代公式对于任意初值 都是收敛的. 事实上，对（4.5）式施行配方手续，易知 以上两式相除得 据此反复递推有 （4.6） 记 对任意 ， 总有 ， 时 ， 解 取初值 ，对 按（4.5）式迭代3次 便得到精度为 的结果 例8 求 . 整理（4.6）式，得 故由上式推知， 即迭代过程恒收敛. 当 （见