nalagr.m

来自「Matlab数值解算法实现代码」· M 代码 · 共 18 行

function yy=nalagr(x,y,xx)
% 用途：Lagrange 插值法数值求解
% 格式：yy=nalagr(x,y,xx), x是节点向量，y是节点上的函数值，xx是插值点(可以是多个)
%       yy返回插值
m=length(x);n=length(y);
if m~=n, error('向量x与y的长度必须一致');end
s=0;
for i=1:n
   t=ones(1,length(xx));
   for j=1:n
      if j~=i, 
         t=t.*(xx-x(j))/(x(i)-x(j));
      end
   end
   s=s+t*y(i);
end
yy=s;