数值分析第八章非线性方程解法讲解.ppt

（局部收敛性）设 f ?C2[a, b]，若 x* 为 f (x) 在[a, b]上的根，且 f ?(x*) ? 0，则存在 x* 的邻域 使得任取初值 Newton’s Method产生的序列{ xk } 收敛到x*，且满足 定理 证明：Newton’s Method 事实上是一种特殊的不动点迭代其中 则 收敛 ? 由 Taylor 展开： 只要 f ?(x*) ? 0, 则令 可得结论. 在单根 附近收敛快 再证： * 其中f (x) 为非线性函数. 如 求 f (x) = 0 的根 第八章 非线性方程求根 例 若干年以前, 美国原子能委员会准备将浓缩的放射性废料装入密封的圆桶内沉至海底. 但是, 当时一些科学家与生态学家都反对这种作法. 科学家用实验测定出圆桶能够承受的最大撞击速度为v=12.2 m/s, 如果圆桶到达海底时的速度超过这个速度, 将会因撞击海底而破裂, 从而引起严重的核污染. 然而原子能委员会却认为不存在这种可能性. 根据圆桶的质量, 体积以及海水的密度与海底的深度, 通过建立数学模型得知圆桶到达海底时的速度v (m/s)满足如下方程: 那么圆桶到达海底时的速度究竟会不会超过12.2 m/s呢? §1 对分区间法(二分法) 原理: 若 f ?C[a, b]，且 f (a) · f (b) 0，则 f 在(a, b) 上必有一根. [a, b] 称为 f (x)=0的有根区间. 下文中, 设 x* 是方程 f (x)=0的根. 取 不妨设 f (a) 0, f (b) 0. ① 若 则 ② 若 取 ③ 若 取 以 作为新的有根区间继续迭代, 得有根区间序列 ? 取 为 的第n个近似值. a b x0 x1 a b 停机准则 或 不能保证 x 的精度 x* ?2 x x* 误差 分析： 第0步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： ①简单; ② 对 f (x) 要求不高(只要连续即可) . ①无法求复根及偶重根 ② 收敛慢 例1 试用二分法求方程 的唯一实根, 要求误差不超过 解 [1, 2]为有根区间; f (x)单调增加, 方程有唯一根. 对分区间次数 取 n=14. 计算结果如下表: 0 1 2 1.5 ?1.6.. 1.75 2.8.. 2 ... ... ... ... ... 12 1.5944825 1.5947266 1.5946046 0.0007... 13 1.5944825 1.5946046 1.5945436 ?0.0003... 14 1.5945436 1.5946046 1.5945741 1.5 2 1 1.5 1.75 1.625 0.54.. f (x) = 0 x = g (x) 等价变换 f (x) 的根 g (x) 的不动点 思路 从一个初值 x0 出发, 计算 x1 = g(x0), x2 = g(x1), …, xk+1 = g(xk), … 若 收敛, 即存在 x* 使得 且 g 连续, 则由 可知 x* = g(x* ),即x* 是 g 的不动点，也就是 f 的根. §2 迭代法 例2 用简单迭代法求方程 的根, 要求精确到0.5?10-6. ? 解法一： 迭代格式: 取初值 计算得 显然此迭代序列发散. 例2 用简单迭代法求方程 的根, 要求精确到0.5?10-6. ? 解法二： 迭代格式: 初值 故 计算得 例2 用简单迭代法求方程 的根, 要求精确到0.5?10-6. ? 解法三： 迭代格式: 初值 计算得 例2 用简单迭代法求方程 的根, 要求精确到0.5?10-6. ? 解法四： 迭代格式: 初值 计算得 上例表明, 对同一方程可构造不同的迭代格式， 产生的迭代序列收敛性也不同. ? 迭代法的收敛性与什么有关？ ? 怎样构造