WOLFRAM

Wolfram Language & System Documentation Center
Wolfram Language Home Page »

Digital Filter Design

Digital Filter Design MethodsOutput Response—Digital Filtering of Signals
Poles and Zeros of Digital Filters
The Wolfram Language provides a comprehensive set of methods for designing digital filters.
Digital Filter Design Methods
LeastSquaresFilterKernel
create an FIR filter using the least-squares method
FrequencySamplingFilterKernel
create an FIR filter using frequency sampling
EquirippleFilterKernel
create an FIR filter using the equiripple method
ToDiscreteTimeModel
create an IIR digital filter from an analog prototype
Digital filter design methods.

Least-Squares Method

This method obtains a finite impulse response (FIR) from a given prototype filter specification in the frequency domain by means of the inverse discrete-time Fourier transform.
Define a zero-phase ideal lowpass filter with a cutoff frequency of 1.2 radians per sample.
H( j ω){
1,|ω| ≤ ωc
0,otherwise
The impulse response is obtained by taking the inverse discrete-time Fourier transform of the filter specification:
In[1]:=1
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-4t8mky
Out[1]=1
Evaluating this expression for integer values of from to will return a filter of length .
Create a lowpass FIR filter of length 15 with a cutoff frequency of 1.2 radians per sample:
In[2]:=2
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-8bxqkz
Out[2]=2
This method minimizes the mean-squared error between the ideal specification and the resulting FIR filter:
In[3]:=3
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-22g3aj
Out[3]=3
The least-squares method is also known as the window-based method. The mean-squared error can be diminished by applying a smoothing window to the FIR returned by LeastSquaresFilterKernel.
Create a Hann window of length 15:
In[4]:=4
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-0wjrma
Out[4]=4
Apply the window to the lowpass FIR filter h:
In[5]:=5
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-yrsp4m
Out[6]=6
Different windows return different attenuations. Select a window type by the degree of the desired attenuation.
Compare attenuation levels of four selected windows
In[6]:=6
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-48z5y7
Out[6]=6

Frequency Sampling Method

Frequency sampling allows the creation of an FIR filter by specifying the filter's amplitude spectrum on the interval 0≤ ω≤ π at a finite number of uniformly spaced frequency locations.
Sample an ideal lowpass filter at for from 0 to 5:
In[7]:=7
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-qx1i21
Out[7]=7
In[8]:=8
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-smnc2x
Out[8]=8
Create an odd-length, even-symmetry lowpass filter:
In[9]:=9
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-5mjg57
Out[10]=10
Show the ideal lowpass filter, amplitude samples, and the resulting filter:
In[11]:=11
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-pnmqw2
Out[11]=11

Equiripple Method

The Parks–McClellan–Rabiner (1979) algorithm is one of the most popular methods of designing linear-phase FIR filters. It minimizes the maximum deviation from the desired ideal frequency response and results in filters with equal ripples in each of the bands of the filter.
Create a multiband filter of length 55:
In[12]:=12
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-348hxd
Out[12]=12
Plot the magnitude response of an equiripple multiband filter:
In[13]:=13
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-foxntk
Out[13]=13

Create a Digital Filter from an Analog Prototype

A popular method of creating digital filters is to transform analog prototypes to their digital equivalents using the bilinear transformation. The result is an infinite impulse response (IIR) filter, represented as a TransferFunctionModel.
ToDiscreteTimeModel
discrete-time approximation of an analog filter
Create an IIR filter from an analog prototype.
Create an analog lowpass filter of Chebyshev type with a cutoff frequency at radians per second:
In[14]:=14
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-myi9yv
Out[14]=14
Transform the analog filter to a digital filter using bilinear transformation:
In[15]:=15
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-yrb4p3
Out[15]=15
Show the frequency response of the resulting IIR filter:
In[16]:=16
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-1ut3li
Out[16]=16
Poles and Zeros of Digital Filters
TransferFunctionPoles
extract poles of analog filters
TransferFunctionZeros
extract zeros of analog filters
Poles and zeros of analog filters.
Compute and show zeros of a lowpass filter:
In[17]:=17
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-dvxna1
Out[18]=18
In[19]:=19
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-5jv6ok
Out[19]=19
In[20]:=20
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-n5f1n9
Out[20]=20
In[21]:=21
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-21kbba
Out[21]=21
Output Response—Digital Filtering of Signals

FIR Filters

ListConvolve
convolve an FIR filter with data
ImageConvolve
convolve an FIR filter with image
DiscreteConvolve
symbolic convolution of two signals
Apply FIR filters to signals.
Create a lowpass equiripple filter:
In[67]:=67
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-x4gkmi
Out[68]=68
Import and filter an audio signal:
In[65]:=65
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-c2c2wb
Out[65]=65
In[69]:=69
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-g1ghga
Out[69]=69
Lowpass filtering of a scanline of the red channel of an image:
In[28]:=28
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-soy8y4
Out[30]=30
Lowpass filtering of an image:
In[31]:=31
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-4ihaml
Out[31]=31

IIR Filters

RecurrenceFilter
apply an IIR filter to a signal
RecurrenceTable
iteratively solve the recurrence equation
OutputResponse
output response of a filter
Apply IIR filters to signals.
Define a transfer function of a digital Butterworth filter:
In[32]:=32
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-rld01i
Compute the impulse response of the IIR filter:
In[33]:=33
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-6xostg
Out[34]=34
Compute the impulse response by using the coefficients of the transfer function model:
In[35]:=35
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-hfwwnz
Out[35]=35
Compute the impulse response using OutputResponse:
In[36]:=36
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-eair8b
Out[36]=36
Compute the impulse response by solving the corresponding recurrence equation model of the filter using RecurrenceTable:
In[37]:=37
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-5lpx81
Out[37]=37
Lowpass filtering of a scanline of the red channel of an image:
In[38]:=38
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-sl5bi
Out[40]=40
Lowpass filtering of an image:
In[41]:=41
✖
https://wolfram.com/xid/06eqt1wysv9dr98be-md9z15
Out[41]=41