泰勒公式在数值分析中的应用讲解

1、2015年度本科生毕业论文（设计）文山学院WENSHANUIMIVE口SITY泰勒公式在数值分析中的应用教学系：数学学院专业：数学与应用数学年级：11级数本（3）班姓名：袁国彦学号：20110701013056导师及职称：程高讲师2015年05月毕业论文（设计）原创性声明本人所呈交的毕业论文（设计）是我在导师的指导下进行的研究工作及取得的研究成果。据我所知，除文中已经注明引用的内容外，本论文（设计）不包含其他个人已经撰写或发表过的研究成果。对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中作了明确说明并表示谢意。作者签名：日期：毕业论文（设计）授权使用说明本论文（设计）作者完全了解文山

2、学院有关保留、使用学生毕业论文（设计）的规定，学校有权保留论文（设计）并向相关部门送交论文（设计）的电子版和纸质版。有权将论文（设计）用于非赢利目的的少量复制并允许论文（设计）进入学校图书馆被查阅。学校可以公布论文（设计）的全部或部分内容。保密的论文（设计）在解密后适用本规定。作者签名：指导教师签名：日期：日期：袁国_毕业论文（设计）答辩委员会（答辩小组）成员名单姓名职称单位备注文山学院本科毕业论文（设计）摘要泰勒公式是微积分中一个重要的公式，它将一些复杂的函数近似的表示为多项式函数，为一些复杂函数的求解带來方便。不仅在数学分析中有着重要的地位，在数值分析中也有着广泛的应用，本文简要介绍了泰勒

3、公式在数值分析中的应用，并讨论泰勒公式在泰勒插值，欧拉方法和牛顿迭代法中的具体应用，在泰勒插值和数值积分中，用泰勒公式展开的多项式去逼近原函数，得出近似解，并分析误差。欧拉方法是通过迭代的方法，求得近似值，通过用不同的步长进行对比，并得到一种通过控制误差来得到步长的方法。牛顿迭代法是求解非线性方程近似解的一种方法，通过程序來得到方程根所在的区间，求出初值，最后控制其误差。泰勒公式需要先取点对原式进行泰勒展开，如何选取，使得泰勒公式展开后，计算的结果在误差的允许范围内，并且计算过程尽量简单，减少计算步骤。关键词：泰勒展开；泰勒插值；数值积分；欧拉方法；牛顿迭代法；数值分析Theapplicati

4、onofTaylorformulainnumericalanalysisABSTRACTTaylorformulaisanlmpoitantformulainCalculus,ItwillbesomefiinctionapproxmiationisexpressedasapolynomialfunctionNotonlyplaysanlmpoitantroleinmathematicalanalysis,anditiswidelyusedinthenumeiicalanalysis,thispaperbneflvintioducestheapplicationofTayloiformulain

5、numeiicalanalysis,anddiscusstheTaylorformulaintheapplicationofTaylorliiteipolation,EulermethodandNewtoniterationmethod,Tayloiinterpolation,polynomialusingtheTavloiexpansiontoapproximatetheonginalfiinction,theappioximatesolutionandeiroranalysis.TheEulermethodisobtainedbyiterativemethod,approximateval

6、ue,comparedtothediffeientstepsize,andamethodtogetlinnstepbycontrollingtheenor.TheNewtoniterativemethodisamethodofappioximatesolutionforsolvingnonlmearequations,tluoughthepiogramtogettherangeofequationroot,andtheenoicontrol.NeedtoselectapouitontheongmalTavloi;howtoselect,theTayloiexpansion,thecalcula

7、tionresultsmtherangeofallowableenoiinthecalculationoftheprocessassimpleaspossible,andtoreducethecomputationalsteps.Keywords:Tavloiexpansion;Tayloimterpolation;Numeiicalintegration;Eulefsmethod;TheNewtoniterationmethod;Numeiicalanalysis目录一、引言1二、泰勒公式的应用32.1泰勒插值-32.2泰勒公式在数值积分中的应用62.3欧拉方法-72.4用泰勒公式求方程根的

