牛顿迭代法111.pptVIP

实验五 牛顿迭代法 * * 【实验目的】 1.了解牛顿迭代法的基本概念. 2.了解牛顿迭代法的收敛性和收敛速度。 3.学习、掌握MATLAB软件有关的命令。 【实验内容】 用牛顿迭代法求解方程 的近似根，误 差不超过10-3。 【实验准备】 1.牛顿迭代法原理 设已知方程 的近似根 ，则在 附近 可用一阶泰勒多项式 近似代替。因此，方程 可近似表示为 。用 近似表示 根差异不大。 设 ，由于 满足 ，解得 重复这以过程，得到迭代格式 这就是著名得牛顿迭代公式，它相应的不动点方程为 牛顿迭代法 2.牛顿迭代法的几何解析 在 处做曲线的切线，切线方程为 令 可得切线与 轴的交点坐标 ，这 就是牛顿迭代法的迭代公式。因此，牛顿法又称“切线法”。 3.牛顿迭代法的收敛性 计算可得 ，设 是 的单根， 有 ， 则 故在 附近，有 。根据不动点原理知牛顿迭 代法收敛。 4.牛顿迭代法的收敛速度 定理（牛顿法收敛定理）设 在区间 上有二 阶连续导数，且满足 ， 在 上不变号， 在[a,b]上不等于0，令 有 ，则对任意 ，牛顿迭代格式收敛于 在 中的唯一实根 ，且 ， 牛顿迭代法为二阶收敛。 对不动点方程 ，它导出的迭代过程有可能发散，也可能收敛得非常缓慢。这时，我们有没有办法改进不动点方程，让迭代过程收敛得快一些呢？ 迭代过程的加速 （1）一个简单办法 注意到 和 都是不动点方程，它们的加权平均 也是不动点方程，而且 和 有完全相同的不动点。适当选取 的值，可以使发散的迭代过程变得收敛，使收敛慢的迭代过程变得收敛迅速。 （2）加速的原因 在下面的实验中我们可以看到， 在不动点 附近的导数值在很大程度上决定了迭代过程的收敛性。 的绝对值越小，收敛性越好。因此，选择 使得 。计算得到理想的 值为 ，相应可计算出 （1）一个简单办法 注意到 和 都是不动点方程，它们的加权平均 也是不动点方程，而且 和 有完全相同的不动点。适当选取 的值，可以使发散的迭代过程变得收敛，使收敛慢的迭代过程变得收敛迅速。 （2）加速的原因 在下面的实验中我们可以看到， 在不动点 附近的导数值在很大程度上决定了迭代过程的收敛性。 的绝对值越小，收敛性越好。因此，选择 使得 。计算得到理想的 值为 ，相应可计算出 （3） 的选取 由于理想的 值为 ，当 变化不大时，可以取 近似计算 。 （4）回到牛顿迭代法的讨论 为求解方程 ，可以使用不动点方程 ，相应的迭代函数为 对 进行加速 所以，牛顿迭代法是对基本迭代格式进行加速的结果。 5.迭代的MATLAB命令 MATLAB中主要用for,while等控制流命令实现迭代。 练习 1 用牛顿迭代法