An Iir-Filter Example: A Butterworth Filter
|
|
- Todd James
- 6 years ago
- Views:
Transcription
1 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 1 2 Analog Butterworth Lowpass-Filters 4 3 Continuous-to-Discrete Transformations Impulse Invariance Transformation Bilinear Transformation 23 4 Discrete-Time Butterworth Filter Example Design Using Impulse-Invariance Transformation Design Using Bilinear Transformation 39 References Iir Butterworth i 2017
2 c Josef Goette, All rights reserved This work may not be translated or copied in whole or in part without the written permission by the author, except for brief excerpts in connection with reviews or scholarly analysis Use in connection with any form of information storage and retrieval, electronic adaptation, computer software is forbidden 33 Iir Butterworth ii 2017
3 1 Introduction In the present document, we use, as usual, the following notation: Continuous-time signal frequencies come in standard characters such as f for the frequency in Hertz ([Hz]), and ω for the radian frequency in [rad/sec], ω = 2πf The corresponding discrete-time signal frequencies use the corresponding symbols with hats: ˆf ˆ= f/fs = ft s and ˆω = ωt s = ω/f s with f s being the sampling frequency in [Hz] and T s the corresponding sampling interval T s = 1/f s in [sec] We abbreviate continuous-time by Ct and discrete-time by Dt Further abbreviations are C2d for continuous-to-discrete, Tf for transfer function, and Lhp for left-half plane We very closely follow [OS75] for the developments in the present document; see our additional remarks on details, which appear after the bibliography at the very end of the document on page Iir Butterworth
4 An Iir-Filter Design Approach Iir (recursive) filters are often design by what might be called analog prototyping design a normalized Ct lowpass prototype filter then: apply frequency-band transformations such as lowpass to lowpass lowpass to highpass lowpass to bandpass lowpass to bandstop then: apply a continuous-to-discrete transformation here: we work-out a Butterworth-filter example The acronym Iir stands for infinite impulse response; theoretically, the impulse response of this kind of filters never completely dies out Of course, stable Iir filters have an impulse response that approaches more and more the zero line as time passes by; if the filter runs on a digital computer with floatingpoint arithmetic, then the impulse response will eventually become smaller than the smallest positive representable number; if it runs on a digital computer with fixed-point arithmetic, then the impulse response might begin to oscillate without end see the discussion to the keyword limit cycles in our document [Goe17] Recall that the acronym Ct means continuous time Matlab supplies the following commands for frequencyband transformations: lp2lp() for the lowpass-to-lowpass transformation; lp2hp() for the lowpass-to-highpass transformation; 33 Iir Butterworth
5 lp2bp() for the lowpass-to-bandpass transformation; and the command lp2bs() for the lowpass-to-bandstop transformation, respectively Matlab has also commands for the continuous-to-discrete transformations, some of which stem from the Control System Toolbox (use help c2d), and some others from the Signal Processing Toolbox (use help impinvar and use help bilinear) 1 As the table of contents reveals, we discuss in the present document the design based on the impulse-invariance transformation and the design based on the bilinear transformation But also note that Matlab supplies higher level commands that integrate the mentioned lower level commands and thus considerably simplify the design of discrete-time and even digital filters 2 from the perspective of the user; you might also want to try out the graphical filter design toolfdatool (use help fdatool to obtain a short description, and type fdatool in the Matlab command window to run the Gui) 1 The mentioned toolboxes just have different names for the commands; the algorithms behind these names are, however, very similar or even identical 2 Recall that we name filters to be discrete time if they filter signals that have the time being discrete but the samples being still real numbers; digital filters filter then signals where not only the time is discrete, but where also the samples are represented be a finite number of bits In the realm of Matlab simulations, we often also call filters and signals to be discrete, if the samples of the signals are represented by floating-point numbers still finitely many bits, but ; and we retain the notion digital for signals with samples that are represented by fixed-point numbers 33 Iir Butterworth
6 2 Analog Butterworth Lowpass-Filters Butterworth Filter Properties magnitude response is maximally flat in passband for a N-th order lowpass, the first (2N 1) derivatives of the squared magnitude function are zero at ω = 0 the approximation to the ideal rectangular lowpass characteristic (brick-wall) is monotonic in passband as well as in stopband squared magnitude function H Ct (jω) 2 1 = ( jω 1 + jω c specified by just 2 parameters the filter order N the 3 db cutoff-frequency ω c ) 2N We also often call these analog filters continuous-time (Ct) filters 33 Iir Butterworth
7 Magnitude-Squared and Magnitude Responses H Ct (jω) 2 1 N = 1, 2, 4, 8 1/2 0 ω c ω H Ct (jω) 1 N = 1, 2, 4, 8 1/ 2 0 ω c ω Recall that the linear gain 1/ 2 corresponds to 3 db; we have 20 log 10 (1/ 2) = Butterworth filters have a frequency magnitude response with no ripple, neither in the passband nor in the stopband As the parameter N in the squared magnitude function on page 4 increases, the filter characteristics become sharper, meaning that if N increases, the characteristic in the passband stays closer to unity whereas the characteristic in the stopband comes close to zero more rapidly Other well known analog filters are the Chebyshev filters, the inverse Chebyshev filters, and the Cauer filters The Chebyshev 33 Iir Butterworth
8 filters have a frequency magnitude response with a ripple in the passband but no ripple in the stopband; inverse Chebyshev filters have a dual characteristic to that of the Chebyshev filters: they have no ripple in the passband but a ripple in the stopband; Cauer filters also called elliptic filters have a ripple in both, in the passband as well as in the stopband Chebyshev filters are often also called Chebyshev type I filters, whereas inverse Chebyshev filters are then called Chebyshev type II filters For a comparable frequency filtering performance, Butterworth filters need the highest order, Chebyshev type I and Chebyshev type II filters need medium orders, and Cauer filters need the lowest order A higher order means a higher number of operations needed to compute an output sample, and the need of a higher number of storage cells 33 Iir Butterworth
9 Transfer Function Poles denote by H Ct (s) the transfer function of the Butterworth filter from squared magnitude function we see that (jω s) 1 H Ct (s)h Ct ( s) = ( s 1 + jω c poles of squared magnitude function are s pk = ( 1) 1/(2N) (jω c ) { ( 1 = ω c exp jπ k )} 2N ) 2N, k = 0, 1,, (2N 1) 33 Iir Butterworth
10 Transfer Function Poles: N = 3 Example for N = 3 2N = 6: s pk = ( 1) 1/6 j ω c ( = e j(π+k2π)) 1/6 e jπ/2 ω c = e jπ/6 } e jk2π/6 {{}} e jπ/2 {{} ω c spacing rotation j 1 ω c 2N = 6-th roots of (-1) corresponding Butterworth poles We see then that the poles of the Butterworth squared magnitude function have an angular spacing of 2π/(2N) = π/n for an N-th order Butterworth filter We further have for these poles that they are distributed around a circle with radius ω c ; that they are distributed symmetrically on either side of the imaginary axis; that there are no poles on the imaginary axis itself; and that there are poles on the real axis if N is odd, but not if N is even 33 Iir Butterworth
11 Transfer Function to determine the Tf H Ct (s) from the Butterworth squaredmagnitude function: perform factorization H Ct (s)h Ct ( s) observe that poles in squared magnitude function appear in pairs if s = s p is a pole then s = s p is also a pole to construct H Ct (s) from squared magnitude function choose one pole of each pair choose stable pole (in Lhp) of each pair p1 poles of N = 3-rd order Butterworth transfer-function: H Ct(s) = p 0 p 1 p 1 (s p 0 )(s p 1 )(s p 1 ) p0= ω c p 1 ω c 33 Iir Butterworth
12 3 Continuous-to-Discrete Transformations C2d-Transformation Procedures there are various procedures to transform an analog prototypefilter design into a discrete-time filter: impulse-invariance transformation procedures based on numerical solution of differential equations first forward-difference (forward Euler) first backward-difference (backward Euler) bilinear transformation We discuss in more detail the impulse-invariance transformation in Subsection 31 below For the procedures based on numerical solution of differential equations we note that the forward Euler transformation is the most simple transformation that has, however, the drawback that an unstable discrete-time system might result from a stable continuous-time (analog) prototype system The transfer function H(z) of the discrete-time filter is obtained here by replacing the variable s in the continuous-time prototype transfer 33 Iir Butterworth
13 function H Ct (s) by (z 1)/T s, where T s denotes the sampling interval, T s = 1/f s with f s being the sampling frequency: forward Euler: s = 1 ) (z 1 = 1 1 z 1 T s T s z 1 The backward Euler transformation has the advantage over the forward Euler transformation that stable continuous-time designs are transformed into stable discrete-time designs The transformation is given by backward Euler: s = 1 T s z 1 z = 1 T s (1 z 1) The bilinear transformation is the most often used transformation in the design of discrete-time filters; 3 therefore, we discuss it in more detail below in Subsection 32 For the interested reader we mention that the bilinear transformation is a member of the larger family of Moebius transformations, which are conformal mappings from complex plane to complex plane; a reference is [Hen74, Chapter 5] 3 In the design of simple discrete-time control algorithms, the Euler approximations are likewise very often used 33 Iir Butterworth
14 31 Impulse Invariance Transformation The Impulse Invariance Transformation given: continuous-time (Ct) prototype (analog prototype) H Ct (s) h Ct (t) ˆ= impulse response of Ct filter find: for the discrete-time filter H(z) h[n] ˆ= impulse response of Dt filter solution: sample the continuous-time impulse response h[n] ˆ= h Ct (t = nt s ) where T s ˆ= sampling time interval The above formulae state, reformulated in words, that the impulse response of the discrete-time filter is obtained from the impulse response of the continuous-time (analog) prototype filter through sampling of the latter 33 Iir Butterworth
15 Impulse Invariance: Transfer Functions it can be shown see the sampling process that the z-transformation of h[n] H(z) is related to the Laplace transformation of h Ct (t) H Ct (s) by H(z) z=e sts = 1 T s k= ( H Ct s + j 2π ) k T s from z = e sts : strips of width 2π/T s is s-plane map onto entire z-plane 33 Iir Butterworth
16 Impulse Invariance: Mapping of Planes 1 j z-plane = 3π T s π T s π T s 3π T s R(s) I(s) s-plane strips of width 2π/T s is s-plane map onto entire z-plane Lhp-part of each strip maps into unit circle stable poles of Ct filter go to stable poles of Dt filter imaginary axis in s-plane maps onto unit circle such that each segment of length 2π/T s is mapped once around the circle We should mention here that the mapping from the s-plane to the z-plane induced by the impulse-invariance procedure is no simple algebraic mapping, as are the mappings that result from the procedures based on numerical solutions of differential equations, an example of which the bilinear transformation we discuss in Subsection Iir Butterworth
17 Impulse Invariance: Frequency Responses the Dt frequency response is expressed in terms of the Ct prototype frequency response as note: H ( z = e jˆω) = 1 T s if and only if k= ( H Ct j ˆω + j 2π ) k T s T s H Ct (jω) 0, for ω π T s then H ( z = e jˆω) = 1 T s H Ct ( j ˆω T s ) else aliasing 33 Iir Butterworth
18 Impulse Invariance: Aliasing 1 H Ct (jω) ω c 2πf s ω 1 T s H(e jˆω ) 2π π 0 ˆω c π 2π ˆω The above figure graphically shows the effects of aliasing: In the upper panel we show the frequency response of a continuoustime filter prototype which is not bandlimited; in the lower panel we show the first three aliasing parts giving, in sum together with all other aliasing parts, the discrete-time frequency response 33 Iir Butterworth
19 Impulse Invariance: Interpretation of Mapping Ct transfer function: partial fraction decomposition H Ct (s) = N k=1 A k s s k, s k ˆ= poles h Ct (t) = N A k e skt u(t) k=1 Dt impulse response and corresponding transfer function h[n] = h Ct (nt s ) = N N A k e s knt s ( u[n] = A k e s k T s ) nu[n] k=1 k=1 H(z) = N k=1 A k 1 (e s kt s)z 1 We denote by u(t) the continuous-time unit-step function and, correspondingly, by u[n] the discrete-time unit-step function Note that the above development is true if all poles have multiplicity 1; for poles with higher multiplicities we must use the corresponding partial fraction decompositions 33 Iir Butterworth
20 Impulse Invariance: Interpretation of Mapping (2) s-plane poles s k transformation of poles: transform to z-plane poles z k = e s kt s coefficients A k are equal for Ct and Dt if Ct filter is stable, then Dt filter is also stable we note however that although poles are mapped by z k = e s kt s the complete planes are not mapped by that relation for example: zeros of Dt transfer function are functions of the coefficients A k and of the poles s k Concerning the stability of the filters we know that the continuous-time filter is stable if and only if all of its poles are in the left-half complex s-plane, real{s k } < 0 By the given mapping of the poles we then have for the discrete-time poles z k that z k = e s k T s < 1, meaning that the discrete-time poles are inside of the unit circle, and in turn, that the discrete-time filter is also stable It is important to note for the impulse-invariant transformation that, although the poles are mapped from the s-plane to the z-plane by the relation z k = e s kt s, the planes themselves are not mapped by that relation For example, zeros are not mapped in that way; the following simple example gives an illustration 33 Iir Butterworth
21 Example Consider the second-order continuous-time transfer function H Ct (s) = = s + a (s + a) 2 + b 2 1/2 s + a + jb + 1/2 s + a jb, having one zero at s = a, the second zero at infinity, and a complex-conjugate pole pair at s = a ± jb The discrete-time transfer function obtained from the impulseinvariance transformation becomes H(z) = = 1/2 1 e ats e jbts z 1 + 1/2 1 e ats e +jbts z 1 1 e ats ( e +jbt s + e jbts ) /2 z 1 (1 e ats e jbts z 1 )(1 e ats e +jbts z 1 ) 1 e ats cos(bt s )z 1 = 1 e ats (e +jbts + e jbts )z 1 + e 2aTs z 2 z ( z e ats cos(bt s ) ) = z 2 2e ats cos(bt s )z + e 2aTs This discrete-time transfer function has one zero at the origin and the other zero at z = e ats cos(bt s ) We thus see that, although the poles are mapped by z k = e s kt s, the zeros are not mapped according to this formula 33 Iir Butterworth
22 Impulse Invariance: Practical Advice if Ct prototype is sufficiently bandlimited, then H ( e jˆω) 1 T s H Ct ( j ˆω T s ) thus: high sampling rates (small T s ) will minimize the aliasing, but the Dt filter may obtain a very high gain advice: instead using use H(z) = H(z) = N k=1 N k=1 the Dt impulse response then is A k 1 (e s kt s) z 1 A k T s 1 (e s kt s) z 1 h[n] = T s h Ct (t = nt s ) 33 Iir Butterworth
23 Impulse Invariance: Generalizations one motivation to use the impulse-invariance procedure if Ct prototype filter is bandlimited then Dt filter frequency response closely approximates the Ct frequency response another motivation to use the procedure: control some aspects of the time response of Dt filter step invariance procedure waveform invariance: extend the concept to preserve the output wave-shape for a variety of inputs final remark: besides aliasing, the impulse-invariance approach transforms the frequency responses linearly The mentioned step-invariance procedure just obtains the step response of the discrete-time filter by sampling the step response of the continuous-time prototype filter The resulting discrete-time filter then might have desired step-response characteristics such as small rise time and low peak overshoot An important feature of the impulse-invariance procedure is that, besides aliasing, the frequency responses of the continuoustime prototype filter and its discrete-time counterpart are linearly related, meaning that the shape of continuous-time frequency response is preserved in the discrete-time filter This result is in contrast to the procedures which use algebraic transformations, an example of which the bilinear transformation we discuss in Subsection 32 below We recall however, that the 33 Iir Butterworth
24 impulse-invariance design technique is only appropriate for essentially bandlimited filters Thus, the design of highpass or bandstop filters requires additional bandlimiting to avoid severe aliasing distortions 33 Iir Butterworth
25 32 Bilinear Transformation Bilinear Transformation Formulae given: Tf of analog prototype H Ct (s) find: Tf H(z) of corresponding discrete-time filter solution: replace s in H Ct (s) by the inverse formula then is s = 2 T s 1 z z 1 z = 1 + (T s/2)s 1 (T s /2)s You might want to verify the given formulae as follows: Start with a first-order continuous-time prototype system and first specify it by its differential equation; second, express the differential equation by its corresponding integral equation formulated for the n-th sampling time instant and its past neighbor (n 1); third, approximate the involved integral by the trapezoidal rule; forth, replace the appearing derivatives by the expression given by the differential equation started with; fifth, rearrange terms to obtain the difference equation of the discretetime system, and from it obtain the transfer function H(z) Generalize to systems of order higher than one 33 Iir Butterworth
26 Bilinear Transformation: Frequency Mapping Ct frequencies ω on imaginary axis in s-plane: s = jω Dt frequencies ˆω on unit circle in z-plane: z = e jˆω unit circle maps to imaginary axis (and vice versa) s ( z = e jˆω) = 2 T s 1 e jˆω 1 + e jˆω = = j 2 T s tan(ˆω/2) = jω therefore, frequency mapping is ω = 2 T s tan (ˆω/2) ˆω = 2 arctan(ωt s /2) To fill-in the dots in the above derivation, first use 1 e jˆω = e j ˆω 2 next use Euler s formulae sin(ˆω/2) = ( ) ( ) e j ˆω 2 e j ˆω 2, 1+e jˆω = e j ˆω 2 e j ˆω 2 + e j ˆω 2 ; ej ˆω 2 e j ˆω 2 2j, cos(ˆω/2) = ej ˆω 2 + e j ˆω 2 to obtain j sin(ˆω/2)/ cos(ˆω/2) which is equal to j tan(ˆω/2) 2, 33 Iir Butterworth
27 Bilinear Transformation: Frequency Mapping ˆω ˆ= Dt frequency ω ˆ= Ct frequency the frequency mapping is ˆω = 2 arctan(ωt s /2) π ˆω ω π 33 Iir Butterworth
28 Bilinear Transformation: Mapping of Planes the mapping s z is z = 1 + (T s/2)s 1 (T s /2)s 1 j image of lefthalf plane image of s = jω z-plane = R(s) I(s) s-plane 33 Iir Butterworth
29 Bilinear Transformation: Mapping of Circles in s-plane: Butterworth poles on a circle with radius ω c, equally spaced in angle bilinear transformation is conformal mapping: Butterworth circle in s-plane maps to a circle in z-plane j 1 1 ω ct s/2 1+ω ct s/2 1+ω ct s/2 1 ω ct s/2 z-plane = ω c s-plane Butterworth N =3,2N =6 Note, however, that the Butterworth circle in the z-plane is neither centered at the origin, nor are the poles equally spaced in angle But the left-half s-plane poles map into poles inside of the unit circle; therefore as is always the case with the bilinear transformation we obtain a stable discrete-time Butterworth filter from a continuous-time Butterworth prototype 33 Iir Butterworth
30 Bilinear Transformation: Frequency Responses ˆω ˆω 0 3dB H(e jˆω )[db] ˆω3dB π π ˆω = 2 arctan(ωt s /2) ω H Ct (jω)[db] 0 3dB ω 3dB ω 33 Iir Butterworth
31 4 Discrete-Time Butterworth Lowpass Filter Example Frequency-Response Specifications here: specifications in the Dt frequency domain requirements: passband magnitude constant to within 1 db passband is 0 ˆω ˆω p ˆ= 02π stopband attenuation > 15 db stopband is 03π ˆ= ˆω s ˆω π thus, if passband magnitude is normalized to 1 20 log 10 H(e j02π ) 1 20 log 10 H(e j03π ) 15 for convenience we may assume that the sampling-time interval is T s ˆ= 1 Note that we must distinguish between two different application scenarios: The first scenario designs a discrete-time filter which is specified in the discrete-time frequency domain This is the situation we consider here and which allows to set the sampling time equal to one, because there exists no real continuoustime filter 33 Iir Butterworth
32 The second application scenario emulates a continuous-time filter with a discrete-time filter Here the specifications are given in the continuous-time frequency domain, and our design has to select an appropriate sampling time interval T s with which we plan to implement the emulating discrete-time filter Obviously, in this scenario the sampling time interval cannot be set to unity 33 Iir Butterworth
33 41 Design Using Impulse-Invariance Transformation Analog Prototype Specifications transform the Dt frequency specifications to corresponding specifications of a Ct prototype recall: impulse invariance design introduces aliasing beside aliasing it linearly maps Ct to Dt frequencies convenient procedure: assume that aliasing is negligible carry out design verify performance of resulting filter thus, mapping of critical frequencies is ˆω p = 02π ω p = ˆω p /T s = ˆω p 20 log 10 H Ct (jω p ) 1 ˆω s = 03π ω s = ˆω s /T s = ˆω s 20 log 10 H Ct (jω s ) 15 Recall that we have assumed that the sampling time interval T s is unity 33 Iir Butterworth
34 Analog Butterworth Parameters the Ct Butterworth squared magnitude function is H Ct (jω) 2 1 = ( ω 1 + ω c ) 2N to do: determine the needed two parameters N and ω c solving given specification inequalities with equality 4 N = 1 2 ( log log ( to meet specifications: select N = 6 ) ) = inserting N = 6 into passband equation gives ω c = Note that because we must round up N to the nearest integer, not both specifications, passband and stopband, can be met exactly If we insert N = 6 into the passband equation, the passband specifications are met exactly and the stopband specifications are exceeded for the Ct filter Such a choice allows some margin for the aliasing that enters in the Dt filter We next supply the detailed steps leading to the above results We start with the two equations in the two unknowns 4 See the derivations in the text below 33 Iir Butterworth
35 N and ω c that we obtain by setting the design specification inequalities to equality: H Ct (jω p ) 2 = H Ct (jω s ) 2 = ( ωp ω c 1 ( ωs ω c ) 2N = 10 01, (1a) 1 + ( ωp ω c ) 2N = (1b) ) 2N = 10 15, (1c) 1 + ( ωs ω c ) 2N = (1d) Taking in (1b) and (1d) logarithms (to any base) we next obtain log log ( ωp ω c ( ωs ω c ) ) = log ( ) 2N = log ( ) 2N Subtracting (2b) from (2a) we obtain = log (ω p ) log (ω c ), (2a) = log (ω s ) log (ω c ) (2b) log (ω p ) log (ω s ) = 1 (log ( ) log ( )) }{{} 2N }{{} = log (ω p /ω s ) ( ) 1 = log Iir Butterworth
36 and in turn ( ) 1 N = 1 log ( ) (3) 2 ωp log ω s Finally, we insert ω p = ˆω p = 02π and ω s = ˆω s = 03π into (3) To numerically obtain the second filter parameter, ω c = 07032, via the passband specifications, we use the parameter N = 6 and ω p = 02π in (2a) Analog Butterworth Poles for N = 6 there are 2N = 12 poles of Butterworth squaredmagnitude function these 12 poles are uniformly distributed in angle on a circle with radius ω c = the Butterworth transfer function uses the N = 6 poles in the Lhp π 6 s-plane 33 Iir Butterworth
37 Analog Butterworth Transfer Function the Lhp poles are p 1 = j p 1 = j p 2 = j p 2 = j p 3 = j p 3 = j the corresponding transfer function becomes H Ct (s) = = p 1 p 1 p 2p 2 p 3p 3 (s p 1 ) (s p 1 )(s p 2)(s p 2 )(s p 3)(s p 3 ) (s s )(s s ) 1 (s s ) We have normalized the above analog transfer function H Ct (s) such that it has a Dc gain of one 33 Iir Butterworth
38 Discrete-Time Butterworth Transfer Function express H Ct (s) as a partial fraction decomposition A k, s k apply to obtain H(z) = N k=1 A k 1 (e s kt s)z 1 ( z 1 ) H(z) = ( z z 2 ) ( ) z 1 + ( z z 2 ) ( ) z 1 + ( z z 2 ) We have obtained the above parallel form of second-order sections by combining the terms of complex-conjugate pole pairs Obviously, we might directly use this parallel form that naturally results from the impulse-invariant design procedure If we desire a cascade form or a direct form, we must combine the separate second-order terms in an appropriate way 33 Iir Butterworth
39 Discrete-Time Butterworth Magnitude-Response H(e jˆω ) ˆ= 1dB ˆ= 15dB 0 ω 02π 03π π H(e jˆω ) db ˆ= 1dB ˆ= 15dB 80 ω 02π 03π π 33 Iir Butterworth
40 Discrete-Time Butterworth Phase-Response π arg(h(e jˆω )) 0 ω 02π 03π π π Final Remarks to Our Impulse-Invariance Design We recall that we have designed the filter to exactly meet with the assumption that we have no aliasing the passband specifications, that is, in the Butterworth case, to have a 1 db attenuation at the passband edge frequency ˆω p = 02π; the design then exceeds the specification at the stopband edge frequency ˆω s = 03π As we observe from the magnitude-response plots on page 37, it is true that at the passband edge the attenuation is slightly below 15 db, indicating that the aliasing is not too strong, or, in other words, that the continuous-time filter is sufficiently bandlimited In an other design this might not be true, such that the resulting discrete-time filter does not meet the specifications If we have that situation, we may try to differently adjust the filter parameters (holding the order fixed), or we may try again with a higher-order filter 33 Iir Butterworth
41 42 Design Using Bilinear Transformation Analog Prototype Specifications Dt frequency specifications must be pre-warped to corresponding analog frequencies, such that critical analog frequencies map to correct Dt frequencies the frequency mapping function of the bilinear transformation is ω = (2/T s )tan(ˆω/2) = 2 tan(ˆω/2) thus, mapping of critical frequencies is ˆω p = 02π ω p = 2 tan(ˆω p /2) = 2 tan(01π) 20 log 10 H Ct (jω p ) 1 ˆω s = 03π ω s = 2 tan(ˆω s /2) = 2 tan(015π) 20 log 10 H Ct (jω s ) 15 Recall that we conveniently have assumed that the sampling time interval is T s = 1 33 Iir Butterworth
42 Analog Butterworth Parameters the Ct Butterworth squared magnitude function is H Ct (jω) 2 1 = ( ω 1 + ω c ) 2N to do: determine the needed two parameters N and ω c solving given specification inequalities with equality 5 N = 1 2 log ( log ) ( tan(01π) tan(015π) ) = to meet specifications: select N = 6 inserting N = 6 into stopband equation gives ω c = To supply detailed steps leading to the above results we may start with (3) on page 34 Using for the bilinear transformation ω p = 2 tan(ˆω p /2) = 2 tan(01π) and ω s = 2 tan(ˆω s /2) = 2 tan(015π), we obtain N = 1 2 ( log log ) ( ωp ω s ) = 1 2 log ( log ) ( tan(01π) tan(015π) 5 See the derivations in the accompanying text ) = Iir Butterworth
43 Because the order N of the Butterworth filter must be an integer, we must select N = 6 in order to meet the specifications If we insert N = 6 into (2b), we can determine ω c to ω c = Note that by using the equation (2b) which comes from the stopband constraint to determine the ω c parameter, we meet the stopband specifications exactly and exceed the passband specifications Such a choice is reasonable, because the bilinear transformation has no aliasing effects Analog Butterworth Poles for N = 6 there are 2N = 12 poles of Butterworth squaredmagnitude function these 12 poles are uniformly distributed in angle on a circle with radius ω c = the Butterworth transfer function uses the N = 6 poles in the Lhp π 6 s-plane 33 Iir Butterworth
44 Analog Butterworth Transfer Function the Lhp poles are p 1 = j p 1 = j p 2 = j p 2 = j p 3 = j p 3 = j the corresponding transfer function becomes H Ct (s) = = p 1 p 1 p 2p 2 p 3p 3 (s p 1 ) (s p 1 )(s p 2)(s p 2 )(s p 3)(s p 3 ) (s s )(s s ) 1 (s s ) Again, we have normalized the above analog transfer function H Ct (s) such that it has a Dc gain of one 33 Iir Butterworth
45 Discrete-Time Butterworth Transfer Function apply the bilinear transformation s = 2 T s 1 z z 1 to the analog Butterworth transfer function H Ct (s) to obtain H(z) = ( 1 + z 1) 6 ( z z 2 ) 1 ( z z 2 ) 1 ( z z 2 ) 33 Iir Butterworth
46 Discrete-Time Butterworth Magnitude-Response H(e jˆω ) ˆ= db ˆ= 15dB 0 ω 02π 03π π H(e jˆω ) db ˆ= db ˆ= 15dB 80 ω 02π 03π π 33 Iir Butterworth
47 Discrete-Time Butterworth Phase-Response π arg(h(e jˆω )) 0 ω 02π 03π π π Final Remarks to Our Bilinear-Transformation Design We recall that we have designed the filter to exactly meet the stopband specifications, that is, in the Butterworth case, to have a 15 db attenuation at the stopband-edge frequency ˆω s = 03π; the design exceeds the specification at the passband-edge frequency ˆω p = 02π As we observe from the magnitude-response plots on page 43, it is true that at the stopband edge the attenuation is at 15 db, whereas at the passband edge the attenuation is only about db, leaving a certain margin to the 1 db required by the specifications Comparison of the two Designs If we compare the bilineartransformation based design on page 43 to the design based on the impulse-invariance transformation on page 37, we see 33 Iir Butterworth
48 that the magnitude function of the bilinear transformation design falls off more rapidly than the magnitude function of the impulse-invariance transformation design This is because the bilinear transformation maps the entire jω axis of the s-plane onto the unit circle, and the continuous-time Butterworth filter of order 6 has a zero of multiplicity 6 at s ; the resulting discrete-time filter then has a zero of multiplicity 6 at z = 1 References [Goe17] Josef Goette Biomedical Signal Processing and Analysis On Fixed-Point Filter Realizations Bern University of Applied Sciences, Script at the Bfh-ti Biel/Bienne, HuCE-microLab, February 2017 [Hen74] Peter Henrici Applied and Computational Complex Analysis, volume 1 of Pure and Applied Mathematics John Wiley & Sons, New York, 1974 [OS75] Alan W Oppenheim and Ronald W Schafer Digital Signal Processing Prentice-Hall Inc, Englewood Cliffs, NJ, 1975 [Ran02] Rangaraj M Rangayyan Biomedical Signal Analysis: A Case-Study Approach IEEE Press, New York, 2002 Bfh-ti Biel/Bienne Library 5708 RANGA In our development, we have mainly followed [OS75, Chapter 5, Sections 51 and 52]; there you also find design examples for Chebyshev and Cauer lowpass filters, as well as examples for using frequency transformations to design filters with highpass, bandpass, and bandstop characteristics On Butterworth filters you might also want to consult [Ran02, Chapter 3, pp 118 ff] 33 Iir Butterworth
Design 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 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 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 informationAn Fir-Filter Example: Hanning Filter
An Fir-Filter Example: Hanning Filter Josef Goette Bern University of Applied Sciences, Biel Institute of Human Centered Engineering - microlab Josef.Goette@bfh.ch February 7, 2018 Contents 1 Mathematical
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 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 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 informationDIGITAL 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 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 informationECE 410 DIGITAL SIGNAL PROCESSING D. Munson University of Illinois Chapter 12
. 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
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 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 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 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 informationLecture 8 - IIR Filters (II)
Lecture 8 - IIR Filters (II) James Barnes (James.Barnes@colostate.edu) Spring 2009 Colorado State University Dept of Electrical and Computer Engineering ECE423 1 / 27 Lecture 8 Outline Introduction Digital
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 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 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 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 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 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. 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 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 informationLecture 8 - IIR Filters (II)
Lecture 8 - IIR Filters (II) James Barnes (James.Barnes@colostate.edu) Spring 24 Colorado State University Dept of Electrical and Computer Engineering ECE423 1 / 29 Lecture 8 Outline Introduction Digital
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 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 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 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 informationDiscrete-Time David Johns and Ken Martin University of Toronto
Discrete-Time David Johns and Ken Martin University of Toronto (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) University of Toronto 1 of 40 Overview of Some Signal Spectra x c () t st () x s () t xn
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 information24 Butterworth Filters
24 Butterworth Filters Recommended Problems P24.1 Do the following for a fifth-order Butterworth filter with cutoff frequency of 1 khz and transfer function B(s). (a) Write the expression for the magnitude
More informationVALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur- 603 203 DEPARTMENT OF INFORMATION TECHNOLOGY Academic Year 2016-2017 QUESTION BANK-ODD SEMESTER NAME OF THE SUBJECT SUBJECT CODE SEMESTER YEAR
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 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 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 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 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 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 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 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 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 informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationLaplace Transform Analysis of Signals and Systems
Laplace Transform Analysis of Signals and Systems Transfer Functions Transfer functions of CT systems can be found from analysis of Differential Equations Block Diagrams Circuit Diagrams 5/10/04 M. J.
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 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 informationChapter 7. Digital Control Systems
Chapter 7 Digital Control Systems 1 1 Introduction In this chapter, we introduce analysis and design of stability, steady-state error, and transient response for computer-controlled systems. Transfer functions,
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 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 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 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 informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
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 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 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 informationMultirate Digital Signal Processing
Multirate Digital Signal Processing Basic Sampling Rate Alteration Devices Up-sampler - Used to increase the sampling rate by an integer factor Down-sampler - Used to decrease the sampling rate by an integer
More informationCast of Characters. Some Symbols, Functions, and Variables Used in the Book
Page 1 of 6 Cast of Characters Some s, Functions, and Variables Used in the Book Digital Signal Processing and the Microcontroller by Dale Grover and John R. Deller ISBN 0-13-081348-6 Prentice Hall, 1998
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 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 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 IIR Filters Using Cascaded Biquads
Design IIR Filters Using Cascaded Biquads This article shows how to implement a Butterworth IIR lowpass filter as a cascade of second-order IIR filters, or biquads. We ll derive how to calculate the coefficients
More informationX. Chen More on Sampling
X. Chen More on Sampling 9 More on Sampling 9.1 Notations denotes the sampling time in second. Ω s = 2π/ and Ω s /2 are, respectively, the sampling frequency and Nyquist frequency in rad/sec. Ω and ω denote,
More informationELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization
ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals 1. Sampling and Reconstruction 2. Quantization 1 1. Sampling & Reconstruction DSP must interact with an analog world: A to D D to A x(t)
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 informationDesign IIR Butterworth Filters Using 12 Lines of Code
db Design IIR Butterworth Filters Using 12 Lines of Code While there are plenty of canned functions to design Butterworth IIR filters [1], it s instructive and not that complicated to design them from
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 informationA.1 THE SAMPLED TIME DOMAIN AND THE Z TRANSFORM. 0 δ(t)dt = 1, (A.1) δ(t)dt =
APPENDIX A THE Z TRANSFORM One of the most useful techniques in engineering or scientific analysis is transforming a problem from the time domain to the frequency domain ( 3). Using a Fourier or Laplace
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #20 Wednesday, October 22, 2003 6.4 The Phase Response and Distortionless Transmission In most filter applications, the magnitude response H(e
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 informationDiscrete Time Signals and Switched Capacitor Circuits (rest of chapter , 10.2)
Discrete Time Signals and Switched Capacitor Circuits (rest of chapter 9 + 0., 0.2) Tuesday 6th of February, 200, 9:5 :45 Snorre Aunet, sa@ifi.uio.no Nanoelectronics Group, Dept. of Informatics Office
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 informationDSP Design Lecture 2. Fredrik Edman.
DSP Design Lecture Number representation, scaling, quantization and round-off Noise Fredrik Edman fredrik.edman@eit.lth.se Representation of Numbers Numbers is a way to use symbols to describe and model
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 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 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 informationDigital Wideband Integrators with Matching Phase and Arbitrarily Accurate Magnitude Response (Extended Version)
Digital Wideband Integrators with Matching Phase and Arbitrarily Accurate Magnitude Response (Extended Version) Ça gatay Candan Department of Electrical Engineering, METU, Ankara, Turkey ccandan@metu.edu.tr
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 informationDiscrete Time Signals and Switched Capacitor Circuits (rest of chapter , 10.2)
Discrete Time Signals and Switched Capacitor Circuits (rest of chapter 9 + 10.1, 10.2) Tuesday 16th of February, 2010, 0, 9:15 11:45 Snorre Aunet, sa@ifi.uio.no Nanoelectronics Group, Dept. of Informatics
More informationFrom Wikipedia, the free encyclopedia
Z-transform - Wikipedia, the free encyclopedia Z-transform Page 1 of 7 From Wikipedia, the free encyclopedia In mathematics and signal processing, the Z-transform converts a discrete time domain signal,
More informationSignal Processing First Lab 11: PeZ - The z, n, and ˆω Domains
Signal Processing First Lab : PeZ - The z, n, and ˆω Domains The lab report/verification will be done by filling in the last page of this handout which addresses a list of observations to be made when
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 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 informationUniversity of Illinois at Chicago Spring ECE 412 Introduction to Filter Synthesis Homework #4 Solutions
Problem 1 A Butterworth lowpass filter is to be designed having the loss specifications given below. The limits of the the design specifications are shown in the brick-wall characteristic shown in Figure
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 information2.161 Signal Processing: Continuous and Discrete Fall 2008
MT OpenCourseWare http://ocw.mit.edu 2.6 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
More informationOptimal Discretization of Analog Filters via Sampled-Data H Control Theory
Optimal Discretization of Analog Filters via Sampled-Data H Control Theory Masaaki Nagahara 1 and Yutaka Yamamoto 1 Abstract In this article, we propose optimal discretization of analog filters or controllers
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 informationAppendix A Butterworth Filtering Transfer Function
Appendix A Butterworth Filtering Transfer Function A.1 Continuous-Time Low-Pass Butterworth Transfer Function In order to obtain the values for the components in a filter, using the circuits transfer function,
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science. Fall Solutions for Problem Set 2
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Issued: Tuesday, September 5. 6.: Discrete-Time Signal Processing Fall 5 Solutions for Problem Set Problem.
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 informationExercises in Digital Signal Processing
Exercises in Digital Signal Processing Ivan W. Selesnick September, 5 Contents The Discrete Fourier Transform The Fast Fourier Transform 8 3 Filters and Review 4 Linear-Phase FIR Digital Filters 5 5 Windows
More informationELEG 5173L Digital Signal Processing Ch. 5 Digital Filters
Department of Electrical Engineering University of Aransas ELEG 573L Digital Signal Processing Ch. 5 Digital Filters Dr. Jingxian Wu wuj@uar.edu OUTLINE 2 FIR and IIR Filters Filter Structures Analog Filters
More informationSIDDHARTH GROUP OF INSTITUTIONS:: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE)
SIDDHARTH GROUP OF INSTITUTIONS:: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : Digital Signal Processing(16EC422) Year & Sem: III-B.Tech & II-Sem Course
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 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 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 informationLectures: Lumped element filters (also applies to low frequency filters) Stub Filters Stepped Impedance Filters Coupled Line Filters
ECE 580/680 Microwave Filter Design Lectures: Lumped element filters (also applies to low frequency filters) Stub Filters Stepped Impedance Filters Coupled Line Filters Lumped Element Filters Text Section
More informationOversampling Converters
Oversampling Converters David Johns and Ken Martin (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) slide 1 of 56 Motivation Popular approach for medium-to-low speed A/D and D/A applications requiring
More informationDigital Signal Processing
Digital Signal Processing Introduction Moslem Amiri, Václav Přenosil Embedded Systems Laboratory Faculty of Informatics, Masaryk University Brno, Czech Republic amiri@mail.muni.cz prenosil@fi.muni.cz February
More informationEE Experiment 11 The Laplace Transform and Control System Characteristics
EE216:11 1 EE 216 - Experiment 11 The Laplace Transform and Control System Characteristics Objectives: To illustrate computer usage in determining inverse Laplace transforms. Also to determine useful signal
More information