International Islamic University Chittagong

Dept. of Electrical and Electronic Engineering

Experiment-6
z-Transform in Matlab

EEE-3604 Digital Signal Processing Sessional

Prepared By
Mohammed Abdul Kader
Assistant Professor, Dept. of EEE, IIUC
Objectives
a) To realize the significance of z-transform
b) To determine the z-transform and inverse z-transform of discrete time signal and systems
in MATLAB
c) To find the pole-zero plot and impulse response of DT system.

2 Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
Significance of 𝒛−𝒏
The z-transform of a discrete-time signal x(n) is defined as the power series

𝑋 𝑧 = 𝑥(𝑛)𝑧 −𝑛
𝑛=−∞
Where z is a complex variable. The relation sometimes called the direct z-transform because
it transforms the time-domain signal x(n) into its complex-plane representation X(z).

Where,
𝐳 = 𝐫𝐞𝐣𝛚
‘r’ is a real number
𝒛−𝒏 = 𝒓−𝒏 𝐞−𝐣𝛚𝐧 𝑒 𝑗𝜔 is Euler's Number
= 𝒓−𝒏 [𝒄𝒐𝒔 𝝎𝒏 − 𝒋𝒔𝒊𝒏(𝝎𝒏)] 𝜔 is the angular frequency in radians per
sample.

3
Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
Significance of 𝒛−𝒏 (Cont.)

n=0:.1:100;
r=1.1; title('cos(wt) when w=pi/4');
w=pi/4; subplot(3,1,3);
y1=r.^(-n); plot(n,y);
y2=cos(w*n); xlabel(' ');
y=y1.*y2; ylabel('amplitude');
subplot(3,1,1); title('Real part of z^-^n= r^-^n *cos(wn)');
plot(n,y1);
xlabel('Sample No');
ylabel('amplitude');
title('r^-^n when r=1');
subplot(3,1,2);
plot(n,y2);
xlabel('Sample No');
ylabel('amplitude');

4
Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
Significance of 𝒛−𝒏 (Cont.)

Plot of 𝑟 −𝑛 𝑐𝑜𝑠 𝜔𝑛 when r<1 and Plot of 𝑟 −𝑛 𝑐𝑜𝑠 𝜔𝑛 when r>1 and
𝜋 𝜋
𝜔 = 4 per samples 𝜔 = 4 per samples
5 Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
Significance of 𝒛−𝒏 (Cont.)

When,
r>1, 𝒛−𝒏 has exponentially decreasing oscillation
r<1, 𝒛−𝒏 has exponentially increasing oscillation
r=1, 𝒛−𝒏 has oscillation of constant amplitude.
So we can say
𝒛−𝒏 represents a set of oscillating
components of constant or increasing or
decreasing amplitude based on value of z

Plot of 𝑟 −𝑛 𝑐𝑜𝑠 𝜔𝑛 when r=1 and

Examle-1 Examle-2
>> syms a n f; >> syms n
>> f=a^n; >> x=1/4^n;
>> ztrans(f) >> xz=ztrans(x)
ans = xz =
-z/(a - z) z/(z - 1/4)
Examle-3

>> syms w n;
>> x=sin(w*n);
>> xz=ztrans(x)

xz =

(z*sin(w))/(z^2 - 2*cos(w)*z + 1)

7 Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
Finding inverse z-transform of a function

>> syms z
>> x=2*z/(2*z-1);
>> xn=iztrans(x)

xn =

(1/2)^n

8
Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
Z-Plane and Impulse response

function []=zplane_impulseResponse(Num,Den)
figure('Name','Z-Plane','position',[179 83 560 420]);
zplane(Num,Den);
set(gca,'Color',[0.89 0.99 0.85]); % gca(get current axes handle), [RGB Triplet ]
title('Poles and Zeros');
figure('Name','Impulse response','position',[751 83 560 420]);
impz(Num,Den);
set(gca,'Color',[0.89 0.92 0.92]);
title('Impulse Response')
end

9
Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
System with complex conjugate poles
𝑧 𝑧
𝐻 𝑧 = 2 =
𝑧 − 0.707𝑧 + 0.2499 𝑧 − 0.3535 + 𝑗0.3535 (𝑧 − (0.3535 − 𝑗0.3535))

zplane_impulseResponse([1 0], [1 -0.707 .2499])

10
Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
System with complex conjugate poles in unit circle
zplane_impulseResponse([1],[1 -1.4144 1]);
𝑧
𝐻 𝑧 = 2
𝑧 − 1.4144𝑧 + 1

11 Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
System with complex conjugate poles outside of unit circle
zplane_impulseResponse([1,0],[1 -2.02 1.14])

12 Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC
 Finding Partial Fraction
>> [R,P,K]=residueZ(num,den);
Where R,P and K forms the partial fraction in the following forms,
𝑅1 𝑧 𝑅2 𝑧
𝑋 𝑧 = + + ⋯⋯+ 𝐾
(𝑧 − 𝑝2 ) (𝑧 − 𝑝2 )

Example: Find the partial fraction of following expression

4 − 7 4 𝑧 −1 + (1 4)𝑧
−2
𝑋 𝑧 =
1 − 3 4 𝑧 −1 + (1 8)𝑧
−2

Code Output
>> num=[4 (-7/4) (1/4)]; r =3 -1
>> den=[1 (-3/4) (1/8)]; p = 0.5000 0.2500
>> [r,p,k]=residuez(num,den) k=2

3z z
𝑋 𝑧 = − +2
(𝑧 − 0.5) (𝑧 − 0.25)

13 Experiment-6, Digital Signal Processing Sessional, Prepared By- Mohammed Abdul Kader, Assistant Prof, Dept. of EEE, IIUC