非线性方程数值解法及其应用.docVIP

非线性方程数值解法及其应用 摘要:数值计算方法主要研究如何运用计算机去获得数学问题的数值解的理论和算法。本文主要介绍非线性方程的数值解法以及它在各个领域的应用。是直接从方程出发，逐步缩小根的存在区间，或逐步将根的近似值精确化，直到满足问题对精度的要求。我将从二分法、Steffensen加速收敛法、Newton迭代法、弦截法来分析非线性方程的解法及应用。 关键字：非线性方程;二分法；Steffensen加速收敛法；代数Newton法；弦截法 前言 随着科技技术的飞速发展，科学计算越来越显示出其重要性。科学计算的应用之广已遍及各行各业，例如气象资料的分析图像，飞机、汽车及轮船的外形设计，高科技研究等都离不开科学计算。因此经常需要求非线性方程 f(x) = O的根。方程f(x) = O 的根叫做函数f(x)的零点。由连续函数的特性知：若f(x)在闭区间[a，b]上连续，且f(a)·f(b)O，则f(x) = O在开区间(a,b)内至少有一个实根。这时称[a,b]为方程f(x) = O的根的存在区间。本文主要是对在区间[1.2]的根的数值解法进行分析，介绍了非线性方程数值解法的四种方法，从而得到在实际问题中遇到非线性方程根的求解问题的解决方法。 非线性方程的数值解法 二分法 二分法的基本思想是将方程根的区间平分为两个小区间，把有根的小区间再平分为两个更小的区间，进一步考察根在哪个更小的区间内。如此继续下去，直到求出满足精度要求的近似值。 设函数f(x)在区间[a,b]上连续，且f(a)·f(b)O，则[a,b]是方程f(x)=O 的根的存在区间，设其内有一实根，记为。取区间[a,b]的中点，并计算，则必有下列三种情况之一成立： = O,就是方程的根； (2)f(a)·f()O，方程的根位于区间[a,]之中，此时令，； (3)f()·f(b)O，方程的根位于区间[,b]之中，此时令。 在(2)、(3)两种情况下，取，并计算，重复上述过程，就可逐次把区间缩短一半，且始终包含根根。当经过k次二分后，根所在的区间的长度为。 若取有根区间的中点作为根的近似值，则在二分过程中，可以获得一个近似根序列，该序列必以根为极限。在实际计算时，不可能完成这种无穷过程，其实也没有这种必要。由于 所以，只要二分足够次(即k足够大)，便有。这里为事先给定的精度，再注意到，所以，在实际计算时，只要某个有限区间的长度小于，就可以停止计算，并取该有限区间的中点作为根的近似值。 二分法的优点是算法简单及近似根序列一定收敛，缺点是收敛速度比较慢。 Steffensen加速收敛法 如图所示，由两点，连接PQ弦有方程： 令y=x,可得点，有解式： 代数Newton法 设n次代数方程 用Newton迭代法求有限区间的实根，则要计算，一般采用秦九韶算法。由Taylor展式? 比较x的同次幂系数得： 则得出代数Newton的公式: 4、弦截法 代数Newton虽然收敛速度快，但需要计算导数，如果函数比较复杂，就会带来一些不便。因此要考虑一种能避开导数运算的迭代公式。 可以用平均变化率：来近似代替。 在牛顿迭代公式中，用来近似代替，就得到迭代公式 按这个公式进行迭代计算的方法就称为弦截法。弦截法的收敛速度比牛顿法慢得多，为了加快收敛速度，改用差商 来代替牛顿迭代公式中的导数，于是得到下列快速弦截法的迭代公式 因为在计算时，用了前面两步的信息和，所以收敛速度就提高了。应该注意，在使用快速弦截法迭代公式时，必须先给出两个初始近似根和。弦截法的优点是收敛速度也相当快。 三、非线性方程数值解法的MATLAB实现 二分法 （1）MATLAB程序： （2）二分法的MATLAB实现及分析： 用二分法计算在区间[1,2]上的根。因为二分法只能求单根，首先搜索函数在区间[1.2]的根的情况。在MATLAB命令窗口输入程序： x=1:0.01:2; y=2*x.^3+2*x.^2-5; plot( x,y) 得到图1。从图1中可以得到函数在区间[1，2]间有唯一的一个大于1.35而小于1.4的单根。建立函数的程序function f=f(x);f=2*x^3+2*x^2-5;在命令窗口输入erfen(‘f ,1,2,0.0001)，得到结果x=1.0929。