《上海交大数值分析课件数值分析7-1二分法》课件.pptVIP

第七章 非线性方程的求根 /* Solutions of Nonlinear Equations */ 求 f (x) = 0 的根 7.1 方程求根与二分法 一、本章解决的问题 二、求根的两个步骤 三、二分法 一、本章解决的问题 在科学计算中常要求解各种方程， 这些方程看似简单，但难于求其精确解。而实际问题：只要能获得满足已定精确度的近似根就可以了。 高次代数方程 超越方程 本章解决的问题： f(x)为非线性函数或高次代数方程，若有数x*使f(x*) = 0成立，则称x*为方程f(x) = 0的根(零点)。 求 f (x) = 0 的根 若f(x)可分解为 m是正整数，且g(x*)≠0 当m = 1，称x*是单根； 当m 1，称x*是m重根. x y 0 y=f(x) 方程根的几何意义 二、求根的两个步骤 (1)确定根的初始近似值(称之为初始近似根) ，一般为一个包含根的区间，称为“有根区间” (2)根的精确化。根据根的初始近似值按某种方法逐步精确化，直至满足预先要求的精度为止。 如何求有根区间呢？ 逐步扫描法 原理：设f(x)在[a, b]连续，且f(a) f(b)0。则由连续函数的性质知f(x)=0在(a, b)内至少有一个根。若f(x)在[a, b]上单调，则f(x)=0在(a, b)上有且仅有一个根。 x y y=f(x) 0 故总假设(a,b)上有唯一根 逐步扫描算法 (1) x0←a； (2) 若 f( x0) f( x0+h)0，则x*必在( x0 , x0+h)中，取 x0或 x0+h作为有根区间，否则转(3); (3) x0← x0+h，转(2)； 例如 考虑方程 解 由于 故方程至少有一个正实根。 设从x=0出发，取h=0.5为步长向右计算，将各个点上的函数值列于下表： x 0 0.5 1.0 1.5 f(x) ― ― ― + 由于 , 且 f(x)在区间[1，1.5]上满足 由此可知在(1，1.5)内有且仅有一个 实根，故可取 作为有根区间。 (1，1.5) 下面将介绍几种常用的数值解法： ——二分法 ——简单迭代法 ——牛顿迭代法 ——弦截法 三、二分法 /* Bisection Method */ 1. 二分法的原理 原理：若 f ?C[a, b]，且 f (a) · f (b) 0，则 f 在 (a, b) 上必有一根 x*。 2. 二分法的实施 将方程根的区间平分为两个小区 间，然后判断根在哪个小区间，舍 去无根的区间，而把有根区间再一 分为二，再判断根属于哪个更小的 区间，如此周而复始，直到求出满 足精度要求的近似根。 a b x1 x2 a b When to stop? 或 不能保证 x 的 精度 x* ?2 x x* 误差分析： 第1步产生的 有误差 第 k 步产生的 xk 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k ： ①简单; ② 对f (x) 要求不高(只要连续即可) . ①无法求复根及偶重根 ② 收敛慢 (b-a) 这里是k+1次 注：用二分法求根，最好先给出 f (x) 草图以确定根的大概位置。或用搜索程序，将[a, b]分为若干小区间，对每一个满足 f (ak)·f (bk) 0 的区间调用二分法程序，可找出区间[a, b]内的多个根，且不必要求 f (a)·f (b) 0 。 例如 求方程 在区间(1, 1.5)内的根,要求用四位小数计算,精确到10-2。 3. 二分法举例 解 这里a=1，b=1.5，我们先估计要二分的次数，按第二种方法，k应满足 求得k=7 k ak bk xk f(xk)的符号 0 1 1.5 1.25 + 1 1.25 1.5 1.375 + 2 1.25 1.3375 1.3125 - 3 1.3125 1.3375 1.3438 + 4 1.3125 1.3438 1.3281 + 5 1.3125 1.3281 1.3203 - 6 1.3203 1.3281 1.3242 - 故原方程在区间(1,1.5)内的根 谢谢！