0.6 Lab 5b - digital filter design (part 2)  (Page 3/6)

To further investigate the Kaiser window, plot the Kaiser windows and their DTFT magnitudes (in dB)for $N=21$ and the following values of $\beta$ :

• $\beta =0$
• $\beta =1$
• $\beta =5$

For each case use at least 512 points in the plot of the DTFT.

Next use a Kaiser window to design a low pass filter, $h\left(n\right)$ , to remove the noise from the signal in nspeech2.mat using equations [link] and [link] . To do this, use equations [link] and [link] to compute the values of $\beta$ and $N$ that will yield the following design specifications:

The low pass filter designed with the Kaiser method will automatically have a cut-offfrequency centered between ${\omega }_p$ and ${\omega }_s$ .

Plot the magnitude of the DTFT of $h\left(n\right)$ for $\left|\omega \right|<\pi$ . Create three plots in the same figure:one that shows the entire frequency response, and ones that zoom in on the passband and stopband ripple, respectively.Mark ${\omega }_p$ , ${\omega }_s$ , ${\delta }_p$ , and ${\delta }_s$ on these plots where appropriate. Note:Since the ripple is measured on a magnitude scale, DO NOT use a decibel scale on this set of plots.

From the Matlab prompt, compute the stopband and passband ripple (do not do this graphically). Record the stopband and passband ripple tothree decimal places.

Filter the noisy speech signal in nspeech2.mat using the filter you have designed. Then compute theDTFT of 400 samples of the filtered signal starting at time $n=20000$ (i.e. 20001:20400 ).Plot the magnitude of the DTFT samples in decibels versus frequency in radians for $\left|\omega \right|<\pi$ . Compare this with the spectrum of the noisy speech signal shownin [link] . Play the noisy and filtered speech signals back using sound and listen to them carefully.

Inlab report

1. Submit the values of $\beta$ and $N$ that you computed.
2. Submit the three plots of the filter's magnitude response. Make sure the plots are labeled.
3. Submit the values of the passband and stopband ripple. Does this filter meet the design specifications?
4. Submit the magnitude plot of the DTFT in dB for the filtered signal.Compare this plot to the plot of [link] .
5. Comment on how the frequency content and the audio quality of the filtered signal have changed afterfiltering.

Fir filter design using parks-mcclellan algorithm

Click here for help on the firpm function for Parks-McClellan filter design. Download the data file nspeech2.mat for the following section.

Kaiser windows are versatile since they allow the design of arbitrary filters which meet specific design constraints.However, filters designed with Kaiser windows still have a number of disadvantages. For example,