8、近似解-92.4.1牛顿迭代法92.4.2扫描法102.4.3误差估计公式10参考文献12致谢13附录14文山学院本科毕业论文（设计）泰勒公式的背景：希腊人在理性数学活动中，已接触到了无限性、联系性等概念，这方面最具有代表性的人物是伊利亚学派的芝诺。他在考虑利用无穷级数求和來得到有限结果时，提出了四个著名的悖论。后來，随着无限小算法的推广，英国的数学家们在大学里教授和研究牛顿的流数术,他们中优秀的代表有泰勒和麦克劳林。泰勒在1715年出版的正的和反的增量方法一书中，陈述了它早在1712年就已获得的著名定理，这就是为人所熟知的泰勒级数。爱丁堡大学教授麦克劳林发现了泰勒级数的特例，称为“麦克劳林级

9、数”。泰勒公式的推导：由导数和微分的概念，如果函错误味找到引用源。在点错误!未找到引用源。可导，则有fCO二f（%）+尸COO咒o）+。（兀一1即在点错误味找到引用源。附近，用一次多项式错误味找到引用源。逼近函数值错误!未找到引用源。时，其误差为错误!未找到引用源。的高阶无穷小量。然而在很多场合,取一次多项式的逼近是不够的，往往需要用二次或高于二次的多项式去逼近，并要求误差为错误!未找到引用源。，其中错误!未找到引用源。往往为多项式的次数，为此，我们考察任-错误!未找到引用源。次多项式错误!未找到引用源。.错误!未找到引用源。逐次求它在点错误!未找到引用源。的各阶导数，得到PnM=%，错误!未

10、找到引用源。，错误味找到引用源。，错误味找到引用源。，错误!未找到引用源。.由此可见，多项式错误味找到引用源。的各项系数由其在点错误味找到引用源。的各阶导数值所唯一确定。对于一般函数错误!未找到引用源。，设它在点错误味找到引用源。存在直到错误!未找到引用源。阶的导数，由这些导数构造一个n次多项式错误!未找到引用源。，错误!未找到引用源。称为函数错误味找到引用源。在点错误味找到引用源。处的泰勒（错误味找到引用源。）多项式，错误!未找到引用源。的各项系数错误!未找到引用源。称为泰勒系数。由上面对多项式系数的讨论，易知错误!未找到引用源。与其泰勒多项式错误!未找到引用源。在点错误味找到引用源。有相同

11、的函数值和直至错误!未找到引用源。阶导数值，即文山学院本科毕业论文(设计)错误!未找到引用源。（3）下面将要证明错误味找到引用源。，即以错误味找到引用源。式所示的泰勒多项式逼近错误!未找到引用源。时，其误差为关于错误!未找到引用源。的高阶无穷小量。定理1.2：若函数错误!未找到引用源。在点错误!未找到引用源。存在直至错误!未找到引用源。阶导数，则有错误味找到引用源。，即fM=f（力）+'J（X一%0）+fJ'（%一咒qF+-（%-1!Z!n证：错误味找到引用源。，错误味找到引用源。,现在只要证由关系式错误味找到引用源。可知,心（）=或（心）=6并易知QnCxo）=Qn（X0）=

12、0$一=0，=汕因为错误味找到引用源。存在，所以点错误味找到引用源。的某邻域错误!未找到引用源。内错误沬找到引用源。存在错误!未找到引用源。阶导函数错误味找到引用源。当错误!未找到引用源。且错误!未找到引用源。时，允许接连使用洛必达法则错误!未找到引用源。次，得乂QnM戈七。Qn（%）力S'limM七o/"a（Q严"。如疋%）&心）n（n1）2（%xQ）1lim7lz心广7（兀）广一巩衍）仗一心）=0其中泰勒公式（4）在错误沬找到引用源。时的特殊形式:凡龙)=代0)+八0)"爭疋户)(0)711X11+o(xn).(6)它也称为(带有佩亚诺余项的)

