yading@10: function y = D(x,N)
yading@10: % Calculate rectangular window transform (Dirichlet kernel)
yading@10: y = sin(N*x/2)./sin(x/2);
yading@10: y(find(y~=y))=N; % avoid NaN if x==0
yading@10: end