第七章 Ch7 非线性方程求根.ppt

第七章 非线性方程求根 Solutions of Nonlinear Equations * §1 二分法 求 f (x) = 0 的根 原理：若 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* §2 Bisection Method 误差 分析： 第1步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： ①简单; ② 对f (x) 要求不高(只要连续即可) . ①无法求复根及偶重根 ② 收敛慢 注：用二分法求根，最好先给出 f (x) 草图以确定根的大概位置。或用搜索程序，将[a, b]分为若干小区间，对每一个满足 f (ak)·f (bk) 0 的区间调用二分法程序，可找出区间[a, b]内的多个根，且不必要求 f (a)·f (b) 0 。 ? + ? + + ? ? 1.25 1.375 1.3125 1.3438 1.3281 1.3203 1.3242 1.5 1.375 1.3438 1.3281 1.0 1.25 1.3125 1.3203 0 1 2 3 4 5 6 f(xk)符号 xk bk ak k §2 迭代法 /* Fixed-Point Iteration */ 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 的根。 So basically we are done! I can’t believe it’s so simple! What’s the problem? Oh yeah? Who tells you that the method is convergent? §3 Fixed-Point Iteration x y y = x x y y = x x y y = x x y y = x x* x* x* x* y=g(x) y=g(x) y=g(x) y=g(x) x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? x0 p0 x1 p1 ? 定义1 设?(x)为定义在区间I上的函数, 且对任何 x?I, 均有?(x)?I,则称?(x)为I到自身上的映射. 定义2 设?(x)为I到自身上的映射,且存在0L1,使对任何x1,x2?I,有|?(x2)-?(x1)|? L|x2-x1|,则称?(x)为I上的压缩映射, L称为Lipschitz常数. 若?(x)为I上的压缩映射,则?(x)在I上连续. 定理2 若?(x)为I到自身上的映射,且?(x)?C1(I), |??(x)|? L1,则?(x)为I上的压缩映射. 证 对任意x1,x2?I,有 |?(x2)-?(x1)|=|??(?)||x2-x1|? L|x2-x1| 定义3 若?(x)为I到自身上的映射,且??I满足,?= ?(?),则称?为?(x)的不动点. 定理3 若?(x)为I上的压缩映射, 则?(x)在I上存 在唯一的一个不动点?,且对任何x0?I,由迭代格式 xk+1=?(xk) , k=0,1,2,… 产生的序列{xk}收敛于?(x)的不动点?. 定理1 证 不妨设I=[a,b],作函数?(x)=?(x)-x,由于x?I时, ?(x)?I,则?(a)=?(a)-a?0,?(b)=?(b)-b?0,由?(x)的连续性, 必存在??I, 使?(?)=?(?)-?=0,即?=?(?),?就是?(x)的不动点. 若?,??I均为?(x)的不动点,则有 |?-?|=|?(?)-?(?)|? L|?-?||?-?| 所以只能?=?,即?(x)在I上仅有一个不动点. 对任意x0?I,有x1=?(x0)?I,递推得{xk}?I,设?是?(x) 的不动点,则 |xk+1-?|=|?(xk)-