13、麦克劳林公错误!未找到引用源。泰勒公式是用一个函数在某点的信息描述其附近取值的公式，如果函数足够光滑的话，在已知函数在某一点的各阶导数值的情况下，泰勒公式可以用这些导数值做系数构建一个多项式來近似函数在这一点的领域中的值，泰勒公式还给出了这个多项式和实际函数值之间的偏差。数值计算中泰勒公式有广泛的应用，泰勒公式的证明与应用方面对于研究者来说一直具有吸引力，其理论方法已经成为研究函数极限和估计误差方面不可或缺的数学工具，在近似计算上有着独特的优势，利用它可以满足很高的精度要求。泰勒公式可以应用于求极限，判断函数极值，求函数在某些点的数值，近似计算等方面。二、泰勒公式的应用2.1泰勒插值实际问题中

14、碰到的函数错误!未找到引用源。是各种各样的，有的表达式很复杂，直接研究函数错误!未找到引用源。可能很困难，面对这种情况，一个很自然的想法是将函数错误味找到引用源。简单化，构造某个简单的函数错误!未找到引用源。作为错误!未找到引用源。的近似函数，通过处理错误!未找到引用源。获得关于错误!未找到引用源。的结果，如果要求近似函数错误!未找到引用源。取给定的离散数据，则称之为错误!未找到引用源。的插值函数。其中泰勒公式展开公式开方法就是一种插值方法，由于代数多项式的结构简单，数值分析方面就相对简单。己知泰勒多项式%CO=代咒(?)+'V(%-x0)+亍(%-x&)2+(%-%0)K.J

15、L!Zj!成立。求作错误!未找到引用源。次多项式错误!未找到引用源。，使其满足条件错误!未找到引用源。，错误!未找到引用源。这里错误!未找到引用源。为一组已给出的数据。容易看出，对于给定的函数错误味找到引用源。，若导数值错误!未找到引用源。已给,则上述泰勒插值的问题的解就是泰勒多项式。运用泰勒公式做近似计算时，一般要用到带有拉格朗日余项的泰勒展开。例2.1.1：求作错误!未找到引用源。在节点错误!未找到引用源。的一阶和二阶泰勒多项式，计算错误!未找到引用源。的近似值，估计误差并与精确值0.723805对比。解：用MATLAB程序求岀错误!未找到引用源。在节点的一阶和二阶泰勒多项式，相关程序见附

16、录1。所以错误!未找到引用源。在错误!未找到引用源。的一次泰勒多项式是：P1W=fM+1!0(兀一兀0)=5+0.05X二次泰勒多项式是：错误!未找到引用源。二错误!未找到引用源。用错误!未找到引用源。作错误味找到引用源。的近似表达式，将错误味找到引用源。代入一次泰勒多项式得：<715=/(x)站血仗)=10.75根据定理1可估计出误差：fM一PiM=匚字(力严>匚辭(x一x0)2=血(x一100)2=-0.028125错误!未找到引用源。与精确值比较，误差为错误!未找到引用源。，具有3位有效数字。二次泰勒多项式的值：卩23)=fOo)+了'J咒一xo)+了(%-x0)2=

17、5+0.05兀-(x-100)2=10.721875这个结果有4位有效数字。我们对兀取的不同取值，通过作图对它们的逼近效果进行对比，程序见附录2。-5-10.910.810.710101010101021105110115图2-1-1图中代表二次泰勒多项式的值，代表一次泰勒多项式的值，“+”代表精确值，可以看出二次的泰勒公式展开逼近的效果更好，且衍越逼近,误差越小。定理错误!未找到引用源。：若函数错误!未找到引用源。在错误!未找到引用源。上存在直至“阶的连续导数，在错误!未找到引用源。内存在错误!未找到引用源。阶导数函数，则对任意给定的错误!未找到引用源。，至少存在一点错误!未找到引用源。，使

