牛顿形式的埃尔米特插值多项式

期末论文课程名称：数值分析院系名称：巢湖学院数学系所在班级：11级数本（2）班学生学号：11020170学生姓名：张秀丽目录【题目】：牛顿形式的埃尔米特插值多项式【摘要】：.........................................................【关键词】：..........................................................【正文】：一、引言二、重节点均差与泰勒插值三、埃尔米特插值典例四、牛顿形式的埃尔米特插值多项式的一些应用领域【结束语】：.........................................................【参考文献】：..........................................................牛顿形式的埃尔米特插值多项式【摘要】：在了解了插值法以后，陆续的又接触和学习到多项式插值、拉格朗日插值、牛顿插值多项式等，但在有些实际问题中，仍需要其它要求，下面又给出有关牛顿的埃尔米特插值的内容。【关键词】：重节点均差、泰勒插值、泰勒插值多项式、埃尔米特插值。【正文】：一、引言插值法是一种古老的数学方法，它来自生产实践。早在一千多年前的隋唐时期制定历法时就应用了二次插值，隋朝刘绰将等距节点二次插值应用于天文计算。但插值理论都是在17世纪微积分产生以后才逐步发展的，牛顿的等距节点插值公式及均差插值公式都是当时的重要成果。近半世纪由于计算机的广泛使用和造船、航空、精密机械加工等实际问题的需要，使插值法在理论上和实践上得到进一步发展，尤其是20世纪40年代后期发展起来的样条插值，更获得广泛应用，成为计算机图形学的基础。在插值法的提出后我们了解了多项式插值；应用各种不同的方法对给定的插值点为求得形如01()...nnPxaaxax=+++的插值多项式我们得到了线性插值与抛物线插值；把线性插值与抛物线插值推广到一般情形，通过讨论如何构造通过n+1个节点01...nxxx的n次插值多项式()nLx，我们定义了n次插值基函数从而得到了拉格朗日插值多项式：()()nnkkkoLxylx==å。利用插值基函数很容易得到拉格朗日插值多项式，公式结构紧凑，在理论分析中甚为重要。但当插值点增减时，计算要全部重新进行，甚为不变，为了计算方便可重新设计一种逐次生成插值多项式的方法，通过一系列的考察与讨论我们利用均差得到了牛顿均差插值多项式001001201()()[,]()[,,]()()...nPxfxfxxxxfxxxxxxx=+-+--++101[,...,]()...()nnfxxxxxx---，随后还涉及了差分形式的牛顿插值公式等。插值多项式要求在插值节点上函数值相等，有的实际问题还要求在节点上倒数值相等，甚至高阶导数值也相等，满足这种要求的插值多项式称为埃尔米特插值多项式。二、重节点均差与泰勒插值先给出一个关于均差的结论。设01[,],,,...,nnfCabxxxÎ为[,]ab上的相异节点，则01[,,...,]nfxxx是其变量的连续函数。如果[,]ab上的节点互异，根据均差定义，若1[,]fCabÎ，则有00'0000()()[,]()limlimxxxxfxfxfxxfxxx-==-.由此定义重节点均差0'0000[,][,]()limxxfxxfxxfx®==.类似地可定义重节点的二阶均差，当10xx¹时，有010000110[,][,][,,]fxxfxxfxxxxx-=-.当10xx®时，有1020''00001201[,,][,,]()2limxxxxfxxxfxxxfx®®==.一般地，可定义n阶重节点的均差，由()01()[,,...,],[,]!nnffxxxabnxx=?则得0()000001()[,,...,][,,...,].!liminnxxfxfxxxfxxxn®==在牛顿均差插值多项式001001201101()()[,]()[,,]()()...[,...,]()...()nnnPxfxfxxxxfxxxxxxxfxxxxxx-=+-+--++--中，若令0(1,2,...,)ixxin?，则由0()000001()[,,...,][,,...,].!liminnxxfxfxxxfxxxn®==可得泰勒多项式()'00000()()()()()...()!nnnfxpxfxfxxxxxn=+-++-.它实际上是在点0x附近逼近()fx的一个带导数的插值多项式，它满足条件()00()(),0,1,...,.kknpxfxkn==称001001201101()()[,]()[,,]()()...[,...,]()...()nnnPxfxfxxxxfxxxxxxxfxxxxxx-=+-+--++--为泰勒插值多项式，它就是一个埃尔米特插值多项式，其余项为(1)10()()(),(,),(1)!nnnfRxxxabnxx++=-?+它与插值余项(1)1()()()()()(1)!nnnnfRxfxLxxnxw++=-=+中0(1,2,...,)ixxin?的结果一致.实际上泰勒插值是牛顿插值的极限形式，是只在一点0x处给出n+1个插值条件()00()(),0,1,...,.kknpxfxkn==得到的n次埃尔米特插值多项式.三、埃尔米特插值典例一般地只要给出m+1个插值条件（含函数值和导数值）就可造出次数不超过m次的埃尔米特插值多项式，由于导数条件各不相同，这里就不给出一般的埃尔米特插值多项式，只讨论两个典型的例子.先考虑满足条件()()(0,1,2)iiPxfxi==及''11()()Pxfx=的插值多项式及其余项表达式.由给定条件，可确定次数不超过3的插值多项式.由于此多项式通过点00(,())xfx，11(,())xfx及22(,())xfx，故其形式为001001201012()()[,]()[,,]()()...()()()nPxfxfxxxxfxxxxxxxAxxxxxx=+-+--++---，其中A为待定常数，可由条件''11()()Pxfx=确定，通过计算可得'101100121012()[,]()[,,]()()fxfxxxxfxxxAxxxx---=--.为了求出余项()()()RxfxPx=-的表达式，可设2012()()()()()()()RxfxPxkxxxxxxx=-=---，其中k(x)为待定函数.构造2012()()()()()()()tftPtkxtxtxtxj=-----，显然()0(0,1,2)jxjj==，且'1()0,()0.xxjj==故()tj在(,)ab内有5个零点（二重根算两个）.假设f具有较好的可微性，反复应用罗尔定理，得(4)()tj在(,)ab内至少有一个零点x，故(4)(4)()()4!()0fxjxx=-=k，于是(41()4!kxfx=）（），余项表达式为(420121()()()()4!Rxfxxxxxxx=---）（），式中x位于012x,,xx和x所界定的范围内.四、牛顿形式的埃尔米特插值多项式的一些应用领域目前，牛顿插值法已经运用到了工程上的各个领域，并解决了许多实际工程中遇到的问题，如物体加热时间的分析、计算；加药量自动标定；智能气体体积分数测量；自动确定支持度阈值；漏磁探测；电力系统采样；凸轮曲线的修正设计等。有时根据实际情况也会使用局部牛顿插值法，如运用局部牛顿插值法提高多狭缝自准直仪准确度。在插值问题中，要求插值多项式通过给定的数据点，但实际上所谓给定的数据本身是有误差的，而且即使插值多项式通过了给定的数据点，在这些给定数据点上的误差很小，但在其他点上的误差可能会很大，这是插值问题的缺点。在实际应用中，可以采用与曲线拟合结合等方法来达到更好的效果其中牛顿形式的埃尔米特插值多项式的应用具体表现在，为求解常微分方程数值解,运用数值积分法,采用埃尔米特插值多项式,推导出三个等距节点的六阶隐式线性多步法公式;并且对所建立公式的精度进行了分析;进一步通过实例运用计算机编程将阿达姆斯外推法等线性多步法和所建立的公式进行了精度比较。结果证明,所建立的隐式线性多步法公式比现有的具有相同节点的线性多步法公式精度更高,求解速度更快,有一定的应用价值。解决了初等数学中的一类较为复杂的求函数值、求范围、作证明的相关问题。通过研究具有任意阶导数信息Hermite插值问题,使用广义差商的一种新的表示方法和构造广义差商表的一种新方法,给出具有任意阶导数信息Hermite插值算法和程序实现,拓展了牛顿差商插值公式和余项公式科学计算中常用的计算方法，其内容包括误差的概念，插值方法，线性代数方程组的解法，非线性方程的求根，数值积分与数值微分，最小二乘法，特征值的计算，常微分方程初值问题的数值解法等五、结束语插值法是函数逼近的一种重要方法，它是数值微分、微分方程数值等数值的基础与工具。由于多项式具有形式简单，计算方便等许多优点，故本文主要介绍插值多项式——牛顿形式的埃尔米特插值多项式，它是插值法中常用和最基本的方法。【参考文献】：《数值计算方法及其应用》，作者：朱长青，出版社：科学出版社，出版日期：2006,01；《数值分析与应用程序》，作者：全惠云，出版社：武汉大学，出版日期：2007,04；《数值分析》第五版，作者：李庆扬、王能超、易大义，出版社：清华大学出版社，《数值计算引论》，作者：白峰杉，出版社：高等教育出版社，出版日期：2004；《计算机方法简明教程》，作者：王能超，出版社：高等教育出版社，出版日期：2004.