8YQJCh4方程式（含案例）求根Q.ppt

4.1 非线性方程的迭代解法 4.1.1 对分法（二分法） 用C[a, b]表示 [a, b] 上一切连续函数所成集合. 连续函数的介值定理 若 f ?C[a, b], 且 f (a) · f (b) 0, 则 f 在 (a, b) 内至少有一实根。 这时称[a,b]为方程 f(x)=0 的有根区间. 对分法求根思想 找有根区间(ak , bk); 用(ak , bk)的中点近似根. 4.1.2 简单迭代法及其收敛性 ( Fixed-Point Iteration ) ?迭代法思想 把 f (x) = 0 转化为等价的方程 x = g(x), 从而， f (s) = 0 = s = g(s). s 叫 g(x)的不动点. ?求 g(x) 不动点的过程 定理4.1 考虑方程 x = g(x), g(x)?C[a, b], 若 (I) 当 x?[a, b] 时， g(x)?[a, b]； (II) ? 0 ? L 1，使得 | g(x) | ? L 1,? x?[a, b]。 则，任取 x0?[a, b]，由 xk+1 = g(xk) 得到的序列 { xk} 收敛于g(x) 在[a, b]上的唯一不动点 s。并且有误差估计式： ?局部收敛性 定义 设 g(x) 有不动点 s , 如果存在 s 的某邻域 | x - s | ?, 使得, 对任意的 x0, 迭代法产生的序列都收敛到 s, 则称迭代法是局部收敛的. 由定理4.1容易得 定理4.2 设 s 为g(x)的不动点, g(x) 在 s 的某个邻域内连续, 且| g(s) | 1, 则迭代法是局部收敛的. 4.1.3 简单迭代法的收敛速度 ?对于一种迭代过程，为了保证它是有效的，需要肯定它的收敛性，同时考察它的收敛速度。 ?收敛阶：在接近收敛的过程中，迭代误差的下降速度。 收敛阶是对迭代法收敛速度的一种度量。 定义 设迭代 xk+1 = g(xk) 收敛到g(x) 的不动点 s. 又 ek = xk ? s，若有常数 C 和 p≥1 , 使得 定理4.3 设函数 g(x)?C[a, b], g(x)?C(a, b), 且 (I) 当 x?[a, b]时, g(x)?[a, b]; 即, (II) 当 x?(a, b) 时, g(x) ≠0, | g(x) | ? L 1, 其中L 是一个常数. 则对任取的 x0?[a, b], 由 xk+1 = g(xk) 得到的序列 { xk} 收敛于 g(x) 在[a, b]上的唯一不动点 s .并且, 当 x0 ≠ s 时, { xk} 是线性收敛的. 证 用定理4.1, 知{ xk} 收敛于 g(x) 在[a, b]上的唯一不动点 s. 再证, 当 x0 ≠ s 时, 序列{ xk} 是线性收敛的. 由Taylor 公式, 有 xk+1=g(xk)= g(s) + g[s+? (xk - s )](xk-s) , 记 ek = xk ? s, 由上式可得, ek+1 = g[s+? ek]ek , 故有 问：会出现 ek = xk ? s= 0 吗? 由于 x?(a, b)， g (x) ≠0， 故当 x0 ≠ s 时， 必有 xk ≠ s ， k = 0,1, 2, ...， 这是因为， x1-s = g(x0) -g(s) = g(?1)(x0-s) ≠0,?1?(a, b), x1≠ s . x2 - s = g(x1) -g(s) = g(?2)(x1-s) ≠0, x2≠ s . ...... 定理4.4 设 s为g(x)的不动点， g(m)(x) 在含 s的某个开区间内连续， (m≥2 )。 如果 g(k)(s) =0，( k =1, 2, ..., m-1)， g(m)(x) ≠ 0， 则存在 ? 0，对于 x0?B?=[s- ?, s+ ?]，使得序列 xk+1= g(xk) 在 B? 内 m 阶收敛于 s. 例 用迭代法求 f (x) = x3 - x -1=0 的正实根. 解: f (1) = -1, f (2)=5， 在 (1,2) 内 f (x) 有一实根. 又在 (1,2) 内 f (x)0， 知 f (x) 有唯一的实根 s. 方案1： 改写原方程为 方案2： 改写原方程为 补例 用简单迭代法求方程 x =e-x 在x=0.5附近的一个根, 要求 解：方程 x =e-x ,迭代函数为 g(x)= e-x , 在区间 [0.5, 0.6] 内有根, 根据定理4.2知,