IIR digital filter design for low pass filter based on impulse invariance and bilinear transformation methods using butterworth analog filter
|
|
- Alannah Little
- 6 years ago
- Views:
Transcription
1 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 08:4 PM Contents Introduction. Filter specifications he design steps. Impulse invariance method bilinear transformation method Summary of analytical derivation method Numerical design examples 9 3. Example using impulse invariance method using = bilinear method Example using impulse invariance method Using bilinear IIR design for minmum order filter 3 4. Impulse invariance method bilinear transformation method References 7 A low pass digital filter IIR is designed based on given specifications. Butterworth analog filter Hs is designed first, then it is converted to digital filter Hz Analytical procedure is illustrated below and simplified to allow one to more easily program the algorithm. 4 numerical examples are used to illustrate the procedure. Introduction his report derives a symbolic procedure to design a low pass IIR digital filter from an analog Butterworth filter using methods: impulse invariance and bilinear transformation. wo numerical examples are then used to illustrate using the symbolic procedure. here are a total of 3 steps. A Mathematica program with GUI is written to enable one to use this design for different parameters. A Matlab script is being written as well.
2 . Filter specifications Filter specifications are 5 parameters. he frequency specifications are analog frequencies, while the attenuation for the passband and the stopband are given in db f c f s δ p δ s he sampling frequency in Hz he passband cutoff frequency in Hz he stopband corner frequency in Hz he attenuation in db at f c he attenuation in db at f s his diagram below illustrates these specifications he frequency specifications are in Hz and they must be first converted to digital frequencies ω where 0 ω π before using the attenuation specifications, he sampling frequency is used to do this conversion since corresponds to π on the digital frequency scale. he design steps. Convert specifications from analog to digital frequencies. Based on design method impulse invariance of bilinear apply the attenuation criteria to determine and N the filter order 3. Using and N find the locations of the poles of the butterworth analog filter Hs. 4. Find Hz from Hs. he method of doing this depends if we are using impulse invariance or bilinear. his step is much simpler for the bilinear method as it does not require performing partial fractions decomposition on Hs Now we begin the analytical design procedure.. Impulse invariance method F We first convert analog specifications to digital specifications: s π = fp ω p, hence ω p = π fp Convert the criteria relative to the digital normalized scale: herefore 0 log H e jωp δ p 0 log H e jωs δ s and ω s = π fs
3 H e jωp δp 0 A H e jω s δs 0 B Notice that we added the scaling factor. Now, Butterworth analog filter squared magnitude Fourier transform is given by H a jω = + jω j hence equations A and B above are now written in terms of the analog Butterworth amplitude frequency response and become herefore the above becomes Ωp + + Ω s δp 0 δs 0 Ωp + + Ω s δp δs Now, for impulse invariance, Ω p = ωp and Ω s = ωs. Change inequalities to equalities and simplify Divide the above equations ωp / + ωs / + ωp / ωs / δp δs = δp = δs δp ωp = ω s δs [log ω p log ω s ] = log δp log δs [ ] [ ] N = log δp log δs log ω p log ω s 3
4 We need to round the above to then nearest integer using the Ceiling function N = N Now for impulse invariance method, use equation above to solve for ωp / log ω p / = δp log ω p / = = log δp ω p / = = log log δp ω p / δp log δp Now that we found N and, next find the poles of H s Since the butterworth magnitude square of the transfer function is H a s = + s j Hence H s poles are found by setting the denominator of the above to zero s + = 0 j s = j s j = e j = e jπ+πk k = 0,,, π+πk s = j e j π+πk = e j π e j π+πk = e j π+k+n We only need to find the LHS poles, which are located at i = 0 N, because these are the stable poles. Hence the i th pole is s i = e j π+i+n For example for i = 0, N = 6 we get s 0 = e j π+n = e j π7 Now we can write the analog filter generated based on the above selected poles, which is, for impulse invariance H a s = K s s i 3 K is found by solving H a 0 = hence k = s i 4
5 Now we need to write poles in non-polar form and plug them into 3 s i = e j π+i+n π + i + N π + i + N = cos + j sin i = 0 N Hence, Where and H a s = N = K s cos π+i+n = ω p / log δp + j sin π+i+n [ ] [ log δp log δs log ω p log ω s ] 4 Now that we have found H s we need to convert it to H z We need to make sure that we multiply poles of complex conjugates with each others to make the result simple to see. Now that we have H a s, we do the A D conversion. I.e. obtain H z from the above H s. When using impulse invariance, we need to perform partial fraction decomposition on 4 above in order to write H s in this form For example, to obtain A j, we write H s = A j = lim s sj H a s = Once we find all the A s, we now write H z as follows H z = s s i i j k s s i exp s i z his completes the design. We can try to convert the above form of H z to a rational expression as H z = Nz Dz. bilinear transformation method F We first convert analog specifications to digital specifications: s π = fp ω p, hence ω p = π fp Convert the criteria relative to the digital normalized scale: 0 log H e jωp δ p 0 log H e jω s δs and ω s = π fs Hence H e jωp δp 0 A H e jωs δs 0 B 5
6 Butterworth analog filter squared magnitude Fourier transform is given by H a jω = + jω j hence equations A and B above are now written in terms of the analog Butterworth amplitude frequency response and become Ωp + Ω + s δp δp 0 = δs 0 = δs Now we assign values for Ω p and Ω s as follows,ω p = tan ω p, Ωs = tan ω s + + tan ω p tan ω s Change inequalities to equalities and simplify Divide the above equations [ log tan ωp log tan ω p tan ω s ωp tan tan ω s tan ωs δp δs = δp = δs δp = δs ] = log δp log δs N = log δp log δs log tan ω p log tan ωs We need to round the above to then nearest integer using the Ceiling function. i.e. N = N 6
7 Now for bilinear transformation we used equation above to find Hence we now solve for + tan ω s log tan ω s = δs log tan ω s = tan ω s = = = log δs log log δs δs tan ω s log δs Now that we found N and we find the poles of H s. Since for bilinear the magnitude square of the transfer function is H a s = + s j Hence H s poles are found by setting the denominator of the above to zero s + = 0 j s = j s j = e j = e jπ+πk k = 0,,, π+πk s = j e j π+πk = e j π e j π+πk = e j π+k+n We only need to find the LHS poles, which are located at i = 0 N, because these are the stable poles. Hence the i th pole is s i = e j π+i+n For example for i = 0, N = 6 we get s 0 = e j π+n = e j π7 For bilinear, Hs is given by H a s = K is found by solving H a 0 =, hence we obtain K s s i 3 k = 7 s i
8 We see that the same expression results for k for both cases. Now we need to write poles in non-polar form and plug them into 3 s i = e j π+i+n π + i + N = cos + j sin π + i + N i = 0 N then Where and H a s = N = K s cos π+i+n = log δp log tan ωp tan ω s log δs log + j sin π+i+n δs log tan ωs 4 Now that we have found H s we need to convert it to H z. After finding Hs as shown above, we simply replace s by z +z. his is much simpler than the impulse invariance method. Before doing this substitution, make sure to multiply poles which are complex conjugate of each others in the denominator of Hs. After this, then do the above substitution 8
9 .3 Summary of analytical derivation method We will now make a table with the derivation equations to follow to design in either bilinear or impulse invariance. Note that the same steps are used in both designs except for step 5, 6, 8, 3. his table make it easy to develop a program. step Impulse invariance common equation bilinear ω p = π fp ω s = π fs 3 α p = δp 4 α s = 5 Ω p = ωp 6 Ω s = ωs 7 N = 8 = Ω p δs log[α p ] log[α s ] logω p logω s log [ αp ] = 9 poles of HS s i = e j π+i+n k = K H a s = s s i do partial fractions: H a s = 3 H z = s s i s i tan ω p tan ω s H a s = Ω s log [αs ] K s s i exp s iz H z = H a s s= z +z 3 Numerical design examples 3. Example Sampling frequency = 0khz, passband frequency f p = khz, stopband frequency f s = 3khz, with δ p db and δ stop 5db 9
10 3.. using impulse invariance method using = step Impulse invariance ω p = π fp ω s = π fs 3 α p = δp 4 α s = 5 Ω p = ωp 6 Ω s = ωs 7 N = 8 = π π π π δs 5 ωp 0.π 0.3π 0.3π log[α p ] log[α s ] logω p logω s Ω p [ ] 0.π log αp log log log 0.π log 0.3π log poles of HS s i = e j π+i+n s i = e j π7+i i = 0 5 s 0 = j , s = j , s = j0.8 s 3 = j0.8, s 4 = j , s 5 = 0.8 j k = s i 0.8 j j j j j j K H a s = s s i 0. 9 s+0.8 j s j s j0.8 s j0.8s j s+0.8 +j multiply complex conjugates s s s s s s s s s s s s+0. 9 A partial fraction H a s = i j j j. 607 s s i s+0.8 j s j s j j j j s j0.8 s j s+0.8 +j j exp s iz exp 0.8 +j z.07 4 j j.607 exp j z + exp j H z = j j.66 8 exp j0.8 z j exp j z + exp 0.8 j z j j j j z j0.90 z j z j j j j z j0.90z j0.5368z H z = +
11 3.. bilinear method Using the above design table, these are the numerical values: = = 0000 step Bilinear ω p = π fp ω s = π fs 3 α p = δp π π π π α s = δs 5 5 Ω p = tan ω p 0000 tan 0.π Ω s = tan ω s 0000 tan 0.3π 038 log[α 7 N = p ] log[α s ] logω p logω s log log log 997 log Ω 8 = s log log [αs ] 9 poles of HS s i = e j π+i+n s i = 535 e j π7+i i = 0 5 s 0 = j4803, s = j836., s = j s 3 = 4803 j3966.4, s 4 = 836 j836., s 5 = j4803 k = s i j j j j j j4803 multiply complex conjugate terms K H a s = s s i s j4803s+836 j836s+4803 j3966.4s+4803+j s+836+j836s j4803 multiply complex conjugate s s s + 67s s s s s s = s s s +67.s s s s s s s s s H z = H a s s= z ODO +z 3. Example Sampling frequency = khz, passband corner frequency f p = khz, stopband corner frequency f s = khz, with criteria δ p 3db and δ stop db
12 3.. using impulse invariance method = step Impulse invariance ω p = π fp ω s = π fs 3 α p = δp 4 α s = 5 Ω p = ωp 6 Ω s = ωs 7 N = 8 = π π π π δs ωp 0.π 0.4π 0.4π log[α p ] log[α s ] logω p logω s Ω p [ ] 0.π log αp log log.0 log 0.π log 0.4π log poles of HS s i = e j π+i+n s i = e j π3+i 4 s 0 = j , s = j k = K H a s = i = 0 s i j j s s i partial fraction H a s = z z. 57 z H z = s j s j s s s s i j s j j s j exp s iz i exp j z i exp j z
13 3.. Using bilinear = 000 step Impulse invariance ω p = π fp ω s = π fs 3 α p = δp π π π π α s =.0 δs 5 Ω p = tan ω p 000 tan 0.π Ω s = tan ω s 000 tan 0.4π 453. log[α 7 N = p ] log[α s ] logω p logω s log log.0 log log Ω 8 = s log.0 log [αs ] 9 poles of HS s i = e j π+i+n s i = e j π3+i 4 s 0 = j593. 3, s = j k = K H a s = i = 0 s i j j s s i s +865.s H z = H a s s= z +z s i593. 3s i z z z z Mathematica GUI program is written to implement the above. It can be downloaded from here Also a Matlab script nma_filter.m.txt was written no GUI to implement this design. his script written on Matlab 007a. his script does not handle the conversion from Hs to Hz well yet, need to work more on this... ofcourse, one can just use Matlab butter function for this. Example output of the above Matlab script is matlab_output.txt 4 IIR design for minmum order filter his is another small note on IIR design for minmum order filter. his document describes how to design an IIR digital filter given a specification in which the filter order is specified. Given the following digram, the specifications for the design will be. Digital filter order N. f pass, the passband corner frequency or the cutoff frequency at 3db. his means A pass = 3db 4. Impulse invariance method We first convert analog specifications to digital specifications: hen the cutoff frequency = ωp for impulse invariance. π = fpass ω p, hence ω p = π fpass 3
14 Next find the poles of H s. Since the butterworth magnitude square of the transfer function is H a s = + s j Hence H s poles are found by setting the denominator of the above to zero s + = 0 j And as I did in the earlier document, the poles of H s are found at s = e j π+k+n We only need to find the LHS poles, which are located at i = 0 N, because these are the stable poles. Hence the i th pole is s i = e j π+i+n Now we can write the analog filter generated based on the above selected poles, which is, for impulse invariance K H a s = 3 K is found by solving H a 0 = hence k = s s i s i Now we need to write poles in non-polar form and plug them into 3 s i = e j π+i+n π + i + N π + i + N = cos + j sin i = 0 N 4
15 Hence, H a s = K s cos π+i+n + j sin π+i+n Now that we have found H s we need to convert it to H z We need to make sure that we multiply poles of complex conjugates with each others to make the result simple to see. Now that we have H a s, we do the A D conversion. I.e. obtain H z from the above H s. When using impulse invariance, we need to perform partial fraction decomposition on 4 above in order to write H s in this form For example, to obtain A j, we write H s = A j = lim s sj H a s = Once we find all the A s, we now write H z as follows H z = s s i i j k s s i exp s i z his completes the design. We can try to convert the above form of H z to a rational expression as H z = Nz Dz 4. bilinear transformation method We first convert analog specifications to digital specifications: π = fpass ω p, hence ω p = π fpass Now = tan ω p Now that we have N and we find the poles of H s. Since for bilinear the magnitude square of the transfer function is H a s = + s j Hence H s poles are found by setting the denominator of the above to zero s + = 0 j Which leads to poles at s i = e j π+i+n. We only need to find the LHS poles, which are located at i = 0 N, because these are the stable poles. For bilinear, Hs is given by H a s = K is found by solving H a 0 =, hence we obtain k = 5 K s s i s i 4 3
16 We see that the same expression results for k for both cases. Now we need to write poles in non-polar form and plug them into 3 s i = e j π+i+n π + i + N π + i + N = cos + j sin i = 0 N then H a s = K s cos π+i+n + j sin π+i+n Now that we have found H s we need to convert it to H z. After finding Hs as shown above, we simply replace s by z +z. his is much simpler than the impulse invariance method. Before doing this substitution, make sure to multiply poles which are complex conjugate of each others in the denominator of Hs. After this, then do the above substitution 4 6
17 5 References. Alan Oppenheim, Ronald Schafer, Digital Signal Processing, Prentice-Hall, inc. 975, Chapter 5.. Mostafa Shiva, Electrical engineering department, California state university, Fullerton, Lecture notes, handout H. 3. John Proakis, Dimitris Manolakis, digital signal processing, 3rd edition, section
2.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 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 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 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 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 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 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 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
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 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 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 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 informationAn Iir-Filter Example: A Butterworth Filter
An Iir-Filter Example: A Butterworth Filter Josef Goette Bern University of Applied Sciences, Biel Institute of Human Centered Engineering - microlab JosefGoette@bfhch February 7, 2017 Contents 1 Introduction
More informationSEL4223 Digital Signal Processing. Inverse Z-Transform. Musa Mohd Mokji
SEL4223 Digital Signal Processing Inverse Z-Transform Musa Mohd Mokji Inverse Z-Transform Transform from z-domain to time-domain x n = 1 2πj c X z z n 1 dz Note that the mathematical operation for the
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 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 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 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 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 informationElectronic Circuits EE359A
Electronic Circuits EE359A Bruce McNair B26 bmcnair@stevens.edu 21-216-5549 Lecture 22 578 Second order LCR resonator-poles V o I 1 1 = = Y 1 1 + sc + sl R s = C 2 s 1 s + + CR LC s = C 2 sω 2 s + + ω
More 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 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 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 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 informationElectronic Circuits EE359A
Electronic Circuits EE359A Bruce McNair B26 bmcnair@stevens.edu 21-216-5549 Lecture 22 569 Second order section Ts () = s as + as+ a 2 2 1 ω + s+ ω Q 2 2 ω 1 p, p = ± 1 Q 4 Q 1 2 2 57 Second order section
More informationSIGNAL PROCESSING. B14 Option 4 lectures. Stephen Roberts
SIGNAL PROCESSING B14 Option 4 lectures Stephen Roberts Recommended texts Lynn. An introduction to the analysis and processing of signals. Macmillan. Oppenhein & Shafer. Digital signal processing. Prentice
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MI OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 008 For information about citing these materials or our erms of Use, visit: http://ocw.mit.edu/terms. Massachusetts
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 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 informationDIGITAL SIGNAL PROCESSING LECTURE 1
DIGITAL SIGNAL PROCESSING LECTURE 1 Fall 2010 2K8-5 th Semester Tahir Muhammad tmuhammad_07@yahoo.com Content and Figures are from Discrete-Time Signal Processing, 2e by Oppenheim, Shafer, and Buck, 1999-2000
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 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 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 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 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 informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More informationLecture 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 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 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 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 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 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 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 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 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 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 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 informationEE 354 Fall 2013 Lecture 10 The Sampling Process and Evaluation of Difference Equations
EE 354 Fall 203 Lecture 0 The Sampling Process and Evaluation of Difference Equations Digital Signal Processing (DSP) is centered around the idea that you can convert an analog signal to a digital signal
More informationDigital Signal Processing. Lecture Notes and Exam Questions DRAFT
Digital Signal Processing Lecture Notes and Exam Questions Convolution Sum January 31, 2006 Convolution Sum of Two Finite Sequences Consider convolution of h(n) and g(n) (M>N); y(n) = h(n), n =0... M 1
More informationEE123 Digital Signal Processing
Optimality H d (e j! ) EE13 Digital Signal Processing Lecture 3 Least Squares: Z! p! s Z Don t Care!care Variation: weighted least-squares H(e j! ) H d (e j! ) d! W (!) H(e j! ) H d (e j! ) d! Design Through
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 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 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 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 informationHilbert Transformator IP Cores
Introduction Hilbert Transformator IP Cores Martin Kumm December 27, 28 The Hilbert Transform is an important component in communication systems, e.g. for single sideband modulation/demodulation, amplitude
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 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 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 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 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 informationPart 4: IIR Filters Optimization Approach. Tutorial ISCAS 2007
Part 4: IIR Filters Optimization Approach Tutorial ISCAS 2007 Copyright 2007 Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 24, 2007 Frame # 1 Slide # 1 A. Antoniou Part4: IIR Filters
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 informationFilter Design Problem
Filter Design Problem Design of frequency-selective filters usually starts with a specification of their frequency response function. Practical filters have passband and stopband ripples, while exhibiting
More 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 informationThere are two main classes of digital lter FIR (Finite impulse response) and IIR (innite impulse reponse).
FIR Filters I There are two main classes of digital lter FIR (Finite impulse response) and IIR (innite impulse reponse). FIR Digital Filters These are described by dierence equations of the type: y[n]
More 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 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 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 informationEE 205 Dr. A. Zidouri. Electric Circuits II. Frequency Selective Circuits (Filters) Low Pass Filter. Lecture #36
EE 05 Dr. A. Zidouri Electric ircuits II Frequency Selective ircuits (Filters) ow Pass Filter ecture #36 - - EE 05 Dr. A. Zidouri The material to be covered in this lecture is as follows: o Introduction
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 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 informationExercise s = 1. cos 60 ± j sin 60 = 0.5 ± j 3/2. = s 2 + s + 1. (s + 1)(s 2 + s + 1) T(jω) = (1 + ω2 )(1 ω 2 ) 2 + ω 2 (1 + ω 2 )
Exercise 7 Ex: 7. A 0 log T [db] T 0.99 0.9 0.8 0.7 0.5 0. 0 A 0 0. 3 6 0 Ex: 7. A max 0 log.05 0 log 0.95 0.9 db [ ] A min 0 log 40 db 0.0 Ex: 7.3 s + js j Ts k s + 3 + j s + 3 j s + 4 k s + s + 4 + 3
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 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 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 informationSIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB
INTRODUCTION SIGNALS AND SYSTEMS LABORATORY 4: Polynomials, Laplace Transforms and Analog Filters in MATLAB Laplace transform pairs are very useful tools for solving ordinary differential equations. Most
More informationExamination with solution suggestions SSY130 Applied Signal Processing
Examination with solution suggestions SSY3 Applied Signal Processing Jan 8, 28 Rules Allowed aids at exam: L. Råde and B. Westergren, Mathematics Handbook (any edition, including the old editions called
More informationSteady State Frequency Response Using Bode Plots
School of Engineering Department of Electrical and Computer Engineering 332:224 Principles of Electrical Engineering II Laboratory Experiment 3 Steady State Frequency Response Using Bode Plots 1 Introduction
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 informationDiscrete Time Systems
Discrete Time Systems Valentina Hubeika, Jan Černocký DCGM FIT BUT Brno, {ihubeika,cernocky}@fit.vutbr.cz 1 LTI systems In this course, we work only with linear and time-invariant systems. We talked about
More informationTime Series Analysis: 4. Linear filters. P. F. Góra
Time Series Analysis: 4. Linear filters P. F. Góra http://th-www.if.uj.edu.pl/zfs/gora/ 2012 Linear filters in the Fourier domain Filtering: Multiplying the transform by a transfer function. g n DFT G
More 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 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 informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
More 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 informationReview of Fundamentals of Digital Signal Processing
Solution Manual for Theory and Applications of Digital Speech Processing by Lawrence Rabiner and Ronald Schafer Click here to Purchase full Solution Manual at http://solutionmanuals.info Link download
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 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 informationDiscrete-Time Signals and Systems. The z-transform and Its Application. The Direct z-transform. Region of Convergence. Reference: Sections
Discrete-Time Signals and Systems The z-transform and Its Application Dr. Deepa Kundur University of Toronto Reference: Sections 3. - 3.4 of John G. Proakis and Dimitris G. Manolakis, Digital Signal Processing:
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 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 informationUNIT - 7: FIR Filter Design
UNIT - 7: FIR Filter Design Dr. Manjunatha. P manjup.jnnce@gmail.com Professor Dept. of ECE J.N.N. College of Engineering, Shimoga October 5, 06 Unit 7 Syllabus Introduction FIR Filter Design:[,, 3, 4]
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 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 information1 Introduction & Objective. 2 Warm-up. Lab P-16: PeZ - The z, n, and O! Domains
DSP First, 2e Signal Processing First Lab P-6: PeZ - The z, n, and O! Domains The lab report/verification will be done by filling in the last page of this handout which addresses a list of observations
More informationINSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad
INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad - 500 043 Title Code Regulation ELECTRONICS AND COMMUNICATION ENGINEERING TUTORIAL QUESTION BANK DIGITAL SIGNAL PROCESSING A60421 R13 Structure
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 informationDiscrete-Time Signals and Systems. Frequency Domain Analysis of LTI Systems. The Frequency Response Function. The Frequency Response Function
Discrete-Time Signals and s Frequency Domain Analysis of LTI s Dr. Deepa Kundur University of Toronto Reference: Sections 5., 5.2-5.5 of John G. Proakis and Dimitris G. Manolakis, Digital Signal Processing:
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 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 information