hermite.m

各种插值拟合的子程序

function y=hermite(X,Y,y1,x)   %Hermit插值
n=length(X);
m=length(x);
for i=1:m
    y0=0;
    for j=1:n
        h=1;
        a=0;
        for k=1:n
            if k~=j
                h=h*((x(i)-X(k))/(X(j)-X(k)))^2;
                a=1/(X(j)-X(k))+a;
            end
        end
        y0=y0+h*(((X(j)-X(i))*(2*a*Y(j)-y1(j))+Y(j)));
    end
    y(i)=y0;
end