>
section 6 of 103 min read

6. Z-Transforms: The Discrete-Time Cousin

The Z-transform is to discrete-time what Laplace is to continuous-time. Same job, different domain.

6.1 Definition

X(z)=n=x[n]znX(z) = \sum_{n=-\infty}^{\infty} x[n]\, z^{-n}

The variable zz is a complex number. You can think of zz as the discrete-time analog of the Laplace ss, with the relationship z=esTsz = e^{sT_s} that comes from sampling a continuous-time signal: a Laplace exponential este^{st} evaluated at t=nTst = nT_s gives esnTs=(esTs)n=zne^{snT_s} = (e^{sT_s})^n = z^n.

That mapping has a beautiful consequence: the imaginary axis in the s-plane (where s=jωs = j\omega, the home of all sinusoids) maps to the unit circle in the z-plane, z=ejωTsz = e^{j\omega T_s}. The frequency response of a discrete-time system is what you get by evaluating H(z)H(z) on the unit circle. The left half of the s-plane (stable region) maps to the inside of the unit circle (stable region for discrete-time).

6.2 Standard Z-transforms

SequenceZ-transformRegion of convergence
δ[n]\delta[n]11All zz
u[n]u[n]1/(1z1)1/(1 - z^{-1})$
anu[n]a^n u[n]1/(1az1)1/(1 - a z^{-1})$
anu[n1]-a^n u[-n-1]1/(1az1)1/(1 - a z^{-1})$
nanu[n]n a^n u[n]az1/(1az1)2a z^{-1}/(1 - a z^{-1})^2$

Two different time signals (the third and fourth rows) have the same algebraic form but different regions of convergence. So the Z-transform is incomplete without specifying the ROC.

6.3 Properties

Just like Laplace and Fourier:

  • Linearity (sums and scalings).
  • Time shift: x[nk]zkX(z)x[n - k] \leftrightarrow z^{-k} X(z). Each one-sample delay multiplies the Z-transform by z1z^{-1}. So z1z^{-1} is the algebraic embodiment of "delay by one sample."
  • Convolution becomes multiplication: x1[n]x2[n]X1(z)X2(z)x_1[n] * x_2[n] \leftrightarrow X_1(z)\, X_2(z).
  • Differentiation in zz: nx[n]zdX(z)/dzn x[n] \leftrightarrow -z\, dX(z)/dz.

6.4 Stability in the z-plane

For a discrete-time LTI system described by transfer function H(z)H(z):

  • All poles inside the unit circle (z<1|z| < 1): stable.
  • Poles on the unit circle: marginally stable.
  • Any pole outside the unit circle: unstable.

This is the discrete-time analog of "left half plane equals stable" for continuous-time, derived directly from the z=esTsz = e^{sT_s} mapping.

6.5 Why Z-transform matters

When you implement a digital filter in software or in an FPGA, the difference equation looks like:

y[n]=b0x[n]+b1x[n1]+a1y[n1]a2y[n2]y[n] = b_0 x[n] + b_1 x[n-1] + \cdots - a_1 y[n-1] - a_2 y[n-2] - \cdots

Take the Z-transform of both sides (each nkn-k shift becomes zkz^{-k}), get an algebraic equation, and solve for H(z)=Y(z)/X(z)H(z) = Y(z)/X(z):

H(z)=b0+b1z1+b2z2+1+a1z1+a2z2+H(z) = \frac{b_0 + b_1 z^{-1} + b_2 z^{-2} + \cdots}{1 + a_1 z^{-1} + a_2 z^{-2} + \cdots}

Now H(z)H(z) tells you everything: stability (pole locations), frequency response (evaluate HH on the unit circle, z=ejωz = e^{j\omega}), transient response. We will use the Z-transform extensively in Chapter 17 (DSP) when we design FIR and IIR digital filters.