matlab拉格朗日插值代码

以下是MATLAB中的拉格朗日插值代码：

% x 和 y 分别为已知数据点的横纵坐标
% xi 为待求点的横坐标
function yi = lagrange_interp(x, y, xi)

n = length(x);
yi = 0;

for i = 1:n
    % 计算拉格朗日基函数
    L = 1;
    for j = 1:n
        if j ~= i
            L = L * (xi - x(j)) / (x(i) - x(j));
        end
    end
    
    % 利用拉格朗日基函数计算插值结果
    yi = yi + y(i) * L;
end

end

使用方法：

将上述代码保存为 `lagrange_interp.m` 文件，放在MATLAB工作路径下。

在命令窗口中输入以下代码即可进行插值计算：

% 定义已知点的横纵坐标
x = [1, 2, 3];
y = [4, 5, 6];

% 定义待求点的横坐标
xi = 2.5;

% 进行拉格朗日插值计算
yi = lagrange_interp(x, y, xi);