Matched Z-transform method

The s-plane poles and zeros of a 5th-order Chebyshev type II lowpass filter to be approximated as a discrete-time filter

The z-plane poles and zeros of the discrete-time Chebyshev filter, as mapped into the z-plane using the matched Z-transform method with T = 1/10 second. The labeled frequency points and band-edge dotted lines have also been mapped through the function z=e^iωT, to show how frequencies along the iω axis in the s-plane map onto the unit circle in the z-plane.

The matched Z-transform method, also called the pole–zero mapping or pole–zero matching method, is a technique for converting a continuous-time filter design to a discrete-time filter (digital filter) design.

The method works by mapping all poles and zeros of the s-plane design to z-plane locations $z=e^{sT}$, for a sample interval $T=1/f_{\mathrm {s} }$. So an analog filter with transfer function:

$H(s)=k_{\mathrm {a} }{\frac {\prod _{i=1}^{M}(s-\xi _{i})}{\prod _{i=1}^{N}(s-p_{i})}}$

is transformed into the digital transfer function

$H(z)=k_{\mathrm {d} }{\frac {\prod _{i=1}^{M}(1-e^{\xi _{i}T}z^{-1})}{\prod _{i=1}^{N}(1-e^{p_{i}T}z^{-1})}}$

The gain $k_{\mathrm {d} }$ must be adjusted to normalize the desired gain, typically set to match the analog filter's gain at DC by setting $s=0$ and $z=1$ and solving for $k_{\mathrm {d} }$.

Since the mapping wraps the s-plane's $j\omega$ axis around the z-plane's unit circle repeatedly, any zeros (or poles) greater than the Nyquist frequency will be mapped to an aliased location.

In the (common) case that the analog transfer function has more poles than zeros, the zeros at $s=\infty$ may optionally be shifted down to the Nyquist frequency by putting them at $z=-1$, dropping off like the BLT.

This transform doesn't preserve time- or frequency-domain properties, and so is rarely used. Alternative methods include the bilinear transform and impulse invariance methods. MZT does provide less high frequency response error than the BLT, however, making it easier to correct by adding additional zeros, which is called the MZTi (for "improved").

Responses of the filter (solid), and its discrete-time approximation (dashed), for nominal cutoff frequency of 1 Hz, sample rate 1/T = 10 Hz. The discrete-time filter does not reproduce the Chebyshev equiripple property in the stopband due to the interference from cyclic copies of the response.