lagran.m

来自「数值分析上机实验」· M 代码 · 共 16 行

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;