MaTLAB 实现 Lagrange插值多项式(简单好理解)

实验题目：Lagrange插值多项式或Newton插值多项式

实验用例： 已知函数y=f(x)的一张表：
x 0 10 20 30 40 50 60 70 80 90 100 110 120
y 5 1 7.5 3 4.5 8.8 15.5 6.5 -5 -10 -2 4.5 7
试验要求（从下面题目中任选一个题目）：
⒈ 利用Lagrange插值多项式求被插值函数f(x)在点x=65处的近似值。建议：画出Lagrange插值多项式的曲线。
%利用Lagrange插值多项式求被插值函数f(x)在点x=65处的近似值，并画出Lagrange插值多项式的曲线。

%利用Lagrange插值多项式求被插值函数f(x)在点x=65处的近似值，并画出Lagrange插值多项式的曲线。
%xx表示x的值 p表示最后的Lagrange插值多项式 Y表示当x等于65的时候f(x)的值
function [p,Y] = lagrange(xx) % MATLAB函数表示法
x=[0,10,20,30,40,50,60,70,80,90,100,110,120];
y=[5 , 1 , 7.5 , 3 , 4.5 , 8.8 , 15.5 , 6.5 , -5 , -10 , -2 , 4.5 , 7];
n=length(x);
syms X;%一个变量x
for i=1:n
    t=1;
    for j=1:n
        if j~=i
            t=t*((X-x(j))/(x(i)-x(j)));
        end
    end
    L(i)=t; %拉格朗日插值基函数
end


u=sum(L.*y); %yk分别与lk相乘，求和  u是最后的表达式
p=simplify(u); %
fprintf('拉格朗日多项式为:\nP(x)=%s\n',p);
Y=double(subs(p,X,xx));  %求出表达式p 当X=65 时 Y的值

xxx = 0:0.1:120; %以0为起始  0.1为递增大小  120为终止数
yyy = subs(p,X,xxx); % 求出所有xxx节点的y值

plot(xxx,yyy); %画图函数 
title('Lagrange插值多项式')
xlabel('X');
ylabel('Y');
xlim([0 120]);

运行注意事项:

由于使用了函数表示方法并带有参数,因此运行需要在命令窗口去调用该函数,在命令窗口敲入命令:

[p,Y] = lagrange(65)
//65就是你想要求被插值函数f(x)在点x=65处的近似值

运行结果截图

总结:

我也只是个新手啦,也是第一次在MATLAB上写实验,有错还望大家多多指教