Mercurial > hg > pmhd

annotate extra/D.m @ 13:844d341cf643 tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Back up before ISMIR
author Yading Song <yading.song@eecs.qmul.ac.uk>
date Thu, 31 Oct 2013 13:17:06 +0000
parents 6840f77b83aa
children
rev   line source
yading@10 1 function y = D(x,N)
yading@10 2 % Calculate rectangular window transform (Dirichlet kernel)
yading@10 3 y = sin(N*x/2)./sin(x/2);
yading@10 4 y(find(y~=y))=N; % avoid NaN if x==0
yading@10 5 end