Butterworth Filter Design: Ω N is the order of the filter
Butterworth Filter Design: Ω N is the order of the filter
Butterworth Filter Design: Ω N is the order of the filter
Lowpass Butterworth filters are characterized by the magnitude-squared frequency response |H()|2 = 1/(1+ (/ c)2N ) where, c is the unnormalized -3dB frequency in rad/sec, N is the order of the filter. Since H(s) H(-s) evaluated at s = j equal to |H()|2 it follows that H(s) H(-s) = 1/(1+(-s2/ c2)N) (2) (1)
The poles of H(s) H(-s) occur on a circle of radius c at equally spaced points, which can be found by using (2) as (-s2/ c2) = (-1)1/N =e j (2k+1) / N , k = 0,1, 2, 3, ....., N-1 (3)
k0
k
s-s
(4)
From (4), we can get the following information N Normalized Denominator Polynomials in Factored Form 1 2 3 4 5 6 7 8 9 10 (1+s) (1+1.414s+s2) (1+s)(1+s+s2) (1+0.765s+s2)(1+1.848s2) (1+s)(1+0.618s+s2)(1+1.618s2) (1+0.518s+s2)(1+1.414s+s2)(1+1.932s+s2) (1+s)(1+0.445s+s2)(1+1.247s+s2)(1+1.802s+s2) (1+0.390s+s2)(1+1.111s+s2)(1+1.663s+s2)(1+1.962s+s2) (1+s)(1+0.347s+s2)(1+s+s2)(1+1.532s+s2)(1+1.879s+s2) (1+0.313s+s2)(1+0.908s+s2)(1+1.414s+s2)(1+1.782s+s2)(1+1.975s+s2)
Example Now let us find the coefficients of IIR Butterworth filter of 1 st order (same procedure can be followed for higher orders) with cutoff frequency of 1 kHz. It is given that the signal being passed through the filter is sampled at 10 kHz. If x(n) represents the input and y(n) represents the filtered output then they are related as
y ( n) = b0 x (n ) + b1 x(n - 1) + a1 y ( n - 1)
where b0 b1 a1 are the filter coefficients to be found out. As N = 1, from (4) we get H(s) = 1/(s+1) ( normalized cutoff frequency is 1 rad/sec) Normalized cutoff frequency in rad / sample () = (1000/10000)2 = 0.2 Unnormalized cutoff frequency in rad / sec ()= tan(/2) = tan(0.1) = 0.3249 For the required low pass filter, analog domain transformation is s s/0.3249. Hence H(s) = 1/(3.0779s+1) ( unnormalized cutoff frequency is 0.3249 rad/sec)
from which we can get b0=0.2495, b1=0.2495, a0=0.51. In the same way we can get the coefficients for higher order filters. If the design requires some other parameters like stop-band edge frequency (s), pass-band tolerance(1), stop-band tolerance (2) then the filter order can be found out using the relation
D log 2 D1 N W log s W p
0.5
where
D1 = 1
(1 - d1 )
-1
D2 =
1 -1 d 22
After finding the filter order the above procedure can be used to find the filter coefficients. If the desired filter is High-pass or Band-pass or Band-stop then the following analog transformations can be applied on low-pass analog filter transfer function.
s
s Wp
s Wp
W pW p ' s
(Low-pass to High-pass)
s 2 + Wl Wu (Low-pass to Band-pass) s ( Wu - W l )
s ( Wu - Wl ) (Low-pass to Band-stop) s 2 + Wl W u
W p ' is new pass band edge frequency and Wl , Wu are lower and upper cutoff frequencies. Characteristics Some important characteristics of Butterworth filters are described below. 1. Lowpass Butterworth filters are all-pole filters. 2. Maximally flat magnitude response and it is monotonic in passband as well as in stopband. 3. Roll-off is slower compared to chebyshev filters. So to get the given transition, order of Butterworth filter required is higher than the order of chebyshev filter. 4. The phase response of a Butterworth filter is more nearly linear than that of a Chebyshev Type II filter of same order. 5. No great differences in the unit impulse and step responses of Butteraworth and Cheyshev Type II filters. References 1. 'Digital signal processing Principles, Algorithms, and Applications' fourth edition, John G. Proakis and Dimitris G. Manolakis. 2. 'Theory and Applications of Digital Signal Processing', L. R. Rabiner and B. Gold. 3. Www.dspguide.com