拉格拉日插值 | MATLAB 源码

---

代码部分如下：

x_train=[0.32 0.34 0.36];
y_train=[0.314567 0.333487 0.352274];
x_test=0.3367;
y_test=myfun(x_test,x_train,y_train);
y_real=sin(0.3367);
disp('拉格朗日差值计算得到的结果为')
disp(y_test)
disp('真实值为')
disp(y_real)
function y_test=myfun(x_test,x_train,y_train)
n = length(x_train);
% 有 n 个数据点，就有 n 个基函数
l = ones(1,n);
for i = 1:n
    for j = 1:n
        if i==j
            continue
        else
            l(i)=l(i)*(x_test-x_train(j))/(x_train(i)-x_train(j));
        end
    end
end
y_test = sum(y_train.*l);
end

---

上述代码可以在大部分的MATLAB版本中运行。如果需要完整代码注释.mlx文件，请关注微信公众号 我不是wc，回复 018 获取。

有问题可以在评论区或者私信问我。

---

2021年9月22日18:56:36