样条插值函数及应用

样条插值函数及应用

摘要

样条函数具有广泛的应用，是现代函数论的一个十分活跃的分支，是计算方法的主要基础和工具之一，由于生产和科学技术向前发展的推动以及电子计算机广泛应用的需要，人们便更多地应用这个工具，也更深刻的认识了它的本质。

在实际问题中所遇到许多函数往往很复杂，有些甚至是很难找到解析表达式的。根据函数已有的数据来计算函数在一些新的点处的函数值，就是插值法所需要解决的问题。

插值法是数值逼近的重要方法之一，它是根据给定的自变量值和函数值,求取未知函数的近似值。早在一千多年前,我国科学家就在研究历法时就用到了线性插值和二次插值。而在实际问题中,有许多插值函数的曲线要求具有较高的光滑性,在整个曲线中,曲线不但不能有拐点,而且曲率也不能有突变。因此，对于插值函数必须二次连续可微且不变号 ,这就需要用到三次样条插值。

关键词三次样条函数；插值法

目录

引言 (1)

第一章三次样条插值 (2)

1.1 样条插值函数简介 (2)

1.2 三次样条函数应用 (3)

第二章AMCM91A 估计水塔水流量 (5)

2.1 理论分析及计算 (6)

2.2运用MATLAB软件计算 (9)

参考文献 (14)

引言

样条函数具有广泛的应用，是现代函数论的一个十分活跃的分支，是计算方法的主要基础和工具之一，由于生产和科学技术向前发展的推动以及电子计算机广泛应用的需要，人们便更多地应用这个工具，也更深刻的认识了它的本质。上世纪四十年代，在研究数据处理的问题中引出了样条函数，例如，在1946年Schoenberg将样条引入数学,即所谓的样条函数，直到五十年代，还多应用于统计数据的处理方面，从六十年代起，在航空、造船、汽车等行业中，开始大量采用样条函数。

在我国，从六十年代末开始，从船体数学放样到飞机外形设计，逐渐出现了一个使用样，逐渐出现了一个使用样条函数的热潮，并推广到数据处理的许多问题中。

在实际生活中有许多计算问题对插值函数的光滑性有较高的要求，例如飞机机翼外形、发动机进、排气口都要求有连续的二阶导数，用三次样条绘制的曲线不仅有很好的光滑度，而且当节点逐渐加密时其函数值整体上能很好地逼近被插函数，相应的导数值也收敛于被插函数的导数值，不会发生“龙格现象”。

现在国内外学者对这方面的研究也越来越重视，根据我们的需要来解决不同的问题，而且函数的形式也在不断地改进，长期以来很多学者致力于样条插值的研究，对三次样条的研究已相当成熟。

第一章 三次样条插值

1.1 样条插值函数简介

在实际问题中所遇到许多函数)(x f 往往很复杂，有些甚至是很难找到解析表达

式的。有时通过实验或者数值计算所得到的也只是一些离散的点()n i x i ...2,10，

=上的函数值，即n i x f i i ...2,1,0),(y ==。

根据函数)(x f 已有的数据来计算函数)(x f 在一些新的点x 处的函数值，就是插值法所需要解决的问题。

插值法的基本思想就是，首先根据已有的函数值来构造一个简单的函数)(y x 作为)(x f 的近似表达式，然后用)(y x 来计算新的点上的函数值作为

)(x f 的近似值。通常可以选多项式函数作为近似函数)(y x ，因为多项式具有各阶导数，

求值也比见方便。常用的有Lagrange 插值、Newton 插值、Hermite 插值和样条插值。线性插值在分段点上仅连续而不可导，三次埃尔米特插值有连续的一阶导数，这样的光滑程度常不能满足物理问题的需要，样条函数可以同时解决这两个问题,使插值函数既是低阶分段函数,又是光滑的函数，并且只需在区间端点提供某些导数信息。

三次样条函数定义：

设在区间[]b ,a 上取1n +个节点b x x x x a n =<<<=...210，函数)(y x f =在各个节点处的函数值为,,...,1,0),(y i n i x f i ==若)(S x 满足：

(1) ;,...1,0,y )(S n i x i i ==

(2) 在区间[]b ,a 上，)(S x 具有连续的二阶导数；

(3) 在区间[]）（1,...,1,0x ,x 1-=-n i i i 上,)(S x 是x 三次的多项式；

则称)(S x 是函数)(y x f =在区间[]b ,a 上的三次样条插值函数。

由以上定义可以看出，虽然每个子区间上的多项式可以各不相同，但在相邻子区间的连接处却是光滑的。因此，样条插值也称为分段光滑插值。

从定义知要求出)(S x ，在每一个小区间[]）（1,...,1,0x ,x 1-=-n i i i 上确定4个待定系

数，共有n 个小区间，故应有4n 个参数。

根据)(S x 在[]b ,a 上二阶导数连续，在节点()1-...2,1n i x i =出满足连续性条件

)0()0(S +=-i i x S x )0()0(S ''+=-i i x S x )0()0(S ''''+=-i i x S x

共有3n-3个条件，再加上)(S x 满足插值条件;,...1,0,y )(S n i x i i ==共有4n-2个条件，因此还需要2个条件才能确定)(S x 。

通常可在区间端点上各加一个条件（称为边界条件），可根据实际问题的要求给定，通常有以下三种：

（1）已知端点的一阶导数值，即

'''

00')(S ,)(S n n f x f x ==

（2）俩端点的二阶导数已知，即

'''''

'00'')(S ,)(S n n f x f x ==

其特殊情况，0)(S )(S ''0''==n x x 称为自然边界条件。

（3）当)(x f 是以0x x n -为周期的函数时，则要求)(S x 也是周期函数。这时边界条件应满足

)0()0(S +=-i i x S x )0()0(S ''+=-i i x S x )0()0(S ''''+=-i i x S x

而此时n y y =0。这样确定的样条函数)(S x ，称为周期函数。

1.2 三次样条函数应用

作函数)2sin()73(42x e x x y x ⨯+-=-在[]1,0取间隔为0.1的点图，用插值进行实验。 使用MATLAB 软件 程序代码如下： %产生原始数据 x=0:0.1:1;

y=(x.^2-3*x+7).*exp(-4*x).*sin(2*x);