%拉格朗日插值多項式擬合
%___________________________________________________________
%調用合驗證：
%調用函數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;