lagran_.m

各种插值拟合的子程序

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