lagrange.m

本次的算法是相位周期展开聚类

function h=lagrange(N,delay)
n=0:N;
h=ones(1,N+1);
for k=0:N
    index=find(n~=k);
    h(index)=h(index)*(delay-k)./(n(index)-k);
end