多項式插值逼近

1.拉格朗日逼近

function [C,L]=lagran(X,Y)
%   X爲橫座標
%   Y爲縱座標
%   L爲拉格朗日多項式係數
%   C=L*Y爲構造的拉格朗日多項式的解C=L1*Y1+......Ln*Yn
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));  
            %poly(X(j)=X-X(j),故conv(V,poly(X(j)))=(X-X1)*....*(X-Xn+1)
        end
    end
    L(k,:)=V;
end
C=Y*L;