资源描述：
《样条插值函数及应用》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、样条插值函数及应用2摘要样条函数具有广泛的应用，是现代函数论的一个十分活跃的分支，是计算方法的主要基础和工具之一，由于生产和科学技术向前发展的推动以及电子计算机广泛应用的需要，人们便更多地应用这个工具，也更深刻的认识了它的本质。在实际问题中所遇到许多函数往往很复杂，有些甚至是很难找到解析表达式的。根据函数已有的数据来计算函数在一些新的点处的函数值，就是插值法所需要解决的问题。插值法是数值逼近的重要方法之一，它是根据给定的自变量值和函数值,求取未知函数的近似值。早在一千多年前,我国科学家就在研究历法时就用到了线性插值和二次插值。而在实际问题中,有许多插值函数的
2、曲线要求具有较高的光滑性,在整个曲线中,曲线不但不能有拐点,而且曲率也不能有突变。因此，对于插值函数必须二次连续可微且不变号,这就需要用到三次样条插值。关键词三次样条函数；插值法2目录引言1第一章 三次样条插值21.1样条插值函数简介21.2三次样条函数应用3第二章 AMCM91A估计水塔水流量52.1理论分析及计算62.2运用MATLAB软件计算9参考文献142引言样条函数具有广泛的应用，是现代函数论的一个十分活跃的分支，是计算方法的主要基础和工具之一，由于生产和科学技术向前发展的推动以及电子计算机广泛应用的需要，人们便更多地应用这个工具，也更深刻的认识了
3、它的本质。上世纪四十年代，在研究数据处理的问题中引出了样条函数，例如，在1946年Schoenberg将样条引入数学,即所谓的样条函数，直到五十年代，还多应用于统计数据的处理方面，从六十年代起，在航空、造船、汽车等行业中，开始大量采用样条函数。在我国，从六十年代末开始，从船体数学放样到飞机外形设计，逐渐出现了一个使用样，逐渐出现了一个使用样条函数的热潮，并推广到数据处理的许多问题中。在实际生活中有许多计算问题对插值函数的光滑性有较高的要求，例如飞机机翼外形、发动机进、排气口都要求有连续的二阶导数，用三次样条绘制的曲线不仅有很好的光滑度，而且当节点逐渐加密时其
4、函数值整体上能很好地逼近被插函数，相应的导数值也收敛于被插函数的导数值，不会发生“龙格现象”。现在国内外学者对这方面的研究也越来越重视，根据我们的需要来解决不同的问题，而且函数的形式也在不断地改进，长期以来很多学者致力于样条插值的研究，对三次样条的研究已相当成熟。14第一章 三次样条插值1.1样条插值函数简介在实际问题中所遇到许多函数往往很复杂，有些甚至是很难找到解析表达式的。有时通过实验或者数值计算所得到的也只是一些离散的点上的函数值，即。根据函数已有的数据来计算函数在一些新的点处的函数值，就是插值法所需要解决的问题。插值法的基本思想就是，首先根据已有的函
5、数值来构造一个简单的函数作为的近似表达式，然后用来计算新的点上的函数值作为的近似值。通常可以选多项式函数作为近似函数，因为多项式具有各阶导数，求值也比见方便。常用的有Lagrange插值、Newton插值、Hermite插值和样条插值。线性插值在分段点上仅连续而不可导，三次埃尔米特插值有连续的一阶导数，这样的光滑程度常不能满足物理问题的需要，样条函数可以同时解决这两个问题,使插值函数既是低阶分段函数,又是光滑的函数，并且只需在区间端点提供某些导数信息。三次样条函数定义：设在区间上取个节点，函数在各个节点处的函数值为若满足：(1)(2)在区间上，具有连续的二阶
6、导数；(3)在区间上,是三次的多项式；则称是函数在区间上的三次样条插值函数。由以上定义可以看出，虽然每个子区间上的多项式可以各不相同，但在相邻子区间的连接处却是光滑的。因此，样条插值也称为分段光滑插值。从定义知要求出，在每一个小区间上确定4个待定系数，共有n个小区间，故应有4n个参数。根据在上二阶导数连续，在节点出满足连续性条件14共有3n-3个条件，再加上满足插值条件共有4n-2个条件，因此还需要2个条件才能确定。通常可在区间端点上各加一个条件（称为边界条件），可根据实际问题的要求给定，通常有以下三种：（1）已知端点的一阶导数值，即（2）俩端点的二阶导数
7、已知，即其特殊情况称为自然边界条件。（3）当是以为周期的函数时，则要求也是周期函数。这时边界条件应满足而此时。这样确定的样条函数，称为周期函数。1.2三次样条函数应用作函数在取间隔为0.1的点图，用插值进行实验。使用MATLAB软件程序代码如下：%产生原始数据x=0:0.1:1;y=(x.^2-3*x+7).*exp(-4*x).*sin(2*x);14%作图subplot(1,2,1);plot(x,y,x,y,'ro')%待求插值点xx=0:0.02:1;yy=interp1(x,y,xx,'spline');%作图subplot(1,2,2)plot(
8、x,y,'ro',xx,yy,'b')运行截图图1.
