复件 复件 lagrange1.m

来自「拉格朗日插值计算对于给定的数据（xj,yi） 求一插值多项式」· M 代码 · 共 16 行

function y=largrange(x0,y0,x)
n=length(x0);m=length(x);
for i=1:m
    z=x(i);
    s=0.0;
    for k=1:n
        p=1.0;
        for j=1:n
            if j~=k
                p=p*(z-x0(j))/(x0(k)-x0(j));
            end
        end
        s=p*y0(k)+s;
    end
    y(i)=s;
end