Scribd
Upload
49 views

Digital Signal Processing Systems, Basic Filtering Types, and Digital Filter Realizations

The document discusses digital signal processing systems and digital filtering. It provides examples of computing the response of discrete-time systems using difference equations and the MATLAB filter function. It also demonstrates using freqz to obtain the frequency response of different filter types, including lowpass, highpass, bandpass and bandstop filters.

Uploaded by

melki
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
49 views

Digital Signal Processing Systems, Basic Filtering Types, and Digital Filter Realizations

The document discusses digital signal processing systems and digital filtering. It provides examples of computing the response of discrete-time systems using difference equations and the MATLAB filter function. It also demonstrates using freqz to obtain the frequency response of different filter types, including lowpass, highpass, bandpass and bandstop filters.

Uploaded by

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

Digital Signal Processing Systems,

Basic Filtering Types, and Digital


Filter Realizations
DSP Labs
The Difference Equation and Digital Filtering
• Given the DSP system

y (n) = 2 x(n) − 4 x(n − 1) − 0.5 y (n − 1) − y (n − 2)

with initial conditions y(-2) = 1, y(-1) = 0, x(-1) = 1, and the input

x(n) = 0.8n u (n)

Compute the system response y(n) for 20 samples!


• MATLAB function filter( ), developed using a direct-form II
realization, can be used to operate digital filtering, and the syntax is
• Plots of the input and system outputs y(n):
Zi = filtic(B,A,Yi,Xi)
y = filter(B,A,x,Zi)

• B and A are vectors for the coefficients bj and aj , whose formats are

A = [1 a1 a2 . . . aN] and B = [b0 b1 b2 . . . bM]

• x and y are the input data vector and the output data vector, respectively.

• Note that the filter function filtic( ) is a MATLAB function used to obtain
initial states required by the MATLAB filter function filter( ) (requirement by
a direct-form II realization) from initial conditions in the difference equation.

• Zi contains initial states required for operating MATLAB function filter( ), that
is,
Zi = [w(-1) w(-2) . . .],
which can be recovered by the MATLAB function, filtic( ). Xi and Yi
are initial conditions with a length of the greater of M or N, given by

Xi = [x(-1) x(-2) . . .] and Yi = [y(-1) y(-2) . . .]

• Especially for zero initial conditions, the syntax is reduced to

y = filter(B,A,x)
Example-1
Compute the system output
y (n) = 0.5 y (n − 2) + x(n − 1)

for the first four samples using the following initial conditions:

a. Initial conditions: y(-2) = 1, y(-1) = 0, x(-1) = -1, and input:


x(n) = 0.5n u (n)

b. Zero initial conditions: y(-2) = 0, y(-1) = 0, x(-1) = 0, and input:


x(n) = 0.5n u (n)
Solution-1
a. The MATLAB codes and results with the nonzero initial conditions are
listed as

>> B = [0 1]; A = [1 0 -0.5];


>> x = [1 0.5 0.25 0.125];
>> Xi = [-1 0]; Yi = [0 1];
>> Zi = filtic(B,A,Yi,Xi);
>> y = filter(B,A,x,Zi)
y =
-0.5000 1.0000 0.2500 0.7500
>>
b. The MATLAB codes and results with the zero initial conditions are listed
as

>> B = [0 1]; A = [1 0 -0.5];


>> x = [1 0.5 0.25 0.125];
>> y = filter(B,A,x)
y =
0 1.0000 0.5000 0.7500
>>
Basic Types of Filtering
• Given a transfer function, the MATLAB function freqz( ) can be used to
determine the frequency response. The syntax is given by

[h,w] = freqz(B,A,N)

whose parameters are defined as:


– h = an output vector containing frequency response
– w = an output vector containing normalized frequency values distributed in the range
from 0 to  radians.
– B = an input vector for numerator coefficients
– A = an input vector for denominator coefficients
– N = the number of normalized frequency points used for calculating the frequency
response
Example-2
Given each of the following digital transfer functions,

z
a. H ( z) =
z − 0.5

b. H ( z ) = 1 − 0.5 z −1

0.5 z 2 − 0.32
c. H ( z) = 2
z − 0.5 z + 0.25
1 − 0.9 z −1 + 0.81z −2
d. H ( z ) =
1 − 0.6 z −1 + 0.36 z −2

1. Plot the poles and zeros on the z-plane.


2. Use MATLAB function freqz() to plot the magnitude frequency response and
phase response for each transfer function.
3. Identify the corresponding filter type, such as lowpass, highpass, bandpass, or
bandstop.
Solution-2
1. The transfer functions of (a) and (c)
need to be converted into the standard
form (delay form) required by the
MATLAB function freqz( ), in which
both numerator and denominator
polynomials have negative powers of z.
Hence, we obtain
z 1
H ( z) = =
z − 0.5 1 − 0.5 z −1

0.5 z 2 − 0.32 0.5 − 0.32 z −2


H ( z) = 2 =
z − 0.5 z + 0.25 1 − 0.5 z −1 + 0.25 z −2

while the transfer functions of (b) and


(d) are already in their standard forms
(delay forms).
2. The MATLAB program for plotting the magnitude frequency response and the
phase response for each case is listed in Program 6.2.
• From the plots of magnitude frequency responses for all cases, we can
conclude that case is a lowpass filter
• From the plots of magnitude frequency responses for all cases, we can
conclude that case is a highpass filter
• From the plots of magnitude frequency responses for all cases, we can
conclude that case is a bandpass filter
• From the plots of magnitude frequency responses for all cases, we can
conclude that case is a bandstop (band reject) filter.
Exercise-1
Given a filter
1 + 2 z −1 + z −2
H ( z) =
1 − 0.5 z −1 + 0.25 z −2

use MATLAB to plot

1. its magnitude frequency response;


2. its phase response.
Exercise-2
Given the difference equation
y (n) = x(n − 1) − 0.75 y (n − 1) − 0.125 y (n − 2)

1. Use the MATLAB functions filter( ) and filtic( ) to


calculate the system response y(n) for n = 0, 1, 2, 3, . . . , 4 with the
input of x(n) = (0.5)n u(n) and initial conditions: x(-1) = -1, y(-2) = 2,
and y(-1) = 1;

2. Use the MATLAB function filter( ) to calculate the system


response y(n) for n = 0, 1, 2, 3, . . . , 4 with the input of x(n) = (0.5)n
u(n) and zero initial conditions: x(-1) = 0, y(-2) = 0, and y(-1) = 0.

You might also like