第5章非线性方程的数值解法.ppt

* 例4 用不同方法求方程 的根 解 这里 ，可改写为各种不同的等价形式 其不动点为 由此构造不同的迭代法： §5.3 迭代法的收敛阶和加速收敛方法 * 取 ，对上述4种迭代法，计算三步所得的结果如下表. 表4 计算结果 * 从计算结果看到迭代法（1）及（2）均不收敛，且它们均不满足定理3中的局部收敛条件. 注意 . 迭代法（3）和（4）均满足局部收敛条件，且迭代法（4）比（3）收敛快，因在迭代法（4）中 . * 则称该迭代过程是p阶收敛的. 定义2 设迭代过程 收敛于方程 的根 ，如果迭代误差 当 时成立下列 渐近关系式 特别地, 时称线性收敛， 时称超线性收敛， 时称平方收敛. 一、迭代法的收敛阶 * 定理4 对于迭代过程 ，如果 在所 求根 的邻近连续，并且 则该迭代过程在点 邻近是p阶收敛的. （2.8） 上述定理说明，迭代过程的收敛速度依赖于迭代函数 的选取. 如果当 时 ，则该迭代过程只可能是线性收敛. * 在例4中，迭代法（3）的 ，故它只是线性 收敛. 而迭代法（4）的 ，而 由定理4知 ，该迭代过程为2阶收敛. * 回顾 二分法：将有根区间逐次减半得到区间序列｛[ak,bk]｝， 收敛于一点，即为方程的根；根据误差 控制迭代次数； 迭代法：根据迭代公式 ，得到迭代序列{xk}; 一般由迭代函数?(x)是否满足条件|?’(x)|1判断收敛; 收敛速度由收敛阶决定 线性收敛：?’(x)?0 二阶收敛： ?’(x)=0，?’’(x)?0 * §5.4 牛顿法（Newton Methods） 一、牛顿法的迭代公式 设已知方程 有近似根 （假定 )， 将函数 在点 展开，有 于是方程 可近似地表示为 （4.1） 牛顿法是一种线性化方法，其基本思想是将非线性方程 逐步归结为某种线性方程来求解. * 这是个线性方程，记其根为 ， 则 的计算公式为 （4.2） 这就是牛顿(Newton)法. * 牛顿法的几何解释(图4): 方程f(x)=0的根x*可解释为曲线y=f(x)与x轴交点的横坐标。 图4 设 是根 的某个近似值，过曲线 上横坐标为 的点 引切线，并将该切线与 轴的交点的横坐标 作为 的新的近似值. 切线方程为 这样求得的值 必满足（4.1）， 从而就是牛顿公式（4.2）的计算结果. 根据几何意义，牛顿法亦称切线法. （4.2） （4.1） * 求导得 牛顿迭代公式的迭代函数为 假定 是 的一个单根，即 ， 则由上式知 于是由定理4可得，牛顿法 在根 的邻近平方收敛. ， 二、牛顿法的收敛性： * 例1 用牛顿法解方程 （4.4） 解 这里牛顿公式为 取迭代初值 ，迭代结果列于表7中. 所给方程（4.4）实际上是方程 的等价形式. 若用不动点迭代到同一精度要迭代17次. 可见牛顿法的收敛速度是很快的. 表7 计算结果 * 牛顿法的计算步骤： 步骤1 准备 选定初始近似值x0，计算f0=f(x0), 步骤2 迭代 按公式 迭代一次得新的近似值x1， 计算 此处 是允许误差，而 其中 是取绝对误差或相对误差的控制常数， 步骤3 控制 如果x1 满足 或 ，则终止迭代，以x1 作为所求的根；否则转步骤4. 一般可取C=1 * 步骤4 修改 如果迭代次数达到预先指定的次数N，或者 ，则方法失败； 否则以 代替 转步骤2继续迭代. * 牛顿法应用举例 对于给定的正数 ，应用牛顿法解二次方程 可导出求开方值C的计算程序 （4.5） 这种迭代公式对于任意初值 都是收敛的. *