研究生数值分析第2章A.docVIP

第二章 方程求根 §2.1 方程求根与二分法 一、引言 在实际问题中经常要求求解方程f(x)=0。满足此方程的解x, 又称为方程的根, 也称x是函数f(x)的零点。如果函数f(x)可以写成: 其中,则当m≥2时,称x*是方程的m重根,或称x*是函数f(x)的m重零点;当m=1时, 称x*是方程的单根，或称x*是函数f(x)的单重零点。 如果f(x)是多项式,称此方程为代数方程,若f(x)是超越函数,就称f(x)=0为超越方程。一般一次方程称为线性方程,而二次以上的代数方程或超越方程称为非线性方程。对于低次代数方程,我们可以找到求根公式,而对于次以上的代数方程和超越方程,我们没有通用的技术来求出精确解,这就需要数值方法来求出方程的近似解。 分离区间：许多方程往往有两个以上的根,在某个区间[a,b]上,如果方程在此区间内只含一个根,我们称此区间为方程的分离区间。 如果f(x)在区间[a,b]上连续,满足f(a)f(b)0,即两个端点值异号,且f(x)在区间[a,b]上严格单调,则利用闭区间上连续函数的性质,可知f(x)在区间[a,b]上存在唯一的零点,其几何意义如下图: 曲线y= f(x)与x轴的交点就是f(x)的零点。 用描述的方法或称搜索的方法,先确定搜索区间的下界a,选择搜索步长h,然后求出遇到异号, 可知区间内有根, 从而可以加密搜索步长进行搜索, 提高精确度。这儿a,b如何确定?h选多大合适?h过大则可能会漏掉根(如图)，(图中，当h过大时, 漏掉了两个根)；h过小则计算量很大, 如要求精度达到0.001,则在长度为1的搜索区间中,要计算1000个函数值, 欲使精度提高一倍, 计算量就要增加10倍。 二、二分法 1、二分法的具体计算过程 设f(x)在[a,b]上连续,现求方程f(x)=0在区间[a,b]上的根。设函数f(x)满足f(a).f(b)0,不妨设f(a)0，f(b)0。 第一步,取区间中点,计算区间中点的函数值, ① 如果, 则就是方程的根; ② 如果, 则在区间上, f(x)在两个端点的函数值异号,于是原方程在区间内有根(图示), 记,下一步在区间[a1,b1]内继续进行; ③ 如果,则在区间上,f(x)在两个端点的函数值异号,于是原方程在区间内有根,记, 下一步在区间内继续进行。 第二步, 求f(x)在区间[a1,b1]的中点的函数值,并检验其正负号, ①如果, 则原方程在区间内有根,并记; ② 如果,则在区间上,原方程有根,记。 于是,我们得到,其区间宽度为: 象这样,继续进行第三步、第四步、...... , 区间宽度每次缩小一半,得到一个区间序列: 此时,f(an). f(bn)0,即原方程在区间[an,bn]内有根,区间宽度为: 当n足够大时,如果此时的区间宽度已达到精度要求,则以区间的中点作为x*的近似值,即; 此时，近似值的误差小于该区间宽度的一半，即。如果精度要求,则要求 ， 两边取自然对数，得: ln(b-a)-(n+1)ln2≤lnε，则 ， 注意到 ，，有． 2、计算流程 根据精度要求可以事先计算出需执行步骤数n． 初态:。 对于n=1,2,...n，做计算，如果，输出；如果，则 ，否则 ，输出． 其几何意义如图: 例1: 求方程f(x)= x3-2x-5=0在区间[2,3]上的近似解，精确到0.001。 解：f(x)= x3-2x-5，ε=0.001，因为f(2)=-10，f(3)=160，故方程在区间[2,3]上有根。又， 取n=9，将计算结果列表如下: 所以x*≈x9=2.0947265，而精确值为 2.0945515...，误差为0 3、二分法的特点 二分法的优点：计算简便，对函数f(x)的要求不高，只要求连续即可，且误差估计容易。二分法的缺点：收敛速度很慢，每计算一步，误差减小一半。 §2.2迭代法及其收敛性 一、不动点迭代法 设方程在区间上有唯一的实根，将方程变形为与其同解方程: 若满足方程，则称为的一个不动点，也是方程的一个根。 如果在区间上连续，则可以在上任取一点作为迭代法的初始值，建立迭代关系(递推关系式): , 称为不动点迭代法，称为迭代函数。 对给定的初值，我们可得到一个数列,如果当时,这个数列收敛到,即 ,则, 则满足方程,即为的不动点。由于方程和是同解方程，所以满足方程。在实际计算中，取足够大，则有，我们把作为原方程的近似解。 计算流程: 选取初值 对 作 如果,跳出循环 否则,置,继续循环 输出 例1:用迭代法求方程的根,精确到0.001。 解:设,,,在区间内有根。 将方程变形为,这里,,在内, 所以迭代是收敛的。 取,则 ， ,迭代结束。 。 几何意义: 取作y轴平行线,交于,作x轴下