United International University

Dept. of Electrical & Electronic Engineering (EEE)

Fall 2022
DIGITAL SIGNAL PROCESSING
LABORATORY EEE 3310
Lab Report – 4

Submitted to:
Dr. Khawza Iftekhar Uddin Ahmed

Submitted by:
GROUP 1

Name ID

Joyanta Debnath 021 182 032

Ranjan Kumer Sarker 021 191 027

Shakib Hasan 021 201 054

Mushfiqul Hoque 021 201 103

Date of Performance: 05th November 2022

Date of Submission: 12th November 2022
 Lab Report - 4
Frequency Domain Analysis using MATLAB
Introduction: From this lab session-04 we get to know about Z-Transformation, IZ-Transformation,
Pole-Zero representation, F-response, Stability & Causality and get some idea about DTFT and finally
we implement some examples using MATLAB on these topics

Theoretical Background: During the class before going every lab session at first our teacher gave
us theoretical knowledge about that particular topic and the examples. This is how we can get to know
about Z-T, IZ-T, Pole-Zero, F-response, Stability & Causality etc. and so on.

In Lab Evaluation:
1. Consider the system having transfer function,
1 + 2 z −1
H (z) =
1 − 0.5 z −1 + z −2
(i) Sketch the impulse response using MATLAB. Comment on the stability of
the system.
(ii) Plot pole-zero diagram and comment on the stability of the system
assuming that the system is causal

Code:
b=[1 2];
a=[1 -.5 1];
subplot(2,1,1);
impz(b,a)
subplot(2,1,2);
zplane(b,a)

Page 2 of 15
 2. Compute and plot the frequency response of a causal LTI discrete-time system
with a transfer function given by,
0.15(1 − z −2 )
(a) H (z) =
1 − 0.5 z −1 + 0.7 z −2

0.15(1 − z −2 )
(b) H (z) =
0.7 − 0.5 z −1 + z −2
What is the difference between the two filters in (a) and (b), respectively? Which
one will you choose for filtering and why?

Code: a
b=[.15 0 -.15];
a=[1 -.5 .7];
[H,w] = freqz(b,a,50); magH = abs(H); phaH = angle(H);
subplot(2,1,1);plot(w/pi,magH);grid
xlabel('frequency in pi units'); ylabel('Magnitude');
title('Magnitude Response')
subplot(2,1,2);plot(w/pi,phaH/pi);grid

Page 3 of 15
xlabel('frequency in pi units'); ylabel('Phase in pi
units');
title('Phase Response')

Code: b
b=[.15 0 -.15];
a=[.7 -.5 1];
[H,w] = freqz(b,a,50); magH = abs(H); phaH = angle(H);
subplot(2,1,1);plot(w/pi,magH);grid
xlabel('frequency in pi units'); ylabel('Magnitude');
title('Magnitude Response')

Page 4 of 15
subplot(2,1,2);plot(w/pi,phaH/pi);grid
xlabel('frequency in pi units'); ylabel('Phase in pi
units');
title('Phase Response')

Homework:
1. Consider the system having transfer function,
1 − 1.2z −1
H( z) =
1 − 0.5z −1
(i) Plot the impulse response using MATLAB. Comment on the stability
of the system.
(ii) Plot pole-zero diagrams and comment on the stability of the system.

Page 5 of 15
Code:
b=[1 -1.2];
a=[1 -0.5];
subplot(2,1,1);
impz(b,a)
subplot(2,1,2);
zplane(b,a)

2. Consider the system having transfer function,

1 + z −1
H ( z) =
1 − 2 z −1 + z −2
(i) Plot the impulse response using MATLAB. Comment on the stability
of the system.
(ii) Plot pole zero diagram and comment on the stability of the system.

Page 6 of 15
Code:
b=[1 1];
a=[1 -2 1];
subplot(2,1,1);
impz(b,a)
subplot(2,1,2);
zplane(b,a)

3. Consider a system which has poles at

p1, 2 = re  j / 4

and a zero at
z1 = 0.707r

Page 7 of 15
 where: (i) r = 0.96 (ii) r = 0.71 and (iii) r = 0.14
(a) Derive analytically the impulse response of the system and show its
dependence on r. Sketch the impulse response for r=0.96
(b) Plot the frequency and phase response for case (i)
(c) Plot the frequency and phase response for case (ii)
(d) Plot the frequency and phase response for case (iii)
(e) Observe the differences in the the three plots.

Code: a

a) P1=.96*e^(j*pi/4) = .68+j.68 . P2 = .68 – j.68. Z1 = 0.707*.96 = .68

b=[.68];
a=[0.68+.68j .68-.68j];
impz(b,a)

Page 8 of 15
Code: b

b=[.68];
a=[0.68+.68j .68-.68j];
[H,w] = freqz(b,a,100);
magH = abs(H);
phaH = angle(H);
subplot(2,1,1);
plot(w/pi,magH);
grid
xlabel('frequency in pi units'); ylabel('Magnitude');
title('Magnitude Response')
subplot(2,1,2);
plot(w/pi,phaH/pi);
grid
xlabel('frequency in pi units'); ylabel('Phase in pi
units');
title('Phase Response')

Page 9 of 15
P1=.71*e^(j*pi/4) = .502+j.502 . P2 = .502 – j.502. Z1 = 0.707*.71 = .502
Code: c
b=[.502];
a=[.502+.502j .502-.502j];
[H,w] = freqz(b,a,100);
magH = abs(H);
phaH = angle(H);
subplot(2,1,1);
plot(w/pi,magH);
grid
xlabel('frequency in pi units'); ylabel('Magnitude');
title('Magnitude Response')
subplot(2,1,2);
plot(w/pi,phaH/pi);
grid
xlabel('frequency in pi units'); ylabel('Phase in pi
units');
title('Phase Response')

Page 10 of 15
P1=.14*e^(j*pi/4) = .098+j.098 . P2 = .098 – j.098. Z1 = 0.707*.14 = .098

Code: c

b=[.098];
a=[.098+.098j .098-.098j];

Page 11 of 15
[H,w] = freqz(b,a,100);
magH = abs(H);
phaH = angle(H);
subplot(2,1,1);
plot(w/pi,magH);
grid
xlabel('frequency in pi units'); ylabel('Magnitude');
title('Magnitude Response')
subplot(2,1,2);
plot(w/pi,phaH/pi);
grid
xlabel('frequency in pi units'); ylabel('Phase in pi
units');
title('Phase Response')

Page 12 of 15
 4. Consider the following system:
z −3 − 1.8z −2 + 1.62z −1 − 0.729
H(z) =
1 − 1.8z −1 + 1.62z −2 − 0.729z −3

(a) Find the poles and zeros of the transfer function.

(b) Plot the frequency and phase response of the system on a linear scale.

(c) Examine the frequency and phase response.

(d) What is the special characteristic of the system.

Code: a

b=[-.729 1.62 -1.8 1];

a=[1 -1.8 1.62 -.729]
zplane(b,a)

Page 13 of 15
Code: b

b=[-.729 1.62 -1.8 1];

a=[1 -1.8 1.62 -.729];
[H,w] = freqz(b,a,200);
magH = abs(H);
phaH = angle(H);
subplot(2,1,1);
plot(w/pi,magH);
grid
xlabel('frequency in pi units'); ylabel('Magnitude');
title('Magnitude Response')
subplot(2,1,2);
plot(w/pi,phaH/pi);
grid
xlabel('frequency in pi units'); ylabel('Phase in pi
units');
title('Phase Response')

Page 14 of 15
Statement: I Have contributed and given my full support in completing this report.

Page 15 of 15