<< Chapter < Page | Chapter >> Page > |
Although not part of the approximation addressed, the phase curve is also very smooth.
An important feature of the Butterworth filter is the closed- form formula for the solution, . The expression for may be determined as
This function has poles evenly spaced around a unit radius circle and zeros at infinity. The determination of is very simple. In order to have a stable filter, is selected to have the left-hand plane poles and zeros at infinity; will necessarily have the right-hand plane poles and the other N zeros at infinity. Thelocation of these poles on the complex plane for , and 4 is shown in [link] .
Because of the geometry of the pole positions, simple formulas are easy to derive for the pole locations. If the realand imaginary parts of the pole location are denoted as
the locations of the poles are given by
for values of where
Because the coefficients of the numerator and denominator polynomials of are real, the roots occur in complex conjugate pairs. The conjugate pairs in [link] , [link] can be combined to be the roots of second-order polynomials so that for even, has the partially factored form of
for . For odd, has a single real pole and, therefore, the form
for
This is a convenient form for the cascade and parallel realizations discussed in elsewhere.
A single formula for the pole locations for both even and odd is
for values of where
One of the important features of the Butterworth filter design formulas is that the pole locations are found by independentcalculations which do not depend on each other or on factoring a polynomial. A FORTRAN program which calculates these values isgiven in the appendix as Program 8. Mathworks has a powerful command for designing analog and digital Butterworth filters.
The classical form of the Butterworth filter given in [link] is discussed in many books [link] , [link] , [link] , [link] , [link] . The less well-known form given in [link] also has many useful applications [link] . If the frequency location of unwanted signals is known, the zeros of the transfer function given by thenumerator can be set to best reject them. It is then possible to choose the pole locations so as to have a passband as flat as theclassical Butterworth filter by using [link] . Unfortunately, there are no formulas for the pole locations; therefore, thedenominator polynomial must be factored.
This section considers the process of going from given specifications to use of the approximation results derived in theprevious section. The Butterworth filter is the simplest of the four classical filters in that all the approximation effort is placed attwo frequencies: and . The transition from passband to stopband occurs at a normalized frequency of . Assuming that this transition frequency or bandedge can later be scaled to any desired frequency, the only parameter tobe chosen in the design process is the order .
The filter specifications that are consistent with what is optimized in the Butterworth filter are the degree of “flatness" at (DC) and at . The higher the order, the flatter the frequency response at these two points. Because ofthe analytic nature of rational functions, the flatter the response is at and , the closer it stays to the desired response throughout the whole passband and stopband. Anindirect consequence of the filter order is the slope of the response at the transition between pass and stopband. The slope ofthe squared-magnitude frequency response at is
The effects of the increased flatness and increased transition slope of the frequency response as N increases are illustrated in Figure 1 from Design of Infinite Impulse Response (IIR) Filters by Frequency Transformations .
In some cases specifications state the response must stay above or below a certainvalue over a given frequency band. Although this type of specification is more compatible with a Chebyshev erroroptimization, it is possible to design a Butterworth filter to meet the requirements. If the magnitude of the frequency response of thefilter over the passband of must remain between unity and , where and , the required order is found by the smallest integer satisfying
This is illustrated in [link] where must remain above 0.9 for up to 0.9, i.e., = 0.9 and = 0.9. These requirements require an order of at least .
If stopband performance is stated in the form of requiring that the response stay below a certain value for frequency above acertain value, i.e., for , the order is determined by the same formula [link] with replaced by .
Note which is called the “half power" frequency because . This frequency is normalized to one for the theory but can be scaled to any value for applications.
To illustrate the calculations, a lowpass Butterworth filter is designed. It is desired that the frequency response stay above0.8 for frequencies up to 0.9. The formula [link] for determining the order gives a value of 2.73; therefore, the order isthree. The analytic function corresponding to the squared-magnitude frequency response in [link] is
The transfer function corresponding to the left-half-plane poles of F'(s) are calculated from [link] to give
The frequency response is obtained by setting which has a plot illustrated in [link] for . The pole locations are the same as shown in [link] c.
Notification Switch
Would you like to follow the 'Digital signal processing and digital filter design (draft)' conversation and receive update notifications?