收敛加速方法.pptVIP

* 简单迭代法 不动点迭代的收敛性 迭代序列的收敛速度 收敛加速的方法 第二章 非线性方程的求根方法 [a,b]称为有根区间. 则 (2) (3) (1) f(ak) f(bk)0 由此可见，如果二分过程无限地进行下去（ ），则有限区间必定缩为一点x*，该点显然就是所求的根。 实际上，我们不可能去完成这种无穷过程，也无必要，只需得到满足一定精度的近似值就可以了。 如果令有根区间[an,bn]的中点 为 x*的近似值，则在二分过程中，得到下列以x*为极限的近似根序列 由于 二分法优点：是方程求根问题的一种直接搜索方法，算法简单、直观、实用，收敛性总能得到保证。 缺点（局限性）：不能求重根；计算速度慢。 思考：为什么不能求重根？ 例2.1 用二分法求方程 在区间 [1, 1.5］内的一个实根，要求误差不超过0.005。 解 由公式估计所要 二分的次数 即只要二分6次，便能达到所要求的精度。 计算结果 - 1.3242 1.3281 1.3203 6 - 1.3203 1.3281 1.3125 5 + 1.3281 1.3438 1.3125 4 + 1.3438 1.375 1.3125 3 - 1.3125 1.375 1.25 2 + 1.375 1.5 1.25 1 - 1.25 1.5 1.0 0 f(xk) xk bk ak k 作业：1、用二分法求方程 在区间[1,2］内的一个实根，要求误差不超过0.005。 将一个计算过程反复进行 一种常见常用的计算技术 构造有效的迭代格式 选取合适的迭代初值 对迭代格式进行收敛性分析 一种圆周率的计算方案: 初值: x0=1 ( n=1,2,3,······ ) 迭代格式: 2.2 迭代法 1 选取初值 把给定的方程 改写成等价形式 f (x)= 0 若存在 x*,使得 ,则称x*为不动点。 在根x*的附近取一点x0作为x*的预测值,也叫迭代初值。 (1) 把x0代入(1)的右端，得 如果 ，则 。 如果 ，把 x1作为根的新的预测值代入(1)，得 如果 ，则 。 如果 ，把 x2作为根的新的预测值代入(1)...... 如此重复上述步骤，则有迭代公式 ( k = 0, 1, 2, ··· ) 2 按迭代格式进行计算 3 判别收敛 其中， ：迭代函数，得到迭代序列 如果迭代序列的极限存在，则迭代过程收敛，显然有 如果迭代序列的极限不存在，则称迭代过程发散。 上述迭代过程也称不动点迭代法。 方程 求根，在几何上就是确定曲线 与直线 的交点 p* 几何意义 x* x2 x1 x0 如果 逐渐逼近p*，－－－迭代过程收敛 y = x y o x 如果 逐渐远离p*，－－－迭代过程发散 （无意义） x2 x1 x0 x* y = x y o x 例2.2 求方程 f (x)=x3 – x – 1 = 0 在x =1.5附近的根 x*。 解 设将方程改写成下列形式 由此得迭代公式 迭代初值取x0 =1.5，计算值用6位数字表示。 迭代结果如下表 1.32494 4 1.32472 8 1.32588 3 1.32472 7 1.33086 2 1.32473 6 1.35721 1 1.32476 5 1.5 0 xk k xk k 从表中可看到 x7与 x8完全相同，这时可认为x8已满足方程， x8 即为所求根的近似值。 上述迭代过程是收敛的。 如果将方程改写成下列形式 据此有迭代公式 迭代初值仍取 x0 =1.5，则有 当 k增大时，xk随之增大而不趋于任何极限，此时迭代过程发散。 通过此例说明，迭代过程只有在一定条件下才可能收