非线性方程的求解毕业论文

非线性方程的求解毕业论文题目(中文):非线性方程的求解(英文):TheSolutionofNonlinearEquations目录绪论................................................................................................................................11非线性方程的简介.....................................................................................................11.1非线性方程的背景...........................................................................................11.2非线性方程的概念...........................................................................................22非线性方程求解的数值方法.......................................................................................32.1二分法..............................................................................................................32.1.1二分法的思想........................................................................................32.1.2二分法的推理........................................................................................32.1.3二分法的应用........................................................................................42.2牛顿迭代法......................................................................................................42.2.1迭代法...................................................................................................42.2.2牛顿迭代法...........................................................................................62.3改进牛顿迭代法............................................................................................102.3.1改进牛顿迭代法的背景......................................................................102.3.2改进的法................................................................................11Newton3牛顿迭代法和改进牛顿迭代法的应用....................................................................123.1牛顿迭代法的应用.........................................................................................123.2改进牛顿迭代法的应用.................................................................................194结束语......................................................................................................................22参考文献......................................................................................................................23致谢..............................................................................................................................24I非线性方程的求解摘要非线性方程在实际问题中经常出现，很多熟悉的线性模型都是在一定的条件下由非线性问题简化得到的;非线性方程在科学与工程计算中的地位越来越重要,因此研究和探讨非线性方程求解的方法是非常有必要的。本文先开始介绍了非线性方程的概念及相关背景，再着重描述了非线性方程的求解的一些常用分法:二分法，迭代法，牛顿迭代法。在这些方法当中，牛顿迭代法是求解非线性方程的一种非常常用并且有效的方法，但是牛顿迭代法有一些应用条件限制，因此提出了改进的牛顿迭代法;针对非线性方程的实例用上面提到的方法进行了数值计算，并且比较了牛顿迭代法和改进牛顿迭代法，最后介绍了牛顿迭代法在实际生活中的应用。【关键字】非线性方程牛顿迭代法数值计算IITheSolutionofNonlinearEquationsAbstractNonlinearequationsappearfrequentlyinpracticalproblems,andmanyofusarefamiliarwiththelinearmodelobtainedbythesimplifiednonlinearproblemsundercertainconditions.Nonlinearequationsarebecomingmoreandmoreimportantinscienceandengineeringcomputing.Therefore,itisnecessarytostudyandexplorewaystosolvenonlinearequations.Firstly,thispaperrecommendsomebasicconceptionsandrelatedbackgroundofnonlinearequations,thendescribesomemethodsofthesolutionofnonlinearequationemphatically,suchas:theprocedureofdichotomy,theiteratemethod,theNewtoniteratemethodandtheimprovedNewtoniteratemethod.ItisveryusefulandeffectivetousetheNewtoniteratemethodforsolvingnonlinearequationsinthosemethods.However,weproposetheimprovedNewtoniteratemethodbecauseofthelimitsoftheNewtoniteratemethod.Also,wehavecarriedontheapproximatecalculationtothenonlinearequationsandhavecomparedtheNewtoniteratemethodwiththeimprovedNewtoniteratemethod,Intheend,weintroducetheapplicationoftheNewtoniteratemethodinthereallife.【Keywords】ThenonlinearequationsTheNewtoniteratemethodNumericalcomputationIII绪论非线性是实际问题中经常出现的，并且在科学与工程计算中的地位越来越重要，很多我们熟悉的线性模型都是在一定的条件下由非线性问题简化得到的，为得到更符合实际的解答，往往需要直接研究非线性模型，从而产生非线性科学，它是21世纪科学技术发展的重要支柱.非线性问题的数学模型有无限维的如微分方程，也有有限维的.从线性到非线性是一个质的变化，方程的性质有本质的不同，求解方法也有很大的差别.非线性方程的数值解法在实际中有广泛的应用，特别是在各种非线性问题的科学计算中更显出它的重要性，而且，随着计算机的广泛应用，有更多的领域涉及到非线性方程的求解问题，例如，动力系统，非线性有限元问题，非线性力学问题，还有非线性最优化与非线性规划问题等，因此，研究性方程的解法就具有重要的实际意义.由于非线性方程的复杂性，在解法上除了极特殊的非线性方程外，直接法几乎是不能使用的，这需借助于二分法，迭代法来求解.从计算的经验来看,Newton迭代法用来求非线性方程一种非常常见的而且是有效的方法，所以我们有必要研究和探讨求解非线性方程的Newton方法.1非线性方程的简介1.1非线性方程的背景非线性科学是一门研究非线性现象共性的基础学科.它是自20世纪六十年代以来，在各门以非线性为特征的分支学科的基础上逐步发展起来的综合性学科，被誉线性科学几乎涉及了自然科学和社会科学为本世纪自然科学的“第三次革命”.非的各个领域，并正在改变人们对现实世界的传统看法.科学界认为:非线性科学的研究不仅具有重大的科学意义，而且对国计民生的决策和人类生存环境的利用也具有实际意义.由非线性科学所引起的对确定论和随机论、有序与无序、偶然性与必然性等范畴和概念的重新认识，形成了一种新的自然观，将深刻地影响人类的思维方法，并涉及现代科学的逻辑体系的根本性问题.1非线性问题的“个性”很强，处理起来十分棘手.历史上曾有过一些解非线性方程的“精品”，但与大量存在的非线性方程相比，只能算是“凤毛麟角”.因此，长期以来，对非线性问题的研究一直分散在自然科学和技术科学的各个领域.20世纪六十年代以来，情况发生了变化.人们几乎同时从非线性系统的两个极端方向取得了突破:一方面从可积系统的一端，即从研究多自由度的非线性偏微分方程的一端获得重大进展.如在浅水波方程中发现了“孤子”，发展起一套系统的数学方法，如反散射法，贝克隆变换等，对一些类型的非线性方程给出了解法;另一方面，从不可积系统的极端，如在天文学、生态学等领域对一些看起来相当简单的不可积系统的研究，都发现了确定性系统中存在着对初值极为敏感的复杂运动.促成这种变化的一个重要原因十计算机的出现和广泛应用.科学家们以计算机为手段，勇敢地探索那些过去不能用解析方法处理的非线性问题，从中发掘出规律性的认识，并打破了原有的学科界限，从共性、普适性方面来探讨非线性系统的行为.在数值计算中，非线性问题也是经常遇到的一类难题，特别是非线性方程组的数值求解问题构成了非线性科学的一个重要组成部分.1.2非线性方程的概念非线性方程，就是因变量与自变量之间的关系不是线性的关系，一般可以表示为.这类方程很多，例如平方关系、对数关系、指数关系、三角函数关系等fx()0,32等.下面这些例子就是常见的非线性方程:，，xxx,，，,10xx,,tan0x.xae,,0非线性方程可分为两类:一类是多项式方程，这类方程可以定义为:2n，.另一类是非多项式方程，fxxxx()0,,，,，,，,,nNC,,,,,,,,01n012n它不能用多项式方程的形式表示，没有固定的形式.求解第一类多项式方程，现在已经有了比较成熟的理论和方法.现在比较常用的一种数值方法是迭代法，能通过迭代次数的增加，从而越来越接近方程的解，求解第二类非多项式方程，是现在数学领域中的一个重点研究方向.一般来说，求解此类方程是采用随机搜索的办法.22非线性方程求解的数值方法2.1二分法2.1.1二分法的思想二分法是区间迭代法的一种.它是重复运用零点存在性定理,每次将区间压缩一半且其中一个区间至少包含一个根,逐步缩短区间,直至最终区间长度满足一定的精度要求为止.2.1.2二分法的推理先考察有根区间,a,b],取中点，将它分成两半，然后进行根的xab,，()/20*搜索，即检查与是否同号，如果确系同号，说明所求的根在x的右fx()fa()x00axb侧，这时令=,=.b101*否则必在x的左侧，这时令a=，b=x，不管出现哪一种情况，新的有xa0110abab根区间[,]的长度仅为,,,的一半.对压缩了有根区间[,]又可施行同样ab1111xabab的过程，即用中点=(+)/2，以将区间[,]再分为两半，然后通过根的搜11111x[]abab索判定所求根在的哪一侧，从而又确定一个新的有根区间，长度是[,]12,211的一半.如此反复二分下去，可得出一系列有根区间，[,][,][][,]abababab,,,,,112,2kkk其中每个区间都是前一个区间的一半，因此[,]ab的长度ba,=，()/2ba,kkkk当时趋向零，就是说，如果二分过程无限地继续下去，这些区间最终必将k,,*收缩于一点x，该点显然就是所求的根.[,]ab每次二分后，设取有根区间的中点xab,，()/2作为根的近似值，则kkkkk*在二分过程中可以获得一个近似根的序列xxxx,,,,，则该序列必以根x为012k极限.不过在实际计算时，不可能完成这个无限过程，其实也没有这个必要，因为数值分析的结果允许带有一定的的误差，由于:3*1k，||()/2()/2xxbaba,,,,,(2.1.1)kkk*只要二分足够多次(即k充分大)便有,，这里为预定的精度.||xx,,,k2.1.3二分法的应用3例1求方程在区间内的一个实根，要求准确到小数fxxx()10,,,,(1.0,1.5)点后第二位.,解这里,而，取的中点x=将区间二fa()(,)abab,,1.0,1.50,()0fb,1.250*等分，由于，即与同号，故所求的根在x右侧，这时应fx()fx()fa()x,0000a=x=,而得到新的有根区间[a,b].1.25,1.5bb,,10111如此反复二分下去，二分过程无需贅述，现在预估所要二分的次数，按误差估计式，只要二分6次，便能达到预定的精度:(6)k,(2.1.1)*||0.005xx,,0二分法计算结果如表1所示表1二分法的计算结果数据表kabxfx()kkkk01.01.51.25-11.25…1.375+2…1.3751.3125-31.3125…1.3438+4…1.34381.3281+5…1.32811.3203-61.3203…1.3242-2.2牛顿迭代法2.2.1迭代法2.2.1.1迭代法的思想迭代法是一种逐步逼近的方法，首先选定方程f(x)=0的一个近似根后，然后使用某个固定公式，反复校正这个根的近似值，使之逐步精确化，一直到满足给定的精度要求为止.42.2.1.2迭代法的推理设方程有根，把方程化为等价方程fx()0,xx,(),(2.2.1)这种方程是隐式的，不能直接得出它的根，但如果给出根的某个猜测值代x0放在的右端，可得，然后，又可取x作为猜测值，进一步得到xx,,()(2.2.1)110，如此反复迭代如果按公式xx,,()21xxk,,,(),0,1,2(2.2.2)kk，1*确定的数列有极限，则称迭代过程式收敛，这时极限值xxx,lim(2.2.2),,kk,,k*显然就是方程的根.这种迭代法又称为不动点迭代法，由迭代过程所产xxx,(),生的数列并不都是收敛于某个数，与迭代方程的选取有关.2.2.1.3迭代法的误差公式假定函数定理1满足下列条件:,()xo1对任意，有xab,[,]axb,,(),o2L,1存在正数，使对任意，有xab,[,]'|()|1,xL,,(2.2.3)*则迭代过程对任意初值均收敛于方程的根，且xx,,()xab,[,]xx,()x,kk，10有如下误差估计式:*k||||/(1)xxLxxL,,,,(2.2.4)k10证明由式有(2.2.3)|||()()|||xxxxLxx,,,,,,,kkkkkk，,,111k,据此反复递推得||xx,，于是对任意正整数，有:,,Lxx||kk，110||||||||xxxxxxxx,,,，,，，,kkkkkkkk，，，,，,，,，,,,,,1121kpkpkk，,，,12，,，，，,,,,()||||/(1)LLLxxLxxL1010*在上式中令,,,，注意到limxx,，即得(2.2.4)，证毕.k，,,,,52.2.1.4迭代法的局部收敛性*'*'*定理2设为方程的根，在的邻近连续且，则迭代xxx,(),()xx,|()|1,x,*过程在邻近具有局部收敛性.xx,,()xkk，1**证明由连续函数的性质，存在的某个邻域R:，使对任意x||xx,,,'成立。此外，对任意，总有，这是因为|()|1,xL,,,()xR,xR,xR,****|()||()()|||||,,,xxxxLxxxx,,,,,,,于是，依据定理1可以断定，迭代过程对任意初值均收敛，xR,xx,,()kk，10证毕.2.2.2牛顿迭代法2.2.2.1牛顿迭代法的背景牛顿迭代法(Newton'smethod)又称为牛顿-拉夫逊方法(Newton-Raphsonmethod)，它是牛顿在17世纪提出的一种在实数域和复数域上近似求解方程的方法.多数方程不存在求根公式，因此求精确根非常困难，甚至不可能，从而寻找方程的近似根就显得特别重要.2.2.2.2牛顿迭代法的推导1通过Taylor进行理论推导xx设是f(x)=0的一个近似根，把f(x)在处泰勒展开:kk'''2fxfxfxxxfxxx()()()()()()/2!,，,，,，kkkkk若取前两项近似代替fx()，则fx()=0的近似线性方程为'fxfxfxxx()()()()0,，,,kkk'xx设fx()0,，设其根为，则的计算公式为k，1k，1kfx()k(0,1,2)k,(2.2.5),,xx，1kk'fx()k这即是牛顿法，称(2.2.5)为牛顿迭代公式，其迭代函数为fx(),,,(2.2.6)()xx'fx()2通过微分中值定理进行推导6*设是根的某个预测值，用迭代公式校正一次得，而由微分中值xxx,,()x010定理有*'*xxxx,,,,,()()10*其中介于与之间.x,x0'假定改变不大，近似地取某个近似值L，则由,()x**xxLxx,,,()101L*得xxx,,1011,,LL可以期望，按上式右端可得1LLxxxxxx,,,，,()210110111,,,LLLx是比更好的近似值.1将每得到一次改进值算作一步，并用和x分别表示第步的校正值和改进值，xkkk则加速迭代计算方案可表述如下:校正xx,(),kk，1L改进(2.2.7),，,()xxxxkkkk，，，1111,L其中xx,()中的,()x可以是多种多样的，当,()x,，xfx()时，相应的迭代公式,是xxfx,，()(2.2.8)kkk，1运用前面的加速技巧，对于迭代过程，其加速公式如下:(2.2.8),xxfx(),，kkk，1,,Lxxxx,，,(),kkkk，，，111,,1LML,,1记，上面两个式子可以合并写成fx()k,,xx，1kkM7这种迭代公式通常称为简化的公式，其相应的的迭代函数是Newtonfx(),,,(2.2.9)()xxM'L需要注意的是，由于是的估计值，而，这里的实际ML,,1,()x,()()xxfx,，''M上是的估计值，如果用代替式中的，则得如下形式的迭代函fx()fx()(2.2.9)数:fx(),,,，()xx'fx()其相应的迭代公式fx()k(2.2.10),,xx，1kk'fx()k这就是著名的公式.Newton2.2.2.3法的几何解释Newton对于方程，如果是线性函数，则对它求根是容易的，法fx()0,fx()Newton实质上是一种线性化方法，其基本思想是将非线性方程逐步归结为某种线fx()0,性方程来求解.图1与轴的交点图fx()x8*方程的解可解释为曲线与轴的交点的横坐标,(见图1)fx()0,yfx,()xx取初值，将在初值处作Taylor展开得:fx()xx00''fx()'20fxfxfxxxxx()()()()(),，,，,，00002!''取线性部分作为的近似值，有:，若，fx()fxfxxx()()()0，,,fx()0,0000则有fx()0xx,,10'fx()0类似，我们也能得到:fx()1xx,,21'fx()1这样一直下去，我们可以得到迭代序列fx()k,,xx，1kk'fx()k由上面图可知，过曲线上横坐标为的点引切线，并将该切线与X轴xPyfx,()00的交点的横坐标坐为新的近似值，类似这样下去，我们可以得一个切线方程:x1'yfxfxxx,，,()()()kkk这样求得的值必满足式，由于这种几何背景，牛顿法亦称切线法.x(2.2.10)k，12.2.2.4法的局部收敛性Newton对于一种迭代过程，为了保证它是有效的，需要肯定它的收敛性，同时考察它的收敛速度.所谓收敛速度，是指在接近收敛过程中迭代误差的下降速度.*1定义设迭代过程xx,,()收敛于方程的根，如果迭代误差xx,,()xkk，1*当时成立下列渐近关系式exx,,k,,kkek，1(0C,为常数),C,ek,,,>1则称该迭代过程是阶收敛的，特别地，=1时称为线性收敛，时称为超线,性收敛，=2时称为平方收敛.(),*xx,,(),()x定理3对于迭代过程，如果在所求根x的邻近连续，并且kk，19'*''*(1)*,,,,,,()()()0xxx,,,,,(2.2.11),()*,()0x,,,,*则该迭代过程在点邻近是阶收敛的.,x'*2证明由于，据定理可以马上断定迭代过程具有局部收,()0x,xx,,()kk，1敛性.*再将在根处展开，利用条件，则有,()xx(2.2.11)k(),,,()**,xxxx,，,()()(),,kk!,**注意到，，xx,,(),()xx,kk，1(),,,()**,由上式可得xxxx,,,()kk，1!,()*,e,()xk，1因此对于迭代误差，有，这表明迭代过程确实是阶收xx,,(),,kk，1,e!,k敛的，证明完毕.由上面定理可知，迭代过程的收敛速度依赖于迭代函数,()x的选取.对于公Newton式，其迭代函数为(2.2.10)''fx()fxfx()()',,,,,，()x()xx''2[()]fxfx()**'*'*假定是fx()的一个单根，即，则由上式知，于是依xfxfx()0,()0,,,()0x,*据定理可以断定，法在根的邻近是平方收敛的.3xNewton2.3改进牛顿迭代法2.3.1改进牛顿迭代法的背景牛顿迭代法是解非线性方程最著名和最有效的方法之一，在单根附近它比一般的迭代法有较快的收敛速度，但要注意它也有缺点:首先，它对迭代初值选取要求'fx()严格，初值选取不好，可能导致不收敛;其次，它迭代一次要计算的值，这k势必会增加计算量，因此在这种情况下，提出改进的法是非常有必要的.Newton102.3.2改进的法Newton2.3.2.1Simpson牛顿法和几何平均牛顿法设是方程的根，是可导函数，显然成立:fx()0,fx(),x'fxfxfxdx()()(),，n,xn若将上面公式的右端积分用数值积分Simpson公式近似代替，并令，则得:x,,,,，,xx'''nn(2.3.1)0()[()4()()],，，，,fxfxffnn62上式中用近似代替，整理得迭代格式(2.3.1)x,n，16()fxn(2.3.2),,xx，1nn，xx'''，1nn，，fxffx()4()()，1nn2中关于是隐式的，这给求解带来很大的麻烦，为了避免隐式求解，我们(2.3.2)xn，1提出了预估校正式:(2.3.2)fx(),nzxn,,,,0.1.2.3nn，1',fx()n,(2.3.3),6()fxnxxn,,,,0,1,2,,，1nnxz，'''nn，1,fxffz()4()()，，nn，1,2式是牛顿迭代法与Simpson公式相结合得到的，我们称它为Simpson牛顿(2.3.3)方法.xx，'''nn，1()4()()fxffx，，nn，1''2若将右端用代替，则得到,|()()|fxfx(2.3.2)nn，16迭代格式:fx()n(2.3.4),,xx，1nn'',|()()|fxfx，1nn''当时，取，当，取,,,1，为了避免隐式求解，同样给出f()0,,,,1f()0,,了式的预估校正式:(2.3.4)11fx(),n,,,zxn,0.1.2.3，1nn',fx()n,(2.3.5),fx()n,,,,xxn,0,1,2,，1nn'',,|()()|fxfz，1nn,的取法同上，我们称为几何平均牛顿方法.(2.3.5),2.3.2.2牛顿下山法*牛顿法的缺点之一是其收敛依赖与初值x的选取，若x偏离所求根较远，x00则牛顿法可能发散，为了防止迭代发散，我们对迭代过程再附加一项条件，既具有单调性:|()||()|fxfx,(2.3.6)kk，1满足这项要求的算法称为下山法.将法与下山法结合起来使用，即我们可在下山法保证函数值稳定下降Newton的前提下，用法加快了收敛速度.为此将法计算结果:NewtonNewtonfx()k,,xx，1kk'fx()kx与前一步的近似值适当加权平均作为新改进值，即kxxx,，,(1),,kkk，，11其中(01),,称为下山因子，在希望挑选下山因子时，希望使单调性条件,,成立.(2.3.6)注意:下山因子的选择是个逐步探索的过程，从开始反复将减半进行试,,,1算，如果能定出值使单调性条件成立，则称下山“成功”，与此相反，如(2.3.6),果在上述过程中找不到使条件成立的下山因子，则称“下山失败”，这时(2.3.6),x需另选初值重算.03牛顿迭代法和改进牛顿迭代法的应用3.1牛顿迭代法的应用ann例1应用牛顿迭代法和，分别导出求的迭()10f(x),x,a,0afx,,,nx12代公式并求nax,，1klim2n,,k()ax,k解方程的牛顿迭代法公式:fx()0,fx()k,,xx，1kk'fx()kfx(),,,()xx'fx()'''2''fxfxfx()()(),fx()''，,,fx()lim()0,x,()1x,,'2'2*xx,fx()fx()''*''''fxfx()()fx()''''',,,,，()()[]xfxlim()x'''**xx,fxfx()()fx()*''*''*xx,,()1()xfx*nk，1lim,,,,,xa*2'*k,,()22()xxfx,kn'1n,''2n,由题意知:fxxa(),,fxnx(),fxnnx()(1),,n,1''fxnxfxn()()1,,1,,[],''nfxxafxx()(),fx()k牛顿迭代法公式,,xx，1kk'fx()k1a1,n=,，(1)xxkknnnax,111nn,,，1klim,,,nn,,kn2axaa,2.ka'(1)(2),，,，nnn再由题意知:()1,(),()(1),,,,，fxfxanxfxannxnxa1,n，1n()fxxx,,,nx'(1),，n()fxanxan'(1),，nfxanxn()1，,,,''(2),，nfxannxx()(1),，牛顿迭代法公式13n，1xfx()kk=(1)，,,,nxxx，1kkk'anfx()k''*nax,fxn()1，，k1lim,,,'*2nn,,k2()fx()2.axa,k结论分析:本题主要应用牛顿迭代法公式和极限还有导数的相关知识，计算量大，着重考察对迭代法思想的深刻理解.32例2方程在附近有根，把方程写成3种不同的等价形x,1.5x,x,1,00式:11(1)，对应迭代格式:x,，x,1，1n，122xxn2323(2)，对应迭代格式:x,1，xx,1，xn，n1112(3)，对应迭代格式:x,x,n，1x,1x,1n讨论这些迭代格式在时的收敛性.若迭代收敛，试估计其收敛速度，选一种x,1.50收敛格式计算出附近的根到4位有效数字.(收敛速度的计算和比较)x,1.50332解，x,[1,]f(x),x,x,12331*，，故方程在[1,]上有根.x，f(1),,1,0f(),,022835539*，故方程在[,]上有根.xf(),,,04246431111149*，故方程在上有根x.[,]f(),,,0828512对于迭代式(1):28102412*3,,,(x),,,2,(),,1，，,,(x),1，(x),,23*3111331xxx2*而,(x),,,0，故该迭代局部收敛，且收敛速度为1阶的.*3x14对于迭代式(2):在上，x,[1,2]2x21/3,，(x),,(x),(1，x),22/33(1，x)*332x2x24*3,,，又,(x),,0,(x),,x,,122/3*2/33333(2x)(1，x)故该迭代在上整体收敛，且收敛速度为一阶的.x,[1,2]1对于迭代式(3):(x),在[1，2]上的值域为，该迭代式不收敛,[1,，,)x,123取迭代式，取初值进行计算，其结果如下:x,1，xx,1.5n，n01，，，x,1.4812x,1.4688x,1.4727x,1.46701324，，，x,1.4662x,1.4659x,1.4657x,1.46565678结论分析:这题主要是分析迭代法的收敛性，以及收敛速度，着重考察对迭代法的收敛性和收敛速度的理解.32例3用Newton迭代法解方程在初值附近的fxxxx()330,，,,,x,1.50根，并用数学工具软件Matlab求解.(保留小数点后6位有效数字)解由题意知初值，由牛顿迭代公式:x,1.50fx()k,,xx，1kk'fx()k代入其中得，x,1.7777781x,1.7333612x,1.7320513x,1.732051432迭代4次后，发现xx与近似相等，因此得到此方程在初fxxxx()330,，,,,34值一个根为1.732051(保留小数点6位有效数字)x,1.50下面用数学工具软件Matlab求解.首先牛顿迭代法在matlab的计算程序如下:15Functionx=newton(fname,dfname,x0,e)%用途:Newton迭代法解非线性方程f(x)=0%格式:x=nanewton(fname,dfname,x0,e)x返回数值解.%fname和dfname分别表示f(x)及其导函数%f(x),x0为迭代初值，e精度要求(默认为1e-4)Ifnargin<4,e=1e-4:%精度默认为1e-4EndX=x0;x0=x+2*e;%使while成立，进入whiler后x0得到赋值Whileabs(x0-x)>eX0=x;X0=x;X=x0-feval(fname,x0)/feval(dfname,x0);End然后，当时，f(x)>0,f(x)>0,即f(x)恒正，所以根在[0，2]，我们先用图解法找x,2初值，在用Newton法程序newton.m求解.Fun=inline('x^3+x^2-3*x-3');Fplot(fun,[0,2]);Gridon;163210-1-2-3-4-500.20.40.60.811.21.41.61.82图2求根图由图2可知方程有唯一正根在[1.6，1.8]之间，我们取初值1.5代入Newton程序之中得:Dfun=inline('3*x^2+2*x-3');Formatlong;Newton(fun,dfun,1.5,1e-4);Formatshort;Ans=1.73205080756888而用Matlab本身的函数fzero求出来的结果为:FormatlongFzero(inline('x^3+x^2-3*x-3'),1.5);FormatshortAns=1.7320508075688817例4住房是居民消费一个主要部分，大部分人选择银行按揭贷款，然后在若干年内逐月分期还款，如果你借了10万，还款额一定超过10万.解设贷款总额为，贷款期限为N个月，采取逐月等额方式偿还本息，若xx0k为第K个月的欠款数，a为月还款，r为月利率，我们得到那些列迭代关系式xrxa,，,(1)kk，1那么xrxa,，,(1)kk,12,，,，，(1)[1(1)]rxark,2kk21,,，,，，，，，，(1)[1(1)(1)(1)]rxarrr0因此得到月还款计算公式:N(1)，rx0a,N(1)1，,r下面是一则报纸在2002年2月12日第二版上一则房产广告:表2房贷数据表建筑面积总价30%首付70%按揭月还款236万10.8万30年1436元85.98m不难算出，你向银行总共借了25.2万，30年内共要还款51.96万，约为当初借款的两倍，这个案例中的贷款年利率的是多少呢,我们根据a=0.1436，=25.2，N=360，由上a的求解公式得到:x036036025.2(1)0.1436[(1)1]0rrr，,，,,我们令360360frrrr()25.2(1)0.1436[(1)1],，,，,则该问题转化为非线性方程求解的问题，令fr()0,求出r我们先用Newton函数求解，在Matlab中输入如下程序:常识1:r应比当时活期存款月利率略高一些，我们用当时的活期存款利率为18fr()0.0198/2作为迭代初值，为剔除r=0这个没有意义的根，我们对稍作变形;Clear;Fun=inline(’25.2*(1+r)^360/0.1436-((1+r)^360-1)/r’,’r’)Fun=Inlinefunction;Fun(r)=25.2*(1+r)^360/0.1436-((1+r)^360-1)/rDfun=inline(’25.2*360*(1+r)^359/0.1436-(360*(1+r)^360-1)/(r^2)’);R=newton(fun,dfun,0.0198/2,1c-4);R=12*r;然后求到结果:R=0.0553由是得到年利率为5.53%.下面我们用Matlab中的fzero函数检验一下:Clear;Fun=inline(’25.2*(1+r)^360-((1+r)^360-1)/r*0.1436’,’r’)Fun=InlinefunctionFun(r)=25.2*(1+r)^360-((1+r)^360-1)/r*0.1436R=fzero(fun,0.0198/2);R=12*rR=0.05533.2改进牛顿迭代法的应用32例1求方程的根，取初值(要求用三种方法)x,1xx，,,1000解(1)由题意知，用牛顿法公式:x,1032xx，,10kkxx,,，kk1232xx，kk将代入，迭代5次得:x,1x,1.8674605(2)由Simpson牛顿法，x,1019fx()