18、得:“、“.厂(龙、2.严司仗°),w/(x)=+_(x-x0)+(x-x0)+(-知)证：作辅助函数：错误!未找到引用源。，G=所要证明的定理为：错误!未找到引用源。或错误!未找到引用源。不妨设错误!未找到引用源。，则错误味找到引用源。与错误!未找到引用源。在错误!未找到引用源。上连续，在错误!未找到引用源。内可导，且巩衍)G(t)'=-(n+1)(%一丁羊0乂因为错误!未找到引用源。，所以由柯西中值定理得：S心）-f（q_y严）G®q）GOq）G&）一G（t）r（-H+1）!其中错误!未找到引用源。，他的余项为：（戈）=f（戈）一為仗）=（料+）!（咒一

19、xo）k+1=%Q+e（x-x0）（Q<e<i）如果错误味找到引用源。，错误味找到引用源。为定数，则取余项不会超过错误!未找到引用源。，从而可以近似地计算某些数值且估计误差。例2丄2：计算错误!未找到引用源。麦克劳林公式展开，并计算错误味找到引用源。的近似值，使其误差不超过错误!未找到引用源。和错误!未找到引用源。解：错误!未找到引用源。当错误味找到引用源。时，错误味找到引用源。取错误!未找到引用源。，便有33爲<一=<10"6910!3628800略取错误!未找到引用源。求得的错误!未找到引用源。的近似值为：111e1+11处2.7182852!3!9!取错

20、误味找到引用源。时，3<«4.81XIO"10<10"9略取错误味找到引用源。求得的错误味找到引用源。的近似值为:111e.1+1+-+_+.+_站2.718281828由上面可以看出如果采用更高次的多项式來逼近错误!未找到引用源。，能在更大范围内满足同一误差，但同时也增大了计算量，所以在计算时应选择适当的阶数。2.2泰勒公式在数值积分中的应用设F（Q为的一个原函数，由牛顿一莱布尼茨公式知，定义在区间心b上的定积分为：文山学院本科毕业论文（设计）f(x)dx=F(x)f-丿ca但有些原函数不能用初等函数表达或者有的原函数十分复杂难以求出或计算，不能用上

21、述公式。理论上定积分是一客观存在的确定的数值，要解决的问题是能否找到其他途径來解决定积分的近似计算。泰勒公式是一件很好的工具，它可实现定积分的近似计算。例221：解定积分的近似值，当取6阶导数值时误差范围是多少？解：定积分的被积函数不可积，可用泰勒公式将其展为幕级数，然后逐项积分，再利用积分后的级数计算。因为：COSA：=1-+-Q-0.23958,2x2：4x4!因为|cos(0x+|n-)|<l,所以此时误差:1R<0.0002316X6!2.3欧拉方法欧拉方法是求解常微分方程初值问题错误!未找到引用源。的重要方法，下面由泰勒公式导出欧拉公式错误味找到引用源。，用泰勒公式将错误

22、味找到引用源。在错误味找到引用源。处展开有兀2y(xn+i)=y(s)+矽牝)+y"(0取右端前两项得错误!未找到引用源。，设用错误!未找到引用源。的近似值错误!未找到引用源。代入上式的右端，记所得结果为错误!未找到引用源。，这样导出的计算公式3.1这就是众所周知的欧拉格式，若初值错误!未找到引用源。是己知的，则可以根据公式逐步算出数值解错误!未找到引用源。，错误!未找到引用源。，错误!未找到引用源。为简化分析，人们常在错误!未找到引用源。的前提下估计误差错误!未找到引用源。，这种误差称为局部截断误差。称一种数值方法的精度是错误!未找到引用源。阶的，则其局部截断误差为错误!未找到引用

23、源。：y(xn+i)-yg)=9"(0这说明欧拉格式仅为一阶方法。例2.3.1：取步长错误味找到引用源。和错误味找到引用源。，用欧拉方法求解初值问题(,2%/、y=y(0<%<1)y(o)=1的解错误!未找到引用源。在点错误味找到引用源。的近似值，即错误味找到引用源。的近似值，并与精确解错误!未找到引用源。进行比较，分析这两种步长的稳定性。解:将程序结果整理成表格，程序见附录3。表2-3-1错误!未找到引用源。的近似值，绝对误差反-fCOI011.10000.0046021.20000.0168021.19180.0086031.27740.0125041.37330.0

