DOCSLIB.ORG

IIR Filter Design 9.1 – IIR Filter

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
IIR Filter Design 9.1 – IIR Filter www.getmyuni.com Digital Signal Processing IIR Filter Design www.getmyuni.com 9.1 – IIR Filter ¾ Difference Equation N N y(n) = ∑ ai x(n − i) +∑bj y(n − j) i=0 j=1 N ¾Transfer Function −i ∑ ai z i=0 H (z) = N − j 1+ ∑bj z j=1 where; N is the filter’s order. ai,bj are the filter’s coefficients. www.getmyuni.com 9.1 – IIR Filter – contd. • IIR filter have infinite-duration impulse responses, hence they can be matched to analog filters, all of which generally have infinitely long impulse responses. • The basic techniques of IIR filter design transform well-known analog filters into digital filters using complex-valued mappings. • The advantage of these techniques lie in the fact that both analog filter design (AFD) tables and the mappings are available extensively in the literature. • The basic techniques are called the A/D filter transformation. • However, thAFDtblthe AFD tables are avail ilblable onl y f or lowpass filters. W e al so want design other frequency-selective filters (highpass, bandpass, bandstop, etc.) • To do this, we need to apply frequency-band transformations to lowpass filters. These transformations are also complex-valued mappings, and they are also available in the literature. www.getmyuni.com 9.2 – Design Methods ¾ The following A/D filter transformation methods are used in calculatinggffff the coefficients of IIR filter: 1. Impulse Invariant Method. 2. Bilinear Z-transform Method. Æ to achieve a diggfital filter that has a certain sppfecification e quivalent to an analogue filter Note ¾ We have no control over the phase characteristics of the IIR filter. Hence IIR filter designs will be treated as magnitude-only designs. www.getmyuni.com 9.3 – Impulse Invariant Method • Here we require that the impulse response of the discrete system (digital filter) be the discrete version of the impulse response of the analogue system (filter), (Hence t he name: impu lse invariant ). H(s) Ł-1 h(t) t = nT h(nT) Z X T H(z) To remove the sampling effect: ¾ Steps: (1/T) δ(t‐nT) 1. Get H(s), continuous time transfer function of an analogue filter that satisfies the prescribed magnitude response. 2. Apply the inverse Laplace transform to get the impulse response h(t). h(n). 3. Obtain a discrete version of h(t) by replacing t by nT Æ h(nT). h(t) 4. Apply the Z-transform to h(nT) to get H(z) and multiply by T. h(nT) n n ∞ www.getmyuni.com Example - 1 ¾ Using the invariant impulse response method, design a digital filter that has the shown pole-zero distribution. jω ω 9 Solution x 0 s + a H (s) = σ [s + (a + jω )][s + (a − jω )] ‐ a s + a 0 0 H (s) = x - ω0 2 2 [(s + a) + ω0 )] −1 −at (1) L {H(s)} = e cos(ω0t) = h(t) −anT (2) h(nT) = e cos(ω0nT) −aT −1 1− e cos(ω0T )z (3) T × Z{h(nT)} = −aT −1 −2aT −2 ×T 1− 2e cos(ω0T )z + e z −aT −1 T − (e T cos(ω0T ))z H (z) = −aT −1 −2aT −2 1− (2e cos(ω0T ))z + e z www.getmyuni.com Example - 1 – cont. −aT −1 T − (e T cos(ω0T ))z H (z) = −aT −1 −2aT −2 1− (2e cos(ω0T ))z + e z a x(n) 0 y(n) X -b1 T a1 X X - b2 T X where, a 0 = T − aT a1 = −e T cos( ω 0T ) − aT b1 = −2e cos( ω 0T ) − 2 aT b2 = e www.getmyuni.com 9.4 – Bilinear Z-Transform (BZT) Method s t z s z s1 Impulse Invariant Method Bilinear Z‐Transform Method • It is the most important method of obtaining IIR filter coefficients. • In the BZT method, the basis operation is to convert an analogue filter H(s) into an equiva lent di g ita l filter H(z) b y usi ng th e bili near approxi mati on. Ha(s) HD(z) sT sT e 2 1 + + ... z = e sT = = 2 1st order bilinear Q − sT approximation 2 1 − sT + ... e 2 1 + sT ∴ z ≅ 2 1 − sT 2 www.getmyuni.com 9.4 – BZT Method – cont. sT 1 + sT sT z ≅ 2 z − z = 1 + 1 − sT 2 2 2 T z - 1 = (z + 1) s 2 2 z −1 ∴s ≅ T z +1 1 H (z) = H (s) 2 z−1 D a s= T z+1 www.getmyuni.com 9.4 – BZT Method – cont. − jΩT 2 z −1 2 e jΩT −1 e 2 s = = × Q jΩT − jΩT T z +1 T e +1 e 2 j ΩT − j ΩT ΩT 2 e 2 − e 2 2 2 j sin ( ) = = 2 j ΩT − j ΩT T 2 2 T 2 cos ΩT e + e ( 2 ) 2 ΩT jω a = j tan ( ) T 2 2 ΩT Pre-Warppgffing effect ω a = tan ( ) T 2 2 −1 ⎛ω T ⎞ Ω = tan ⎜ a Warppgffing effect T ⎝ 2 ⎠ Digital frequency www.getmyuni.com 9.5 – Frequency Warping Effect (BZT) Analogue frequency • The effect of the bilinear Z-transform is called ωa frequency warping: every feature appears in the frequency response o f th e conti nuous-time filt er appears as it is in the frequency response of the discrete-time filter but at different frequency. jω when − ∞ ≤ ω a ≤ ∞ S domain tan(θ) π π Ω − ≤ Ω ≤ σ Digital frequency T T ‐π/2 jΩ θ Ω s Ω s Ωs/2 π/2 − ≤ Ω ≤ S domain 2 2 1 - Ωs/2 • To compensate for this, every frequency specification the designer has a control over (ωc,ωs,…) has to be ‘prewarped’ by setting it with; 2 ω = tan (ΩT ) 2 T 2 www.getmyuni.com 9.6 - IIR Filter Design Procedure Using BZT 1 Given specification in digital domain 2 Convert it into analog filter specification(prewarping) 3 Design analog filter (Butterworth , Chebyshov , elliptic):H(s) 4 Apply bilinear transform to get H(z) out of H(s) ω ω ωs 3 2 ΩT 2 ω a = tan ( ) T 2 ωp | H(jω)| Ω 1 1 1 | H(ejΩ)| 2 1 1+ε A 1 1 + ε 2 4 1 1 H (z) = H (s) 2 1− z −1 A s= ⋅ Ω −1 Ω π T 1+ z p Ωs www.getmyuni.com 9.7 – IIR Filter Design Steps Using BZT 1. Prewarp any critical frequency in the digital filter specifications(ωc,ωp,ωs,…) using; 1 ω = tan Ω T ;T = ( 2 ) Fs 2. Use the diggfital filter sp pfecifications to f ind a suitable normalized p rototyp e analogue LPF, H(s), e.g., for butterworth; 1 • 1st order: H (s) = s +1 1 • 2nd order: H (s) = s2 + 2s +1 For butterworth prototype filter 1 • 3rd order: H (s) = (s +1)(s2 + s +1) Where ⎛ 10 0.1 As − 1 ⎞ log ⎜ ⎟ 10 ⎜ 0.1 A p A the passband ripple in dB ⎝ 10 − 1 ⎠ p N ≥ As the stopband attenuation in dB ⎛ ω ⎞ 2.log ⎜ s ⎟ 10 ⎜ ⎟ ⎝ ω p ⎠ www.getmyuni.com 9.7 – Design Steps – cont. 2. Denormalization according to filter type; (LPF, HPF, BPF, BSF). • LPF s s ω c • HPF s ω c s 2 2 2 • BPF s (s + ω 0 ) sω ;where, ω0 = ω1ω2 ω0 = ω1ω2 2 2 • BSF ω = ω2 −ω1 s sω (s + ω 0 ) 3. MfMap from s-didomain to z-didomain; 1− z −1 s = 1+ z −1 4. Realize the IIR filter . www.getmyuni.com Example – 2 st ¾ Design a 1 order Butterworth HPF with Ωc = 1 rad/sec, Fs = 1Hz. 9 Solution 1 H (s) = s ω s s +1 c 1 s H (s) = = ωc s +1 s +ωc ⎛ Ω cT ⎞ ⎛ 1 ⎞ ω c = tan ⎜ = tan ⎜ ⎝ 2 ⎠ ⎝ 2 ⎠ (1− z −1) (1+ z −1) 1− z −1 H (z) = = (1− z −1) 1− z −1 +ω +ω z −1 +ω c c (1+ z −1) c www.getmyuni.com Example - 2 – cont. 1 1− z −1 H (z) = 1 ω −1 k = 1+ωc c −1 1+ z 1+ωc ωc +1 k a x(n) 0 y(n) X X - b1 T a1 X X where, a 0 = 1 a1 = −1 ω c − 1 b1 = ω c + 1 www.getmyuni.com Example - 3 ¾ Design a 3rd order Butterworth LPF to have a cutoff frequency at 4 kHz using the BZT method & assuming a sampling frequency of 10 kHz. 9 Solution 1 −4 fc = 4 kHz & T = =10 sec Fs 3 −4 ΩcT = 2π × 4×10 ×10 = 2.513 rad Ω T prewarping : ω = tan( c ) = 3.0762 rad c 2 sec 1 H (s) = (s +1)(s2 + s +1) 1 H(z) = 2 2 ⎛ 1 ⎛1− z −1 ⎞ ⎞⎛⎛ 1 ⎞ ⎛1− z−1 ⎞ 1 ⎛1− z−1 ⎞ ⎞ ⎜ ⎜ +1 ⎜⎜ ⎜ + ⎜ +1 ⎜ω ⎜1+ z −1 ⎜⎜ω ⎜1+ z−1 ω ⎜1+ z−1 ⎟ ⎝ c ⎝ ⎠ ⎠⎝⎝ c ⎠ ⎝ ⎠ c ⎝ ⎠ ⎠ www.getmyuni.com Example - 3 – cont. (1+z−1) (1+z−1)2 H(z)= −1 −1 2 −2 −1 2 (1.3249 + 0.6751 z )(0.105 (1−z ) + 0.3249 (1−z )+(1+z ) ) (1 + z −1 ) (1 + 2 z −1 + z −2 ) = 1 . 3249 (1 + 0.5095 z −1 ) ()1.4305 + 1.789 z −1 + 0.7806 z − 2 (1 + z −1 ) (1 + 2 z −1 + z −2 ) = × 1.3249 (1 + 0.5095 z −1 ) 1.4305 ()1 + 1.2506 z −1 + 0.5457 z − 2 (1 + z −1 ) (1 + 2 z −1 + z −2 ) = 0.52763 × × (1 + 0.5095 z −1 ) (1 + 1.2506 z −1 + 0.5457 z − 2 ) www.getmyuni.com Example - 3 – cont.
Load more

