Scribd
Upload
19 views10 pages

Lab 10 SS

Uploaded by

obaid ur rehman
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
19 views10 pages

Lab 10 SS

Uploaded by

obaid ur rehman
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 10

LAB #10

Z-transform and Inverse Z-transform Analysis

Objective:

To study the Z-transform and Inverse Z-transform practically using MATLAB.

Tool Used: MATLAB

Description:

In mathematics and signal processing, the Z-transform converts a discrete time-domain


signal, which is a sequence of real or complex numbers, into a complex frequency-domain
representation.

The Z-transform, like many other integral transforms, can be defined as either a one-sided or
two-sided transform.

Bilateral Z-transform

The bilateral or two-sided Z-transform of a discrete-time signal x[n] is the function X(z)
defined as

Unilateral Z-transform

Alternatively, in cases where x[n] is defined only for n ≥ 0, the single-sided or unilateral Z-
transform is defined as
In signal processing, this definition is used when the signal is causal.

As analog filters are designed using the Laplace transform, recursive digital filters are
developed with a parallel technique called the z-transform. The overall strategy of these
two transforms is the same: probe the impulse response with sinusoids and exponentials
to find the system's poles and zeros. The Laplace transforms deals with differential
equations, the s-domain, and the s-plane. Correspondingly, the z-transform deals with
difference equations, the z-domain, and the z-plane. However, the two techniques are not
a mirror image of each other; the s-plane is arranged in a rectangular coordinate system,
while the z-plane uses a polar format. Recursive digital filters are often designed by
starting with one of the classic analog filters, such as the Butterworth, Chebyshev, or
elliptic. A series of mathematical conversions are then used to obtain the desired digital
filter. The Z transform of a discrete time system X[n] is defined as Power Series.

Rational Z-transform to factored Z-transform:


Example:
Let the given transfer function be in the rational form,
2z4+16z3+44z2+56z+32
G(z)=
3z4+3z3-15z2+18z-12
It is required to convert it into factored form, so that we can find the poles and zeros
mathematically by applying quadratic equation.
Matlab command required for converting rational form to factored form be
‘Zp2sos’
The factored form of G(z) as evaluated by ‘zp2sos’ be,

G(z)=( 0.6667 + 0.4z-1 + 0.5333 z-2) (1.000 + 2.000 z-1 +2.000 z-2)
(1.000 + 2.000z-1 -4.000z-2 )(1.000 - 1.000 z-1 + 1.000 z-2)

Factored Z-transform / zeros,poles to rational Z-transform:


It is the inverse of the above case, when the transfer function is given in factored form and it
is required to convert in rational form then a single ‘matlab’ command can serve the purpose.
Example:
Lets use the above result i-e;transfer function in factored for,

G(z)=( 0.6667 + 0.4z-1 + 0.5333 z-2) (1.000 + 2.000 z-1 +2.000 z-2)
(1.000 + 2.000z-1 -4.000z-2 )(1.000 - 1.000 z-1 + 1.000 z-2)

For building up transfer function in rational form we find the poles and zers of above system
simply by using MATLAB ‘root’ command or by hand. Or simply we have poles and zeros
of the given system we can find the transfer function in factored form.

Matlab command that converts poles and zeros of the system in to transfer function is ‘zp2tf’
.

Rational Z-transform to partial fraction form:

This technique is usually used, while taking the inverse Z-transform and when the order
‘H(z)’ is high so that it is quite difficult to solve it mathematically.
Example:
Consider the transfer function in the rational form i-e;
18z3
G(z)=
18z3+3z2-4z-1
We can evaluate the partial fraction form of the above system using MATLAB command.
The partial fraction form be,
G(z)= 0.36 + 0.24 + _0.4

1 – 0.5z-1 1+0.33 z-1 (1+0.33 z-1)

Matlab command that converts rational z-transform in to partial fraction form is


‘residuez’.
Partial fraction form to Z-transform:
This technique is used when it is required to convert partial fraction expression in to rational
Z-transform.
Example:

Take the partial fraction form of above,


The partial fraction form be,
G(z)= 0.36 + 0.24 + _0.4

1 – 0.5z-1 1+0.33 z-1 (1+0.33 z-1)

Matlab command that converts partial fraction form into rational z-transform is
‘residuez’
Zplane:
Zero-pole plot
zplane(b,a)
This function displays the poles and zeros of discrete-time systems.

MATLAB:
syms z n
a=ztrans(1/16^n)

Inverse Z-Transform:

MATLAB:
syms Z n
iztrans(3*Z/(Z+1))

Pole Zero Diagrams For A Function In Z Domain:


Z plane command computes and display the pole-zero diagram of Z function.
The Command is
Zplane(b,a)
To display the pole value, use root(a)
To display the zero value, use root(b)
Matlab Code:

b=[0 1 1 ]
a= [1 -2 +3]
roots(a)
roots(b)
zplane(b,a);

ans =
1.0000 + 1.4142i
1.0000 - 1.4142i

ans=
-1
Frequency Response:

The Freqz function computes and display the frequency response of given Z- Transform of
the function
freqz(b,a,Fs)
b= Coeff. Of Numerator
a= Coeff. Of Denominator
Fs= Sampling Frequency

Matlab Code:

b=[2 5 9 5 3]
a= [5 45 2 1 1]
freqz(b,a);

Example:

Plot the magnitude and phase of the frequency response of the given digital filter
Using freqz function:
Y(n) = 0.2x(n) + 0.52y(n-1) – 0.68(y(n-2)

Matlab Code:

b = [0.2];
a= [1, -0.52, 0.68];
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');
POST LAB TASK:

Task#1:
Express the following z-transform in factored form , plot its poles and zeros,and then
determine its ROCs.
2z4+16z3+44z2+56z+32
G(z)=
3z4+3z3-15z2+18z-12

Task#2:

Determine the partial fraction expansion of the z-transform G(z) given by

18z3
G(z)=
18z3+3z2-4z-1
Conclusion:

In this lab, I learnt about the Z-transformation. Z-transform converts a discrete time-
domain signal, which is a sequence of real or complex numbers, into a complex
frequency-domain representation.
Then I performed the Z – Transform and inverse Z transform of various equation and
calculated my results in MATLAB.

You might also like