Next  |  Prev  |  Up  |  Top  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

FFT of Truncated(2048) Ideal Impulse Response

\epsfig{file=eps/AnalyticHilbertIdealFR.eps,width=6in}
h = [h(N/2+1:N),h(1:N/2)]; % zero-centered form (almost)
hw = wzp .* h; % Apply window to ideal impulse response
Hw = fft(hw);  % Frequency response we really get

% Compute total stopband attenuation:
ierr = norm(Hw(N/2+2:N))/norm(Hw)
     = 0.0378 = 3.8 percent
For spectral displays, rotate buffer so negative frequencies are on the left and dc is at $ k=N/2$ :
Hwp = [Hw(N/2+2:N), Hw(1:N/2+1)]; % Neg. freqs on left
Hwpn = abs(Hwp); Hwpn = Hwpn/max(Hwpn);
plot(f,20*log10(Hwpn)); grid;
...


Next  |  Prev  |  Up  |  Top  |  JOS Index  |  JOS Pubs  |  JOS Home  |  Search

[Comment on this page via email]

``The Window Method for FIR Digital Filter Design}'', by Julius O. Smith III, (From Lecture Overheads, Music 421).
Copyright © 2020-06-27 by Julius O. Smith III
Center for Computer Research in Music and Acoustics (CCRMA),   Stanford University
CCRMA  [Automatic-links disclaimer]