Convolution and Filtering
The mathematical foundation of filtering is convolution. For a finite impulse response (FIR) filter, the output y(k) of a filtering operation is the convolution of the input signal x(k) with the impulse response h(k):
If the input signal is also of finite length, you can implement the filtering operation
using the MATLAB®
conv function. For example, to filter a five-sample random vector with a
third-order averaging filter, you can store x(k) in a vector
x, h(k) in a vector
h, and convolve the
x = randn(5,1); h = [1 1 1 1]/4; % A third-order filter has length 4 y = conv(h,x)
y = -0.3375 0.4213 0.6026 0.5868 1.1030 0.3443 0.1629 0.1787
yis one less than the sum of the lengths of
Filters and Transfer Functions
The transfer function of a filter is the Z-transform of its impulse response. For an FIR filter, the Z-transform of the output y, Y(z), is the product of the transfer function and X(z), the Z-transform of the input x:
The polynomial coefficients h(1), h(2), …, h(n + 1) correspond to the coefficients of the impulse response of an nth-order filter.
The filter coefficient indices run from 1 to (n + 1), rather than from 0 to n. This reflects the standard indexing scheme used for MATLAB vectors.
FIR filters are also called all-zero, nonrecursive, or moving-average (MA) filters.
For an infinite impulse response (IIR) filter, the transfer function is not a polynomial, but a rational function. The Z-transforms of the input and output signals are related by
where b(i) and a(i) are the filter coefficients. In this case, the order of the filter is the maximum of n and m. IIR filters with n = 0 are also called all-pole, recursive, or autoregressive (AR) filters. IIR filters with both n and m greater than zero are also called pole-zero, recursive, or autoregressive moving-average (ARMA) filters. The acronyms AR, MA, and ARMA are usually applied to filters associated with filtered stochastic processes.
Filtering with the
For IIR filters, the filtering operation is described not by a simple convolution, but by a difference equation that can be found from the transfer-function relation. Assume that a(1) = 1, move the denominator to the left side, and take the inverse Z-transform to obtain
In terms of current and past inputs, and past outputs, y(k) is
which is the standard time-domain representation of a digital filter. Starting with y(1) and assuming a causal system with zero initial conditions, the representation is equivalent to
To implement this filtering operation, you can use the MATLAB
filter stores the coefficients in two
row vectors, one for the numerator and one for the denominator. For example, to solve the
you can use
b = 1; a = [1 -0.9]; y = filter(b,a,x);
filtergives you as many output samples as there are input samples, that is, the length of
yis the same as the length of
x. If the first element of a is not 1, then
filterdivides the coefficients by a(1) before implementing the difference equation.