Recommended publications
  • Virtual Analog (VA) Filter Implementation and Comparisons Copyright © 2013 Will Pirkle
    App Note 4"Virtual Analog (VA) Filter Implementation and Comparisons "Copyright © 2013 Will Pirkle Virtual Analog (VA) Filter Implementation and Comparisons Will Pirkle I have had several requests from readers to do a Virtual Analog (VA) Filter Implementation plug-in. The source of these designs is a book electronically published in June 2012 named The Art of VA Filter Design by Vadim Zavalishin. This awesome piece of work is free and available from many sources including my own site www.willpirkle.com. This short book is an excellent introduction to basic analog filtering theory as well as digital transformations. It is so concise in this respect that I am considering using it as part of the text materials for an Advanced Analog Circuits class I teach. I highly recommend this excellent book - you will need to understand its content to use this App Note; for the most part I use the same variable names as the book so you will want to use it as a reference. Zavalishinʼs derivations and descriptions are so well thought out and so well written that it makes no sense for me to repeat them here and I donʼt think it can be simplified any more that it already is in his book. Another reason that I enjoyed this book is that the author followed a similar derivation to reach the bilinear transform as my DSP professor (Claude Lindquist) did when I was in graduate school; in fact I use part of that same derivation in my classes and book. Also similar was the use of an integrator as a prototype filter to generate the discrete time transforms.
    [Show full text]
  • ECGR4124 Digital Signal Processing Exam 2 Spring 2017 Name
    ECGR4124 Digital Signal Processing Exam 2 Spring 2017 Name: _____________________________________ LAST 4 NUMBERS of Student Number: _____ Do NOT begin until told to do so Make sure that you have all pages before starting NO TEXTBOOK, NO CALCULATOR, NO CELL PHONES/WIRELESS DEVICES Open handouts, 2 sheet front/back notes, NO problem handouts, NO exams, NO quizzes DO ALL WORK IN THE SPACE GIVEN Do NOT use the back of the pages, do NOT turn in extra sheets of work/paper Multiple-choice answers should be within 5% of correct value Show ALL work, even for multiple choice ACADEMIC INTEGRITY: Students have the responsibility to know and observe the requirements of The UNCC Code of Student Academic Integrity. This code forbids cheating, fabrication or falsification of information, multiple submission of academic work, plagiarism, abuse of academic materials, and complicity in academic dishonesty. Unless otherwise noted: F{} denotes Discrete time Fourier transform {DTFT, DFT, or Continuous, as implied in problem} F-1{} denotes inverse Fourier transform ω denotes frequency in rad/sample, Ω denotes frequency in rad/second ∗ denotes linear convolution, N denotes circular convolution x*(t) denotes the conjugate of x(t) Useful constants, etc: e ≈ 2.72 π ≈ 3.14 e2 ≈ 7.39 e4 ≈ 54.6 e-0.5 ≈ 0.607 e-0.25 ≈ 0.779 1/e ≈ 0.37 √2 ≈ 1.41 e-2 ≈ 0.135 √3 ≈ 1.73 e-4 ≈ 0.0183 √5 ≈ 2.22 √7 ≈ 2.64 √10 ≈ 3.16 ln( 2 ) ≈ 0.69 ln( 4 ) ≈ 1.38 log10( 2 ) ≈ 0.30 log10( 3 ) ≈ 0.48 log10( 10 ) ≈ 1.0 log10( 0.1 ) ≈ -1 1/π ≈ 0.318 sin(0.1) ≈ 0.1 tan(1/9) ≈ 1/9 cos(π / 4) ≈ 0.71 cos( A ) cos ( B ) = 0.5 cos(A - B) + 0.5 cos(A + B) ejθ = cos(θ) + j sin(θ) 1/10 5 Points Each, Circle the Best Answer 1.
    [Show full text]
  • Digital Signal Processing I Exam 2 Fall 1999 Session 17 Live: 21 Oct
    Digital Signal Processing I Exam 2 Fall 1999 Session 17 Live: 21 Oct. 1999 Cover Sheet Test Duration: 75 minutes. Open Book but Closed Notes. Calculators NOT allowed. This test contains four problems. All work should be done in the blue books provided. Do not return this test sheet, just return the blue books. Prob. No. Topic of Problem Points 1. Digital Upsampling 35 2. Digital Subbanding 25 3. Multi-Stage Upsampling/Interpolation 20 4. IIR Filter Design Via Bilinear Transform 20 1 Digital Signal Processing I Exam 2 Fall 1999 Session 17 Live: 21 Oct. 1999 Problem 1. [35 points] X (F) H (ω) a LP 1/4W 2 F ω W W π 3π π π 3π π 4 4 4 4 x (t) h [n] a Ideal A/D x [n] w [n] Lowpass Filter y [n] 2 π LP 3π F = 4W ω = ω = s p 4 s 4 gain =2 Figure 1. The analog signal xa(t) with CTFT Xa(F ) shown above is input to the system above, where x[n]=xa(n/Fs)withFs =4W ,and sin( π n) cos( π n) h [n]= 2 4 , −∞ <n<∞, LP π n − n2 2 1 4 | |≤ π 3π ≤| |≤ such that HLP (ω)=2for ω 4 , HLP (ω)=0for 4 ω π,andHLP (ω) has a cosine π 3π roll-off from 1 at ωp = 4 to 0 at ωs = 4 . Finally, the zero inserts may be mathematically described as ( x( n ),neven w[n]= 2 0,nodd (a) Plot the magnitude of the DTFT of the output y[n], Y (ω), over −π<ω<π.
    [Show full text]
  • IIR Filters (II)
    Lecture 8 - IIR Filters (II) James Barnes ([email protected]) Spring 2009 Colorado State University Dept of Electrical and Computer Engineering ECE423 – 1 / 27 Lecture 8 Outline ● Introduction ● Digital Filter Design by Analog → Digital Conversion ● (Probably next lecture) ”All Digital” Design Algorithms ● (Next lecture) Conversion of Filter Types by Frequency Transformation Colorado State University Dept of Electrical and Computer Engineering ECE423 – 2 / 27 ❖ Lecture 8 Outline Introduction ❖ IIR Filter Design Overview Method: Impulse Invariance for IIR FIlters Approximation of Derivatives Bilinear Transform Matched Z-Transform Introduction Colorado State University Dept of Electrical and Computer Engineering ECE423 – 3 / 27 IIR Filter Design Overview ● Methods which start from analog design ✦ Impulse Invariance ✦ Approximation of Derivatives ✦ Bilinear Transform ✦ Matched Z-transform All are different methods of mapping the s-plane onto the z-plane ● Methods which are ”all digital” ✦ Least-squares ✦ McClellan-Parks Colorado State University Dept of Electrical and Computer Engineering ECE423 – 4 / 27 ❖ Lecture 8 Outline Introduction Method: Impulse Invariance for IIR FIlters ❖ Impulse Invariance ❖ Impulse Invariance (2) ❖ Impulse Invariance (3) ❖ Impulse Invariance (5) ❖ Impulse Invariance Procedure ❖ Impulse Invariance Example Method: Impulse Invariance for IIR FIlters ❖ Impulse Invariance Example (2) Approximation of Derivatives Bilinear Transform Matched Z-Transform Colorado State University Dept of Electrical and Computer Engineering ECE423 – 5 / 27 Impulse Invariance We start by sampling the impulse response of the analog filter: ha(t) h[n]= ha(nt0) t0 Sampling Theorem gives relation between Fourier Transform of sampled and continuous ”signals”: ∞ 1 ω 2πk H(z)|z=ejω = Ha(j − j ), (1) t0 t0 t0 k=X−∞ where ω =Ωt0 = 2πf/fs and f is the analog frequency in Hz.
    [Show full text]
  • Signal Approximation Using the Bilinear Transform
    SIGNAL APPROXIMATION USING THE BILINEAR TRANSFORM Archana Venkataraman, Alan V. Oppenheim MIT Digital Signal Processing Group 77 Massachusetts Avenue, Cambridge, MA 02139 [email protected], [email protected] ABSTRACT The analysis presented in this paper has application in contexts This paper explores the approximation properties of a unique basis where only a fixed number of DT values can be used to represent a expansion, which realizes a bilinear frequency warping between a CT signal. For example, in a binary detection problem, one might continuous-time signal and its discrete-time representation. We in- want to limit the number of digital multiplies used to compute the vestigate the role that certain parameters and signal characteristics inner product of two CT signals. Numerical simulations of this sce- have on these approximations, and we extend the analysis to a win- nario suggest that the bilinear expansion achieves a better detection dowed representation, which increases the overall time resolution. performance than Nyquist sampling for certain signal classes. Approximations derived from the bilinear representation and from Nyquist sampling are compared in the context of a binary detection 2. THE BILINEAR REPRESENTATION problem. Simulation results indicate that, for many types of signals, the bilinear approximations achieve a better detection performance. As derived in [1], the network shown in Fig. 1 realizes a one-to-one Index Terms— Signal Representations, Approximation Meth- frequency warping between the Laplace and Z-transform domains ods, Bilinear Transformations, Signal Detection according to the bilinear transform. Specifically, the Laplace trans- form F (s) of the signal f(t) and the Z-transform F (z) of the se- f[n] 1.
    [Show full text]
  • Visual and Intuitive Approach to Explaining Digitized Controllers
    Paper ID #16436 Visual and Intuitive Approach to Explaining Digitized Controllers Dr. Daniel Raviv, Florida Atlantic University Dr. Raviv is a Professor of Computer & Electrical Engineering and Computer Science at Florida Atlantic University. In December 2009 he was named Assistant Provost for Innovation and Entrepreneurship. With more than 25 years of combined experience in the high-tech industry, government and academia Dr. Raviv developed fundamentally different approaches to ”out-of-the-box” thinking and a breakthrough methodology known as ”Eight Keys to Innovation.” He has been sharing his contributions with profession- als in businesses, academia and institutes nationally and internationally. Most recently he was a visiting professor at the University of Maryland (at Mtech, Maryland Technology Enterprise Institute) and at Johns Hopkins University (at the Center for Leadership Education) where he researched and delivered processes for creative & innovative problem solving. For his unique contributions he received the prestigious Distinguished Teacher of the Year Award, the Faculty Talon Award, the University Researcher of the Year AEA Abacus Award, and the President’s Leadership Award. Dr. Raviv has published in the areas of vision-based driverless cars, green innovation, and innovative thinking. He is a co-holder of a Guinness World Record. His new book is titled: ”Everyone Loves Speed Bumps, Don’t You? A Guide to Innovative Thinking.” Dr. Daniel Raviv received his Ph.D. degree from Case Western Reserve University in 1987 and M.Sc. and B.Sc. degrees from the Technion, Israel Institute of Technology in 1982 and 1980, respectively. Paul Benedict Caballo Reyes, Florida Atlantic University Paul Benedict Reyes is an Electrical Engineering major in Florida Atlantic University who expects to graduate Spring 2016.
    [Show full text]
  • The Z-Transform
    CHAPTER 33 The z-Transform Just as analog filters are designed using the Laplace transform, recursive digital filters are developed with a parallel technique called the z-transform. The overall strategy of these two transforms is the same: probe the impulse response with sinusoids and exponentials to find the system's poles and zeros. The Laplace transform deals with differential equations, the s-domain, and the s-plane. Correspondingly, the z-transform deals with difference equations, the z-domain, and the z-plane. However, the two techniques are not a mirror image of each other; the s-plane is arranged in a rectangular coordinate system, while the z-plane uses a polar format. Recursive digital filters are often designed by starting with one of the classic analog filters, such as the Butterworth, Chebyshev, or elliptic. A series of mathematical conversions are then used to obtain the desired digital filter. The z-transform provides the framework for this mathematics. The Chebyshev filter design program presented in Chapter 20 uses this approach, and is discussed in detail in this chapter. The Nature of the z-Domain To reinforce that the Laplace and z-transforms are parallel techniques, we will start with the Laplace transform and show how it can be changed into the z- transform. From the last chapter, the Laplace transform is defined by the relationship between the time domain and s-domain signals: 4 X(s) ' x(t) e &st dt m t '&4 where x(t) and X(s) are the time domain and s-domain representation of the signal, respectively.
    [Show full text]
  • IIR Filters (II)
    Lecture 8 - IIR Filters (II) James Barnes ([email protected]) Spring 2014 Colorado State University Dept of Electrical and Computer Engineering ECE423 – 1 / 29 Lecture 8 Outline ● Introduction ● Digital Filter Design by Analog → Digital Conversion ● (Probably next lecture) ”All Digital” Design Algorithms ● (Next lecture) Conversion of Filter Types by Frequency Transformation Colorado State University Dept of Electrical and Computer Engineering ECE423 – 2 / 29 ❖ Lecture 8 Outline Introduction ❖ IIR Filter Design Overview Method: Impulse Invariance for IIR FIlters Approximation of Derivatives Bilinear Transform Matched Z-Transform Introduction Colorado State University Dept of Electrical and Computer Engineering ECE423 – 3 / 29 IIR Filter Design Overview ● Methods which start from analog design ✦ Impulse Invariance ✦ Approximation of Derivatives ✦ Bilinear Transform ✦ Matched Z-transform All are different methods of mapping the s-plane onto the z-plane ● Methods which are ”all digital” ✦ Least-squares ✦ McClellan-Parks Colorado State University Dept of Electrical and Computer Engineering ECE423 – 4 / 29 ❖ Lecture 8 Outline Introduction Method: Impulse Invariance for IIR FIlters ❖ Impulse Invariance ❖ Impulse Invariance (2) ❖ Impulse Invariance (3) ❖ Impulse Invariance (4) ❖ Impulse Invariance Procedure ❖ Impulse Invariance Example Method: Impulse Invariance for IIR FIlters ❖ Impulse Invariance Example (2) Approximation of Derivatives Bilinear Transform Matched Z-Transform Colorado State University Dept of Electrical and Computer Engineering ECE423 – 5 / 29 Impulse Invariance We start by sampling the impulse response of the analog filter: ha(t) h[n]= ha(nt0) t0 Sampling Theorem gives relation between Fourier Transform of sampled and continuous ”signals”: ∞ 1 ω 2πk H(z)|z=ejω = Ha(j − j ), (1) t0 t0 t0 k=X−∞ where ω =Ωt0 = 2πf/fs and f is the analog frequency in Hz.
    [Show full text]
  • Al-Alaoui Operator and the New Transformation Polynomials for Discretization of Analogue Systems
    Electr Eng (2008) 90:455–467 DOI 10.1007/s00202-007-0092-0 ORIGINAL PAPER Al-Alaoui operator and the new transformation polynomials for discretization of analogue systems Mohamad Adnan Al-Alaoui Received: 22 August 2007 / Accepted: 1 December 2007 / Published online: 8 January 2008 © Springer-Verlag 2007 Abstract The “new transformation polynomials for discret- The bilinear transform meets the above requirements. ization of analogue systems” was recently introduced. The However, it introduces a warping effect due to its nonlinear- work proposes that the discretization of 1/sn should be done ity, albeit it can be ameliorated somewhat by a pre-warping independently rather than by raising the discrete representa- technique. tion of 1/s to the power n. Several examples are given in to The backward difference transform satisfies the first con- back this idea. In this paper it is shown that the “new transfor- dition, but the second condition is not completely satisfied, mation polynomials for discretization of analogue systems” since the imaginary axis of the s-plane maps onto the cir- is exactly the same as the parameterized Al-Alaoui opera- cumference in the z-plane centered at z = 1/2 and having a tor. In the following sections, we will show that the same radius of 1/2. The mapping meets condition 2 rather closely results could be obtained with the parameterized Al-Alaoui for low frequencies [1–6]. operator. Other transforms were introduced in attempts to obtain better approximations [7–11]. In particular, in [7,8]the Keywords Al-Alaoui operator · Digital filters · approach interpolates the rectangular integration rules and Discretization · s-to-z transforms the trapezoidal integration rule.
    [Show full text]
  • Digital Signal Processing
    Chapter 5. Digital Signal Processing Digital Signal Processing Academic and Research Staff Professor Alan V. Oppenheim, Professor Arthur B. Baggeroer, Dr. Charles E. Rohrs. Visiting Scientists and Research Affiliates Dr. Dan E. Dudgeon1, Dr. Yonina Eldar2, Dr. Ehud Weinstein3, Dr. Maya R. Said4 Graduate Students Thomas Baran, Sourav Dey, Zahi Karam, Alaa Kharbouch, Jon Paul Kitchens, Melanie Rudoy, Joseph Sikora III, Archana Venkataraman, Dennis Wei, Matthew Willsey. UROP Student Jingdong Chen Technical and Support Staff Eric Strattman Introduction The Digital Signal Processing Group develops signal processing algorithms that span a wide variety of application areas including speech and image processing, sensor networks, communications, radar and sonar. Our primary focus is on algorithm development in general, with the applications serving as motivating contexts. Our approach to new algorithms includes some unconventional directions, such as algorithms based on fractal signals, chaotic behavior in nonlinear dynamical systems, quantum mechanics and biology in addition to the more conventional areas of signal modeling, quantization, parameter estimation, sampling and signal representation. When developing new algorithms, we often look to nature for inspiration and as a metaphor for new signal processing directions. Falling into this category to a certain extent, is our previous work on fractals, chaos, and solitons. 1 BAE Systems IEWS, Senior Principal Systems Engineer, Nashua, New Hampshire. 2 Department of Electrical Engineering, Faculty of Engineering, Technion-Israel Institute of Technology, Israel. 3 Department of Electrical Engineering, Systems Division, Faculty of Engineering, Tel-Aviv University, Israel; adjunct scientist, Department of Applied Ocean Physics and Engineering, Woods Hole Oceanographic Institution, Woods Hole, Massachusetts. 4 Visiting Scientist, Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts.
    [Show full text]
  • Design of IIR Filters.Pdf
    Design of IIR Filters Elena Punskaya www-sigproc.eng.cam.ac.uk/~op205 Some material adapted from courses by Prof. Simon Godsill, Dr. Arnaud Doucet, Dr. Malcolm Macleod and Prof. Peter Rayner 1 IIR vs FIR Filters 2 IIR as a class of LTI Filters Difference equation: Transfer function: To give an Infinite Impulse Response (IIR), a filter must be recursive, that is, incorporate feedback N ≠ 0, M ≠ 0 the recursive (previous output) terms feed back energy into the filter input and keep it going. (Although recursive filters are not necessarily IIR) 3 IIR Filters Design from an Analogue Prototype • Given filter specifications, direct determination of filter coefficients is too complex • Well-developed design methods exist for analogue low-pass filters • Almost all methods rely on converting an analogue filter to a digital one 4 Analogue filter Rational Transfer Function () 5 Analogue to Digital Conversion Im (z) Re (z) 6 Impulse Invariant method 7 Impulse Invariant method: Steps 1. Compute the Inverse Laplace transform to get impulse response of the analogue filter 2. Sample the impulse response (quickly enough to avoid aliasing problem) 3. Compute z-transform of resulting sequence 8 Example 1 – Impulse Invariant Method Consider first order analogue filter = 1 - Corresponding impulse response is δ - v The presence of delta term prevents sampling of impulse response which thus cannot be defined Fundamental problem: high-pass and band-stop filters have functions with numerator and denominator polynomials of the same degree and thus cannot be designed using this method 9 Example 2 – Impulse Invariant Method Consider an analogue filter Step 1.
    [Show full text]
  • (73) Assignee: Kind of Loud Technologies, LLC, E.R.E
    US007831645B1 (12) United States Patent (10) Patent No.: US 7,831,645 B1 Berners et al. (45) Date of Patent: Nov. 9, 2010 (54) DIGITAL RESONANT SHELF FILTER OTHER PUBLICATIONS (75)75 Inventors: E.id P. s "E. S. R. SEGs) Azizi.Parametric Seyed-Ali. and Graphic “A New Equalizer Concept Banks'.of Interference Audio CompensationEngineering Soc., for Onanan S. ADel, Menlo Park, Convention Paper 5482, Sep. 21-24, 2001, New York, NY. pp. 1-8. (73) Assignee: Kind of Loud Technologies, LLC, Azizi,E.R.E. Seyed-Ali, “A New C t of Interf C tion f Santa Cruz, CA (US) Convention Paper 5629, May 10-13, 2002, Munich, Germany, pp. 1-8. (*) Notice: Subject to any disclaimer, the term of this Berners, David P. et al., “Discrete-Time Shelf Filter Design for patent is extended or adjusted under 35 Analog Modeling”. Audio Engineering Soc., Convention Paper, Oct. U.S.C. 154(b) by 1274 days. 11-13, 2003, New York, NY. pp. 1-5. Blesser, Barry A., et al., “Audio Dynamic Range Compression for (21) Appl. No.: 11/249,162 Minimum Perceived Distortion', IEEE Transactions on Audio and Electroacoustics, vol. AU-17, No. 1, Mar. 1969, pp. 22-32. (22) Filed: Oct. 11, 2005 (Continued) Primary Examiner Tan V Mai Related U.S. Application Data (74) Attorney, Agent, or Firm Pillsbury Winthrop Shaw (60) Provisional application No. 60/617,343, filed on Oct. Pittman LLP 8, 2004. s (57) ABSTRACT (51) Int. Cl. A method and system for designing a discrete-time filter G06 L/00 (2006.01) having a transfer function which approximates that of an G06F 7/10 (2006.01) - 0 52) U.S.
    [Show full text]