小熊的数值分析.docVIP

插值法的实际应用 ——信科132 姓名 熊玉玮 摘要:插值法是计算数学中的一种重要的方法，而且计算问题可以说是现代社会各个领域普遍存在的共同问题，无论哪一行哪一业都有许多数据需要处理，插值法正在科学技术中发挥越来越大的作用. 关键字：插值 分段线性插值 三次样条插值 温度预测 插值法是函数逼近的一种重要方法，是数值计算的基本课题在许多实际问题及科学研究中，因素之间往往存在着函数关系，然而，这种关系经常很难有明显的解析表达，通常只是由观察与测试得到一些离散数值f（x）a=x0x1…xn=b上的函数值f0,f1,…,fn，记hk=xk+1-xk,h=max hk，求一折线函数Ih(x)满足 （1） （2）Ih(xk)=fk(k=0,1,…，n) （3）Ih(x)在每个小区间 [xk,xk+1]上是线性函数. 则称Ih(x)为分段线性插值函数. 由定义可知Ih(x)在每个小区间 [xk,xk+1]上可表示为 在整个区间[a,b]上定义分段线性插值函数为 其中 分段线性插值函数的误差估计可利用插值余项（2.23）得到 或 其中，由此可得到，在 [a,b]上一致成立.故Ih(x)在 [a,b]上一致收敛到f(x). 二、三次样条插值 定义：设在区间[a,b]上给定一组节点a=x0x1…xn=b上的函数值y0,y1,…，yn，其函数S(x)满足： （1）S(x)在每个子区间[xk-1,xk](k=1,2,…,n)上都是次数不超过3的多项式； （2）S(xk)=yk,k=0,1,…，n; （3）S(x)在[a,b]上的二阶导数连续。则称S(x)为三次样条插值函数。 从定义知要求出S(x)，在每个区间[xk-1,xk]上要确定4个待定系数，共有n个小区间，故应确定4n个常数，根据S(x)在[a,b]上二阶导数连续，在节点xk(k=1,2,…,n-1)处应满足连续性条件 共有3n-3个条件，再加上S(x)满足插值条件S(xk)=yk，共有4n-2个条件，因此还需要2个条件才能确定S(x)，通常可在区间[a,b]的端点a=x0,b=xn上各加一个条件（称为边界条件），可根据实际问题的要求给定。常见的有以下3种： （1）已知两端的一阶导数值，即 （2）已知两端的二阶导数值，即 其特殊情况为 称为自然边界条件。 （3）当f(x)是以x-x0为周期的周期函数时，则要求S(x)也是周期函数，这时边界条件应满足 而此时有y0=yn，这样确定的样条函数S(x)称为周期样条函数。 四、温度预测上的实际应用 某班级组织一次实践活动，为了是同学了解温度的变化，在12h内，每隔1h测量一次温度，温度依次为：5，8，9，15，25，29，31，30，22，25，27，24.（单位：oC），试分别用分段线性插值、三次样条插值方法估计在3.2h,6.5h,7.1h,11.7h的温度值，每隔1/10 h估计一次温度值并画出其图形. ①用分段线性插值、三次样条插值方法估计在3.2h,6.5h,7.1h,11.7h的温度值 应用分段差值多项式 进行编程。 程序如下： hours=1:12; temps=[5,8,9,15,25,29,31,30,22,25,27,24]; t=interp1(hours,temps,[3.2,6.5,7.1,11.7]) T=interp1(hours,temps,[3.2,6.5,7.1,11.7],spline) 结果为 t=10.2000 30.0000 30.9000 24.9000 T=9.6734 30.0427 31.1755 25.3820 ②每隔1/10 h估计一次温度值并画出其图形，如图一（程序为附件1） 图一 由上可以很明显的看出，由于计算机的加入，使其非常复杂的运算变得非常简单，因而计算机对计算机对插值法的作用真的很大. 参考文献 [1] 徐萃薇编.《计算方法引论》（第二版）[M].高等教育出版社,1985:13-147页 [2] 张韵华、奚梅成编.《数值计算方法与算法》（第二版）[M].科学出版社,2006:10-40页 [3] 朱长青编著.《数值计算方法及其应用》（第二版）[M].科学教育出版社,2000:23-58页 附件 附件1: hours=1:12; temps=[5,8,9,15,25,29,31,30,22,25,27,24]; h=1:0.1:12; t=interp1(hours,temps,h,spline); plot(hours,temps,+,h,t,hours,temps,r:) xlabel(时间),ylabel(温度) 1