lagrange.m

这是我从一本基于MATLAB实现的数值分析算法的书中自带的光盘程序

%1-1 拉格朗日插值
%x,y:数组形势的数据表
%xi:待计算y值的横坐标数组
%yi:用Lagrange插值算出的y值数组
function fi=Lagrange(x,f,xi)
fi=zeros(size(xi));
np1=length(f);
for i=1:np1
    z=ones(size(xi));
    for j=1:np1
        if i~=j,
           z=z.*(xi-x(j))/(x(i)-x(j));
       end
    end
    fi=fi+z*f(i);
 end