拉格朗日插值多项式matlab程序

下面是拉格朗日插值多项式的 Matlab 程序：

function y = lagrange_interpolation(x, y, xi)
% x: 插值节点的 x 坐标
% y: 插值节点的 y 坐标
% xi: 需要插值的 x 坐标

n = length(x);
m = length(xi);
L = ones(n, m);

for i = 1:n
    for j = 1:n
        if i ~= j
            L(i,:) = L(i,:) .* (xi - x(j)) / (x(i) - x(j));
        end
    end
end

y = 0;
for i = 1:n
    y = y + y(i) * L(i,:);
end
end

这个程序有三个输入参数，分别是插值节点的 x 坐标、插值节点的 y 坐标和需要插值的 x 坐标。输出结果是在需要插值的点 xi 处的插值结果。