2 Low-Pass Filter Design
Total Page:16
File Type:pdf, Size:1020Kb
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING 2.161 Signal Processing – Continuous and Discrete 1 Introduction to Continuous Time Filter Design 1 Introduction The design techniques described here are based on the creation of a prototype low-pass filter, with subsequent transformation to other filter forms (high-pass, band-pass, band-stop) if necessary. 2 Low-Pass Filter Design The prototype low-pass filter is based upon the magnitude-squared of the frequency response 2 function |H(jΩ)| , or the frequency response power function. The phase response of the filter is not considered. We begin by defining tolerance regions on the power frequency response plot, as shown in Fig. 1. |H( jW ) |2 1 1 1 + e 2 1 1 + l 2 0 W W W ( 0 c r ra d / sec ) pa s s ba n d tran s ition ba n d stop ba n d Figure 1: Tolerance regions in the frequency response plot. The filter specifications are that ≥|H j |2 > 1 | | < 1 ( Ω) 2 for Ω Ωc (1) 1+ 1D. Rowell, October 2, 2008 1 2 1 and |H(jΩ)| < for |Ω| > Ωr, (2) 2 1+ λ where Ωc is the cut-off frequency, Ωr is the rejection frequency, and and λ are design parameters that select the filter attenuation at the two critical frequencies. For example, if 2 =1, atΩc the power response response |H(jΩc)| =0.5, the -3 dB response frequency. In general we expect the response function to be monotonically decreasing in the transition band. The filter functions examined in this document will be of the form |H j |2 1 . ( Ω) = 2 (3) 1+ f (Ω) where f(Ω) → 0asΩ→ 0, and f(Ω) →∞as Ω →∞to generate a low-pass filter action. 2.1 The Butterworth Filter The Butterworth filter is defined by the power gain 2 1 |H(jΩ)| = 2N (4) 2 1+ (Ω/Ωc) where N is a positive integer defining the filter order. Note that λ does not appear in this formulation, but clearly N and λ are interrelated, since at Ω = Ωr 1 1 ≥ (5) 2 2 2N 1+ λ 1+ (Ωr/Ωc) which may be solved to show log(λ/) N ≥ (6) log(Ωr/Ωc) The power response function of Butterworth filters for N =1... 5 is shown in Fig. 2. Butterworth filters are also known as “maximally flat” filters because the response has the maximum number of vanishing derivatives at Ω = 0 and Ω = ∞ for filters of the form of Eq. 3. 2.1.1 The Transfer Function of the Butterworth Filter The poles of the power gain transfer function may be found from the characteristic equation � � s 2N 2 1+ = 0 (7) jΩc which yields 2N roots (poles) that lie on a circle: −1/N jπ(2n+N −1)/2N pn =Ωc e n =1... 2N (8) −1/N with radius r =Ωc , and angular separation of π/N rad. Notice that if N is odd a pair −1/N of poles will lie on the real axis at s = ±Ωc , while if N is even the roots will form complex conjugate pairs. 2 |H(jω )| 2 1 0.9 0.8 0.7 Power response 0.6 0.5 0.4 0.3 N increasing 0.2 N=1 0.1 N=5 0 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 Normalized frequency ω/ω c Figure 2: Power response of Butterworth filters for N =1... 5( = 1). 2 Figure 3a shows the six poles of |H(s)| for a third-order (N = 3) Butterworth filter. For a causal system we must have |H j |2 H j H j H s H −s | ( Ω) = ( Ω) ( Ω) = ( ) ( ) s=jΩ 2 which allows us to take take the N poles of |H(jΩ)| in the left half-plane as the poles of the filter H(s), that is the poles specified by (n =1...N ) in Eq. (8) above, −1/N jπ(2n+N −1)/2N pn =Ωc e n =1...N (9) as is shown in Fig. 3b. If the filter is to have unity gain at low frequencies lim |H(jΩ)| =1 Ω→0 we require the complete Butterworth transfer function to be (−p1)(−p2) ...(−pN ) H(s)= (s − p1)(s − p2) ...(s − pN ) N (−1) p1p2 ...pN = , (s − p1)(s − p2) ...(s − pN ) where only the N stable left-half-plane poles are included. 3 W W j j s - plane s - pla n e n= 1 n= 1 n= 6 W e - 1 /N x x W e-1/ N x r = c r = c 2p / 6 p/ 3 s n=2 x x x s n= 2 n= 5 x x n= 3 x n= 3 n=4 (a) (b) 2 Figure 3: The poles of (a) |H(s)| , and (b) H(s) for a third-order (N = 3) Butterworth filter. jW W r = c x s - plane p/ 4 s x Figure 4: Poles of a second-order Butterworth filter (N =2,= 1). 4 Example Consider a second-order (N =2, = 1) Butterworth filter with cut-off frequency Ωc. Equation (9) generates a pole pair as shown shown in the pole-zero of Fig. 4. The transfer function is that of a second-order system with damping ratio ζ =0.707 and undamped natural frequency Ωn =Ωc. Show that the power frequency response satisfies the Butterworth specification of Eq. (4). The transfer function for the second-order filter is 2 Ωc H(s)= √ (10) s2 s 2 + 2Ωc +Ωc and the frequency response is 2 Ωc √ H(jΩ) = H(s)|s =jΩ = (11) 2 − 2 j (Ωc Ω )+ 2ΩcΩ 4 2 Ωc 1 | H(jΩ)| = H(jΩ)H(jΩ) = = (12) 4 4 / 4 (Ωc +Ω ) 1+(Ω Ωc) which is of the Butterworth form of Eq. (4). Example Design a Butterworth low-pass filter to meet the power gain specifications shown in Fig. 5. Comparing Figs. 1 and 5 1 =0.9 −→ =0.3333 1+ 2 1 =0.05 −→ λ =4.358 1+ λ2 Then log(λ/) N ≥ =3.70 log(Ωr/Ωc) we therefore select N=4. The 4 poles (p1 ...p4) lie on a circle of radius r = −1/N Ωc =13.16 and are given by | pn| =13.16 pn = π(2n +3)/8 for n =1... 4, giving a pair of complex conjugate pole pairs p1,4 = − 5.04 ± j12.16 p2,3 = − 12.16 ± j5.04 5 |H( jW ) |2 1 0. 9 0. 0 5 0 0 10 20 W pa s s ba n d tran s ition ba n d stop ba n d Figure 5: Filter specification for Butterworth design. Bode Diagram 50 0 −50 Magnitude (dB) −100 −1500 −90 −180 Phase (deg) −270 −360 −1 0 1 2 3 10 10 10 10 10 Frequency (rad/sec) Figure 6: Bode plots for fourth-order Butterworth design example. 6 The transfer function, normalized to unity gain, is 29993 H(s) = (s2 + 10:07s + 173:2)(s2 + 24:32s + 173:2) and the filter Bode plots are shown in Fig. 6. 2.2 Chebyshev Filters The order of a filter required to met a low-pass specification may often be reduced by relaxing the requirement of a monotonically decreasing power gain with frequency, and allowing “ripple” to occur in either the pass-band or the stop-band. The Chebyshev filters allow these conditions: 2 1 Type 1 jH(jΩ)j = 2 2 (13) 1 + ² TN (Ω=Ωc) 1 Type 2 jH(jΩ)j2 = (14) 2 2 2 1 + ² (TN (Ωr=Ωc)=TN (Ωr=Ω)) Where TN (x) is the Chebyshev polynomial of degree N. Note the similarity of the form of the Type 1 power gain (Eq. (13)) to that of the Butterworth filter, where the function N TN (Ω=Ωc) has replaced (Ω=Ωc) . The Type 1 filter produces an all-pole design with slightly different pole placement from the Butterworth filters, allowing resonant peaks in the pass band to introduce ripple, while the Type 2 filter introduces a set of zeros on the imaginary axis above Ωr, causing a ripple in the stop-band. The Chebyshev polynomials are defined recursively as follows T0(x) = 1 T1(x) = x 2 T2(x) = 2x ¡ 1 3 T3(x) = 4x ¡ 3x . TN (x) = 2xTN¡1(x) ¡ TN¡2(x); N > 1 (15) with alternate definitions ¡1 TN (x) = cos(N cos (x)) (16) ¡1 = cosh(N cosh (x)) (17) The Chebyshev polynomials have the min-max property: Of all polynomials of degree N with leading coe±cient equal to one, the polynomial N¡1 TN (x)=2 has the smallest magnitude in the interval jxj · 1. This \minimum maximum" amplitude is 21¡N . 7 In low-pass filters given by Eqs. (13) and (14), this property translates to the following characteristics: Filter Pass-Band Characteristic Stop-Band Characteristic Butterworth Maximally flat Maximally flat 2 Chebyshev Type 1 Ripple between 1 and 1=(1 + ² ) Maximally flat 2 Chebyshev Type 2 Maximally flat Ripple between 1 and 1=(1 + ¸ ) 2.2.1 The Chebyshev Type 1 Filter With the power response from Eq. (13) 2 1 jH(jΩ)j = 2 2 1 + ² TN (Ω=Ωc) and the filter specification from Fig. 1, the required filter order may be found as follows. At the edge of the stop-band Ω = Ωr 1 1 jH(jΩ j2 = · r 2 2 2 1 + ² TN (Ωr=Ωc) 1 + ¸ so that ³ ´ ¡1 ¸ · ²TN (Ωr=Ωc) = ² cosh N cosh (Ωr=Ωc) and solving for N cosh¡1 (¸=²) N ¸ ¡1 (18) cosh (Ωr=Ωc) The characteristic equation of the power transfer function is à ! à ! 2 2 s s j 1 + ² TN = 0 or TN = § jΩc jΩc ² ¡1 Now TN (x) = cos(N cos (x)), so that à à !! ¡1 s j cos N cos = § (19) jΩc ² ¡1 If we write cos (s=jΩc) = γ + j®, then s = Ωc (j cos (γ + j®)) = Ωc (sinh ® sin γ + j cosh ® cos γ) (20) which defines an ellipse of width 2Ωc sinh(®) and height 2Ωc cosh(®) in the s-plane.