第7章 非线性方程的数值解法.pptVIP

第7章 非线性方程与方程组的数值解法/* Numerical Solutions of Nonlinear Equations*/ 7.1 方程求根与二分法 7.2 不动点迭代法及其收敛性 7.2.2 不动点的存在性与迭代法的收敛性 7.2.3 局部收敛性与收敛阶 7.3 迭代收敛的加速方法/* Accelerating Method*/ 7.3.2 斯蒂芬森迭代法 7.4 牛顿法 7.4.3 简化牛顿法与牛顿下山法 –牛顿法的改进 7.4.4 重根情形 7.5 弦截法与抛物线法 作业 (1) P238 1,2,3 (2) P238 5,8 (3) P238 7,9,12 7.3.1 埃特金加速收敛方法 设 是根 的某个近似值，用迭代公式迭代一次得 由微分中值定理 其中 介于 与 之间. 假定 改变不大， （3.1） 若将校正值 再迭代一次，又得 有 一般情形 （3.2） 埃特金（Aitken） 加速方法 可以证明 它表明序列 的收敛速度比 的收敛速度快. 假定 改变不大， 有 Aitken 加速： x y y = x y = (x) x* x0 P(x0, x1) x1 x2 P(x1, x2) 比 收敛得略快。 将 视为新的初值， 重复上述步骤 ? Steffensen 加速： 斯蒂芬森(Steffensen)迭代法 （3.3） 要求 的根 不动点迭代 （3.4） 其中 （3.5） 定理5 若 为 定义的迭代函数 的不动点，则 为 的不动点. 反之，若 为 的不动点，设 存在， ， 则 是 的不动点，且斯蒂芬森迭代法是2阶收敛的. 迭代 取 . 例8 用斯蒂芬森迭代法求解方程 计算结果如下表. 解 是发散的. 用斯蒂芬森迭代法 斯蒂芬森加速可使原本不 收敛的迭代改进到收敛. 例9 求方程 在 中的解. 由方程得等价形式 ，取对数得 构造迭代法 且当 时， ， 由于 , 根据定理2此迭代法是收敛的. 若 为 阶收敛，则 为 阶收敛. 解 若取 迭代16次得 ，有六位有效数 字. 若用斯蒂芬森加速，计算结果如下 ： 这里计算2步结果与 相同， 说明用迭代法（3.3）的收敛速度比迭代法（2.2）快得多. 注意：对本来就是P(1)阶收敛的方法，改用Stefensen迭代方法优点不多。 7.4.1 牛顿法及其收敛性 设已知方程 有近似根 （假定 )， 将函数 在点 展开，有 于是方程 可近似地表示为 （4.1） 这是个线性方程，记其根为 ， 则 的计算公式为 （4.2） 牛顿法 几何意义 切线法 收敛性 （4.2） x y x* x0 只要 f ?C1，每一步迭代都有f ’( xk ) ? 0， 而 ，则 x*就是 f 的根。 迭代函数 假定 是 的一个单根，即 ， 则由上式知 局部平方收敛 （4.3） 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 §6.1.4 牛顿迭代法 例10 用牛顿法解方程 （4.4） 解 取迭代初值 ，迭代结果列于表7-7中. 所给方程（4.4）实际上是方程 的等价形式. 牛顿公式为 若用不动点迭代到同一精度 要迭代17次. 可见牛顿法的收敛速度是很快的. 牛顿法的计算步骤： 步骤1 准备 选定初始近似值 ，计算 步骤2 迭代 按公式 迭代一次， 得新的近似值 ， 计算 步骤3 控制 如果 满足