第2章非线性方程(组)的数值解法2解析.ppt

* 图2-6 简化牛顿法几何示意 * 例2-14 在 附近的一个根 . 解 1） 设取迭代初值 ，牛顿法公式为 迭代3次得到的结果 有6位有效数字. 作为迭代初值，则依牛顿法公式（2.23）迭代一次得 .这个结果反而比 更偏离了所求的根 计算得 2） 如果改用 用牛顿法求解方程 * 通常,牛顿迭代法的收敛性依赖于初始值的选取,如果初始值偏离所求的根比较远,则牛顿法可能发散。为了防止迭代发散,我们对牛顿迭代法的迭代过程再附加一项要求,即具有单调性 满足这项要求的算法称下山法。 2、牛顿下山法 牛顿下山法原理：将牛顿迭代法与下山法结合起来使用,即在下山法保证函数值下降的前提下,用牛顿迭代法加快收敛速度。把这一算法称为牛顿下山法。即 其中λ 0 λ 1 为下山因子 牛顿法计算的结果： 与前一步的近似值xk适当加权平均作为新的改进值 下山因子的选择是个逐步探索的过程，设从λ 1开始反复将λ减半进行试算, 即逐次取λ为 下山因子的选择： 从中挑选下山因子，直至找到其中某个λ使单调性条件 成立，则称“下山成功”，否则“下山失败”， 这时需另选初值重算。 * 例14中取 逐次取半进行计算， 时可求得 ，此时 而 ，显然 由 计算 时 均能使下降条件成立，计算结果如下： ， ， 即为 的近似. 通过 当 由 计算 ，它不满足下降条件。 * 一般情况只要使下降条件成立，则可 从而使 牛顿下山法是目前方程求根的一个 速度虽然没有牛顿法快，但它对初 值的选取范围放宽很多. 收敛. 有效方法，当 时， 得到 它的收敛 * 2.5.3 割线法 牛顿迭代法的收敛速度很快，但是每迭代一次，都要计算 的值，如果出现 接近于零，可能 不为零，但当 时，导数的计算工作量也较大.通常希望避免计算导数， 而导出一个类似牛顿迭代的公式. 导致溢出.既使 较复杂 割线法（Secant Method）原理： 时，由前面算出的 两点的割线斜率代替导数，可得 在牛顿迭代中，为避免计算导数，用经过两点的割线来代替切线，即求 * 代入牛顿迭代公式 中，有割线法迭代公式 . 其割线方程为 ， 割线与 的交点记为 则得到割线法迭代公式.割线法又称为弦截法、弦线法. 过两点 割线法的几何意义： 迭代公式的导数是曲线 的割线（弦线）的斜率， * 图2-7 割线法几何示意 * 如果 在零点 附近有 ，且初始值 充分接近 ，则割线法的迭代过程是收敛的，其收敛速度为 其中 分别是第 次和第 【注】 割线法具有超线性收敛速度，其收敛阶为 连续的二阶导数， 次迭代近似值. 割线法的收敛性： * 用割线法求方程 解 割线法迭代公式为 在1.5附近的一个根.保留四位有效数字。 取初始值 迭代计算，得到 . 例2-15 * 用割线法求方程 在0.5附近. 用割线法迭代计算，得到 比较例2-13牛顿法的计算结果可以看出， 割线法的收敛速度也是相当快的. 解 取初始值 的根 例2-16 例2-13中 * 2.5.4 求重根的修正牛顿法 ，整数 ， 则 为方程 的 重根，此时有 只要 仍可用牛顿法计算，此时迭代函数 的导数为 且 ，所以牛顿法求重根只是线性收敛. 设 * 一、求重根的修正牛顿法1 ，则 因此利用迭代法 , 求 重根，具有2阶收敛.应用该方法的 的重数 . 选取 （2.28） 前提是：要知道 * 二、求重根的修正牛顿法2 若 是 的 重根，则 ，（2.29） 构造求重根的迭代法，还可令 故 是 的单根.对它用牛顿法，其迭代函数为 从而可构造迭代法 ， 它是二阶收敛的. * 方法1） 方法2） 方法3） 1 1.座机电话号码3 1.座机电话号码7 1.座机电话号码6 2 1.座机电话号码3 1.座机电话号码6 1.座机电话号码8 3 1.座机电话号码9 1.座机电话号码2 1.座机电话号码2 表2-8　三种方法数值结果 计算三步，方法2）及方法3）均达到10位有效数字，而用牛顿法只有线性收敛，要达到同样精度需要迭代30次. * 【注】 求重根的修正牛顿法1需要已知根的重数，因此不实用.求重根的修正牛顿法2需要求函数的二阶导数，并且当所求根为单根时，不能改善本来已经二次收敛的牛顿法.对于实际问题，往往事先并不知道所求根是否是重根，需要通过试算来判断，如当牛顿法收敛很慢时通常为重根. 小结： 一、牛顿迭代法： 1、牛顿迭代函数： 2、牛顿迭代公式： 3、牛顿迭代收敛性： 定理7 设 是方程 的单根, 且f x 在 的某邻域内有连续的二阶导数, 则牛顿法在 附近局部收敛, 且至少二阶收敛, 有 1、简化牛顿法： 迭代公式： 迭代函数： 即取 若 在根 附近成立，则 迭代法（2.22）局部收敛. 在（2.22）中取 ，则称为简化牛顿法， 二、简化牛顿法与牛顿下山法 * 例17