五种插值法的对比研究 开题报告

五种插值法的对比研究

1. 选题依据

1.1 选题背景

插值法是一种古老的数学方法，插值法历史悠久。据考证，在公元六世纪时， 我国刘焯(zhuo) 已经把等距二次插值法应用于天文计算。十七世纪时，Newton 和 Gregory(格雷格里) 建立了等距节点上的一般插值公式，十八世纪时，Lagrange(拉格朗日) 给出了更一般的非等距节点插值公式。 而它的基本理论是在微积分产生以后逐渐完善的，它的实际应用也日益增多，特别是在计算机工程中。许多库函数的计算实际上归结于对逼近函数的计算。

1.2 研究的目的和意义

插值法是数值分析中最基本的方法之一。 在实际问题中碰到的函数是各种各样的，有的甚至给不出表达式，只提供了一些离散数据，例如，在查对数表时， 要查的数据在表中找不到，就先找出它相邻的数，再从旁边找出它的修正值， 按一定关系把相邻的数加以修正，从而找出要找的数，这种修正关系实际上就是一种插值。 在实际应用中选用不同类型的插值函数，逼近的效果也不同。在数值计算方法中，我们学习过五种基本的插值方法，即Lagrange 插值、Newton 插值、分段线性插值、分段三次Hermite 插值、样条插值函数。所以通过从这五种插值法的基本思想、特征、性质和具体实例入手，探讨五种插值法的优缺点和适用范围，让学习者能够迅速而准确的解决实际问题，掌握插值法的应用。

2. 研究的方法

从具体实例入手并结合Matlab 在科学计算中的优势，通过实验对它们的精度和效率进行比较分析。

3. 论文结构

3.1 论文的总体结构

第一部分 导言

主要介绍选题的背景、目的及意义、研究现状、文献综述等。

第二部分 五种插值法的基本思想、性质及特点

在数值计算方法中，插值法是计算方法的基础，数值微分、数值积分和微分方程数值解都建立在此基础上。

插值问题的提法是：已知f(x)(可能未知或非常复杂函数)在彼此不同的n+1 个实点0x ,1x ,…n x 处的函数值是f(0x )，f(1x )，…,f(n x )，这时我们简单的说f(x)有n+1 个

离散数据对0n i i )}y ,{(x i .要估算f(x)在其它点x 处的函数值，最常见的一种办法就是插

值，即寻找一个相对简单的函数y(x)，使其满足下列插值条件：y(i x )=f(i x )，i=0,1,…,n.，并以y(x)作为f(x)的近似值.其中y(x)称为插值函数，f(x)称为被插函数。

多项式插值是最常见的一种函数插值.在一般插值问题中，由插值条件可以唯一确定一个次数不超过n 的插值多项式满足上述条件.从几何上看可以理解为：已知平面上n+1 个不同点，要寻找一条次数不超过n 的多项式曲线通过这些点.插值多项式一般有两种常见的表达形式，一个是拉格朗日（Lagrange ）插值多项式，另一个是牛顿（Newton ）插值多项式. 且Lagrange 插值公式恒等于Newton 插值公式.

分段线性插值与样条插值可以避免高次插值可能出现的大幅度波动现象，在实际应用中通常采用分段低次插值来提高近似程度，比如可用分段线性插值或分段三次埃尔米特插值来逼近已知函数，但它们的总体光滑性较差.为了克服这一缺点，一种全局化的分段插值方法———三次样条插值成为比较理想的工具.

（1）拉格朗日插值

Lagrange 插值是n 次多项式插值，其成功地利用构造插值基函数的方法解决了求n 次多项式插值函数问题。对Lagrange n 次插值多项式，首先构造n+1个插值点0x 1x ,....,n x 上的n 次插值基函数)(x l i ))...()()...(()

)...()()...(110110n i i i i i i n i i x x x x x x x x x x x x x x x x --------=+-+-（，)...,2,1,0(n i =

有了这n+1个n 次插值基函数，n 次Lagrange 插值多项式就容易写出来了，具体表达式为)

()()(0x l x f x Ln i n i i ∑==。

表1 插值数值表

Lagrange 插值的方法是：对给定的n 个插值节点,0x 1x ,....,n x 及对应的函数值n y y y y ,......,,,210，利用n 次Lagrange 插值多项式，则对插值区间任意的x 的函数值y 可以通过下式Ln （x ）来求解。

表（1）中的n 次Lagrange 插值多项式Ln （x ）的数学公式为：)

()()(0x l x f x Ln i n i i ∑==。

其中，)(x l i （i=0,1,2,3...,n ）是插值基函数，且∏

=--=n j j i j i x x x x x l 0)(。

Lagrange 插值多项式的余项为R(x)=)()()!1(1)()()1(x f n x L x f n n ωξ++=

-,其中))()(()(10n x x x x x x x ---=ω。

（2）牛顿插值

Newton 插值也是n 次多项式插值，它提出另一种构造插值多项式的方法，与Lagrange 插值相比，具有承袭性和易于变动节点的特点。

Newton 插值的方法：由表（1）构造的牛顿插值多项式为

],...,,[))...((...],,[))((],[)()()(1010210101000n n x x x f x x x x x x x f x x x x x x f x x x f x N ---++--+-+=用它插值时，首先要计算各阶差商，而各阶差商的计算可归结为一阶差商的逐次计算，一般的

111022010)

,...,,(),,...,,(),...,,(-----=k k k k k n x x x x x f x x x x f x x x f

其余项为：)

,...,,()()()(10n x x x f x N x f x Rn =-=。 （3）分段线性插值

分段线性插值函数,记为y(x),y(x)具有下列性质:

①y(x) 可以分段表示,在每个小区间],[1i i x x -上,它是线性函数`)(x y i ;

②)(x y i i i f x f ==)(,(i=0,1,2,3...,n ）.

③ 在整个区间[a,b]上,y(x) 连续.

作分段线性插值的目的在于克服Lagrange 插值方法可能发生的不收敛性缺点.所谓分段线性插值就是利用每两个相邻插值基点作线性插值，即可得如下分段线性插值函数：11)()()(+++=i i i i f x l f x l x y ，],[1+∈i i x x x ，i=0,1,...n.其中11)(++--=

i i i i x x x x x l ，

i i i

i x x x x x l --=++11)(.

特点：插值函数序列具有一致收敛性，克服了高次Lagrange 插值方法的缺点，故可通过增加插值基点的方法提高其插值精度. 但存在基点处不光滑、插值精度低的缺点.从几何上看所谓分段线性插值就是通过插值基点用折线段连接起来逼近原曲线，这也是计算机绘制图形的基本原理.