计算方法 6 非线性方程迭代法.ppt

第六章 解非线性方程的迭代法 非线性方程（常见的有：代数方程 [多项式]；超越方程 [三角、指数、对数方程]）求根问题： 关于方程求根，需考虑如下几个问题： 1. 如何选取初始点？ 2. 如何构造迭代序列（迭代格式）？ 3. 迭代序列是否收敛？在什么条件下收敛？ 4. 若收敛，收敛速度如何？并给出定量的刻画. 5. 讨论近似解的误差估计及算法的复杂性. 确定方程的有根区间; 计算根的近似值. ? ? 1. 方程是否有根（存在性）？ 2. 单根？重根? 多根？有多少个根？ 零点定理 构造迭代法 1. 二分法 – 一种逐步搜索法 1. 绘图法. 2. 逐步搜索法 a b x0:= x1 (a+b)/2 x* a b x0 x1 x* a0 b0 优点：算法简洁；方法可靠； 只要求f(x)连续 缺点：不能求偶数重根和复根， 且收敛较慢. 2. 简单迭代法 显然， 同解构造 注意， 反之亦然. 等价 转化！ 是解线性方程迭代法的推广！ x y y = x x* y=φ(x) p0 p1 ? x y y = x x* y= φ(x) p0 p1 ? 迭代法的基本步骤如下： 1、给出方程在有根区间上局部等价形式 2、选取适当的初始点 ，产生迭代序列 3、求极限 易知，该值即为方程的根. 一定收敛吗？ 注：迭代法的收敛与发散，依赖于迭代函数的构造（构造的方法很多！）.迭代函数须满足什么条件，迭代法才能收敛？ 称为全局Lipschitz条件 称为Lipschitz常数 不动点存在的条件？充分条件！ 充分性条件 收敛速度？误差估计？ 称序列是适定的，它表明 迭代法算出的每个点是有 意义的！ 它表明由迭代法产生的序列收敛！ 事前误差估计 事后误差估计 迭代法的全局收敛性 例2.1 在区间[1, 2]内讨论迭代格式的敛散性. 例2.2 在用迭代法求方程 在区间[0, 1]内的一个根. 表明：在大的有根区间上，条件不一定都成立！所以，使用迭代法时往往在根的附近进行！ 0 迭代法的局部收敛性 收敛阶的概念 衡量迭代法收敛快慢的重要标准. p的大小反映了序列{xk}收敛速度的快慢. 迭代格式的收敛速度 与迭代函数有关！ ① ② 3. Aitken（埃特金）加速迭代法 ① ② 4. Newton（牛顿）迭代法，又称为牛顿-拉弗森方法（Newton-Raphson method） 1. 牛顿法的基本思想 将非线性方程逐步线性化，以线性方程的解逼近非线性方程的解。 将非线性方程逐步线性化，如何实现？？ 取 xk ? x*，将 f (x) 在 xk 处做一阶Taylor展开: ，? 在 xk 和 x 之间 2. 牛顿迭代法的原理 可将 (x* ? xk)2 看成高阶小量，则有： 取 忽略高阶项 x y x* 令 牛顿法最初由艾萨克·牛顿，1671年完成，在牛顿死后的1736年 公开发表。约瑟夫·拉弗森也曾于1690年提出此方法 称上式称为Newton迭代格式.用Newton迭代格式求方程根的方法称为Newton法. 2、牛顿迭代法的局部收敛性定理 其中 ，则 收敛 ? 证明： 牛顿迭代法事实上是一种特殊的不动点迭代 由泰勒展开： 在单根附近收敛快！ 只要 ，则令 可得结论. 在 和 之间 牛顿迭代法的改进 ? 重根 问题1: 若　　　 ，牛顿迭代法是否仍收敛？ 设 x* 是 f 的 m 重根，则： 且 。 从而有 答2: 若根的重数未知，可将 f 的重根转化为另一函数的单根。 答1: 有局部线性收敛性，但重数 m 越高，收敛越慢。 问题2: 如何加速重根的收敛速度 ? 令 　　　　　 ，则 f 的重根是 ? 的单根，且 对 构造出相应的牛顿迭代格式，迭代函数为 从而可构造出相应的迭代法格式为 若已知根的重数为 n，可将迭代格式改为， 所以格式是平方收敛的 ① 收敛速度快，稳定性好； ② 精度高。 ① 在重根附近收敛速度会降阶； ② 每次都要计算函数及其导数值，计算量大。 优点 缺点 注：牛顿法是局部收敛的，所以要求初值点 选在解 的附近，实际计算时，常先用简单迭代法算几步，估计出一 个质量较好的初值！！ 主要缺陷！！ 5. 弦截（