lagrange.m

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

function yi=Lagrange(x,y,xi)
%Lagrange为插值多项式
%x为向量，全部插值节点
%y为向量，插值节点处函数值
%xi为标量，被估计函数自变量
n=length(x);m=length(y);
if n~=m
error('The length of X and Y must be equal!');
return;
end
p=zeros(1,n);
for k=1:n
    t=ones(1,n);
    for j=1:n;
        if j~=k
            %输入的插值节点必须互异
            if abs(x(k)-x(j))<eps
                error('the DATA is error')
                return;
            end
            t(j)=(xi-x(j))/(x(k)-x(j));
        end
    end
    p(k)=prod(t);    
end
yi=sum(y.*p);