tc_hermite.m

％ Nevill等插值算法地实现 不会花费大家的时间的 不希望其他人没有帐号自由下载此源码

function H=tc_hermite(x,y,dy,xx)
% Hermite插值     07年12月6日
% x表示输入函数变量 y表示函数变量对应的值  dy为对应的导数
%x中不应该含有相同的数值
m=length(x);
if(isempty(dy)==1)
    dy=gradient(y,x);
end  %if         如果此处导数缺省采用均差代替 

for j=1:m
    a=1;b=0;
for i=1:m
    if i~=j
        a=a*((xx-x(i))/(x(j)-x(i)))^2;  %  表示L(x)^2的相乘
        b=b+1/(x(j)-x(i));                                                        end
end
h(j)=a*(1-2*(xx-x(j))*b);   % 表示函数f（x）的基函数
g(j)=(xx-x(j))*a;           % 表示函数导数的基函数
end
H=y*h'+dy*g';
%检验x=[1 4 9];y=[1 2 3];dy=[1/2 1/4 1/6];xx=5;H=tc_hermite(x,y,dy,xx)