lagran_.m

来自「matlab常用插值算法源代码,内有本人编写的和下载的程序」· M 代码 · 共 13 行

function fi=Lagran_(x,y,xi) %根据给出的数据向量x和y，构造拉格朗日插值多项式，%并计算在xi处的值。
fi=zeros(size(xi));  % 全零阵
np1=length(y);
for i=1:np1
    z=ones(size(xi)); %全1阵
    for j=1:np1
        if i~=j
            z=z.*(xi-x(j))/(x(i)-x(j)); %插值基函数
        end
    end
    fi=fi+z*y(i);
end
return