清华大学附属教育集团常见插值法.docVIP

常见插值法 【摘 要】插值方法在数值分析中起着非常重要的作用。在此介绍一些常见的插值方法及 其应用范例。 【关键字】数值分析；插值方法；应用； 插值法定义 ....... ....... 插值法又称“内插法”，是利用函数f (x)在某区间中插入若干点的函数值，作出适当的特定函数，在这些 表(1) 插值点 点上取已知值，在区间的其他点上用这特定函数的值作为函数f (x)的近似值，这种方法称为插值法。如果这特定函数是多项式，就称它为插值多项式。 2.常见的插值法及其构造 ? Lagrange插值法 （a).公式推导： 表（1）的Lagrange插值的插值多项式 ,(j=0,1,2....n)。 其中插值基函数是 ,(i,j=0,1 2...n) 。 其插值余项为 其中， (b).matlab实现方法： Matlab没有直接求解的相关函数，现编译如下： function yi = Lagarange_chazhi(x,y,xi) % 求拉格朗日插值,并返回一个输入为xi时的函数值 % x 为插值点向量，至少有三项 % y 为插值点值的向量，项数与x相同 m = length(x); %求插值个数 m1 = length(y); if m=2 error(项数不足!); end if m~=m1 error(!!!y的项数应与x相同!!!); end %对参数的判断 lag_hanshu = 0; syms X; for (l = 1:m) %构造插值基函数 la = y(l); for a = (1:l-1) la = la*(X-x(a))/(x(l)-x(a)); end for a = (l+1:m) la = la*(X-x(a))/(x(l)-x(a)); end format long lag_hanshu = lag_hanshu+la; %求解出插值函数 end yi = subs( lag_hanshu,X,xi); %返回插值函数输入为xi时的值 End (c).方法缺陷：当插值点个数时，将产生龙格现象： 经典例子，对进行拉格朗日插值图（1）中从左到右，从上到下，n分别为 图(1) Lagarange插值法的龙格现象 4,5...11，可以看出，当后，它的\插值函数在两个端点处发生剧烈的波动，造成较大的误差。所以拉格朗日插值法一般不适用于高次插值。 ? 分段线性插值法 (a).插值原理： 为避免高次产生的龙格现象，采用分段插值，相邻两节点间的函数为一次线性函数, 图形为线段，在[a,b]间为折线，如图（2）。 图(2) 分段线性插值 (b).matlab实现方法： Matlab中没有现成的函数，现编译如下： function Fenduan_liner(x,y) %对x,y进行分段线性插值 %用虚线画出插值后的函数 m = length(x); %求插值个数 m1 = length(y); if m~=m1 error(!!!y的项数应与x相同!!!); end %对参数的判断 hold on; for ii = 1:m-1 plot([x(ii) x(ii+1)],[y(ii) y(ii+1)]); %画出线段 end plot(x,y,o); %画出插值点 end (c).方法缺陷：在节点处曲线不平滑。 ? 保形插值（Hermite插值） .插值原理： 针对分段插值法不光滑的问题，Hermite插值引入插值点的导数。 表(2) Hermite插值法 由于若且已知函数表及导数表，则存在唯一不超过次多项式满足插值条件 则，通过求解方程，可得出插值函数 其插值余项为 且与x有关。 (b)