§8函数方程的近似求解 、数值方法 求方程的解主要方法有两种 解析方法和数值方法。 数值方法是一种求近似解的方法。是设法构造 个可实际计算的过程,并通过运行这个过程产生 方程的精确解的一系列近似值。这些近似值(精 度较高)理论上将收敛(在一定的条件下)于方 程的精确解,称其为数值解或近似解。 数值方法是用数学工具解决实际问题过程中 的一个重要过程
1 §8 函数方程的近似求解 一、数值方法 1、求方程的解主要方法有两种: 解析方法和数值方法。 数值方法是一种求近似解的方法。是设法构造一 个可实际计算的过程,并通过运行这个过程产生 方程的精确解的一系列近似值。这些近似值(精 度较高)理论上将收敛(在一定的条件下)于方 程的精确解,称其为数值解或近似解。 数值方法是用数学工具解决实际问题过程中 的一个重要过程
2、迭代法 数值计算中常用的求近似值的方法是迭代法。 f(x)=0化为等价形式x=F(x) 即x是方程∫(x)=0的解,则成立 x=F(x), 反之亦然。F(x)称为迭代函数。 取初始值x,按 k+1=F(xk)k=0,1 产生序列{xk}(设每个xk都属于F(x)的定义域), 上述的计算过程称为迭代
2 2、迭代法 数值计算中常用的求近似值的方法是迭代法。 f x( ) 0 化为等价形式 x F x ( ) 即 x 是方程 f (x) = 0 的解,则成立 x F x( ) , 反之亦然。 F (x) 称为迭代函数。 取初始值 x0 ,按 1 ( ) 0,1, 2, k k x F x k 产生序列 {xk }(设每个 xk 都属于 F (x) 的定义域), 上述的计算过程称为迭代
3、 Newton迭代法 求函数方程∫(x)=0的解: 是求曲线y=f(x)与x轴的交点的横坐标x, 若曲线在x=xk处的切线方程为个y=) y-f(k=f(x(x-xR 与x轴的交点的横坐标为 f(xk) +1 k+1-kk-1 即 Newton迭代法实质上是通过一系列的切线 与x轴的交点的横坐标,来逼近曲线与x轴的三 交点的横坐标,所以又称为 Newton切线法。 3
3 3、Newton 迭代法 求函数方程 f (x) = 0 的解: 是求曲线 y = f (x) 与 x 轴的交点的横坐标 xk , 若曲线在 x = xk 处的切线方程为 ( ) ( )( ) k k k y f x f x x x 与 x 轴的交点的横坐标为 ( ) ( ) k k k f x x x f x k 1 x x y 0 y = f (x) x * xk+1 xk xk-1 即 Newton 迭代法实质上是通过一系列的切线 与 x 轴的交点的横坐标,来逼近曲线与 x 轴的 交点的横坐标,所以又称为 Newton 切线法
定理:设f(x)在a,b中有二阶连续导数,且满足 (1)f(x)∫(b<0; (2)f(x)在(a,b)保号; (3)∫(x)在(a,b)保号; 又设x是a和b中满足 f(x0)·f"(x0)>0 的点,则以x为初值的 Newton迭代过程 ∫(x) k+1=k k=0,1,2, 产生序列{xk}单调收敛于方程∫(x)=0在{a,b 中的唯一解
4 定理:设 f (x) 在 [a, b] 中有二阶连续导数,且满足 (1) f (x) · f (b) < 0 ; (2) f x ( ) 在(a, b) 保号 ; (3) f x ( ) 在(a, b) 保号 ; 又设x0 是 a 和 b 中满足 0 0 f x f x ( ) ( ) 0 的点,则以 x0 为初值的 Newton 迭代过程 1 ( ) 0,1, 2, ( ) k k k k f x x x k f x 产生序列 {xk }单调收敛于方程 f (x) = 0 在 [a, b] 中的唯一解
例:解方程Ix=sinx 解:记f(x)=nx-sinx,考虑区间[,el, 则有∫6|=ln3-sin4<0,f(e)=le-sne>0, 2 2 而∫(x) 元 C0sx>0x∈[,引, 元 cte+sinx> sine 2>0x∈[,el, 元 2 所以符合定理的全部条件; 因为f(e)"(e)>0 所以初值取为x0=e,利用 +1 f(x)k=0,1,2 f(xk)
5 例:解方程 ln sin . x x 解:记 f x x x ( ) ln sin ,考虑区间 [ , ] , 2 e 则有 ln sin 0 , 2 2 2 f f e e e ( ) ln sin 0 , 而 1 ( ) cos 0 [ , ] , 2 f x x x e x 2 1 f x x ( ) sin x 2 4 sin 0 [ , ] , 2 e x e 所以符合定理的全部条件; 因为 f e f e ( ) ( ) 0 , 所以初值取为 x0 = e , 利用 1 ( ) 0,1, 2, ( ) k k k k f x x x k f x