lagrangen.m

MatLab的拉格朗日插值法

%1.此程序为n个节点及函数值的n-1次Lagrange插值多项式求x处的近似值程序
%2.lagrange insert
%3.x0为n维插值节点向量，y0为n维被插函数值向量
%4.x为插值点（也可以是m维向量）
%5.y为插值点x处的函数值（近似值，也可以是与x同维的向量）
%6.本程序于2001.7.3.编写? 
function y=lagrangen(x0,y0,x)
n=length(x0);
%m=length(x);
%for i=1:m
   z=x;%(i);
   s=0;
   for k=1:n
      L=1;
      for j=1:n
         if j~=k
            L=L.*(z-x0(j))/(x0(k)-x0(j));
         end
      end
      s=s+L*y0(k);
   end
   y=s;%y(i)=s;
%end