第三章迭代法s1s2.pptVIP

第三章迭代法 §3.1 二分法 §3.2 迭代法原理 §3.3 Newton迭代法和迭代加速 §3.4 解线性方程组的迭代法 §3.1 二分法 根的估计 二分法 非线性方程的根 根的估计 引理3.1（连续函数的介值定理） 设f(x)在[a,b]上连续，且f(a)f(b)0，则存在x*?(a,b)使f(x*)=0。 例3.1 证明x3?3x?1 = 0 有且仅有3个实根，并确定根的大致位置使误差不超过? =0.5。 解: 单调性分析和解的位置 选步长h=2?, 扫描节点函数值 异号区间内有根 第一节 二分法 算法 误差分析 例3.2 x3?x?1 = 0, [1,2] 两位有效数字 ?=0.5*10-1 , k?(ln 20/ ln 2)-1, 取k=4 §3.2 迭代法原理 迭代法的思想 不动点原理 局部收敛性 收敛性的阶 第二节 迭代法原理 例 x3?x?1 = 0, [1,2], 取 x0=1.5 迭代公式1： 迭代公式2： 计算结果： 不动点原理 收敛性分析 全局收敛与局部收敛 例 用不同方法求 x2?3 = 0的根 , 取 x0=2. 迭代公式1： 迭代公式2： 迭代公式3： 迭代公式4： 计算结果： 收敛性的阶 p 阶收敛 p 阶收敛 * 求 f (x) = 0 的根 代数方程： f (x) = a0 + a1x + . . . + anxn 超越方程： f (x) 含超越函数，如 sin(x), ex, lnx 等 实根与复根 根的重数 f (x) = ( x – x*)m · g(x) 且 g(x*) ? 0, 则 x* 为 f (x) 的 m 重根 有根区间：[a, b] 上存在 f (x) = 0 的一个实根 在有根的前提下求出方程的近似根。 研究 内容: 若 f ? C[a, b]，且 f (a) · f (b) 0，则 f 在 (a, b) 上必有一根。 基本原理： 具体方法： 通过二等分不断缩小有根区间的长度，直到满足精度为止。 a b x0 x1 x* 何时终止? 或 不能保证 x 的精度 算法 3.1 (二分法) 给定有根区间 [a, b] ( f(a) · f(b) 0) 和 精度要求 ? 1. 令 x = (a+b)/2 2. 如果 b – a = 2?, 停机，输出 x ,否则执行第3步 3. 如果 f (a) f (x) 0 , 则令 b = x，否则令 a = x, 返回第1步 P50. Matlab源程序：nabisect.m 用二分法求根，通常先给出 f (x) 草图以确定根的大概位置。 记 a0 = a, b0 = b, 第 k 步的有根区间为 [ak, bk] 对于给定的精度 ? ,可估计二分法所需的步数 k ： 取 简单易用 无法求复根及偶重根 对 f (x) 要求不高，只要连续即可 收敛速度慢 ? (x) 的不动点x* f (x) = 0 x = ? (x) 称为迭代函数 等价变换 基本思想 从一个给定的初值 x0 出发，计算 x1 = ? (x0), x2 = ? (x1), … 若 收敛，即存在 x* 使得 ，则由 ? 的连续性和 可得 x* = ? (x*)，即 x* 是 ? 的不动点，也就是 f (x) 的零点。 具体做法： 不动点迭代 f (x) 的零点x* 几何含义： 求曲线 y = ?(x) 与直线 y = x 的交点 xk+1 = ?(xk) x y y = x x y y = x x y y = x x y y = x x* x* x* x* y= ?(x) y= ?(x) y= ?(x) y=?(x) x0 p0 x1 p1 x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x2 ? 公式1： 公式2： 怎么判断迭代公式收敛或发散呢？ 定理3.1 设 ?(x)?C[a, b] 且一阶导数连续，若 (2) ? 0 ? L 1，使得 | ?’(x) | ? L 对 ? x?[a, b] 成立 则函数 f (x) = x - ?(x) 在 [a, b] 中有唯一的零点 x*。 (压缩映像定理，不动点原理) x* 称为 ?(x) 的不动点 ?(x*) = x* (1) a ? ?(x) ? b 对一切 x?[a, b] 都成立 简证： f(a) = a - ? (a)? 0 , f(b) = b - ? (b) ? 0 f(x) 在[a, b] 上有零点。 唯一性：反