24、317041.35820.0166051.43510.02090.61.53150.0483061.50900.0257071.58030.0311081.64980.0373081.68110.0686091.71780.0445101.82690.0949101.78480.0527由表格可以看出当步长错误味找到引用源。时，欧拉方法计算错误味找到引用源。的近似值为1.7848,与精确解的误差绝对值为0.0527o当错误!未找到引用源。时，欧拉方法计算错误!未找到引用源。的近似值为1.8269,与精确解的误差绝对值为0.0949。所以步长错误!未找到引用源。时计算更稳定。但从每一步来看，步长

25、越小，截断误差就越小；但随着步长的缩小，在一定求解范围内所要完成的步数就增加了。步数的增加不但引起计算量的增大，而且可能导致舍入误差的严重积累。所以微分方程的数值解法也需要选择步长。先以错误!未找到引用源。为步长，求出一个近似解，记错误!未找到引用源。，然后将步长取半，即取错误味找到引用源。为步长从错误味找到引用源。跨两步到错误！未找到引用源。在求得一个近似解错误!未找到引用源。这样可以通过检查步长折半前后两次计算结果的偏差來判断所选取的步长是否合适，若要求求初步误差为0.01,则程序见附录4,取其中的儿项进行对比：表2-3-2错误!未找到引用源。的近似值，绝对误差ly«.-yOJ0

26、.06251.06250.00180.51.42670.013411.76650.0344得到步长错误!未找到引用源。，并且初次迭代的误差控制到0.0018c2.4用泰勒公式求方程根的近似解2.4.1牛顿迭代法牛顿公式的导出:设方程错误味找到引用源。的近似根为错误味找到引用源。，则错误味找到引用源。在点错误!未找到引用源。附近可用一阶泰勒多项式p（%）=f（%+f0（咒一代替，故方程错误!未找到引用源。可以近似的表示为错误!未找到引用源。，后者是个线性方程，它的求根是容易的，我们取错误!未找到引用源。的根作为错误!未找到引用源。的新近似根，记错误!未找到引用源。，则有：这就是著名的牛顿公式，相

27、应的迭代方程为错误味找到引用源。，其中仏）牛顿法是一种逐步线性化方法，这种方法的基本思想是将非线性方程错误!未找到引用源。的求根问题归结为计算一系列线性方程错误!未找到引用源。十错误!未找到引用源。=0的根。牛顿法有明显的儿何解释，方程错误!未找到引用源。的根错误!未找到引用源。在儿何上解释为曲线错误味找到引用源。与错误味找到引用源。轴的交点的横坐标。设错误味找到引用源。是根错误!未找到引用源。的某个近似值，对曲线错误味找到引用源。上横坐标为错误味找到引用源。的点错误味找到引用源。引切线，设该切线与错误!未找到引用源。轴的交点的横坐标记为错误!未找到引用源。（见图），则这样获得的错误!未找到引

28、用源。即为按牛顿法求得的近似根，由于这种儿何背景，牛顿法亦称切线法。-11-图2-4-12.4.2扫描法扫描法是一种在计算机上较实用的方法，简单的来说就是将有根区间分为若干个子区间，然后从有根区间的左端点开始，一个一个小区间检查是否是隔根区间。对于代数方程fM=aoxn+的咒九7+ax+a”=0（a0工0）记错误味找到引用源。，则其根的上、下界分别为错误味找到引用源。和错误味找到引用源。，由此即可确定其有根区间错误！未找到引用源。，其中二），b=1+2.43误差估计公式估计以错误味找到引用源。作为错误味找到引用源。的近似值的误差，由中值定错误!未找到引用源。错误!未找到引用源。，错误!未找到引

