Scribd
Upload
97 views

DSP LAB of Filter

This document discusses digital filter design in MATLAB. It describes estimating the order of IIR filters using functions like buttord based on filter type (Butterworth, Chebyshev 1/2, elliptic) and specifications. It also discusses designing the actual IIR filter transfer function using functions like butter. FIR filter design using windowing functions fir1 and fir2 to reduce Gibbs phenomenon is also covered. The lab tasks involve designing various filter types in MATLAB to meet given specifications and analyzing the results.

Uploaded by

Mujtaba Farrukh
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
97 views

DSP LAB of Filter

This document discusses digital filter design in MATLAB. It describes estimating the order of IIR filters using functions like buttord based on filter type (Butterworth, Chebyshev 1/2, elliptic) and specifications. It also discusses designing the actual IIR filter transfer function using functions like butter. FIR filter design using windowing functions fir1 and fir2 to reduce Gibbs phenomenon is also covered. The lab tasks involve designing various filter types in MATLAB to meet given specifications and analyzing the results.

Uploaded by

Mujtaba Farrukh
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Lab 12: Digital Filter Design

Objective:
The process of deriving the transfer function H(z) whose frequency response H(ejw) approximates
the given frequency response specifications is called digital filter design. After H(z) has been
obtained, it is then realized in the form of a suitable filter structure. In the previous laboratory
exercise, the realizations of FIR and IIR transfer functions have been considered. In this laboratory
exercise you will learn how to design an IIR or FIR digital filter to meet a specified magnitude or
gain response.

12.1 IIR Filter


The most common method of IIR filter design is based on the bilinear transformation of a prototype
analog transfer function. The analog transfer function is usually one of the following types: Butterworth,
Type 1 Chebyshev, Type 2 Chebyshev, and elliptic transfer functions. The difference between these
filter types can be explained by considering the analog lowpass filter. The Butterworth lowpass transfer
function has a maximally-flat magnitude response at dc, that is, ω = 0, and a monotonically decreasing
magnitude response with increasing frequency.
The Type 1 Chebyshev lowpass transfer function has an equiripple magnitude response in the pass band
and a monotonically decreasing magnitude response with increasing frequency outside the pass band.
The Type 2 Chebyshev lowpass transfer function has a monotonically decreasing magnitude response
in the pass band with increasing frequency and an equiripple magnitude response in the stop band.
Finally, the elliptic low pass transfer function has equiripple magnitude responses both in the pass band
and in the stop band.

12.2 Estimation of Order of IIR Filter


The first step in the filter design process is to choose the type of filter approximation to be employed
and then to estimate the order of the transfer function from the filter specifications. The MATLAB
command for estimating the order of a Butterworth filter is

where the input parameters are the normalized pass band edge frequency Wp, the normalized stop band
edge frequency Ws, the pass band ripple Rp in dB, and the minimum stop band attenuation Rs in dB.
Both Wp and Ws must be a number between 0 and 1 with the sampling frequency assumed to be 2 Hz.
The output data are the lowest order N meeting the specifications and the normalized cutoff frequency
Wn. If Rp = 3 dB, then Wn = Wp. buttord can also be used to estimate the order of a high pass, a band
pass, and a band stop Butterworth filter. For a highpass filter design, Wp > Ws. For bandpass and
bandstop filter designs, Wp and Ws are two-element vectors specifying both edge frequencies, with the
lower edge frequency being the first element of the vector. In the latter cases, Wn is also a two-element
vector.
For estimating the order of a Type 1 Chebyshev filter, the MATLAB command is

For designing a Type 2 Chebyshev filter, the MATLAB command for estimating the order is

Finally, in the case of an elliptic filter design, the command is


As before, Wp and Ws are the pass band and stop band edge frequencies with values between 0 and 1.
Likewise, Rp and Rs are the pass band ripple and the minimum stop band attenuation in dB. N contains
the estimated lowest order and Wn is the cutoff frequency. It should be noted that for band pass and
band stop filter designs, the actual order of the transfer function obtained using the appropriate filter
design command is 2N.

