ECE 410 DIGITAL SIGNAL PROCESSING D. Munson University of Illinois Chapter 12
|
|
- Marjorie Newton
- 5 years ago
- Views:
Transcription
1 . ECE 40 DIGITAL SIGNAL PROCESSING D. Munson University of Illinois Chapter IIR Filter Design ) Based on Analog Prototype a) Impulse invariant design b) Bilinear transformation ( ) ~ widely used ) Computer-Aided Optimization Designs in category ) proceed by first designing an analog filter having a frequency response with the desired shape, and then transforming it to a digital filter. To use these design methods, we must first learn just a bit about analog filter design. Elements of Analog Filter Design Notation: H L (s) h a (t) e -st dt (Laplace transform) H a (Ω) H L (jω) (Fourier transform) Consider only lowpass Butterworth, Chebyshev, and elliptic (Cauer) filters. For each of these types of filters, H L (s) is found indirectly from a specified H a (Ω). We need H L (s) because, later, this is what will be transformed into H(z). For Butterworth, Chebyshev, and elliptic filters, H a (Ω) has the form: H a (Ω) M(Ω ) + F(Ω ) rational with real coeffs.
2 . How do we get H L (s) from H a (Ω)? That is, how do we find H L (s) satisfying Answer: H L (jω) H a (Ω) M(Ω )? (*) ) First find poles and zeros of M( s ) where s is a complex variable. Since M has real coefficients and is a function of s, the poles and zeros will have symmetry around both the real and imaginary axii. ) Take H L (s) to be the left-half-plane pole factors (for stability) and left-half-plane zero factors (for smallest delay, called minimum phase ). But, does this work? Need to show H L (jω) M(Ω ) () Have: H L (s) H L ( s) M( s ) RHP factors LHP factors which implies So, () will be true if H L (jω) H L ( jω) M(Ω ) H L ( jω) H (jω) () L This follows, though, because the poles and zeros of H L (s) are symmetric around the real axis, and therefore occur in complex-conjugate pairs. For any pole pair or zero pair (s p) (s p*) in H L (s), we have (s p) (s p ) s jω ( jω p) ( jω p*) [(jω p*) (jω p)]* (s p)(s p ) s jω which proves (), and therefore ().
3 .3 Example Find H L (s). H a (Ω) M(Ω ) +Ω We have M( s ) s ( s) ( + s) LHP factor H L (s) s + Let s check to see if H a (Ω) above M(Ω ): H a (Ω) H L (jω) jω+ jω+ Ω + M(Ω ) Example Suppose M(Ω ) +Ω + Ω 4 (not for a B, C, or E filter!) Then M(-s ) s +s 4. ( s)( + s) (s +γ)(s+γ*)(s γ)(s γ*)
4 .4 γ e jπ 4 + j Pole-zero diagram: Take H L (s) + s (s +γ)(s+γ ) (LHP factors) Can check that H L (jω) M(Ω ). M( Ω ) for B, C, and E Filters Butterworth M(Ω ) Result: + F(Ω ) with F(Ω ) Ω n for nth-order filter M(Ω ) H a (Ω).0.5 monotone ~ no ripple.0 Ω M(Ω ).5 at Ω Ω c.0. For a general cutoff frequency Ω c, use F(Ω ) Ω Ωc n.
5 .5 Here, we are defining the cutoff frequency to be the value of Ω where H a (Ω) is reduced to onehalf its maximum height, or correspondingly, H a (Ω) reaches times its maximum value. This definition of cutoff frequency is common, particularly for smooth frequency responses that contain little or no ripple. Optimality: This M(Ω ) has maximum # of derivatives 0 at the origin for its order. Thus, the response is very flat across lower frequencies. Can show poles of M( s ) H L (s) H L ( s) are equally spaced on the unit circle. This fact helps in factoring M( s ). Chebyshev F(Ω ) ε C n (Ω) where ε is a real constant chosen by the designer and C n ( ) is the nth-order Chebyshev polynomial: C n (Ω) ( ) Ω ( ) Ω > cos n cos (Ω) cosh ncosh (Ω) with cosh t et + e t Can show: C 0 (Ω), C (Ω) Ω, C (Ω) Ω, and that there is a recursion relation: C n+ (Ω) Ω C n (Ω) C n (Ω) Result:
6 .6 M(Ω ) H a (Ω).0 + ε.0 Ω For this type of filter, the cutoff frequency is defined to be the value of Ω where H a (Ω) first drops below /(+ ε ) or, correspondingly, H a (Ω) first drops below +ε. This is a Type Chebyshev filter. Its response is equiripple in the passband and monotone decreasing in the stopband. It has a narrower transition band than a Butterworth filter. Tradeoff: Smaller ε gives smaller passband ripple but a wider transition band. Poles of M( s ) H L (s) H L ( s) lie on an ellipse. There is also a Type Chebyshev filter that has a monotone response in its passband and equiripple behavior in its stopband. Elliptic F(Ω ) ε J n (Ω) where J n is the Jacobi elliptic function. The defining formulas for J n are so cumbersome that they are not presented here. Result:
7 .7 M(Ω ) H a (Ω).0 δ p Ω p Ω s δ s Ω The response is equiripple in both the passband and stopband. Elliptic filters are optimal in the sense that for a given n, δ p, δ s, Ω p, the transition bandwidth Ω s - Ω p is the smallest possible. H a (Ω) for B, C, and E filters is reasonably linear till you get near the edge of the passband, where it can be quite nonlinear. The phase response is closest to linear for B, then C. Elliptic is worst. All-pass filters are sometimes cascaded onto elliptic filters to compensate for the nonlinear phase of elliptic filters. All-pass filters have H a (Ω) constant and the coefficients are chosen to shape H a (Ω) in a desired way. Bilinear Transformation Start with analog prototype H L (s). Take H(z) H L (s) sα z +z α is a real positive constant that we will be able to choose to control the cutoff frequency of the digital filter.
8 .8 s α z is a bilinear transformation (BLT) of the z-plane to the s-plane. For + z example, the point z o maps to the point s o as shown below. Im(z) z-plane α Im(s) s-plane z o π 4 so jα tan π 8 Re(z) Re(s) z o e jπ 4 To see this, note: s o α z o + z o α e jπ 4 + e jπ 4 α e jπ 8 e jπ 8 e j π 8 e jπ 8 e jπ 8 + e jπ 8 α jsin π 8 cos π 8 j α tan π 8 π So, if we design a digital filter using the BLT, then H(z o ) H d 4 H L (s o ) H a α tan π 8. will have the same value as We must have a much broader understanding than this, however. Questions: i) Stable H L (s) stable H(z)? ii) How is H d (ω) H( e jω )related to H L (s)? Answer i) by considering a point s s o and determining what z it gets mapped to. BLT mapping is s α z + z s ( + z ) α ( z )
9 .9 z ( s +α) α s z α+s α s So, a point s o σ o + j Ω o. gets mapped to: z o α+σ o + jω o α σ o jω o (*) i.e., H(z o ) H L (s o ). From (*) z o ( ) +Ω o ( ) +Ω o α+σ o α σ o / So: < σ o < 0 z o σ ο 0 > σ ο > 0 giving: Note: a) Left-half s-plane is mapped inside the unit circle in the z-plane. b) Right-half s-plane is mapped outside the unit circle in the z-plane. c) Imaginary axis in s-plane is mapped onto the unit circle in the z-plane. [ a) stable H L (s) results in stable H(z) ] [ c) H d (ω) H(e jω ) depends only on H L (jω)h a (Ω)] What is the relationship between H d (ω) and H a (Ω)? H d (ω) is given by
10 .0 Note: H d (ω) H(e jω ) H L (s) sα e jω +e jω α e jω + e jω α e jω/ e j ω/ e jω / e jω / e j ω/ + e j ω / α jsin ω cos ω j α tan ω So, H d (ω) H L (s) s jα tan ω Since H a (Ω) H L (jω) we have H d (ω) H a αtan ω ( ) This equation tells us exactly how, when using the bilinear transformation design method, the frequency response of the designed digital filter will depend on the frequency response of the analog prototype.
11 . Pictorial interpretation of ( ): α tan ω Ω o H d (ω o ) H a (Ω o ) -π ω o π ω So, H d (ω) takes on exactly the same set of values as H a (Ω), but there is a squashing of the analog frequency axis, according to the above curve. This mapping is nonlinear, and it has to be, since the infinite-length analog frequency axis < Ω < is mapped onto the finite-length interval π < ω < π. Because of this, H d (ω) won't look like H a (Ω) in general. It will be a warped version. This effect is sometimes called frequency warping. Example Applying the BLT to an analog prototype having frequency response: H a (Ω) Ω
12 . results in: H d (ω) π ω In general, to design H d (ω) having a desired shape, we would need to design H a (Ω) so that after application of the BLT, the frequency warping produces the desired H d (ω). Thus, we would need to prewarp. Let the desired H d (ω) be D(ω). ( ) want H a αtan ω H a (Ω) D tan Ω α D(ω) ( ) If we could design H L (s) to satisfy ( ), the BLT, would then give H d (ω) D(ω). This is problematic, however. The design of H L with a general shape will require computer optimization. Thus, we may as well design H d directly, using computer optimization. Fortunately, for LPF s, BPF s and HPF s, there is no such problem. For these kinds of filters, frequency warping just affects the cutoff frequencies, e.g., H a (Ω) H d (ω) BLT From ( ): Ω c Ω π ω π ω c Ω c αtan ω c ( )
13 .3 So, we can simply pick Ω c to give the desired ω c. Equivalently, if H a (Ω) is normalized to Ω c, we can choose α in the BLT to give the desired ω c. This suggests two alternative, but equivalent, design procedures:. a) Choose α arbitrarily, say α. b) Then design the analog prototype to have cutoff Ω c given by ( ). c) Apply the BLT.. a) Use an analog prototype with Ω c. b) Choose α to give the desired ω c. From ( ), α Ω c cot ω c so that the BLT method becomes H(z) H L (s) s Ω c cot ω c z +z Analog prototype filters are usually designed with normalized cutoff frequencies Ω c.0. In this case, the bilinear transformation method of design reduces to H(z) H L (s) s cot ω c z +z (BLT) In this course, we will use the second option, with Eq. (BLT), to perform the bilinear transformation method of design. Example
14 .4 Design a first-order Butterworth digital filter with ω c π 4, using the BLT method Solution: First, find the analog prototype: H a (Ω) B(Ω ) +Ω B( s ) s s + s Pole locations: Take LHP factor for H L (s): H L (s) s + Apply (BLT): H(z) H(z) s + scot π z 8 +z.44 z z + z.44 ( z ) + + z + z z
15 z.44 z Now, since H(z) resulted from a BLT design using an analog prototype filter with H a (0) and H d ( ) 0, we know H d (0) H a (0) and H a (π) H a ( ) 0 Furthermore, since the digital cutoff is π, we should have 4 π H d 4 Thus, H d (ω) should look like H d (ω) π 4 The correctness of the values of H d (ω) at ω 0, π can be easily verified from the transfer function H(z). Note that H d (0) H(e j0 ) H() H d (π) H(e jπ ) H( ) For ω 0 or π we can always evaluate the magnitude response via the lengthier calculation H d (ω) e jω.44 e jω π ω ( cosω) + (.99sinω) (.44 cosω) + (.44sinω)
16 .6 Finally, a plot of H d (ω) would show that the phase is nearly linear for ω < π 4 and becomes more nonlinear for ω π 4 and larger. Example Similar to previous example, but now design a second-order Butterworth digital filter with ω c π 4. To find analog prototype, note H a (Ω) M(Ω ) +Ω 4 B( s ) + s 4 e j 3π 4 j π e 4 e j 3π 4 e j π 4 Left-half-plane poles are used for H L (s): H L (s) s e j 3π 4 s e j3π 4 s + s + H(z) H L (s) scot π 8 z +z ( ) ( ) (.44) z z z + z + ( + z ) (.44) ( z + z ) +.44 ( z ) + ( + z + z )
17 .7 + z + z z + 3.4z z z z z Now, once again, from the shape of the analog prototype H a (Ω) we know π H d (0), H d 4, H d(π) 0 However, for the second-order filter, the frequency response makes a sharper transition around ω c π and looks like H d (ω) π 4 π ω The phase H d (ω) will again be quite linear across the passband and more nonlinear across the stopband. Let s try to get a better feel for the BLT mapping by considering one further example. Example Suppose that an analog prototype filter H L (s) has the frequency response H a (Ω) 3 3 Ω
18 .8 and that the bilinear transformation is used to produce a digital filter with H(z) H L (s) z s. Sketch H d (ω). Label all critical frequencies and amplitudes. +z Solution H d (ω) is a squashed version of H a (Ω), described by H d (ω) H a αtan ω H a tan ω Thus, the value of H d (ω) equals the value of H a (Ω) at Ω tan ω. So a feature (e.g., jump) in H a (Ω) that occurs at Ω Ω 0 will occur in H d (ω) at ω 0 tan Ω 0 The interesting features in H a (Ω) occur at Ω 0,, and 3. In addition let s also consider Ω 0.5. The corresponding ω 0 are given in the table below. Ω ω π 0.66 π π π Thus, H d (ω) looks like H d (ω) π.795 π.66 π.705 π π π.66 π.705 π.795 π π ω Notice that H d (ω) takes on exactly the same set of values taken on by H a (Ω). It just does so at different frequencies, which causes a change in shape.. Computer-Aided Optimization
19 .9 Used for IIR designs with general magnitude and/or phase specifications (i.e., H d (ω) not LPF, HPF, or BPF). Digital Frequency Transformation By making a change of variable in H L (s) or H(z), we can transform an analog or digital LPF to a LPF having a different cutoff frequency, or to a HPF or BPF. We will consider only digital transformations (transformations of H(z)). Procedure: Let H(z) be the transfer function of a low-pass digital filter. Then simply substitute for z in H(z), using the expressions below, to produce filters having the described characteristics. Lowpass Lowpass where z z β β z ω c cutoff of new filter [ ] [ ] β sin (ω c ω c )/ sin (ω c +ω c )/ Lowpass Highpass z z β βz ω c cutoff of new filter [ ] [ ] β cos (ω c + ω c )/ cos (ω c ω c )/ Lowpass Bandpass z z β z +β β z β z +
20 .0 ω lower cutoff of BPF ω u upper cutoff of BPF β γ K/(K + ) β (K )/(K + ) [ ] [ ] γ cos (ω u +ω )/ cos (ω u ω )/ K cot ω u ω tan ω c Relationship Between Pole and Zero Locations and Frequency Response In digital filter design we choose the coefficients of H(z) a 0 + a z + + a N z N + b z + + b N z N to shape the frequency response H d (ω) H(e jω ) in a desired way. Since H(z) can also be written in terms of its poles and zeros as H(z) a 0 N i z z i z p i this provides an alternative parameterization of H(z). Choosing the pole and zero locations of the filter is basically equivalent to choosing the {a i } and {b i }. In this connection, it is worth exploring how the pole and zero locations affect the shape of H d (ω). The general shape of H d (ω) often can be visualized from knowledge of the pole and zero locations of H(z). This is especially true for situations where poles are near the unit circle and zeros are either on or near the unit circle. Example Consider a causal, stable all-pole filter with H(z) z αcosω ( p )z +α z αe jω p ( ) z αe jω p ( )
21 . and 0 < α <. H(z) approaches at the pole locations and approaches zero for z large. Thus, we expect H(z) to look somewhat like a two-pole circus tent: Im(z) ω p ω p Re(z) Here, the pole locations are marked by and are at a distance α from the origin. Since H(z) is causal, its ROC is {z : z > α} and the above tent covers only this set of z. (The tent is undefined elsewhere.) Since α <, ROC H includes the unit circle. Thus, the frequency response H d (ω) H(e jω ) is well defined and is a circular slice of the circus tent, around the unit circle. Now, since H(z) is infinite at z α e ±jω p, we expect that H(z) will be large for z near the poles. If α is nearly one then the poles are close to the unit circle and H d (ω) will be large for ω such that e jω is close to the poles. This suggests that H d (ω) will look like H d (ω) where the peaks occur near ω ω p. Example ω p ω p ω
22 . H(z) z ( cosω p )z + z αcosω ( p )z +α z e jω p ( ) z e jω p ( ) z αe jω p ( ) z αe jω p ( ) As before, H(z) approaches at the pole locations z α e ±jω p. H(z) 0 at the zero locations z e ±jω p. Thus, H(z) is similar to the previous two-pole circus tent except it touches the ground at the zero locations z e ±jω p as shown below, where poles are indicated by and zeros are represented by. Im(z) ω p ω p Re(z) Since the pole and zero locations are close together, it may seem difficult to determine what H d (ω) {H(z) for z on the unit circle} will look like. A simple observation helps greatly, though. Suppose so that H(z) ( z z ) z z ( z p ) z p ( ) ( ) ( ) H d (ω) H(e jω ) ejω z e jω z e jω p e jω p Then, the value of H d (ω) at some specific frequency ω* is
23 .3 e jω z e jω z e jω p e jω p. The factor e jω z i is the distance in the complex plane between e jω* and z i. Likewise, the factor e jω p i is the distance in the complex plane between e jω* and p i. Thus, H d (ω ) is the product of the distances between e jω* and z, and between e jω* and z, divided by the product of the distances between e jω* and p, and between e jω* and p. In our example, these distances can be visualized for three different frequencies ω by examining the figure below. Im(z) z p e jω * e jω 3 * ω p ω p e jω * Re(z) p z We see that the distances from e jω * to the zero z at e jω p and to the neighboring pole p are nearly the same. Likewise, the distances from e jω * to the zero z at e jω p and its neighboring pole p are nearly the same. Thus, H d ω ( ) e jω z e jω p e jω z e jω p ()() The same result holds for any ω that is close to zero. Similarly, in the case of e jω 3* we see that distances to all poles and zeros are nearly equal. Thus. H d ( ω 3 ),
24 .4 which holds for any ω 3 roughly satisfying π ω 3 π. The distance between e jω * and the zero at e jω p approaches zero as ω ω p. Thus H d ( ω ) 0 for ω close enough to ω p. These considerations lead to the H d (ω) sketched below. H d (ω) π ω p ω p π ω Here, we can find the precise values of H d (0) and H d (π) by using H d (0) H() cosω p +α αcosω p H d (π) H( ) + cosω p +α + αcosω p. The above frequency response is that of a crude notch filter where signal components near ω ω p are greatly attenuated and signal components at other frequencies are passed with nearly unit amplitude. The nulls (notches) in H d (ω) at ω ω p are caused by the zeros of H(z) at z e ±jω p. It is only slightly more difficult to gain a rough idea of what H d (ω) will look like. From ( ) we have H d (ω) (e jω z ) + (e jω z ) (e jω p ) (e jω p ) Each term (e jω z i ) or (e jω p i ) is a vector in the complex plane. (e jω z i ) and (e jω p i ) are simply the angles of these vectors with respect to the positive real axis. This interpretation can be helpful when trying to visualize H d (ω) for low-order filters. In general, however, Matlab should be used to plot both H d (ω) and H d (ω) for higher-order filters Example
25 .5 For an FIR lowpass filter, use Matlab to find the zero locations of H(z). (All poles are at z 0.) You will find that zeros in H d (ω) within the stopband are caused by zeros of H(z) on the unit circle. Other zeros of H(z) are strategically placed off the unit circle to give a flat response in the passband. Example For a Butterworth lowpass filter you will find that some poles are located near the unit circle e jω for ω corresponding to the cutoff frequency. Example For an elliptic lowpass filter you will find poles near e jω for ω cutoff frequency, and some zeros on the unit circle at locations corresponding to the stopband.
DIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS. 3.6 Design of Digital Filter using Digital to Digital
DIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS Contents: 3.1 Introduction IIR Filters 3.2 Transformation Function Derivation 3.3 Review of Analog IIR Filters 3.3.1 Butterworth
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #24 Tuesday, November 4, 2003 6.8 IIR Filter Design Properties of IIR Filters: IIR filters may be unstable Causal IIR filters with rational system
More informationDigital Control & Digital Filters. Lectures 21 & 22
Digital Controls & Digital Filters Lectures 2 & 22, Professor Department of Electrical and Computer Engineering Colorado State University Spring 205 Review of Analog Filters-Cont. Types of Analog Filters:
More informationDigital Signal Processing IIR Filter Design via Bilinear Transform
Digital Signal Processing IIR Filter Design via Bilinear Transform D. Richard Brown III D. Richard Brown III 1 / 12 Basic Procedure We assume here that we ve already decided to use an IIR filter. The basic
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 05 IIR Design 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationPS403 - Digital Signal processing
PS403 - Digital Signal processing 6. DSP - Recursive (IIR) Digital Filters Key Text: Digital Signal Processing with Computer Applications (2 nd Ed.) Paul A Lynn and Wolfgang Fuerst, (Publisher: John Wiley
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA November 1, 2009 1 / 27 1 The z-transform 2 Linear Time-Invariant System 3 Filter Design IIR Filters FIR Filters
More informationFilter Analysis and Design
Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved
More informationDesign of IIR filters
Design of IIR filters Standard methods of design of digital infinite impulse response (IIR) filters usually consist of three steps, namely: 1 design of a continuous-time (CT) prototype low-pass filter;
More informationDIGITAL SIGNAL PROCESSING. Chapter 6 IIR Filter Design
DIGITAL SIGNAL PROCESSING Chapter 6 IIR Filter Design OER Digital Signal Processing by Dr. Norizam Sulaiman work is under licensed Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International
More informationFilters and Tuned Amplifiers
Filters and Tuned Amplifiers Essential building block in many systems, particularly in communication and instrumentation systems Typically implemented in one of three technologies: passive LC filters,
More informationLike bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform.
Inversion of the z-transform Focus on rational z-transform of z 1. Apply partial fraction expansion. Like bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform. Let X(z)
More informationMITOCW watch?v=jtj3v Rx7E
MITOCW watch?v=jtj3v Rx7E The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free. To
More informationDigital Signal Processing Lecture 8 - Filter Design - IIR
Digital Signal Processing - Filter Design - IIR Electrical Engineering and Computer Science University of Tennessee, Knoxville October 20, 2015 Overview 1 2 3 4 5 6 Roadmap Discrete-time signals and systems
More informationLecture 7 Discrete Systems
Lecture 7 Discrete Systems EE 52: Instrumentation and Measurements Lecture Notes Update on November, 29 Aly El-Osery, Electrical Engineering Dept., New Mexico Tech 7. Contents The z-transform 2 Linear
More informationLecture 9 Infinite Impulse Response Filters
Lecture 9 Infinite Impulse Response Filters Outline 9 Infinite Impulse Response Filters 9 First-Order Low-Pass Filter 93 IIR Filter Design 5 93 CT Butterworth filter design 5 93 Bilinear transform 7 9
More informationLecture 16: Filter Design: Impulse Invariance and Bilinear Transform
EE58 Digital Signal Processing University of Washington Autumn 2 Dept. of Electrical Engineering Lecture 6: Filter Design: Impulse Invariance and Bilinear Transform Nov 26, 2 Prof: J. Bilmes
More informationINFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters 4. THE BUTTERWORTH ANALOG FILTER
INFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters. INTRODUCTION 2. IIR FILTER DESIGN 3. ANALOG FILTERS 4. THE BUTTERWORTH ANALOG FILTER 5. THE CHEBYSHEV-I
More informationChapter 7: IIR Filter Design Techniques
IUST-EE Chapter 7: IIR Filter Design Techniques Contents Performance Specifications Pole-Zero Placement Method Impulse Invariant Method Bilinear Transformation Classical Analog Filters DSP-Shokouhi Advantages
More informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More informationLecture 3 - Design of Digital Filters
Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can
More informationECE503: Digital Signal Processing Lecture 5
ECE53: Digital Signal Processing Lecture 5 D. Richard Brown III WPI 3-February-22 WPI D. Richard Brown III 3-February-22 / 32 Lecture 5 Topics. Magnitude and phase characterization of transfer functions
More informationUNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter?
UNIT - III PART A. Mention the important features of the IIR filters? i) The physically realizable IIR filters does not have linear phase. ii) The IIR filter specification includes the desired characteristics
More informationSignals and Systems. Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI
Signals and Systems Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE LONDON Effect on poles and zeros on frequency response
More informationDigital Signal Processing
Digital Signal Proceing IIR Filter Deign Manar Mohaien Office: F8 Email: manar.ubhi@kut.ac.kr School of IT Engineering Review of the Precedent Lecture Propertie of FIR Filter Application of FIR Filter
More informationElectronic Circuits EE359A
Electronic Circuits EE359A Bruce McNair B26 bmcnair@stevens.edu 21-216-5549 Lecture 22 578 Second order LCR resonator-poles V o I 1 1 = = Y 1 1 + sc + sl R s = C 2 s 1 s + + CR LC s = C 2 sω 2 s + + ω
More informationDigital Signal Processing Lecture 9 - Design of Digital Filters - FIR
Digital Signal Processing - Design of Digital Filters - FIR Electrical Engineering and Computer Science University of Tennessee, Knoxville November 3, 2015 Overview 1 2 3 4 Roadmap Introduction Discrete-time
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even
More informationElectronic Circuits EE359A
Electronic Circuits EE359A Bruce McNair B26 bmcnair@stevens.edu 21-216-5549 Lecture 22 569 Second order section Ts () = s as + as+ a 2 2 1 ω + s+ ω Q 2 2 ω 1 p, p = ± 1 Q 4 Q 1 2 2 57 Second order section
More informationAn Iir-Filter Example: A Butterworth Filter
An Iir-Filter Example: A Butterworth Filter Josef Goette Bern University of Applied Sciences, Biel Institute of Human Centered Engineering - microlab JosefGoette@bfhch February 7, 2017 Contents 1 Introduction
More informationINF3440/INF4440. Design of digital filters
Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 INF3440/INF4440. Design of digital filters October 2004 Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 Last lectures:
More information1 1.27z z 2. 1 z H 2
E481 Digital Signal Processing Exam Date: Thursday -1-1 16:15 18:45 Final Exam - Solutions Dan Ellis 1. (a) In this direct-form II second-order-section filter, the first stage has
More informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
More informationResponses of Digital Filters Chapter Intended Learning Outcomes:
Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing
More informationChapter 7: Filter Design 7.1 Practical Filter Terminology
hapter 7: Filter Design 7. Practical Filter Terminology Analog and digital filters and their designs constitute one of the major emphasis areas in signal processing and communication systems. This is due
More informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture : Design of Digital IIR Filters (Part I) Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner (Univ.
More informationV. IIR Digital Filters
Digital Signal Processing 5 March 5, V. IIR Digital Filters (Deleted in 7 Syllabus). (dded in 7 Syllabus). 7 Syllabus: nalog filter approximations Butterworth and Chebyshev, Design of IIR digital filters
More informationFilter Design Problem
Filter Design Problem Design of frequency-selective filters usually starts with a specification of their frequency response function. Practical filters have passband and stopband ripples, while exhibiting
More informationDSP-CIS. Chapter-4: FIR & IIR Filter Design. Marc Moonen
DSP-CIS Chapter-4: FIR & IIR Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ PART-II : Filter Design/Realization Step-1 : Define
More informationAll-Pole Recursive Digital Filters Design Based on Ultraspherical Polynomials
RADIOENGINEERING, VOL. 3, NO. 3, SEPTEMBER 4 949 All-Pole Recursive Digital Filters Design Based on Ultraspherical Polynomials Nikola STOJANOVIĆ, Negovan STAMENKOVIĆ, Vidosav STOJANOVIĆ University of Niš,
More informationUNIVERSITY OF OSLO. Faculty of mathematics and natural sciences. Forslag til fasit, versjon-01: Problem 1 Signals and systems.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 1th, 016 Examination hours: 14:30 18.30 This problem set
More informationCoefficients of Recursive Linear Time-Invariant First-Order Low-Pass and High-Pass Filters (v0.1)
Coefficients of Recursive Linear Time-Invariant First-Order Low-Pass and High-Pass Filters (v0. Cliff Sparks www.arpchord.com The following is a quick overview of recursive linear time-invariant first-order
More informationDHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A
DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A Classification of systems : Continuous and Discrete
More informationLecture 2 OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE
OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE EEE 43 DIGITAL SIGNAL PROCESSING (DSP) 2 DIFFERENCE EQUATIONS AND THE Z- TRANSFORM FALL 22 Yrd. Doç. Dr. Didem Kivanc Tureli didemk@ieee.org didem.kivanc@okan.edu.tr
More informationQuadrature-Mirror Filter Bank
Quadrature-Mirror Filter Bank In many applications, a discrete-time signal x[n] is split into a number of subband signals { v k [ n]} by means of an analysis filter bank The subband signals are then processed
More informationAPPLIED SIGNAL PROCESSING
APPLIED SIGNAL PROCESSING DIGITAL FILTERS Digital filters are discrete-time linear systems { x[n] } G { y[n] } Impulse response: y[n] = h[0]x[n] + h[1]x[n 1] + 2 DIGITAL FILTER TYPES FIR (Finite Impulse
More informationThere are two main classes of digital lter FIR (Finite impulse response) and IIR (innite impulse reponse).
FIR Filters I There are two main classes of digital lter FIR (Finite impulse response) and IIR (innite impulse reponse). FIR Digital Filters These are described by dierence equations of the type: y[n]
More informationBasic Design Approaches
(Classic) IIR filter design: Basic Design Approaches. Convert the digital filter specifications into an analog prototype lowpass filter specifications. Determine the analog lowpass filter transfer function
More informationLINEAR-PHASE FIR FILTERS DESIGN
LINEAR-PHASE FIR FILTERS DESIGN Prof. Siripong Potisuk inimum-phase Filters A digital filter is a minimum-phase filter if and only if all of its zeros lie inside or on the unit circle; otherwise, it is
More informationButterworth Filter Properties
OpenStax-CNX module: m693 Butterworth Filter Properties C. Sidney Burrus This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 3. This section develops the properties
More informationOptimal Design of Real and Complex Minimum Phase Digital FIR Filters
Optimal Design of Real and Complex Minimum Phase Digital FIR Filters Niranjan Damera-Venkata and Brian L. Evans Embedded Signal Processing Laboratory Dept. of Electrical and Computer Engineering The University
More informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set
More information-Digital Signal Processing- FIR Filter Design. Lecture May-16
-Digital Signal Processing- FIR Filter Design Lecture-17 24-May-16 FIR Filter Design! FIR filters can also be designed from a frequency response specification.! The equivalent sampled impulse response
More informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 3 Signals & Systems Prof. ark Fowler Note Set #28 D-T Systems: DT Filters Ideal & Practical /4 Ideal D-T Filters Just as in the CT case we can specify filters. We looked at the ideal filter for the
More informationThe Approximation Problem
EE 508 Lecture 3 The Approximation Problem Classical Approximating Functions - Thompson and Bessel Approximations Review from Last Time Elliptic Filters Can be thought of as an extension of the CC approach
More informationOptimum Ordering and Pole-Zero Pairing of the Cascade Form IIR. Digital Filter
Optimum Ordering and Pole-Zero Pairing of the Cascade Form IIR Digital Filter There are many possible cascade realiations of a higher order IIR transfer function obtained by different pole-ero pairings
More informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More informationLecture 14: Windowing
Lecture 14: Windowing ECE 401: Signal and Image Analysis University of Illinois 3/29/2017 1 DTFT Review 2 Windowing 3 Practical Windows Outline 1 DTFT Review 2 Windowing 3 Practical Windows DTFT Review
More informationTransform Analysis of Linear Time-Invariant Systems
Transform Analysis of Linear Time-Invariant Systems Discrete-Time Signal Processing Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-Sen University Kaohsiung, Taiwan ROC Transform
More informationFourier Series Representation of
Fourier Series Representation of Periodic Signals Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline The response of LIT system
More informationMultimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2
Multimedia Signals and Systems - Audio and Video Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Kunio Takaya Electrical and Computer Engineering University of Saskatchewan December
More informationProblem Set 9 Solutions
Problem Set 9 Solutions EE23: Digital Signal Processing. From Figure below, we see that the DTFT of the windowed sequence approaches the actual DTFT as the window size increases. Gibb s phenomenon is absent
More informationTime Series Analysis: 4. Linear filters. P. F. Góra
Time Series Analysis: 4. Linear filters P. F. Góra http://th-www.if.uj.edu.pl/zfs/gora/ 2012 Linear filters in the Fourier domain Filtering: Multiplying the transform by a transfer function. g n DFT G
More informationNotes on L (Optimal) Filters
Notes on L (Optimal) Filters by C. Bond, 20 Background In 959 A. Papoulis published a paper which completed the description of a new class of filters with optimal properties, originally reported in 958.
More informationDSP. Chapter-3 : Filter Design. Marc Moonen. Dept. E.E./ESAT-STADIUS, KU Leuven
DSP Chapter-3 : Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ Filter Design Process Step-1 : Define filter specs Pass-band, stop-band,
More informationAnalog and Digital Filter Design
Analog and Digital Filter Design by Jens Hee http://jenshee.dk October 208 Change log 28. september 208. Document started.. october 208. Figures added. 6. october 208. Bilinear transform chapter extended.
More informationDigital Filters Ying Sun
Digital Filters Ying Sun Digital filters Finite impulse response (FIR filter: h[n] has a finite numbers of terms. Infinite impulse response (IIR filter: h[n] has infinite numbers of terms. Causal filter:
More informationComputer Engineering 4TL4: Digital Signal Processing
Computer Engineering 4TL4: Digital Signal Processing Day Class Instructor: Dr. I. C. BRUCE Duration of Examination: 3 Hours McMaster University Final Examination December, 2003 This examination paper includes
More informationLABORATORY 3 FINITE IMPULSE RESPONSE FILTERS
LABORATORY 3 FINITE IMPULSE RESPONSE FILTERS 3.. Introduction A digital filter is a discrete system, used with the purpose of changing the amplitude and/or phase spectrum of a signal. The systems (filters)
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 2011) Final Examination December 19, 2011 Name: Kerberos Username: Please circle your section number: Section Time 2 11 am 1 pm 4 2 pm Grades will be determined by the correctness of your answers
More informationMatched Second Order Digital Filters
Matched Second Order Digital Filters Martin Vicanek 14. February 016 1 Introduction Second order sections are universal building blocks for digital filters. They are characterized by five coefficients,
More informationDesign of Stable IIR filters with prescribed flatness and approximately linear phase
Design of Stable IIR filters with prescribed flatness and approximately linear phase YASUNORI SUGITA Nagaoka University of Technology Dept. of Electrical Engineering Nagaoka city, Niigata-pref., JAPAN
More informationIIR digital filter design for low pass filter based on impulse invariance and bilinear transformation methods using butterworth analog filter
IIR digital filter design for low pass filter based on impulse invariance and bilinear transformation methods using butterworth analog filter Nasser M. Abbasi May 5, 0 compiled on hursday January, 07 at
More informationPerfect Reconstruction Two- Channel FIR Filter Banks
Perfect Reconstruction Two- Channel FIR Filter Banks A perfect reconstruction two-channel FIR filter bank with linear-phase FIR filters can be designed if the power-complementary requirement e jω + e jω
More informationIT DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A
DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING IT6502 - DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A 1. What is a continuous and discrete time signal? Continuous
More informationAnalysis of Finite Wordlength Effects
Analysis of Finite Wordlength Effects Ideally, the system parameters along with the signal variables have infinite precision taing any value between and In practice, they can tae only discrete values within
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 20) Final Examination December 9, 20 Name: Kerberos Username: Please circle your section number: Section Time 2 am pm 4 2 pm Grades will be determined by the correctness of your answers (explanations
More informationDISCRETE-TIME SIGNAL PROCESSING
THIRD EDITION DISCRETE-TIME SIGNAL PROCESSING ALAN V. OPPENHEIM MASSACHUSETTS INSTITUTE OF TECHNOLOGY RONALD W. SCHÄFER HEWLETT-PACKARD LABORATORIES Upper Saddle River Boston Columbus San Francisco New
More informationHomework Assignment 11
Homework Assignment Question State and then explain in 2 3 sentences, the advantage of switched capacitor filters compared to continuous-time active filters. (3 points) Continuous time filters use resistors
More information(Refer Slide Time: 01:28 03:51 min)
Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture 40 FIR Design by Windowing This is the 40 th lecture and our topic for
More informationUniversity of Illinois at Urbana-Champaign ECE 310: Digital Signal Processing
University of Illinois at Urbana-Champaign ECE 0: Digital Signal Processing Chandra Radhakrishnan PROBLEM SET : SOLUTIONS Peter Kairouz Problem. Hz z 7 z +/9, causal ROC z > contains the unit circle BIBO
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
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. I Reading:
More informationECSE 512 Digital Signal Processing I Fall 2010 FINAL EXAMINATION
FINAL EXAMINATION 9:00 am 12:00 pm, December 20, 2010 Duration: 180 minutes Examiner: Prof. M. Vu Assoc. Examiner: Prof. B. Champagne There are 6 questions for a total of 120 points. This is a closed book
More informationw n = c k v n k (1.226) w n = c k v n k + d k w n k (1.227) Clearly non-recursive filters are a special case of recursive filters where M=0.
Random Data 79 1.13 Digital Filters There are two fundamental types of digital filters Non-recursive N w n = c k v n k (1.226) k= N and recursive N M w n = c k v n k + d k w n k (1.227) k= N k=1 Clearly
More informationUNIVERSITI SAINS MALAYSIA. EEE 512/4 Advanced Digital Signal and Image Processing
-1- [EEE 512/4] UNIVERSITI SAINS MALAYSIA First Semester Examination 2013/2014 Academic Session December 2013 / January 2014 EEE 512/4 Advanced Digital Signal and Image Processing Duration : 3 hours Please
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
More informationDESIGN OF CMOS ANALOG INTEGRATED CIRCUITS
DESIGN OF CMOS ANALOG INEGRAED CIRCUIS Franco Maloberti Integrated Microsistems Laboratory University of Pavia Discrete ime Signal Processing F. Maloberti: Design of CMOS Analog Integrated Circuits Discrete
More informationOptimum Ordering and Pole-Zero Pairing. Optimum Ordering and Pole-Zero Pairing Consider the scaled cascade structure shown below
Pole-Zero Pairing of the Cascade Form IIR Digital Filter There are many possible cascade realiations of a higher order IIR transfer function obtained by different pole-ero pairings and ordering Each one
More informationThe Approximation Problem
EE 508 Lecture The Approximation Problem Classical Approximating Functions - Elliptic Approximations - Thompson and Bessel Approximations Review from Last Time Chebyshev Approximations T Type II Chebyshev
More informationDiscrete Time Systems
Discrete Time Systems Valentina Hubeika, Jan Černocký DCGM FIT BUT Brno, {ihubeika,cernocky}@fit.vutbr.cz 1 LTI systems In this course, we work only with linear and time-invariant systems. We talked about
More informationVel Tech High Tech Dr.Ranagarajan Dr.Sakunthala Engineering College Department of ECE
Subject Code: EC6502 Course Code:C302 Course Name: PRINCIPLES OF DIGITAL SIGNAL PROCESSING L-3 : T-1 : P-0 : Credits 4 COURSE OBJECTIVES: 1. To learn discrete Fourier transform and its properties 2. To
More informationMultirate signal processing
Multirate signal processing Discrete-time systems with different sampling rates at various parts of the system are called multirate systems. The need for such systems arises in many applications, including
More informationFilter structures ELEC-E5410
Filter structures ELEC-E5410 Contents FIR filter basics Ideal impulse responses Polyphase decomposition Fractional delay by polyphase structure Nyquist filters Half-band filters Gibbs phenomenon Discrete-time
More informationToday. ESE 531: Digital Signal Processing. IIR Filter Design. Impulse Invariance. Impulse Invariance. Impulse Invariance. ω < π.
Today ESE 53: Digital Signal Processing! IIR Filter Design " Lec 8: March 30, 207 IIR Filters and Adaptive Filters " Bilinear Transformation! Transformation of DT Filters! Adaptive Filters! LMS Algorithm
More informationSignals and Systems. Problem Set: The z-transform and DT Fourier Transform
Signals and Systems Problem Set: The z-transform and DT Fourier Transform Updated: October 9, 7 Problem Set Problem - Transfer functions in MATLAB A discrete-time, causal LTI system is described by the
More informationLet H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 )
Review: Poles and Zeros of Fractional Form Let H() = P()/Q() be the system function of a rational form. Let us represent both P() and Q() as polynomials of (not - ) Then Poles: the roots of Q()=0 Zeros:
More informationTime Series Analysis: 4. Digital Linear Filters. P. F. Góra
Time Series Analysis: 4. Digital Linear Filters P. F. Góra http://th-www.if.uj.edu.pl/zfs/gora/ 2018 Linear filters Filtering in Fourier domain is very easy: multiply the DFT of the input by a transfer
More informationSpeaker: Arthur Williams Chief Scientist Telebyte Inc. Thursday November 20 th 2008 INTRODUCTION TO ACTIVE AND PASSIVE ANALOG
INTRODUCTION TO ACTIVE AND PASSIVE ANALOG FILTER DESIGN INCLUDING SOME INTERESTING AND UNIQUE CONFIGURATIONS Speaker: Arthur Williams Chief Scientist Telebyte Inc. Thursday November 20 th 2008 TOPICS Introduction
More informationCh. 7: Z-transform Reading
c J. Fessler, June 9, 3, 6:3 (student version) 7. Ch. 7: Z-transform Definition Properties linearity / superposition time shift convolution: y[n] =h[n] x[n] Y (z) =H(z) X(z) Inverse z-transform by coefficient
More informationThe Approximation Problem
EE 508 Lecture The Approximation Problem Classical Approximating Functions - Elliptic Approximations - Thompson and Bessel Approximations Review from Last Time Chebyshev Approximations T Type II Chebyshev
More information