lagran.m

来自「余向阳博士Matlab教程的源代码。 包括：1.非线性方程的解法；2.线性方程」· M 代码 · 共 17 行

function[c,l]=lagran(x,y)  
% 这里x为n个节点的横坐标所组成的向量，y为纵坐标所组成的向量。
% c为所得插值函数的系数组成的向量。
w=length(x);
n=w-1;
l=zeros(w,w);
for k=1:n+1
  v=1;
     for j=1:n+1
         if k~=j
             v=conv(v,poly(x(j)))/(x(k)-x(j));
    end
end
  l(k,:)=v;
end
c=y*l;