第2章非线性方程求解21.pptVIP

作业 习题二 1; * ? 2009, Henan Polytechnic University 第二章非线性方程求根 * 计算方法 在科学技术中,许多问题可归结为求非线性方程 　　　　 f (x)=0 (2.1) 的根,即函数 f (x) 的零点。 第二章 非线性方程求根 如果f(x)是多项式函数，则称为代数方程，否则称为超越方程（三角方程，指数、对数方程等）。 §2.1 方程求根与二分法 一、引言 通常方程根的数值解法大致分为三个步骤进行 ①?判定根的存在性。即方程有没有根？如果有 根，有几个根？ ② 确定根的分布范围。即将每一个根用区间隔 离开来，这个过程实际上是获得方程各根的 初始近似值。 ③ 根的精确化。将根的初始近似值按某种方法 逐步精确化，直到满足预先要求的精度为止。 为了确定根的初值，首先必须圈定根所在的范围，称为圈定根或根的隔离。 在上述基础上，采取适当的数值方法确定具有 一定精度要求的初值。 由高等数学的知识可知,设 f (x)为区间[a, b]上的连续函数,如果f (a)· f (b)0,则[a, b]中至少有一个实根。如果f (x)在[a, b]上还是单调递增或递减，则仅有一个实根。 y=f(x) a b y x 对于求方程根的问题，几何上就是求曲线 y=f (x)与x轴交点的横坐标。 求隔根区间的方法有： (1) 画图法 (2) 逐步搜索法 若方程在区间[a, b]内只有一个根，则称区间 [a, b]为隔根区间。 (1) 画图法 画出y= f (x)的草图，从而看出曲线与x轴交点的大致位置。 也可将f (x) = 0分解为?1(x)= ?2(x)的形式，?1(x)与?2(x)两曲线交点的横坐标所在的子区间即为隔根区间。 0 y x 例如： xlgx-1= 0，可以改写为lgx=1/x 画出对数曲线y=lgx 与曲线y= 1/x, 得到它们交点的横坐标所在区间。 (2) 逐步搜索法 对于给定的f (x), 设有根区间为[a, b],从x0=a出发,以步长h=(b-a)/n (n是正整数),在[a, b]内取点: xi=x0+ih (i=0,1,2,…,n), 从左至右检查f (xi)的符号,如发现xi与端点x0的函数值异号,则得到f (x)一个隔根区间[xi-1, xi]。 逐步搜索法： 例1 确定方程 f (x)=x3-x-1=0 的隔根区间。 解：用试凑法，不难发现 f(0)0 f(2)0 在区间(0,2)内至少有一个实根 设从x=0出发,取h=0.5为步长向右进行根的 搜索,列表如下 x f(x) 0 0.5 1.0 1.5 2 – – – + + 可以看出，在[1.0, 1.5]内必有一根。 二、二分法的基本思想 首先确定隔根区间,并将隔根区间二等分,通过 判断f (x)的符号,逐步将隔根区间缩小,直至隔根区 间足够地小,以便求出满足精度要求的近似根。 0 x y a b f(x) x0 ① 取隔根区间[a, b]的中点,将它对分两半,分点 ,这样就可缩小隔根区间。 具体过程如下： ② 对缩小了的隔根区间 ，取 0 x y a b f(x) x0 a1 b1 x1 然后再确定有根区间 ,其长度是 的二分之一。 a2 b2 得到新的隔根区间 ，施行同样的手法,即取 中点 ,将区间 再分为两半, 上述每个区间都是前一个区间的一半，因此， 的长度 ③ 如此反复下去,若不出现 ,即可得出一系列有根区间序列： 当n→∞时趋于零,这些区间最终收敛于一点 ，即为所求的根 。 二分法区间变换（单击播放） 换（二分法中点变单击播放） 三、二分法的收敛性分析 只要二分足够多次(即n足够大),便有 这里ε为给定精度,由于 ,则 　 每次二分后,取隔根区间 的中点 作为根的近似值，得到一个近似根的序列 该序列以根 为极限。 当给定精度ε＞0后,要想 成立, 只要满足 即可，亦即当: 做到第n次二分,计算得到的 就是满足精度要求的近似根 。 1 （先验估计） 在程序中通常用相邻的 与