Butterworth filter

From Wikipedia, the free encyclopedia

Linear analog electronic filters
Butterworth filter
Chebyshev filter
Elliptic (Cauer) filter
Bessel filter
Gaussian filter
Optimum "L" (Legendre) filter
edit

The Butterworth filter is one type of electronic filter design. It is designed to have a frequency response which is as flat as mathematically possible in the passband. Another name for them is 'maximally flat magnitude' filters.

The Butterworth type filter was first described by the British engineer Stephen Butterworth in his paper "On the Theory of Filter Amplifiers", Wireless Engineer (also called Experimental Wireless and the Radio Engineer), vol. 7, 1930, pp. 536-541.

1 Overview
2 A Simple Example
3 The transfer function
4 Filter Design
- 4.1 Cauer topology
- 4.2 Sallen-Key topology
5 Comparison with other linear filters

[edit] Overview

The frequency response of the Butterworth filter is maximally flat (has no ripples) in the passband, and rolls off towards zero in the stopband. When viewed on a logarithmic Bode plot, the response slopes off linearly towards negative infinity. For a first-order filter, the response rolls off at −6 dB per octave (−20 dB per decade). (All first-order filters, regardless of name, are actually identical and so have the same frequency response.) For a second-order Butterworth filter, the response decreases at −12 dB per octave, a third-order at −18 dB, and so on. Butterworth filters have a monotonically decreasing magnitude function with ω. The Butterworth is the only filter that maintains this same shape for higher orders (but with a steeper decline in the stopband) whereas other varieties of filters (Bessel, Chebyshev, elliptic) have different shapes at higher orders.

Compared with a Chebyshev Type I/Type II filter or an elliptic filter, the Butterworth filter has a slower roll-off, and thus will require a higher order to implement a particular stopband specification. However, Butterworth filter will have a more linear phase response in the passband than the Chebyshev Type I/Type II and elliptic filters.

[edit] A Simple Example

A third order low pass filter (Cauer topology). The filter becomes a Butterworth filter with cutoff frequency ωc=1 when (for example) C2=4/3 farad, R4=1 ohm, L1=1/2 and L3=3/2 henry.

A third order low pass filter (Cauer topology). The filter becomes a Butterworth filter with cutoff frequency ω_c=1 when (for example) C₂=4/3 farad, R₄=1 ohm, L₁=1/2 and L₃=3/2 henry.

Log density plot of the transfer function H(s) in complex frequency space for the third order Butterworth filter with ω_c=1. Note the three poles which lie on a circle of unit radius in the left half plane.

A simple example of a Butterworth filter is the 3rd order low-pass design shown in the figure on the right, with $C 2 = 4 / 3$ farad, $R 4 = 1$ ohm, $L 1 = 1 / 2$ and $L 3 = 3 / 2$ henry. Taking the impedance of the capacitors C to be 1/Cs and the impedance of the inductors L to be Ls, where $s = σ + j ω$ is the complex frequency, the circuit equations yields the transfer function for this device:

$H(s)=\frac{V_o(s)}{V_i(s)}=\frac{1}{1+2s+2s^2+s^3}$

The gain $G (ω)$ as a function of angular frequency $ω$ is given by:

$G^2(\omega)=|H(j\omega)|^2=\frac{1}{1+\omega^6}\,$

and the phase is given by:

$\Phi(\omega)=\arg(H(j\omega))\,$

The group delay is defined as the derivative of the phase with respect to angular frequency and is a measure of the distortion in the signal introduced by phase differences for different frequencies. The $log 10$ (gain) and the delay for this filter are plotted in the graph on the left. It can be seen that there are no ripples in the gain curve in either the passband or the stop band.

The log of the absolute value of the transfer function H(s) is plotted in complex frequency space in the second graph on the right. The function is defined by the three poles in the left half of the complex frequency plane. These are arranged on a circle of radius unity, symmetrical about the real s axis. The gain function will have three more poles on the right half plane to complete the circle.

By replacing each inductor with a capacitor and each capacitor with an inductor, a high-pass Butterworth filter is obtained. If we change each capacitor and inductor into a resonant capacitor and inductor in parallel, with the proper choice of component values, a band-pass Butterworth filter is obtained.

log 10

(gain) and group delay of the third order Butterworth filter with ω_c=1

[edit] The transfer function

Plot of the gain of Butterworth low-pass filters of orders 1 through 5. Note that the slope is 20n dB/decade where n is the filter order.

Like all filters, the typical prototype is the low-pass filter, which can be modified into a high-pass filter, or placed in series with others to form band-pass and band-stop filters, and higher order versions of these.

The gain $G (ω)$ of an n-order Butterworth low pass filter is given in terms of the transfer function H(s) as:

$G^2(\omega)=\left |H(j\omega)\right|^2 = \frac {G_0^2}{1+\left(\frac{\omega}{\omega_c}\right)^{2n}}$

where

n = order of filter
ω_c = cutoff frequency (approximately the -3dB frequency)
$G 0$ is the DC gain (gain at zero frequency)

It can be seen that as n approaches infinity, the gain becomes a rectangle function and frequencies below ω_c will be passed with gain $G 0$ , while frequencies above ω_c will be suppressed. For smaller values of n, the cutoff will be less sharp.

