非线性方程的求根方法.pptVIP

取x0=2，对上述四种方法，计算三步所得结果如下：k xk (1) (2) (3) (4)0 x0 2 2 2 2x1 3 1.5 1.75 1.75x2 9 2 1.73475 1.732143x3 87 1.5 1.732361 1.732051注：x*=1.7320508……*第30页，共50页，星期日，2025年，2月5日例用迭代法求方程在内的实根。取解：对方程进行如下三种变形：建立迭代格式：这是一个发散的迭代格式。*第31页，共50页，星期日，2025年，2月5日建立迭代格式：该迭代格式收敛。建立迭代格式：该迭代格式收敛。*第32页，共50页，星期日，2025年，2月5日例.用迭代法求方程在内的实根。取解：对方程进行如下三种变形：理论分析：由上述定理知，迭代格式发散，和计算结果吻合。理论分析：由定理知，迭代格式收敛，和计算结果吻合。*第33页，共50页，星期日，2025年，2月5日理论分析：由定理知，迭代格式收敛，和计算结果吻合。而且，，②和③都收敛，但③收敛的效果比②好。*第34页，共50页，星期日，2025年，2月5日例求在[0，1]内的一个实根.将方程化为等价方程因为又所以对于[0，1]中任意初值，迭代序列收敛，计算结果如下表，取*第35页，共50页，星期日，2025年，2月5日非线性方程的求根方法*第1页，共50页，星期日，2025年，2月5日BisectionMethod方程求根的二分法Fixed-PointIteration迭代法及其收敛性NewtonMethodofNonlinearEquationsNewton迭代法*第2页，共50页，星期日，2025年，2月5日在实际应用中有许多非线性方程的例子，例如（1）在光的衍射理论（thetheoryofdiffractionoflight)中,需要求x-tanx=0的根（2）在行星轨道（planetaryorbits）的计算中,对任意的a和b，需要求x-asinx=b的根（3）在数学中，需要求n次多项式xn+a1xn-1+...+an-1x+an＝0的根*第3页，共50页，星期日，2025年，2月5日方程求解是科学计算中一个重要的研究对象;几百年前就已经找到了代数方程中二次至五次方程的求解公式;但是,对于更高次数的代数方程目前仍无有效的精确解法;对于无规律的非代数方程的求解也无精确解法;因此,研究非线性方程的数值解法成为必然。*第4页，共50页，星期日，2025年，2月5日非线性方程的一般形式：f(x)=0代数方程：f(x)=a0+a1x+……+anxn（an?0）超越方程：f(x)中含三角函数、指数函数、或其他超越函数。用数值方法求解非线性方程的步骤：（1）找出有根区间；（只含一个实根的区间称隔根区间）（2）近似根的精确化。从隔根区间内的一个或多个点出发，逐次逼近，寻求满足精度的根的近似值。*第5页，共50页，星期日，2025年，2月5日BisectionMethod二分法

二分法的基本思想： 假定f(x)=0在[a,b]内有唯一单实根x*，考察有根区间[a,b]，取中点x0=(a+b)/2，若f(x0)=0，则x*=x0，否则，（1）若f(x0)f(a)0，则x*在x0右侧，令a1=x0，b1=b;（2）若f(x0)f(a)0，则x*在x0左侧，令a1=a，b1=x0。介值定理设函数f(x)在区间[a,b]连续，且f(a)f(b)0，则方程f(x)=0在区间(a,b)内至少有一个根。*第6页，共50页，星期日，2025年，2月5日以[a1,b1]为新的隔根区间，且仅为[a,b]的一半，对[a1,b1]重复前过程，得新的隔根区间[a2,b2]，如此二分下去，得一系列隔根区间：[a,b]?[a1,b1]?[a2,b2]?……?[ak,bk]?……其中每个区间都是前一区间的一半，故[ak,bk]的长度：当k趋于无穷时趋于0。即若二分过程无限继续下去，这些区间最后必收敛于一点x*，即方程的根。*第7页，共50页，星期日，2025年，2月5日性质：f(an)·f(bn)0;bn–an=(b–a)/2nBisectionMethod 每次二分后，取有根区间的中点xk=(ak+