Scribd
Upload
8 views

DSP Lab Report 5

Uploaded by

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

DSP Lab Report 5

Uploaded by

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

Z-Transform and Inverse Z-

Transform

3 / 28 / 24

05
EE-20-A

Muhammad Raza Madni 210401034


Riyan Abdullah 210401013
Lab # 05
Objectives:
By the end of this lab students will be able to:

• Become familiar with the computation of the z-transform of various signals and with the
determination of ROC.
• Learn and understand the properties of the z-transform and how these properties can be used to
simplify computations.
• Learn the process of inverting the z-transform by using method of partial fraction expansion.
• Understand how LTI systems are represented in the z-domain and the relationship to the frequency
response.
• Learn how to solve difference equations that describe LTI systems with initial conditions.
• Understand the difference between bilateral and unilateral z-transform.

Equipment’s:
• Laptop
• MATLAB (software)

Introduction to MATLAB
MATLAB is a high-level programming language that has been used extensively to solve complex
engineering problems.
Task # 1
Using MATLAB, determine the partial-fraction expansion of the z-transform X(z) given by
18z3
X(z) =
18z3 + 3z2 – 4z -1

Answers:
Following is the MATLAB function code.
syms z
b = [18 0 0 0];
a = [18 3 -4 -1];
[r, p ,k] = residue(b,a);
sum=[];
for i=1:length(a)-1
sum = [sum ; (k + r(i)/1-p(i)*z^(-i))];
end
disp(sum);
Task # 2
Determine the inverse z-transform of the 11 terms of following equation.
1- 2z-1
X(z) =
(1-0.2z-1) (1+0.6z-1)

Answers:
Following is the MATLAB function code.

Task 2
b = [1 -2 0];
a = [1 0.4 -0.12];
impz(b,a,11);

Task # 3
Determine the Inverse z-transform of the following sequences, using partial fraction expansion method.
x(z) = __1 – z-1 – 4z-2 + 4z-3_______
1 – 2.75z-1 + 1.625z-2 – 0.25z-3

Answers:
Following is the MATLAB function code.
syms z
b = [1 -1 -4 4];
a = [1 -2.75 1.625 -0.25];
[r, p ,k] = residue(b,a);
sum=[];
for i=1:length(a)-1
sum = [sum ; (k + r(i)/1-p(i)*z^(-i))];
end
invs=iztrans(sum);
disp(sum);
disp(invs);

Task # 4
A digital filter is described by the difference equation,

Y(n) = x(n) + x(n-1) + 0.9y(n-1) – 0.81y(n-2)


Answers:
Following is the MATLAB function code.
b = [1 1 0];
a = [1 -0.9 0.81];
w = (0:1:500)*pi/500;
H=freqz(b,a,w);
magH = abs(H);
phaH = angle(H)*180/pi;
subplot(2,1,1);
plot(w/pi,magH);
title('Magnitude Response');
xlabel('frequency in pi units');
ylabel('│H│');
subplot(2,1,2);
plot(w/pi,phaH);
title('Phase Response');
xlabel('frequency in pi units');
ylabel('Degrees');

…the end !

You might also like