關於部落格
  • 122904

    累積人氣

  • 31

    今日人氣

    0

    追蹤人氣

Lagrange Base Code for Octave/Matlab

function M=Lag_Base(l,xi,x)
% Lag_Base - The Lagrange Basis
%               (x(1)-x)(x(2)-x)...(x(l-1)-x)(x(l+1)-x)...(x(k)-x)
%   L_l(x)=-----------------------------------------------------------------
%           (x(1)-x(l))(x(2)-x(l))...(x(l-1)-x(l))(x(l+1)-x(l))...(x(k)-x(l))
%   INPUT ARGUMENTS:
%     l - the index of Lagrange Base
%     xi - the points which we want to interpolat
%     x - the varlue
if length(x)==1
    n=length(xi);
    %Main file
    M=1;
    for j=1:n
        if j~=l
            %(xi(j)-x)/(xi(j)-xi(l))
            M=M*(xi(j)-x)/(xi(j)-xi(l));
        end
    end
else
    for iter=1:length(x)
        M(:,iter)=Lag_Base(l,xi,x(iter));
    end
end
Download this code:

相簿設定
標籤設定
相簿狀態