几种插值方法比较与应用2.docxVIP

几种插值方法的比较与应用 摘要：本文是对学过的插值方法进行了总结、比较，使我们在进行工程计算的过程中更清楚的知道 哪一种方法适合哪一种类型，了解哪种方法在已知条件下可以得到更优的结果以满足计算要求。 关键词：数值分析，插值，多项式 1前言 在许多实际问题及科学研究中，因素之间往往存在着函数关系，然而，这种关系经常很难有明 显的解析表达，通常只是由观察与测试得到一些离散数值。有时，即使给出了解析表达式，但却由 于表达式过于复杂，不仅使用不方便，而且不易于进行计算与理论分析。解决这类问题的方法有两 种：一种是插值法，另一种是拟合法。插值法是一种古老的数学方法，它来自生产实践。早在一千 多年前，我国科学家在研究历法上就应用了线性插值与二次插值，但它的基本理论确实在微积分产 生之后才逐渐完善的，其应用也日益增多，特别是在计算机软件中，许多库函数，如 sinx,cosx,ex等 的计算实际上归结于它的逼近函数的计算。逼近函数一般为只含有算术运算的简单函数，如多项式、 有理分式(即多项式的商)。在工程实际问题当中，我们也经常会碰到诸如此类的函数值计算问题。 被计算的函数有时不容易直接计算，如表达式过于复杂或者只能通过某种手段获得该函数在某些点 处的函数值信息或者导数值信息等。因此，我们希望能用一个“简单函数”逼近被计算函数，然后 用该简单函数的函数值近似替代被计算函数的函数值。这种方法就叫插值逼近或者插值法。插值法 要求给出函数f x的一个函数表，然后选定一种简单的函数形式，比如多项式、分段线性函数及三 角多项式等，通过已知的函数表来确定一个简单的函数 x作为f x的近似，概括地说，就是用 简单函数为离散数组建立连续模型。 2插值法的基本概念 2.1插值法的定义 设函数y = f x在区间［a,b ］上有定义，且已知在点a _兀一…_ xn _b上得值 f x = yi(i =0,1/ , n)，若存在一个简单函数 x，使得 ;：K 二 yMi =0,1/ ,n) 成立，就称 x为f x的插值函数，点xji =°，1，…，n)为插值节点，包括插值节点的区间 La,bl成 为插值区间，求插值函数 ‘ x的方法成为插值法。 若「x为次数不超过n的代数多项式 x 二a。盼 anXn 其中的a (i =0,1/ ,n)为实数，就称；x为插值多项式，相应的插值法称为多项式插值。 若■； x为分段的多项式，就称为分段插值。 22截断误差 截断误差(余项)：若在la,b ］上用：n x近似f x，则 Rn X］=f X?-申 n x 称为插值多项式的截断误差，又称为插值多项式的余项。 代数插值法有 Lagra nge插值法、逐次线性插值法、 Newt on插值法、Hermite插值法、分段插值 法等。其基本思想都是用高次代数多项式或分段的低次代数多项式作为被插函数 f x的近似表达 式。 3几种常见的代数插值法 3.1 Lagrange 插值法 设函数y=f x在区间la,b ］上有定义，且在 la,b ］上n 1个不同点乂0,为，…，xn上得函数值 y^f k (i=o,l「n )，若存在一个至少 n次的插值多项式 它们的次数不超过 n,且满足；x 二a。 a/ a 它们的次数不超过 n,且满足 其中ai(i =0,1; ,n)为实数。先构造函数li x (i =0,1,…，n), 然后以对应点处的函数值为系数作线性组合，即得所要求的多项式。由多项式l 然后以对应点处的函数值为系数作线性组合，即得所要求的多项式。 由多项式 li x有n个根 Xj( j =0,1厂，i - ij ? 1, n,)故它必有如下形式 n X _ XjL（i =0,1； ,n）j =0 xi _ Xjj=Ex） _（X -怡） （X - Xi」）（x - Xj 1）（X - Xn） （Xi — X。）…（xi — Xi」）（Xj — Xj n X _ Xj L（i =0,1； ,n） j =0 xi _ Xj j=E 这些函数称为 Lagrange插值基函数，而 n x 是至多 n次多项式，且满足 n ：：n(Xk) =7 y』(xQ 二 yk (k =0,1,2，…，n)称其为 n 次 Lagrange 插值多项式。 i ￡ 3.2 Newton插值法 f (x?)一f (x■) 设有函数f(x),心人公2,…，为一系列互不相等的点，称 f［X,Xj］ 1 j (i = j)为f X 备_Xj 关于点x「Xj的一阶差商。一般的称 fix。/!,,兀］二f［xo,xi; ,xk」］ f［xo,xi，X］为f x关于 XoM , x的k阶差商。 表1差商表 X f(x) 一阶差商 二阶差商 三阶差商 Xo Xi X2 X3 f(Xo) f (Xi) f(X2) f(X3