插值法的matlab实现（拉格朗日、牛顿、样条）

拉格朗日

function yy = lagronge()
x=[25,40,50,60];
y=[95,75,63,54];
xx=70;
a = length(x);
f = 0;
for i = 1:a
    t = 1;
    for j = 1:a
        if j~=i
            t=t*(xx-x(j))/(x(i)-x(j));
        end
    end
    f = f+t*y(i);
end
fprintf('%.1f',f);
end

牛顿插值法（先计算出差商表再进行计算）

n = 10;
f = zeros(n,n);
xi = [1994,1995,1996,1997,1998,1999,2000,2001,2002,2003];
x = 2010;
yi = [67.052,68.008,69.803,72.024,73.400,72.063,74.669,74.487,74.065,76.777];
for k = 1:n
    f(k) = yi(k);
end
for i=2:n       
    for k=i:n
        f(k,i)=(f(k,i-1)-f(k-1,i-1))/(xi(k)-xi(k+1-i));