代码部分如下:
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