拉格朗日插值法matlab程序

好的，以下是拉格朗日插值法的 MATLAB 程序：

function y = lagrange(x, xdata, ydata)
% x: 插值点
% xdata: 已知点的横坐标
% ydata: 已知点的纵坐标

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

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