29、用源。，因而记错误味找到引用源。，则例2.4.1：:方程错误!未找到引用源。有儿个根，并求出最小正根的近似解，是误差不超过错误!未找到引用源。泰勒公式在数值分析中的应用解：（1）首先用扫描法求出根错误!未找到引用源。属于有根区间错误!未找到引用源。，错误!未找到引用源。和错误味找到引用源。，所以方程错误!未找到引用源。有3个根，且在错误!未找到引用源。上有最小正根，见附录5-1。（2）在错误味找到引用源。上求出根的初值，即0和1谁是根的初值，程序见附录5-2,则求的初值为错误!未找到引用源。（3）错误!未找到引用源。在错误!未找到引用源。上的最小值为错误!未找到引用源。，则可以跟据误差公式和牛

30、顿迭代法求出满足误差条件的近似根，程序见附录5-3。结果为：xO=031250000000000xO=0.25210736468500xO=0.25121933490266该程序能求出根所在的区间及初值，并对误差进行控制，从而得到误差允许范围内的近似根，且第三个根的符合误差的精度要求时程序停止运算。-#-文山学院本科毕业论文（设计）参考文献1 华东师范人学数学系数学分析M.北京：高等教育出版社，2001,134-1582 张春红.泰勒公式在近似计算及数值积分中的应用J.黑龙江科学，2014,第5卷第七期，1373 孔珊珊.泰勒公式在数值计算中的应用J.济宁学院学报，2011,第32卷第三期，7

31、0-724 王能超.数值分析简明教程M.北京:高等教育出版社，2003,135-1365 孙凤芝.数值计算方法与实验M.哈尔滨：黑龙江出版社，2013,25-26-19-致谢本论文是在我的导师程老师的亲切关怀下完成的，老师渊博的专业知识，精益求精的工作态度对我影响深远。从课题的选择到论文的最终完成，每一步都是在程老师的悉心指导下完成的，倾注了导师大量的心血，在此我向他表示衷心的谢意。这四年來感谢文山学院数学学院的老师对我的培养，他们在学业上的细心指导为我打下了良好的基础，在这里我要向诸位老师深深鞠上一躬。时光匆匆流去，转眼毕业在即，从开始进入课题到论文的完成都离不开老师、同学给我的热情帮助，很

32、庆幸我遇到了如此多的良师益友，无论在学习上、生都给予了我无私的帮助，谨以最朴实的话语致以敬意。附录1:输入程序：clear;symsx;f=xA(1/2);tl=taylor(f,2,100,x)t2=taylor(fz3,100zx)结果：tl=5+l/20*xt2=5+l/20*x-l/8000*(x-100)A2附录2：输入程序：clear;x=100:1:115;f=sqrt(x);t2=5+0.05*x;t3=5+005*x-(1/8000)*(x-100)A2;plot(x,f,'+o')结果：附录3：输入程序：先建立euler.m文件functionxOzyO=e

33、uler(fun,xO,xlzyOzn);h=(xl-xO)/n;fori=l:n;yl=yO+h大(fun(xOzyO)y=(l+2*(xO+h)A(l/2)jd=abs(yl-y)yO=yl;xO=xO+h;end然后取步长错误味找到引用源。和错误!未找到引用源。，即错误!未找到引用源。和错误!未找到引用源。(1)错误!未找到引用源。时输入程序：fun=inline(1yO-2大xO/yO1);xOzyO=euler(funz0,lzlz10)结果：yl=1.1000y=1.0954jd=0.0046yl=1.1918y=11832jd=0.0086yl=1.2774y=12649jd=0

34、.0125yl=1.3582y=1.3416jd=0.0166yl=1.4351y=1.4142jd=00209yl=1.5090y=14832jd=0.0257yl=1.5803y=1.5492jd=0.0311yl=1.6498y=1.6125jd=0.0373yl=1.7178y=1.6733jd=0.0445yl=1.7848y=17321jd=0.0527x0=1.0000yO=1.7848（2）错误!未找到引用源。时fun=inline(1yO-2大xO/yO1);xO,yO=euler(fun,0z1,1,5)结果：yl=1.2000yl=1.3733yl=1.5315yl=1.6811yl=1.8269xO=1y=1.1832y=1.3416y=1.4832y=