480 likes | 756 Views
CE 40763 Digital Signal Processing Fall 1992 Design of digital FIR filters using the Windowing Technique. Hossein Sameti Department of Computer Engineering Sharif University of Technology. Design of Digital Filters. LTI Systems h(n). FIR. IIR. Determine coefficients of
E N D
CE 40763Digital Signal ProcessingFall 1992Design of digital FIR filtersusing the Windowing Technique HosseinSameti Department of Computer Engineering Sharif University of Technology
Design of Digital Filters LTI Systems h(n) FIR IIR Determine coefficients of h(n) [or P(z) and Q(z)] With rational transfer function No rational transfer function Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Design of digital filters • Design Stages • Specifications Application dependent • Design h(n) Determine coefficients of h(n) • Realization Direct form I,II, cascade and parallel • Implementation Programming in Matlab/C, DSP, ASIC,… • Design of FIR filters • Windowing Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Motivation: impulse response of ideal-low-pass filter • IDTFT of ideal low-pass filter: Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Motivation: impulse response of ideal low-pass filter Multiply by a rectangular window • It can be shown that if we have a linear-phase ideal filter and we multiply it by a symmetric window function, we end up with a linear-phase FIR filter. Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Incorporation of Generalized Linear Phase • Windows are designed with linear phase in mind • Symmetric around M/2 • So their Fourier transform are of the form • Will keep symmetry properties of the desired impulse response • Assume symmetric desired response • With symmetric window • Periodic convolution of real functions Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Design of FIR filters using windows • The steps in the design of FIR filters using windows are as follows: • Start with the desired frequency response results in the sinc function in time domain • Compute • Determine the appropriate window function w(n) • Calculate A finite-length window function Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Desired frequency response • Two properties should be considered: • 1) The amplitude is unity in the pass band and it is zero in the stop band: • 2) The phase is linear: Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Example: Design of a high-pass FIR filter • First, we have to decide on the type of the filter. • Assume Type I filter • (linear-phase) Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Example: Design of a high-pass FIR filter IIR filter Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Example: Design of a high-pass FIR filter • It is a high-pass FIR filter with 7 taps that approximates the high-pass IIR filter. • How can we quickly check that the resulting FIR filter has the desired properties that we were looking for? (i.e., it is a high-pass linear-phase filter)? Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Reminder: DTFT Pairs HosseinSameti, ECE, UBC, Summer 2012 Originally Prepared by: MehrdadFatourechi,
Windowing in frequency domain • What condition should we impose on W(ω) so that H(ω) looks like Hd(ω) ? • Impulse function in the frequency domain, means an infinitely-long constant in the time-domain • Larger window means more computation Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Windowing in Frequency Domain • Windowed frequency response • The windowed version is smeared version of desired response • If w[n]=1 for all n, then W(ej) is pulse train with 2 period Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Rationale for the shape of the filter Ideal filter RectangularWindowfunction (Oppenheim and Schaffer, 2009) Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Filter Specifications Pass-band: Stop-band: Pass-band ripple: Stop-band ripple: Transition width: (Oppenheim and Schaffer, 2009) • What is the ideal situation? Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Filter Specifications Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Observations Width of transition is not sharp! • The width of transition depends on the width of the main lobe of the window. • Ripples in the passband / stopband are proportional to the peaks of side lobes of the window. Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Controlling the width of the main lobe • Q: How can we control the transition width (size of the main lobe)? • A1: using the size of the window Uncertainty principle Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Controlling the width of the main lobe • Q: How can we control the size of transition width (size of the main lobe)? • A2: Shape of the window; in other words, windows with a fixed size that have different shapes can have different main lobe width. • Rectangular window Smallest; and Blackman largest main lobe width Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Controlling the peak of the side lobe • Q: How can we control the peak of the side lobes so that we can get a good ripple behavior in the FIR filter? • A: using the shape of the window Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Controlling the peak of the side lobe • Q: Can we control the peak of the side lobes by changing the size of the window? • A: It can be shown that changes are not significant. Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Demonstration using Kaiser window Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Properties of Windows • Prefer windows that concentrate around DC in frequency • Less smearing, closer approximation • Prefer window that has minimal span in time • Less coefficient in designed filter, computationally efficient • So we want concentration in time and in frequency • Contradictory requirements • Example: Rectangular window Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Rectangular Window • Narrowest main lobe • 4/(M+1) • Sharpest transitions at discontinuities in frequency • Large side lobes • -13 dB • Large oscillation around discontinuities • Simplest window possible Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Bartlett (Triangular) Window • Medium main lobe • 8/M • Side lobes • -25 dB • Hamming window performs better • Simple equation Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Hanning Window • Medium main lobe • 8/M • Side lobes • -31 dB • Hamming window performs better • Same complexity as Hamming Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Hamming Window • Medium main lobe • 8/M • Good side lobes • -41 dB • Simpler than Blackman Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Blackman Window • Large main lobe • 12/M • Very good side lobes • -57 dB • Complex equation Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Frequency response of some popular windows (M=50) rectangular Bartlett Hanning Hamming Blackman Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Peak Approximation Error • Approximation Error, defined in passband and stopband. • Peake Approximation Error is the maximum value of Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Comparison of different windows Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Good design strategy Main lobe Shape of the window width of the window Main lobe Side lobe Good design strategy: 1) Use shape to control the behavior of the side lobe. 2) Use width to control the behavior of the main lobe. Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Kaiser window Zeroth order modified Bessel function of the first kind Number of taps Parameter to control the shape of the Kaiser window and thus the trade-off between the width of the main lobe and the peak of the side lobe. Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Demonstration of Kaiser window M=20 Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Demonstration of Kaiser window Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Comparison with popular windows Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Design Guidelines using Kaiser window • Calculate the transition bandwidth • Calculate • Choose • Choose Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Example: Design of LPF using Kaiser window Specs: Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Example: Design of LPF using Kaiser window Specs: Type II filter Use Bessel equation to get w(n) Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Example: Design of LPF using Kaiser window Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Example: Design of LPF using Kaiser window Q: Does it satisfy the specs? Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology
Summary • Windowing method is a fast and efficient solution to design FIR filters. • Using Kaiser windows, the window can be chosen automatically. Hossein Sameti, Dept. of Computer Eng., Sharif University of Technology