拉格拉日插值 | MATLAB 源码

拉格朗日插值法MATLAB实现

最新推荐文章于 2024-05-07 14:34:43 发布

原创最新推荐文章于 2024-05-07 14:34:43 发布 · 534 阅读

3 ·

CC 4.0 BY-SA版权

抓住一只白嫖怪(｀・ω・´)

文章标签：

#matlab

MATLAB 专栏收录该内容

20 篇文章

订阅专栏

本文提供了一个使用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