nalagr.m

来自「matlab具有强大的数值计算能力」· M 代码 · 共 17 行

function yy=nalagr(x,y,xx)
%用途： Lagrange插值法数值求解
%格式： yy=nalagr(x,y,xx) x是节点向量，y是节点上的函数插值，xx是插值点
% （可以是多个），yy返回插值
m=length(x);n=length(y);
if m~=n, error('x!=y');end
s=0;
for i=1:n
    t=ones(1,length(xx));
    for j=1:n
        if j~=i,
            t=t.*(xx-x(j))/(x(i)-x(j));
        end
    end
    s=s+t*y(i);
end
yy=s;