lagrange.m

来自「数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序」· M 代码 · 共 14 行

function yy=lagrange(x,y,xi)   %建立lagrange插值; 
m=length(x);n=length(y);       %m为向量x的长度,n为向量y的长度;
if m~=n,error('向量x与y的长度必须一致');end   %若m与n的值不相等，则输出提示语句;
s=0;
for i=1:n
    z=ones(1,length(xi));
    for j=1:n
        if j~=i
            z=z.*(xi-x(j))/(x(i)-x(j));      %建立lagrange插值的基本格式;
        end
    end
    s=s+z*y(i);
end
yy=s;                          %输出yy的值;