We wish to determine the transfer function H(s) where $s = σ + j ω$ . Since H(s)H(-s) evaluated at s = jω is simply equal to |H(ω)|², it follows that:

$H(s)H(-s) = \frac {G_0^2}{1+\left (\frac{-s^2}{\omega_c^2}\right)^n}$

The poles of this expression occur on a circle of radius ω_c at equally spaced points. The transfer function itself will be specified by just the poles in the negative real half-plane of s. The k-th pole is specified by:

$-\frac{s_k^2}{\omega_c^2} = (-1)^{\frac{1}{n}} = e^{\frac{j(2k-1)\pi}{n}} \qquad\mathrm{k = 1,2,3, \ldots, n}$

and hence,

$s_k = \omega_c e^{\frac{j(2k+n-1)\pi}{2n}}\qquad\mathrm{k = 1,2,3, \ldots, n}$

The transfer function may be written in terms of these poles as:

$H(s)=\frac{G_0}{\prod_{k=1}^n (s-s_k)/\omega_c}$

The denominator is a Butterworth polynomial in s.

[edit] Normalized Butterworth polynomials

The Butterworth polynomials may be written in complex form as above, but are usually written with real coefficients by multiplying pole pairs which are complex conjugates, such as $s 1$ and $s n$ . The polynomials are normalized by setting $ω c = 1$ . The normalized Butterworth polynomials then have the general form:

$B_n(s)=\prod_{k=1}^{\frac{n}{2}} \left[s^2-2s\cos\left(\frac{2k+n-1}{2n}\,\pi\right)+1\right]$ for n even

$B_n(s)=(s+1)\prod_{k=1}^{\frac{n-1}{2}} \left[s^2-2s\cos\left(\frac{2k+n-1}{2n}\,\pi\right)+1\right]$ for n odd

To four decimal places, they are:

n	Factors of Polynomial $B n (s)$
1	$(s + 1)$
2	$s 2 + 1.4142 s + 1$
3	$(s + 1)(s 2 + s + 1)$
4	$(s 2 + 0.7654 s + 1)(s 2 + 1.8478 s + 1)$
5	$(s + 1)(s 2 + 0.6180 s + 1)(s 2 + 1.6180 s + 1)$
6	$(s 2 + 0.5176 s + 1)(s 2 + 1.4142 s + 1)(s 2 + 1.9319 s + 1)$
7	$(s + 1)(s 2 + 0.4450 s + 1)(s 2 + 1.2470 s + 1)(s 2 + 1.8019 s + 1)$
8	$(s 2 + 0.3902 s + 1)(s 2 + 1.1111 s + 1)(s 2 + 1.6629 s + 1)(s 2 + 1.9616 s + 1)$

[edit] Maximal flatness

Assuming $ω c = 1$ and $G 0 = 1$ , the derivative of the gain with respect to frequency can be shown to be:

$\frac{dG}{d\omega}=-nG^3\omega^{2n-1}$

which is monotonically decreasing for all $ω$ since the gain G is always positive. The gain function of the Butterworth filter therefore has no ripple. Furthermore, the series expansion of the gain is given by:

$G(\omega)=1 - \frac{1}{2}\omega^{2n}+\frac{3}{8}\omega^{4n}+\ldots$

In other words, all derivatives of the gain up to but not including the 2n-th derivative are zero, resulting in "maximal flatness".

[edit] High Frequency Roll Off

Again assuming $ω c = 1$ , the slope of the log of the gain for large ω is:

$\lim_{\omega\rightarrow\infty}\frac{d\log(G)}{d\log(\omega)}=-n$

In decibels, the high frequency roll off is therefore 20n dB/decade. (The factor of 20 is used because the power is proportional to the square of the voltage gain.)

[edit] Filter Design

There are a number of different filter topologies available to implement a linear analog filter. These circuits differ only in the values of the components, but not in their connections.

[edit] Cauer topology

The Cauer topology uses passive components (shunt capacitors and series inductors) to implement a linear analog filter. The Butterworth filter having a given transfer function can be realised using a Cauer 1-form. The k^th element is given by:

$C_k = 2 \sin \left [\frac {(2k-1)}{2n} \pi \right ]$ ; k = odd

$L_k = 2 \sin \left [\frac {(2k-1)}{2n} \pi \right ]$ ; k = even

[edit] Sallen-Key topology

The Sallen-Key topology uses active and passive components (op amps and capacitors) to implement a linear analog filter. Each Sallen-Key stage implements a conjugate pair of poles; the overall filter is implemented by cascading all stages in series. If there is a real pole (in the case where $n$ is odd), this must be implemented separately, usually as an RC circuit, and cascaded with the op-amp stages.

The Sallen-Key transfer function is given by

$H(s)=\frac{1}{1+C_2(R_1+R_2)s+C_1C_2R_1R_2s^2}$

We wish the denominator to be one of the quadratic terms in a Butterworth polynomial. Assuming that $ω c = 1$ , this will mean that

$C_1C_2R_1R_2=1\,$

and

$C_2(R_1+R_2)=2\cos\left(\frac{2k+n-1}{2n}\right)$

This leaves two component values undefined, which may be chosen at will.

[edit] Comparison with other linear filters

Here is an image showing the gain of a discrete-time Butterworth filter next to other common kind of filters obtained with the same number of coefficients: