lagelangri.m

来自「牛顿迭代法」· M 代码 · 共 27 行

function yi=lagelangri(x, y, xi)
% 拉格朗日插值多项式，其中
%   x --- 向量，全部的插值节点
%   y --- 向量，插值节点处的函数值
%   xi --- 标量，自变量x
%   yi --- xi 处的函数估计值
    
n=length(x); m=length(y);
if n~=m
    error('The lengths 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);