Peaking and Shelving Filter Properties
|
|
- Benjamin Walters
- 6 years ago
- Views:
Transcription
1 MUS424: Signal Processing Techniques for Digital Audio Effects Handout #22 Jonathan Abel, David Berners May 18, 24 Lecture #14: May 18, 24 Lecture Notes 12 Peaking and Shelving Filter Properties Shelf Filter Specification A shelf filter h(ω) takes on a gain of l at DC, a gain of l π at high frequencies, and a gain of (l l π ) 1 2 at a transition frequency ϕ. A low shelf filter has the high-frequency gain l π fixed at one, with l free to vary. In a high shelf filter, the high-frequency gain l π is varied with the DC gain l fixed at one.
2 2 MUS424: Handout #22 First-Order Shelf Filter Analog Prototype Consider the following first-order analog prototype filter, h(s) = l π s/ρ + l. s/ρ + 1 Note that the analog prototype filter takes on the desired DC and high-frequency gains, h(ω = ) = l, h(ω ) = l π, with the factor ρ controlling the transition frequency. The analog prototype filter is designed to have a transition frequency of one. Setting h(j 1) 2 = l l π, the factor ρ may be specified. h(j 1) 2 = l l π = l2 π/ρ 2 + l 2 1/ρ 2 + 1, which implies ρ = ( lπ l )1 2. Substituting into the expression for h(s), we have h(s) = (l l π ) 1 2 s + (l /l π ) 1 2 (l /l π ) 1 2s + 1.
3 MUS424: Handout #22 3 First-Order Shelf Filter Analog Prototype Properties The analog prototype shelf filter has transfer function h(s) = (l l π ) 1 2 s + (l /l π ) 1 2 (l /l π ) 1 2s + 1. It has a real pole and real zero at reciprical frequencies about the transition frequency s = j 1, s p = (l π /l ) 1 2, s z = (l /l π ) 1 2. Note that on a db scale the normalized analog shelf filter is antisymmetric in log s about the transition frequency, h(1/s) 1/s + (l /l π ) 1 2 (l /l π ) 1 2/s + 1 = 1 + (l /l π ) 1 2s (l /l π ) s 1/h(s).
4 4 MUS424: Handout #22 First-Order Digital Shelf Filter via Bilinear Transform The digital shelf filter is formed via bilinear transform on the analog prototype, We have and, h(z) = h(s(z)), s = 2 T 1 z z 1. h(s) = (l l π ) 1 2 s + 1/ρ s/ρ + 1, ρ = (l π/l ) 1 2, h(z) = (l l π ) 1 2 ρ(1 z 1 ) + T 2 (1 + z 1 ) (1 z 1 ) + ρ T 2 (1 + z 1 ), = (l l π ) 1 ( T ρ) + (T 2 ρ)z 1 (ρ T 2 + 1) + (ρt 2 1)z 1, [ ] ( ) = (l l π ) 1 ρ + T 1 1 T/2ρ z T/2ρ 1 + ρ T. 2 1 )z 1 ( 1 ρt/2 1+ρT/2 Note that the pole at ρ and zero at 1/ρ have been transformed according to the bilinear transform.
5 MUS424: Handout #22 5 Selecting T The sampling period T is chosen to put the transition frequency of the analog prototype, Ω = 1, at the desired frequency ϕ. Evaluating the bilinear transform on the unit circle, z 1 = e jω, s = 2 T 1 z z 1 = 2 T 1 e jω 1 + e jω, = 2 T ejω/2 e jω/2 e jω/2 + e jω/2 = 2 T j sin(ω/2) cos(ω/2), we see that the frequency ω on the unit circle corresponds to the frequency Ω = 2 T tan(ω/2) on the s-plane imaginary axis. (Note that for small frequencies ω 1, tan(ω/2) ω/2, and Ω 1 T. To map the analog prototype filter transition frequency of 1 to the desired transition frequency ϕ, we set T 2 = tan(ϕ/2).
6 6 MUS424: Handout #22 Example Shelf Filters 2 second-order shelf filter transfer function magnitudes, ft = 125 * 2.^[:6] second-order shelf filter transfer function magnitudes, gain = [2:2:2]
7 MUS424: Handout #22 7 Peak Filter Specification A peak (or notch) filter p(ω) takes on a gain of 1 at DC and high frequencies, and achieves a maximum (or minimum) gain of l at some point between the transition frequencies ϕ ±, at which the gain is l. The frequency of the magnitude extremum is called the center frequency, ϕ c.
8 8 MUS424: Handout #22 Digital Peak Filter The second-order digital filter with coefficients given by p(z) = b + b 1 z 1 + b 2 z a 1 z 1 + a 2 z 2 a 2 = 2Q sin ϕ c, 2Q + sinϕ c a 1 = b 1 = (1 + a 2 ) cos ϕ c, b = 1 2 (1 + a 2) (1 a 2)l, b 2 = 1 2 (1 + a 2) 1 2 (1 a 2)l, implements a peak (or notch) filter with maximum (or minimum) gain l at a center frequency ϕ c between the specified transition frequencies ϕ ±, at which the filter takes on magnitude l. The center frequency ϕ c and the inverse bandwidth Q may be written in terms of the transition frequencies ϕ ± and db peak gain λ, { = acos κ sign{κ} ( } κ 2 1 )1 2, ϕ c κ = 1 + cos ϕ cos ϕ + cos ϕ + cosϕ + Q = 1 [ l sin 2 ϕ c (cos ϕ + cos ϕ + ) 2 2 cos ϕ c cosϕ cos ϕ + In the case that ϕ + + ϕ = π, we have ] 1 2. ϕ c = π/2, l Q = 2 cotδ, δ = 1 2 (ϕ ϕ + ).
9 MUS424: Handout #22 9 Example Peak Filters 2 peak filter transfer function magnitudes, ft = 125 * 2.^[:6] peak filter transfer function magnitudes, gain = [2:2:2]
10 1 MUS424: Handout #22 Graphic Equalizer 6 graphic equalizer shelf and peak filters A graphic equalizer is implemented as a cascade of peak and shelf filters having a prioi specified bandwidths and user-controlled gains.
11 MUS424: Handout #22 11 Graphic Equalizer Behavior 1 traditional graphic equalizer with controls set to [ 5 5 ] db By adjusting the gains of the peak and shelf filters in the cascade, a wide variety of transfer function magnitudes may be produced. Although the idea is that the transfer function magnitude should smoothly interpolate the specified gains, it doesn t always work out that way.
12 12 MUS424: Handout #22 Graphic Equalizer Behavior 1 traditional graphic equalizer with controls set to 5*ones(1,8) db Contributions from adjacent bands cause the filter to overshoot the desired gains.
13 MUS424: Handout #22 13 Graphic Equalizer Behavior 1 traditional graphic equalizer with narrow-band peak filters Traditional graphic equalizers don t smoothly interpolate the given gains. Small filter bandwidths lead to rippled transfer functions. Large filter bandwidths overshoot desired gains.
14 14 MUS424: Handout #22 Alternative Gain Computation 1 novel graphic equalizer The idea is to find a set of peak and shelf filter gains which account for the overlap between bands, so that the resulting peak and shelf filter cascade interpolates the specified band gains.
15 MUS424: Handout #22 15 Peak Filter Self Similarity 2 peak filter transfer function magnitudes, gain = [2:2:2] normalized peak filter transfer function magnitudes, gain = [2:2:2] Peaking filters parameterized by a maximum db gain λ, achieved somewhere between transition frequencies ϕ and ϕ +, at which the db gain is λ/2, are approximately self similar on a log magnitude scale. α log p(ω; λ,ϕ ± ) log p(ω;α λ,ϕ ± ).
16 16 MUS424: Handout #22 Shelf Filter Self Similarity 2 second-order shelf filter transfer function magnitudes, gain = [2:2:2] normalized shelf filter transfer function magnitudes, gain = [2:2:2] Low shelf filters specified by a low-frequency db gain λ, a db gain λ/2 at a transition frequency ϕ, and a high-frequency gain of one are also approximately self similar on a log magnitude scale. α log h(ω;λ, ϕ) log h(ω;α λ,ϕ).
17 MUS424: Handout #22 17 Filter Design Approach Consider a cascade g(ω; θ) of K peak and shelf filters having gains and λ k, k = 1,...,K and transition frequencies ϕ k, k = 1,...,K 1 stacked in the column θ, g(ω; θ) = h(ω; λ 1,ϕ 1 ) h(ω;λ K,ϕ K 1 ) K 1 p(ω; λ k,ϕ k 1, ϕ k ). k=2 Because of the self similarity property, the db magnitude of the cascade, denoted by γ(ω;θ), γ(ω;θ) = def 2log 1 {g(ω;θ)}, = σ(ω; λ 1,ϕ 1 ) + σ(ω; λ K,ϕ K 1 ) + K 1 π(ω;λ k,ϕ k 1,ϕ k ), k=2 is approximately linear in the filter gains, γ(ω; θ) λ 1 σ(ω; 1.,ϕ 1 ) + λ K σ(ω; 1.,ϕ K 1 ) + K 1 λ k π(ω; 1., ϕ k 1,ϕ k ). k=2
18 18 MUS424: Handout #22 Filter Design Approach At a particular frequency ω i, we have γ(ω i ; θ) [ σ 1 (ω i ) π 2 (ω i ) π K 1 (ω i ) σ K (ω i ) ] λ, where λ is the stack of db band gains, λ = [ λ 1 λ K ], and where σ k (ω i ) and π(ω i ) are the transfer function db magnitudes of shelf and peak filters with specified transition frequencies and 1. db nominal gains, evaluated at ω i, and σ k (ω i ) = 2 log 1 {h(ω i ; 1.dB,ϕ k )} π(ω i ) = 2 log 1 {p(ω i ; 1.dB,ϕ k 1,ϕ k )}. Stacking instances of γ(ω; θ) evaluated at a set of frequencies ω i to form the column γ, we have γ Bλ, B = [ σ 1 π 2 π K 1 σ K ], where σ k and π k are columns of db magnitudes evaluated at ω i of shelf and peak filters having gains of 1. db and specified transition frequencies.
19 MUS424: Handout #22 19 Graphic Equalizer Design The self similarity of second-order peaking and shelving filters leads to an approximate linear relationship between the db filter gains and the db gain of the cascade, γ Bλ, Therefore, given a set of shelf and peak filters having specified transition frequencies, and positive definite weighting matrix W, the gains ˆλ = (B WB) 1 B Wγ will approximately minimize the weighted square difference between a desired db magnitude response γ and the shelf and peak filter cascade db magnitude at frequencies ω i, γ. For a graphic equalizer with K 1 fixed band edges, the frequencies ω i can be chosen as the K band centers, and the gains ˆλ simply computed as the control gains γ scaled by the basis inverse, ˆλ = B 1 γ.
20 2 MUS424: Handout #22 Graphic Equalizer Design Example traditional and novel graphic equalizers For a graphic equalizer with K 1 fixed band edges, the frequencies ω i can be chosen as the K band centers, and the gains ˆλ simply computed as the control gains γ scaled by the basis inverse, ˆλ = B 1 γ. Note that for a grphic equalizer with fixed band edges, B 1 may be computed a priori. To account for discrepancies in the self similarity property, ˆλ may be computed iteratively, forming B using the gains from the previous solution.
21 MUS424: Handout #22 21 Iterated Solution graphic equalizer with iterated band gains Iterative solution of the equation for the corrected band gains can correct for discrepancies in filter self similarity, λ n+1 = B(λ n ) 1 γ. It seems likely that the iteration converges since the sensitivity of the basis to changes in the gain are small, and changes in the basis B are monotonic with changes in the gains λ.
22 22 MUS424: Handout #22 Transfer Function Modeling measured transfer function magnitude, selected extrema A cascade of peak and shelf filters can be fit to an arbitrary transfer function by selecting a set of transition frequencies, and fitting the gains. Transition frequency selection method. Tabulate extrema or inflection point frequencies. Pick transition frequencies as geometric means of significant extrema or at significant inflection points.
23 MUS424: Handout #22 23 Transfer Function Modeling measured (solid) and modeled (dashed) transfer function magnitudes Band gain selection method. Pick a dense sampling of frequencies ω i, say Bark or ERB spaced. Form basis, and compute gains, ˆλ = (B WB) 1 B Wγ. The resulting transfer function will not interpolate the specified gains γ, but rather approximate the desired magnitude, minimizing the weighted mean square db difference.
24 24 MUS424: Handout #22 Critical Bandwidth and Critical Band Smoothing 1 original and critical-band smoothed power spectra When applied to broad band signals, details of a filter smaller than a critical bandwidth are irrelevant only the total power matters.
25 MUS424: Handout #22 25 Critical Bandwidth The Bark frequency scale divides the human hearing range into about 25 ciritcal bands. Frequency f in khz and critical-band rate b in Bark are related via the following expressions, b = 1.3 [ log(1 + f 2 ) ] 1 2, f = [ exp { (b/1.3) 2} 1 ] 1 2.
26 26 MUS424: Handout #22 Critical Bandwidth Critical bandwidth as a function of frequency is approximately w = f 3 2.
27 MUS424: Handout #22 27 Critical Band Smoothing 3 original and smoothed power spectra By smoothing the filter power over bandwidths proportional to critical bandwidth, filter complexity is reduced while retaining psychoacoustically relevant cues. Critical-band smoothing may be done via a running mean over frequency, or via windowing of a frequencywarped impulse response.
28 28 MUS424: Handout #22 Critical Band Smoothing 7 original and smoothed impulse responses time - milliseconds As expected, impulse responses are shortened by criticalband smoothing.
29 MUS424: Handout #22 29 Filter Phase 3 raw, minimum phase and linear phase impulse responses time - milliseconds Humans are relatively insensitive to filter phase. A linear phase filter delays all frequency components the same amount τ, H(ω) = e jτω H(ω). A minimum phase filter delays input energy the least of all filters having the same magnitude response, m m h 2 min (n) h 2 (n), m. n= n=
PS403 - 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 informationReverberation Impulse Response Analysis
MUS424: Signal Processing Techniques for Digital Audio Effects Handout #18 Jonathan Abel, David Berners April 27, 24 Lecture #9: April 27, 24 Lecture Notes 7 Reverberation Impulse Response Analysis 1 CCRMA
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 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 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 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 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 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 informationPoles and Zeros in z-plane
M58 Mixed Signal Processors page of 6 Poles and Zeros in z-plane z-plane Response of discrete-time system (i.e. digital filter at a particular frequency ω is determined by the distance between its poles
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 informationConvention Paper Presented at the 125th Convention 2008 October 2 5 San Francisco, CA, USA
Audio Engineering Society Convention Paper Presented at the 125th Convention 2008 October 2 5 San Francisco, CA, USA The papers at this Convention have been selected on the basis of a submitted abstract
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 information2.161 Signal Processing: Continuous and Discrete
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete Fall 28 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
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 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 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 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 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 information2.161 Signal Processing: Continuous and Discrete Fall 2008
IT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete all 2008 or information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. assachusetts
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 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 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 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 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 informationCITY UNIVERSITY LONDON. MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746
No: CITY UNIVERSITY LONDON MSc in Information Engineering DIGITAL SIGNAL PROCESSING EPM746 Date: 19 May 2004 Time: 09:00-11:00 Attempt Three out of FIVE questions, at least One question from PART B PART
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 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 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 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 informationLecture 04: Discrete Frequency Domain Analysis (z-transform)
Lecture 04: Discrete Frequency Domain Analysis (z-transform) John Chiverton School of Information Technology Mae Fah Luang University 1st Semester 2009/ 2552 Outline Overview Lecture Contents Introduction
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
MIT OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. M MASSACHUSETTS
More informationEE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley
University of California Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences Professors : N.Morgan / B.Gold EE225D Digital Filters Spring,1999 Lecture 7 N.MORGAN
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 informationComputer-Aided Design of Digital Filters. Digital Filters. Digital Filters. Digital Filters. Design of Equiripple Linear-Phase FIR Filters
Computer-Aided Design of Digital Filters The FIR filter design techniques discussed so far can be easily implemented on a computer In addition, there are a number of FIR filter design algorithms that rely
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 informationCourse and Wavelets and Filter Banks. Filter Banks (contd.): perfect reconstruction; halfband filters and possible factorizations.
Course 18.327 and 1.130 Wavelets and Filter Banks Filter Banks (contd.): perfect reconstruction; halfband filters and possible factorizations. Product Filter Example: Product filter of degree 6 P 0 (z)
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 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 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 informationCMPT 889: Lecture 9 Wind Instrument Modelling
CMPT 889: Lecture 9 Wind Instrument Modelling Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University November 20, 2006 1 Scattering Scattering is a phenomenon in which the
More informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
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 informationAutomatic Control (MSc in Mechanical Engineering) Lecturer: Andrea Zanchettin Date: Student ID number... Signature...
Automatic Control (MSc in Mechanical Engineering) Lecturer: Andrea Zanchettin Date: 29..23 Given and family names......................solutions...................... Student ID number..........................
More informationSignals, Instruments, and Systems W5. Introduction to Signal Processing Sampling, Reconstruction, and Filters
Signals, Instruments, and Systems W5 Introduction to Signal Processing Sampling, Reconstruction, and Filters Acknowledgments Recapitulation of Key Concepts from the Last Lecture Dirac delta function (
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 information/ (2π) X(e jω ) dω. 4. An 8 point sequence is given by x(n) = {2,2,2,2,1,1,1,1}. Compute 8 point DFT of x(n) by
Code No: RR320402 Set No. 1 III B.Tech II Semester Regular Examinations, Apr/May 2006 DIGITAL SIGNAL PROCESSING ( Common to Electronics & Communication Engineering, Electronics & Instrumentation Engineering,
More informationLecture 8 Finite Impulse Response Filters
Lecture 8 Finite Impulse Response Filters Outline 8. Finite Impulse Response Filters.......................... 8. oving Average Filter............................... 8.. Phase response...............................
More informationLectures on APPLICATIONS
APP0 University of alifornia Berkeley ollege of Engineering Department of Electrical Engineering and omputer Science obert W. Brodersen EES40 Analog ircuit Design t ISE Lectures on APPLIATINS t FALL.0V
More informationLecture 14: Minimum Phase Systems and Linear Phase
EE518 Digital Signal Processing University of Washington Autumn 2001 Dept. of Electrical Engineering Lecture 14: Minimum Phase Systems and Linear Phase Nov 19, 2001 Prof: J. Bilmes
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 information( ) John A. Quinn Lecture. ESE 531: Digital Signal Processing. Lecture Outline. Frequency Response of LTI System. Example: Zero on Real Axis
John A. Quinn Lecture ESE 531: Digital Signal Processing Lec 15: March 21, 2017 Review, Generalized Linear Phase Systems Penn ESE 531 Spring 2017 Khanna Lecture Outline!!! 2 Frequency Response of LTI System
More informationDigital Signal Processing. Midterm 1 Solution
EE 123 University of California, Berkeley Anant Sahai February 15, 27 Digital Signal Processing Instructions Midterm 1 Solution Total time allowed for the exam is 8 minutes Some useful formulas: Discrete
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. Massachusetts
More information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
More informationRational Implementation of Distributed Delay Using Extended Bilinear Transformations
Rational Implementation of Distributed Delay Using Extended Bilinear Transformations Qing-Chang Zhong zhongqc@ieee.org, http://come.to/zhongqc School of Electronics University of Glamorgan United Kingdom
More informationLecture 9 Time-domain properties of convolution systems
EE 12 spring 21-22 Handout #18 Lecture 9 Time-domain properties of convolution systems impulse response step response fading memory DC gain peak gain stability 9 1 Impulse response if u = δ we have y(t)
More informationLinear Prediction 1 / 41
Linear Prediction 1 / 41 A map of speech signal processing Natural signals Models Artificial signals Inference Speech synthesis Hidden Markov Inference Homomorphic processing Dereverberation, Deconvolution
More informationBridge between continuous time and discrete time signals
6 Sampling Bridge between continuous time and discrete time signals Sampling theorem complete representation of a continuous time signal by its samples Samplingandreconstruction implementcontinuous timesystems
More information2A1H Time-Frequency Analysis II
2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.11: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 21 QUESTION BOOKLET
More informationECGR4124 Digital Signal Processing Exam 2 Spring 2017
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,
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 informationNAME: ht () 1 2π. Hj0 ( ) dω Find the value of BW for the system having the following impulse response.
University of California at Berkeley Department of Electrical Engineering and Computer Sciences Professor J. M. Kahn, EECS 120, Fall 1998 Final Examination, Wednesday, December 16, 1998, 5-8 pm NAME: 1.
More informationReview: transient and steady-state response; DC gain and the FVT Today s topic: system-modeling diagrams; prototype 2nd-order system
Plan of the Lecture Review: transient and steady-state response; DC gain and the FVT Today s topic: system-modeling diagrams; prototype 2nd-order system Plan of the Lecture Review: transient and steady-state
More informationSolutions. Number of Problems: 10
Final Exam February 4th, 01 Signals & Systems (151-0575-01) Prof. R. D Andrea Solutions Exam Duration: 150 minutes Number of Problems: 10 Permitted aids: One double-sided A4 sheet. Questions can be answered
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
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 informationE : Lecture 1 Introduction
E85.2607: Lecture 1 Introduction 1 Administrivia 2 DSP review 3 Fun with Matlab E85.2607: Lecture 1 Introduction 2010-01-21 1 / 24 Course overview Advanced Digital Signal Theory Design, analysis, and implementation
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 informationLOPE3202: Communication Systems 10/18/2017 2
By Lecturer Ahmed Wael Academic Year 2017-2018 LOPE3202: Communication Systems 10/18/2017 We need tools to build any communication system. Mathematics is our premium tool to do work with signals and systems.
More informationSolutions: Homework Set # 5
Signal Processing for Communications EPFL Winter Semester 2007/2008 Prof. Suhas Diggavi Handout # 22, Tuesday, November, 2007 Solutions: Homework Set # 5 Problem (a) Since h [n] = 0, we have (b) We can
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 informationChapter 9. Linear Predictive Analysis of Speech Signals 语音信号的线性预测分析
Chapter 9 Linear Predictive Analysis of Speech Signals 语音信号的线性预测分析 1 LPC Methods LPC methods are the most widely used in speech coding, speech synthesis, speech recognition, speaker recognition and verification
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 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 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 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 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 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 informationThis homework will not be collected or graded. It is intended to help you practice for the final exam. Solutions will be posted.
6.003 Homework #14 This homework will not be collected or graded. It is intended to help you practice for the final exam. Solutions will be posted. Problems 1. Neural signals The following figure illustrates
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 informationDiscrete Time Systems
1 Discrete Time Systems {x[0], x[1], x[2], } H {y[0], y[1], y[2], } Example: y[n] = 2x[n] + 3x[n-1] + 4x[n-2] 2 FIR and IIR Systems FIR: Finite Impulse Response -- non-recursive y[n] = 2x[n] + 3x[n-1]
More informationModeling Buck Converter by Using Fourier Analysis
PIERS ONLINE, VOL. 6, NO. 8, 2010 705 Modeling Buck Converter by Using Fourier Analysis Mao Zhang 1, Weiping Zhang 2, and Zheng Zhang 2 1 School of Computing, Engineering and Physical Sciences, University
More informationCMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation
CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is
More informationProblem Set 2: Solution Due on Wed. 25th Sept. Fall 2013
EE 561: Digital Control Systems Problem Set 2: Solution Due on Wed 25th Sept Fall 2013 Problem 1 Check the following for (internal) stability [Hint: Analyze the characteristic equation] (a) u k = 05u k
More informationDamped Oscillators (revisited)
Damped Oscillators (revisited) We saw that damped oscillators can be modeled using a recursive filter with two coefficients and no feedforward components: Y(k) = - a(1)*y(k-1) - a(2)*y(k-2) We derived
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 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 informationECE 8440 Unit 17. Parks- McClellan Algorithm
Parks- McClellan Algorithm ECE 8440 Unit 17 The Parks- McClellan Algorithm is a computer method to find the unit sample response h(n for an op=mum FIR filter that sa=sfies the condi=ons of the Alterna=on
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 informationDiscrete Time Signals and Systems Time-frequency Analysis. Gloria Menegaz
Discrete Time Signals and Systems Time-frequency Analysis Gloria Menegaz Time-frequency Analysis Fourier transform (1D and 2D) Reference textbook: Discrete time signal processing, A.W. Oppenheim and R.W.
More informationDigital Signal Processing 2/ Advanced Digital Signal Processing, Audio/Video Signal Processing Lecture 10, Frequency Warping, Example
Digital Signal Processing 2/ Advanced Digital Signal Processing, Audio/Video Signal Processing Lecture 10, Gerald Schuller, TU Ilmenau Frequency Warping, Example Example: Design a warped low pass filter
More informationFourier Series Example
Fourier Series Example Let us compute the Fourier series for the function on the interval [ π,π]. f(x) = x f is an odd function, so the a n are zero, and thus the Fourier series will be of the form f(x)
More informationOp-Amp Circuits: Part 3
Op-Amp Circuits: Part 3 M. B. Patil mbpatil@ee.iitb.ac.in www.ee.iitb.ac.in/~sequel Department of Electrical Engineering Indian Institute of Technology Bombay Introduction to filters Consider v(t) = v
More informationCHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME
CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME Shri Mata Vaishno Devi University, (SMVDU), 2013 Page 13 CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME When characterizing or modeling a random variable, estimates
More informationECE 598: The Speech Chain. Lecture 5: Room Acoustics; Filters
ECE 598: The Speech Chain Lecture 5: Room Acoustics; Filters Today Room = A Source of Echoes Echo = Delayed, Scaled Copy Addition and Subtraction of Scaled Cosines Frequency Response Impulse Response Filter
More informationReview of Fundamentals of Digital Signal Processing
Chapter 2 Review of Fundamentals of Digital Signal Processing 2.1 (a) This system is not linear (the constant term makes it non linear) but is shift-invariant (b) This system is linear but not shift-invariant
More informationRecursive, Infinite Impulse Response (IIR) Digital Filters:
Recursive, Infinite Impulse Response (IIR) Digital Filters: Filters defined by Laplace Domain transfer functions (analog devices) can be easily converted to Z domain transfer functions (digital, sampled
More informationLecture 5 - Assembly Programming(II), Intro to Digital Filters
GoBack Lecture 5 - Assembly Programming(II), Intro to Digital Filters James Barnes (James.Barnes@colostate.edu) Spring 2009 Colorado State University Dept of Electrical and Computer Engineering ECE423
More information