12.3 IIR Filter Design


After the filter type has been selected and its order estimated the next step is to determine the transfer
function of the filter. To this end MATLAB provides functions for all four types of filters. For
designing Butterworth digital low pass or band pass filters, the command is

where the input parameters N and Wn are determined through the use of the function buttord, and
the output is the vectors num and den containing, respectively, the coefficients of the numerator
and denominator polynomials of the transfer function in ascending powers of z-1 . If Wn is a scalar,
butter returns a low pass transfer function of order N, and if Wn is a two-element vector, it returns
a band pass transfer function of order 2N. For designing a Butterworth digital high pass filter of
order N, the command is

The command

returns the transfer function of a Butterworth band stop filter of order 2N provided W n is a two
element vector. For designing a Type 1 Chebyshev digital filter, the commands are

For designing a Type 2 Chebyshev digital filter, the commands are

A low pass transfer function of order N is returned in each case if W n is a scalar, and a band pass
transfer function of order 2N is returned if Wn is a two-element vector. In each of the above
commands, filter type is high for designing a high pass filter with Wn being a scalar, and filter type
is stop for designing a bandstop filter with Wn being a two-element vector.

12.4 FIR Filter Design


Conceptually the simplest approach to FIR filter design is to simply truncate to a finite number of
terms the infinite-length impulse response coefficients obtained by computing the inverse discrete
time Fourier transform of the desired ideal frequency response. However, a simple truncation
results in an oscillatory behavior in the respective magnitude response of the FIR filter, which is
more commonly referred to as the Gibbs phenomenon.
The Gibbs phenomenon can be reduced by windowing the infinite-length impulse response
coefficients by an appropriate finite-length window function. The functions fir1 and fir2 can be
employed to design windowed FIR digital filters in MATLAB. Both functions yield a linear-phase
design. The function fir1 can be used to design conventional low pass, high pass, band pass, and
band stop linear-phase FIR filters. The command
returns in vector b the impulse response coefficients, arranged in ascending powers of z-1 , of a
lowpass or a bandpass filter of order N for an assumed sampling frequency of 2 Hz. For low pass
design, the normalized cutoff frequency is specified by a scalar Wn, a number between 0 and 1.
For band pass design, Wn is a two-element vector [Wn1, Wn2] containing the specified pass band
edges where 0 <Wn1 < Wn2 < 1. The command

with N an even integer, is used for designing a high pass filter. The command

with Wn a two-element vector, is employed for designing a band stop FIR filter.

Lab Task
1. Determine the lowest order of a digital IIR low pass filter of all four types in MATLAB.
The specifications are as follows: sampling rate of 40 kHz, pass band edge frequency of 4
kHz, stop band edge frequency of 8 kHz, pass band ripple of 0.5 dB, and a minimum stop
band attenuation of 40 dB. Comment on your results.
2. Determine the lowest order of a digital IIR high pass filter of all four types in MATLAB.
The specifications are as follows: sampling rate of 3,500 Hz, pass band edge frequency of
1,050 Hz, stopband edge frequency of 600 Hz, pass band ripple of 1 dB, and a minimum
stop band attenuation of 50 dB. Comment on your results.
3. Determine the lowest order of a digital IIR band pass filter of all four types in MATLAB.
The specifications are as follows: sampling rate of 7 kHz, pass band edge frequencies at 1.4
kHz and 2.1 kHz, stop band edge frequencies at 1.05 kHz and 2.45 kHz, pass band ripple of
0.4 dB, and a minimum stop band attenuation of 50 dB. Comment on your results.
4. Given: sampling rate of 20 kHz, δp = 0.002 and δS = 0.002, and stop band edge = 2.3 kHz.
Using the function fir1, design a linear-phase FIR low pass filter meeting the specifications
given above and plot its gain and phase responses. Estimate the order of filter. Does your
design meet the specifications? If it does not, adjust the filter order until the design meets the
specifications. What is the order of the filter meeting the specifications?
5. Design a Butterworth Bandstop Digital Filter.

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy