hmt.m

HERMITE插值函数

%hermite 插值函数  其中n为节点数，x为已知节点数组，f为函数值，f1为一阶导数值，xi为欲插值的数组。
function y=hmt(n,x,f,f1,xi)
y=0.0;
for i=1:n
     h=1.0;
     a=0.0;
     for j=1:n
        if j ~= i
          h=h*((xi-x(j))/(x(i)-x(j)))^2 ;
          a=1/(x(i)-x(j))+a;
        end
     end
     y=y+h*(x(i)-xi)*(2*a*f(i)-f1(i)+f(i));
end