newton_int.m

来自「为配套王能超版本数值分析简明教程编写,方便学生教师演示!」· M 代码 · 共 28 行

function yi=Newton_int(x,y,xi)
%newton插值公式
n=length(x);m=length(y);
if n~=m
error('The length of X and Y must be equal!');
return;
end
%计算均差表Y
Y=zeros(n);Y(:,1)=y';
for k=1:n-1;
    for i=1:n-k
        %%输入的插值节点必须互异
        if abs(x(i+k)-x(i))<eps
                error('the DATA is error')
                return;
            end
            Y(i,k+1)=(Y(i+1,k)-Y(i,k))/(x(i+k)-x(i));
        end
    end
    %计算牛顿插值公式
    yi=0;
    for i=1:n
        z=1;
        for k=1:i-1
            z=z*(xi-x(k));
        end
        yi=yi+Y(1,i)*z;
    end