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

几种插值方法的比较与应用 摘要：本文是对学过的插值方法进行了总结、比较，使我们在进行工程计算的过程中更清楚的知道哪一种方法适合哪一种类型，了解哪种方法在已知条件下可以得到更优的结果以满足计算要求。 关键词：数值分析，插值，多项式 1 前言 在许多实际问题及科学研究中，因素之间往往存在着函数关系，然而，这种关系经常很难有明显的解析表达，通常只是由观察与测试得到一些离散数值。有时，即使给出了解析表达式，但却由于表达式过于复杂，不仅使用不方便，而且不易于进行计算与理论分析。解决这类问题的方法有两种：一种是插值法，另一种是拟合法。插值法是一种古老的数学方法，它来自生产实践。早在一千多年前，我国科学家在研究历法上就应用了线性插值与二次插值，但它的基本理论确实在微积分产生之后才逐渐完善的，其应用也日益增多，特别是在计算机软件中，许多库函数，如等的计算实际上归结于它的逼近函数的计算。逼近函数一般为只含有算术运算的简单函数，如多项式、有理分式（即多项式的商）。在工程实际问题当中，我们也经常会碰到诸如此类的函数值计算问题。被计算的函数有时不容易直接计算，如表达式过于复杂或者只能通过某种手段获得该函数在某些点处的函数值信息或者导数值信息等。因此，我们希望能用一个“简单函数”逼近被计算函数，然后用该简单函数的函数值近似替代被计算函数的函数值。这种方法就叫插值逼近或者插值法。插值法要求给出函数的一个函数表，然后选定一种简单的函数形式，比如多项式、分段线性函数及三角多项式等，通过已知的函数表来确定一个简单的函数作为的近似，概括地说，就是用简单函数为离散数组建立连续模型。 2 插值法的基本概念 2.1 插值法的定义 设函数在区间上有定义，且已知在点上得值，若存在一个简单函数，使得 成立，就称为的插值函数，点为插值节点，包括插值节点的区间成为插值区间，求插值函数的方法成为插值法。 若为次数不超过n的代数多项式 其中的为实数，就称为插值多项式，相应的插值法称为多项式插值。若为分段的多项式，就称为分段插值。 2.2 截断误差 截断误差（余项）：若在上用近似，则 称为插值多项式的截断误差，又称为插值多项式的余项。 代数插值法有Lagrange插值法、逐次线性插值法、Newton插值法、Hermite插值法、分段插值法等。其基本思想都是用高次代数多项式或分段的低次代数多项式作为被插函数的近似表达式。 3 几种常见的代数插值法 3.1 Lagrange插值法 设函数在区间上有定义，且在上个不同点上得函数值，若存在一个至少n次的插值多项式 其中为实数。先构造函数，它们的次数不超过n，且满足 然后以对应点处的函数值为系数作线性组合，即得所要求的多项式。由多项式有n个根，故它必有如下形式 这些函数称为Lagrange插值基函数，而是至多n次多项式，且满足称其为n次Lagrange插值多项式。 3.2 Newton插值法 设有函数为一系列互不相等的点，称为关于点的一阶差商。一般的称为关于的k阶差商。 表1 差商表 一阶差商 二阶差商 三阶差商 其中 显然，是满足插值条件的至多n次多项式。可得。 3.3 Hermite插值法 设，已知互异点，，…，及所对应的函数值为，，…，，导数值为，，…，，则满足条件 的次Hermite插值多项式为 其中 称为Hermite插值基函数，是Lagrange插值基函数，若，插值误差为 ， 4 插值法的应用 拉格朗日插值法和牛顿插值法是两种常用的简便的插值法。但牛顿插值法则更为简便，与拉格朗日插值多项式相比较，它不仅克服了“增加一个节点时整个计算工作必须重新开始”的缺点，而且可以节省乘、除法运算次数。同时，在牛顿插值多项式中用到的差分与差商等概念，又与数值计算的其他方面有着密切的关系。 例1 已知函数表如下 0.1 0.2 0.3 0.4 0.5 0.6 0.09983 0.19867 0.29552 0.38942 0.47943 0.56464 计算的值。 解：方法一 利用拉格朗日插值法计算。 (1) 因为，取插值点 使用的拉格朗日插值法（线性插值），得 因 所以 故有： 可见，用作为的近似值，可保证有三位有效数字。 （2）取插值点 使用的拉格朗日插值法（线性插值），得 因 所以 故有： 可见，用作为的近似值，可保证有四位有效数字。 方法二