样条插值函数及应用..doc

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