lagrange.m

来自「该算法实现了拉格朗日插值并比较了龙格现象」· M 代码 · 共 15 行

function fi=lagrange(x,f,xi) %lagrange插值算法
%x自变量
%f函数值
%xi插值点
%fi返回插值
 fi=zeros(size(xi));
 npl=length(f);
 for i=1:npl
     z=ones(size(xi));
     for j=1:npl
         if i~=j,z=z.*(xi-x(j))/(x(i)-x(j));end
     end
     fi=fi+z*f(i);
 end 
 return