lagranl.m

来自「数值计算实践是我们必备的一门学科！这是我在学习之余利用MATLAB编写的程序」· M 代码 · 共 18 行

function [C,L,L1,l]=lagranl(X,Y)
%L为拉格朗日插值多项式;
%C为拉格朗日插值多项式的系数向量;
%X为节点的横坐标向量;
%Y为节点的纵坐标向量;
%L1为基函数的系数矩阵;
m=length(X);L=ones(m,m);
for k=1:m
    V=1;
    for i=1:m
        if k~=i
            V=conv(V,poly(X(i)))/(X(k)-X(i));
        end
    end
    L1(k,:)=V;
    l(k,:)=poly2sym(V)
end
C=Y*L1;L=Y*l;