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

一、使用两个迭代值的组合方法： §2.4.6 迭代收敛的加速方法 /* Accelerating Method*/ 本节讨论迭代法加速收敛问题，常用于线性收敛迭代法 将 x = g(x) 等价地改造为 当 和 时，有 相应的迭代公式为 或者 选取特殊的 ，有可能使迭代加速。 x y y = x y = g(x) x* 如： 迭代公式为 几何意义如图示 注： (1)这种迭代对原迭代公式(*)的各近似值在根 的两侧往复地趋于 时较为有效； 中点 (2)只有 且 较大时，加速效果才明显。 又如： 新的迭代函数为 当 时 根据定理2.5知，迭代法至少是二阶的. 但由于 不知道，故 也得不到，因此可取其 的近似值，即 从而有 二、Steffensen(斯蒂芬森)加速迭代法：（三个迭代值组合） x y y = x y = g(x) x* x0 从初值 出发，计算出 在曲线 上得到 两个点 用直线连接 、 两点它与 的交点设为 点的坐标为 将 视为新的初值，重复上述步骤 一般地， 由 组合得到迭代式 或者 这个方法称为斯蒂芬森(Steffensen)迭代方法 若令 则得到斯蒂芬森(Steffensen)迭代法的另一种形式： Steffensen迭代法的优点：可以改进收敛速度，有时也能把不收敛的迭代法改进为收敛的二阶方法. 艾特肯（Aitken）加速方法： 例4： 已知方程 在 上有一个根（正根） 下面选取3种迭代格式： 1、 即 2、 即 3、 即 显示计算结果 设不动点迭代 的迭代函数 在其不动点 的某邻域内具有二阶连续导数， 则斯蒂芬森(Steffensen)的迭代技术是二阶收敛的，而 且其极限仍为 。 证明： 设斯蒂芬森(Steffensen)的迭代为 其中 首先证明有 相同的不动点 设 则 ? 反之，设 型 洛比塔法则 ? * 第二章 非线性方程的数值解法 /* Numerical Solutions of Nonlinear Equations*/ 本章主要内容： 1、二分法(重点) 2、不动点迭代的构造及其收敛性判定 3、Newton（重点）和Steffensen迭代 4、割线法 5、非线性方程组的迭代解法 历史背景 代数方程的求根问题是一个古老的数学问题。理论上， 次代数方程在复数域内一定有 个根(考虑重数)。早在16世纪就找到了三次、四次方程的求根公式，但直到19世纪才证明大于等于5次的一般代数方程式不能用代数公式求解，而对于超越方程就复杂的多，如果有解，其解可能是一个或几个，也可能是无穷多个。一般也不存在根的解析表达式。因此需要研究数值方法求得满足一定精度要求的根的近似解。 求方程 几何意义 基本定理 如果函数 在 上连续，且 则至少有一个数 使得 ，若同时 的一阶导数 在 内存在且保持定号，即 (或 )则这样的 在 内唯一。 a b x* §1 二分法 /* Bisection Method */ 原理：若 f ?C[a, b]，且 f (a) · f (b) 0，则 f 在 (a, b) 上至少有一实根。 基本思想：逐步将区间分半，通过判别区间端点函数值的符号，进一步搜索有根区间，将有根区间缩小到充分小，从而求 出满足给定精度的根 的近似值。 以此类推 终止法则？ a b x1 x2 a b When to stop? 或 不能保证 x 的精度 x* ?2 x x* 二分法算法 给定区间[a,b] ，求f(x)=0 在该区间上的根x. 输入: a和b; 容许误差 TOL; 最大对分次数 Nmax. 输出: 近似根 x. Step 1 Set k = 1; Step 2 Compute x=f((a+b)/2); Step 3 While ( k ? Nmax) do steps 4-6 Step 4 If |x| TOL , STOP; Output the solution x. Step 5 If x*f