分段线性插值

摘要 用函数来表示变量间的数量关系广泛应用于各学科领域，但是在实际问题中，往往是通过实验、观测以及计算等方法，得到的是函数在一些点上的函数值。如何通过这些离散数据找到函数的一个满足精度要求且便于使用的近似表达式，是经常遇到的问题。对于这类问题我们解决的方法为插值法，而最常用也最简单的插值方法就是多项式插值。当然用插值法得到的近似表达式必须满足插值条件即假设给定了n+1个点的自变量的值以及函数值，近似函数必须要过这n+1个点。多项式插值，从几何角度看，就是寻求n次代数曲线y=Pn（x）通过n+1个点作为f（x）的近似。但是随着插值节点个数的增加，高次插值多项式的近似效果并不理想。根据大量实验得出，在进行高次多项式插值时，会出现龙格现象。龙格（Runge）现象即当n趋于无穷大时，x在某一邻域内,f(x)收敛，而在这个区域外f(x)发散。因此，为了解决这样的一个问题，我们可以通过缩小插值区间的办法达到减小误差的目的，所以本实验将针对低次分段插值多项式来做具体的讨论和学习。关键词：龙格现象 分段差值1、实验目的1) 通过对分段线性插值算法程序的编写，提高自己编写程序的能力2) 体会分段线性插值是如何消除龙格现象的。3) 用实验报告的形式展现，提高自己在写论文方面的能力2、算法理论设y=f(x)在节点a=x0x1xn=b处的函数值为yi=f(xi),i=0,1, ,n。为了提高近似程度，可以考虑用分段线性插值来逼近原函数，这时的插值函数为分段函数：Sx=s1x,xx0,x1 s2x,xx1,x2 snx,xxn-1,xn在区间xi-1,xi上的线性函数为Six=yi-1x-xixi-1-xi+yix-xi-1xi-xi-1,i=1,2,n误差为：Rix=fx-Six=fi2!x-xi-1x-xi,xi-1ixi易见，S(x)是平面上以点xi,yi(i=0,1,2,n)为节点的折线，有如下的特点：1. S(x)在xi-1,xi上为次数不超过一次的多项式；2. S(x)Ca,b；3. S(x)C1xi-1,xi；如果f(x)C2a,b,由线性插值的误差公式得到 Ri(x)=fi2(x-xi-1)(x-xi)12xxi-1,ximax fxxxi-1,ximax x-xi-1x-xi 令Mi=maxxi-1xxifx,则有 Ri(x)=fi2(x-xi-1)(x-xi)12xxi-1,ximax fxxxi-1,ximax x-xi-1x-xi=Mi214xi-xi-12=Mi8hi2,hi=xi-xi-1关于整体误差：Rx=fx-S(x)可以按如下方式考虑，若记M=1inmax Mi ,h=1inmax hi ,则对任一xa,b都有Rx=fx-S(x)1inmax Ri181inmax Mi1inmax hi2M8h2于是，当h0,时R(x)0，说明分段线性插值S(x)收敛于fx。3、数值算例已知点坐标如下表所示：xi0.50.70.91.11.3yi0.45790.6440.7830.8910.964用分段线性插值法，求解当x为0.8时，对应y的值解：具体程序如下所示：#include stdafx.hfloat Fdline(float x,float y,float x1,int len)int i=0;float s=0;for(i=0;i=xi & x1xi+1)break;s=(x1-xi)/(xi-1-xi)*yi-1+(x1-xi-1)/(xi-xi-1)*yi;return s;float Fdline(float x,float y,float x1,int len);void main()float x=0.5,0.7,0.9,1.1,1.3;float y=0.479,0.644,0.783,0.891,0.964;int len=sizeof(x)/sizeof(x0);float x1=0;float s=0;printf(请输入要求解的x1的值:n);scanf(%f,&x1);s=Fdline(x,y,x1,len);printf(经过分段三次Hermite插值的结果为:n);printf(%fn,s);运行结果：5、对结果进行分析根据分段线性插值的原理，可以看出分段线性插值虽然有很好的收敛性质，但却不是光滑的，所以线性插