拉格朗日插值多项式拟合.m

来自「拉格朗日插值多项式拟合,牛顿插值多项式,欧拉方程解偏微分方程,使用极限微分求解导」· M 代码 · 共 27 行

%拉格朗日插值多项式拟合
%___________________________________________________________
%调用合验证：
%调用函数lagran(X,Y)即得到拟合多项式C
%然后用polyval(C,z0)去验证，即求的是当x=zo的时候，多项式的值
%___________________________________________________________
% X是自变量，Y是对应的函数值
%拟合问题即是给定两组向量，分别是x量合y量
%然后用相应的多项式去拟合
%___________________________________________________________
function [C,L]=lagran(X,Y)
N=length(X)
n=N-1;
L=zeros(N,N);
for k=1:N
    V=1;
    for j=1:N
        if k~=j
            %其中poly(z0)得到的是一个多根为z0的多项式，比如poly(2),得到的是【1 －2】即x－2
            V=conv(V,poly(X(j)))/(X(k)-X(j));
        end
    end
        L(k,:)=V;
    end
    C=Y*L;