lagrange.m

来自「本次的算法是相位周期展开聚类」· M 代码 · 共 7 行

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