Further Design Parameters:
f1 = 530; % lower passband limit (Hz) N = 2^(nextpow2(8*fs/f1)) % FFT size from xition width if N<8*M, N = 8*(M-1); end; % spectral interpolation
Bandlimited ``Ideal'' Hilbert transformer frequency response (FR), including transition bands matched to FFT size such that :
H = [ ([0:k1-2]/(k1-1)).^8,ones(1,k2-k1+1), ... ([k1-2:-1:0]/(k1-1)).^8,zeros(1,N/2-1)]; Hp = [H(N/2+2:N), H(1:N/2+1)]; // pos. freqs only f = [-0.5 + 1/N : 1/N : 0.5]; normalized freq axis plot(f,Hp); grid; ...