<< Chapter < Page Chapter >> Page >

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

  • β = 0
  • β = 1
  • β = 5

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

To create theKaiser windows, use the Matlab command kaiser(N,beta) command where N is the length of the filter and betais the shape parameter β . To insure at least 512 points in the plot use the command DTFT(window,512) when computing the DTFT.
Submit the plots of the 3 Kaiser windows and the magnitude of their DTFT's in decibels.Comment on how the value β affects the shape of the window and the sidelobes of the DTFT.
DTFT of a section of noisy speech.

Next use a Kaiser window to design a low pass filter, h ( n ) , 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 β and N that will yield the following design specifications:

ω p = 1 . 8 ω c = 2 . 0 ω s = 2 . 2 δ p = 0 . 05 δ s = 0 . 005

The low pass filter designed with the Kaiser method will automatically have a cut-offfrequency centered between ω p and ω s .

ω c = ω p + ω s 2

Plot the magnitude of the DTFT of h ( n ) for | ω | < π . 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 ω p , ω s , δ p , and δ 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.

To compute the passband ripple, find the value of the DTFT at frequencies corresponding to the passband using the command H(abs(w)<=1.8) where H is the DTFT of h ( n ) and w is the corresponding vector of frequencies. Then use this vector to compute the passband ripple.Use a similar procedure for the stopband ripple.

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 | ω | < π . 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

Do the following:
  1. Submit the values of β 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,

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Purdue digital signal processing labs (ece 438). OpenStax CNX. Sep 14, 2009 Download for free at http://cnx.org/content/col10593/1.4
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Purdue digital signal processing labs (ece 438)' conversation and receive update notifications?

Ask