S14 插值及拟合.pptVIP

* 插值就是根据被插函数给出的函数表计算所要点的函数值。 x x0 x1 x2 …… xn y y0 Y1 y2 …… yn 插值与拟合 求任意x(x≠xi)时y的值 满足 则称P(x)为f(x)的n次插值多项式。这种插值法通常称为代数插值法。 代数插值法 * 是x的n次函数，且 由以上条件知, 都是n次函数 的零点,故可设 拉格朗日插值多项式 * 代入得 称 为关于基点 的n次插值基函数(i=0,1,…,n) * 当n=1时可得两点插值公式： 当n=2时可得三点插值公式（抛物线）： * y: 求解函数值 p： * * SUBROUTINE ENLGR(X,Y,N,T,Z) DIMENSION X(N),Y(N) DOUBLE PRECISION X,Y,T,Z,P Z=0.0 IF (N.EQ.2) THEN !！两点插值公式 Z=(Y(1)*(T-X(2))-Y(2)*(T-X(1)))/(X(1)-X(2)) END IF I=1 DO J=1,N IF (X(J).LT.T)I=I+1 !确定插值点T的位置 ENDDO K=I-4 !确定k 的值,左端点 IF (K.LT.1) K=1