0.6 Lab 5b - digital filter design (part 2)

Questions or comments concerning this laboratory should be directedto Prof. Charles A. Bouman, School of Electrical and Computer Engineering, Purdue University, West Lafayette IN 47907;(765) 494-0340; bouman@ecn.purdue.edu

Introduction

This is the second part of a two week laboratory in digital filter design.The first week of the laboratory covered some basic examplesof FIR and IIR filters, and then introduced the concepts of filter design.In this week we will cover more systematic methods of designing both FIR and IIR filters.

Filter design using standard windows

Download DTFT.m for the following section.

We can generalize the idea of truncation by using different windowing functions to truncate anideal filter's impulse response. Note that by simply truncating the ideal filter's impulse response, we are actuallymultiplying (or “windowing”) the impulse response by a shifted $rect\left(\right)$ function. This particular type of window is called a rectangular window. In general, the impulse reponse $h\left(n\right)$ of the designed filter is related to the impulse response $h_{ideal}\left(n\right)$ of the ideal filter by the relation

where $w\left(n\right)$ is an $N$ -point window. We assume that

where ${\omega }_c$ is the cutoff frequency and $N$ is the desired window length.

The rectangular window is defined as

The DTFT of $w\left(n\right)$ for $N=21$ is shown in [link] . The rectangular window is usually not preferredbecause it leads to the large stopband and passband ripple as shown in [link] .

More desirable frequency characteristics can be obtained by making a better selection for the window, $w\left(n\right)$ . In fact, a variety of raised cosine windowsare widely used for this purpose. Some popular windows are listed below.

1. Hanning window (as defined in Matlab, command hann(N) ):
   $$w\left(n\right)=\left\{\begin{array}{cc}0.5-0.5cos\frac{2\pi n}{N-1}\hfill & n=0,1,...,N-1\hfill \\ 0\hfill & \text{otherwise}\hfill \end{array}\right)$$
2. Hamming window
   $$w\left(n\right)=\left\{\begin{array}{cc}0.54-0.46cos\frac{2\pi n}{N-1}\hfill & n=0,1,...,N-1\hfill \\ 0\hfill & \text{otherwise}\hfill \end{array}\right)$$
3. Blackman window
   $$w\left(n\right)=\left\{\begin{array}{cc}0.42-0.5cos\frac{2\pi n}{N-1}+0.08cos\frac{4\pi n}{N-1}\hfill & n=0,1,...,N-1\hfill \\ 0\hfill & \text{otherwise}\hfill \end{array}\right)$$

In filter design using different truncation windows, there are two key frequency domain effects that are important to the design: the transitionband roll-off , and the passband and stopband ripple (see [link] below). There are two corresponding parameters in the spectrum of each type ofwindow that influence these filter parameters. The filter's roll-off is related tothe width of center lobe of the window's magnitude spectrum. The ripple is influenced by the ratio of the mainlobe amplitude to the first sidelobe amplitude (or difference if using a dB scale).These two window spectrum parameters are not independent, and you should see a trend as you examine the spectra for different windows.The theoretical values for the mainlobe width and the peak-to-sidelobe amplitude are shown in [link] .