Inverse Z Transform in Matlab Programming

In Digital Signal Processing (DSP) and control systems engineering, the Z-Transform transforms discrete-time signals into the complex frequency domain (Z-domain) to simplify stability analysis and filter design. However, to implement physical systems or analyze time-domain responses, you must convert the frequency representation back to a discrete sequence using the Inverse Z-Transform.

MATLAB provides built-in tools within the Symbolic Math Toolbox, primarily the iztrans() function, to compute exact symbolic inverse Z-transforms.

What Is the Inverse Z-Transform?

The Inverse Z-Transform converts a Z-domain transfer function (X(z)) back to its corresponding discrete-time sequence (x[n]):

This conversion is vital for determining impulse responses, evaluating difference equations, and inspecting transient behaviors in digital filters.

MATLAB Function Syntax: iztrans()

The standard syntax for computing inverse Z-transforms in MATLAB is:

x = iztrans(F) x = iztrans(F, z, n)

• F: Symbolic expression in the Z-domain.
• z: Z-domain variable.
• n: Discrete-time index variable.

Step-by-Step Examples

Example 1: Basic First-Order System

Find the inverse Z-transform of the first-order transfer function:

syms z n X = z / (z - 0.5);
x = iztrans(X, z, n)

Output:

x = (1/2)^n

This result represents a geometrically decaying discrete-time sequence characteristic of stable first-order systems.

Example 2: Second-Order Polynomial Transfer Function

Find the inverse Z-transform of:

syms z n X = z^2 / (z^2 - 1.5*z + 0.5);
x = iztrans(X, z, n)

Output:

x = 1^n + (1/2)^n

The response combines a unit step component and a decaying exponential, matching second-order system dynamics.

Example 3: Generalized Symbolic Constant

Compute the inverse Z-transform using an arbitrary constant a:

syms z n a X = z / (z - a);
x = iztrans(X, z, n)

Output:

x = a^n

Visualizing Time-Domain Sequences in MATLAB

To inspect discrete-time sequence values, evaluate the output over a time index array and plot it with stem():

n_vals = 0:10;
x_vals = (0.5).^n_vals;
stem(n_vals, x_vals, 'filled', 'MarkerFaceColor', 'b');
title('Discrete-Time Sequence x[n] = (0.5)^n');
xlabel('Sample index n');
ylabel('Amplitude x[n]');
grid on;

Summary

• Use iztrans() to solve inverse Z-transforms symbolically in MATLAB.
• Always define syms z n before declaring Z-domain transfer functions.
• Pair your symbolic analytical results with stem() for discrete-time visualization.