计算方法论文浅谈拉格朗日插值法.docVIP

盐城师范学院课程考查论文 课程名称：《计算方法》 学 院： 专 业： 班 级： 姓 名： 学 号： 论文题目：浅谈拉格朗日插值法 成绩： 浅谈拉格朗日插值法 插值法又称“内插法”，是利用函数f (x)在某区间中若干点的函数值，作出适当的特定函数，在这些点上取已知值，在区间的其他点上用这特定函数的值作为函数f (x)的近似值，这种方法称为插值法。如果这特定函数是多项式，就称它为插值多项式。 下面将主要谈谈拉格朗日插值法。 一、 问题的背景在实际问题中常遇到这样的函数y=f(x),其在某个区间[a,b]上是存在的。但是,通过观察或测量或试验只能得到在区间[a,b]上有限个离散点x0 ,x1 ,…,xn 上的函数值yi =f(xi ) ,(i=0,1,…,n) 。 或者f(x)的函数f(x)表达式是已知的,但却很复杂而不便于计算;希望用一个既能反映函数f(x)的特性，又便于计算的简单函数来描述它。 ??二、插值问题的数学提法:已知函数在n+1个点x0 ,x1 ,…,xn 上的函数值?????? yi =f(xi ), (i=0,1,…,n) 求一个简单函数y=P(x),使其满足:??????? P(xi )=yi ,(i=0,1,…,n) 。即要求该简单函数的曲线要经过y=f(x)上已知的这个n+1个点： ????? (x0 ,y0 ),(x1 ,y1 ),…,(xn ,yn ),同时在其它x∈[a,b]上要估计误差:???????????? R(x) = f(x) - P(x) 其中P(x)为f(x)的插值函数，x0 ,x1 ,…,xn 称为插值节点，包含插值节点的区间[a,b] 称为插值区间，求插值函数P(x)的方法称为插值法。若P(x)是次数不超过n的代数多项式，就称P(x)为插值多项式，相应的插值法称为多项式插值。若P(x)是分段的多项式，就是分段插值。若P(x)是三角多项式，就称三角插值。 下面是对拉格朗日插值法的介绍： ?已知函数y=f(x)在n+1个不同的点x0?,x1?,…,x2?上的函数值分别为： y0?,y1?,…,yn?,求一个次数不超过n的多项式Pn?(x),使其满足:?Pn?(xi?)=yi?, (i=0,1,…,n),即n+1个不同的点可以唯一决定一个n次多项式。 （1）. 插值基函数过n+1个不同的点分别决定n+1个n次插值基函数??????????l0?(x),l1?(x),…,ln?(X)每个插值基本多项式li?(x)满足：? (1)?li?(x)是n次多项式; (2)?li?(xi?)=1,而在其它n个li?(xk?)=0 ,(k≠i)。由于li?(xk?)=0 ,(k≠i), 故有因子:???? (x-x0?)…(x-xi-1?)(x-xi+1?)…(x-xn?)因其已经是n多项式，故而仅相差一个常数因子。令:????????? li?(x)=a(x-x0?)…(x-xi-1?)(x-xi+1?)…(x-xn?)由li?(xi?)=1,可以定出a, 进而得到： （x）= （2）.?n次拉格朗日型插值多项式Pn?(x)Pn?(x)是n+1个n次插值基本多项式l0?(x),l1?(x),…,ln?(X)的线性组合,相应的组合系数是y0?,y1?,…,yn?。即:?????????Pn?(x)=y0?l0?(x)+y1?l1?(x)+…+yn?ln?(x)?,从而Pn?(x)是一个次数不超过n的多项式,且满足 ??????Pn?(xi?)=yi?, (i=0,1,2,…,n). 应用程序设计: 一方面利用计算机求拉格朗日插值多项式的值编程比较容易，另一方面，许多应用程序中可以利用插值的方法“试探”一下，看能否得到更精确的结果，所以通常的应用只是对已经得到的一些点处得函数值再利用插值法计算1个或少量几个点处得插值的结果。在这种情况下，有一个便于应用的程序比起节省计算量来说显得更为重要，所以在是实际应用中，应用拉格朗日插值多项式进行计算的机会也相当多。下面是拉格朗日插值法用VB编的程序。 Private Sub Form_Load()Dim X(5) As Double, Y(5) As DoubleDim X0 As Double X(1) = 5: X(2) = 7: X(3) = 10X(4) = 14: X(5) = 19 Y(1) = 35: Y(2) = 47: Y(3) = 56Y(4) = 78: Y(5) = 87 X0 = 15 PrintPrintPrint ??五个点 (5,35)，(7,47)，(10,56)，(14,78)，(19,87)插值X0=15 Print ????