Signal Processing. Lecture 10: FIR Filter Design. Ahmet Taha Koru, Ph. D. Yildiz Technical University Fall
|
|
- Sibyl Brown
- 6 years ago
- Views:
Transcription
1 Signal Processing Lecture 10: FIR Filter Design Ahmet Taha Koru, Ph. D. Yildiz Technical University Fall ATK (YTU) Signal Processing Fall 1 / 47
2 Introduction Introduction ATK (YTU) Signal Processing Fall 2 / 47
3 Introduction What is FIR filter? FIR = Finite impulse response Impulse response has finite duration Settles to zero in finite time. x[n] n Figure: A FIR filter impulse response. ATK (YTU) Signal Processing Fall 3 / 47
4 Introduction Definition For a causal discrete-time FIR filter of order N, the output is y[n] = a 0 x[n] + a 1 x[n 1] + a 2 x[n 2] + + a N x[n N] N y[n] = a k x[n k] k=0 The impulse response of the filter is N h[n] = a k δ[n k] k=0 ATK (YTU) Signal Processing Fall 4 / 47
5 Introduction Structure of FIR Filter x[n] z 1 z 1 z 1 h[0] h[1] h[2] h[3] y[n] Figure: Block diagram of a FIR filter ATK (YTU) Signal Processing Fall 5 / 47
6 Introduction FIR Properties Require no feedback of the output. Simpler implementation. Always stable. For x B x y N a k B x k=0 Can easily be designed to be linear phase by making the coefficient sequence symmetric. Disadvantage: Requires more taps and computational power is required compared to an IIR filter with similar sharpness or selectivity, especially for low frequency cut-offs. ATK (YTU) Signal Processing Fall 6 / 47
7 Introduction Design Methods Window design method Frequency sampling method Weighted least squares design Linear programming Parks-McClellan method Equiripple FIR filters design Etc... ATK (YTU) Signal Processing Fall 7 / 47
8 Introduction Filtering Specifications H(e jω ) 1 + δ 1 1 δ ω Transition band δ Ripple [ω s, π] Stop band [0, ω p] Pass band ω δ ω p ω c ω s π ω ATK (YTU) Signal Processing Fall 8 / 47
9 Window Design Method Window Design Method ATK (YTU) Signal Processing Fall 9 / 47
10 Window Design Method What is window design method? Ideal low-pass filters have infinite length of duration. Windowing is truncating the ideal low pass impulse response. Windowed impulse response is not ideal but an approximation There are different types of windows. ATK (YTU) Signal Processing Fall 10 / 47
11 Window Design Method An example: Low pass filter h i [n] n w[n] n h FIR [n] = h i [n] w[n] n ATK (YTU) Signal Processing Fall 11 / 47
12 Window Design Method Ideal Filters Window Design Method: Ideal Filters ATK (YTU) Signal Processing Fall 12 / 47
13 Window Design Method Ideal Filters Low-pass filter H LP (e jω ) H LP (e jω ) = { 1, ω ωc 0, ω c < ω π π ω c ω c π n h LP [n] h LP [n] = sin(ω cn), πn Matlab implementation: >> hlp = sinc(wc/pi*n)*wc/pi; n ATK (YTU) Signal Processing Fall 13 / 47
14 Window Design Method Ideal Filters High-pass filter H HP (e jω ) H HP (e jω ) = { 0, ω ωc 1, ω c < ω π π ω c ω c π n h HP [n] h HP [n] = δ[n] sin(ω cn) πn Matlab implementation: >> hhp = [dirac(n2)>1] - sinc(wc/pi*n2)*wc/pi; n ATK (YTU) Signal Processing Fall 14 / 47
15 Window Design Method Ideal Filters Band-pass filter 0, ω ω c1 H BP (e jω ) = 1, ω c1 < ω ω c2 0, ω c2 < ω π H BP (e jω ) h BP [n] = sin(ω c 2 n) πn sin(ω c 1 n) πn ω c2 -ω c1 ω c1 ω c2 π n Matlab implementation: >> hbp = sinc(wc2/pi*n)*wc2/pi - sinc(wc1/pi*n)*wc1/pi; ATK (YTU) Signal Processing Fall 15 / 47
16 Window Design Method Ideal Filters Band-stop filter 1, ω ω c1 H BS (e jω ) = 0, ω c1 < ω ω c2 1, ω c2 < ω π H BS (e jω ) ω c2 -ω c1 ω c1 ω c2 π n h BS [n] = sin(ω c 1 n) +δ[n] sin(ω c 2 n) πn πn Matlab implementation: >> hbs = sinc(wc1/pi*n)*wc1/pi + [dirac(n)>1]... - sinc(wc2/pi*n)*wc2/pi; ATK (YTU) Signal Processing Fall 16 / 47
17 Window Design Method Ideal Filters Multi-level Frequency Filter H ML (e jω ) = A k, for ω k 1 ω ω k H ML (e jω ) h ML [n] = L l=1 with A L+1 = 0. [ (A l A l+1 ) sin(ω ] ln) πn ω 1 ω 2 ω 3 π n Matlab implementation: >> n = -100:100; % Length of filter >> A = [3, 1, 0, 2, 0]; % List magnitudes >> w = [0.3*pi, 0.55*pi, 0.8*pi,1*pi]; % List frequencies >> h = zeros(size(n)); >> for k = 1:4 >> h = h + (A(k) - A(k+1))*sinc(w(k)/pi*n)*w(k)/pi; >> end ATK (YTU) Signal Processing Fall 17 / 47
18 Window Design Method Ideal Filters Hilbert Transformer { H HT (e jω j, π ω < 0 ) = j, 0 ω π { 0, for n even h HT [n] = for n odd 2 πn, Matlab implementation: >> h = mod(n,2)*2/pi./n; ATK (YTU) Signal Processing Fall 18 / 47
19 Window Design Method Ideal Filters Differantiator H DIFF (e jω ) = jω { 0, n = 0 h DIFF [n] = cos(πn) n, n > 0 Matlab implementation: >> h = cos(pi*n)./n; >> h(find(isinf(h))) = 0; ATK (YTU) Signal Processing Fall 19 / 47
20 Window Design Method Window Types Window Design Method: Window Types ATK (YTU) Signal Processing Fall 20 / 47
21 Window Design Method Window Types Window Types: Rectangular w[n] = u[n + M] u[n M 1] w[n] 1 M M n ATK (YTU) Signal Processing Fall 21 / 47
22 Window Design Method Window Types Window Types: Bartlett w[n] = 1 n M + 1, w[n] 1 M n M M 1 M + 1 n ATK (YTU) Signal Processing Fall 22 / 47
23 Window Design Method Window Types Window Types: Hann w[n] = 1 2 [ ( )] 2πn 1 + cos, M n M 2M + 1 w[n] 1 M M n ATK (YTU) Signal Processing Fall 23 / 47
24 Window Design Method Window Types Window Types: Hamming ( ) 2πn w[n] = cos, M n M 2M + 1 w[n] 1 M M n ATK (YTU) Signal Processing Fall 24 / 47
25 Window Design Method Window Types Window Types: Blackman ( ) ( ) 2πn 4πn w[n] = cos cos, M n M 2M + 1 2M + 1 w[n] 1 M M n ATK (YTU) Signal Processing Fall 25 / 47
26 Window Design Method Window Types Properties of some fixed window functions Type of Window Main Lobe Width Relative Sideslobe Level Minimum Stopband Attenuation Transition Bandwidth ω Rectangular 4π/(2M + 1) 13.3 db 20.9 db 0.92 π/m Bartlett 4π/(2M + 1) 26.5 db - - Hann 8π/(2M + 1) 31.5 db 43.9 db 3.11 π/m Hamming 8π/(2M + 1) 42.7 db 54.5 db 3.32 π/m Blackman 12π/(2M + 1) 58.1 db 75.3 db 5.56 π/m ATK (YTU) Signal Processing Fall 26 / 47
27 Parks-McClellan Method Parks-McClellan Method ATK (YTU) Signal Processing Fall 27 / 47
28 Parks-McClellan Method What is Parks-McClellan Algorithm Is an iterative algorithm for finding the optimal Chebyshev FIR filters. Powerful tool for design and implement efficient and optimal FIR filters. Can easily be calculated with help of a computer. ATK (YTU) Signal Processing Fall 28 / 47
29 Parks-McClellan Method What does algorithm do? A desired filter frequency response Tries to minimize error the difference with desired response Matlab Code: firpm ATK (YTU) Signal Processing Fall 29 / 47
30 Parks-McClellan Method Example: Low-pass Filter H d (e jω ) 1 Desired frequency response H(e jω ) ω p ω s π Parks-McClellan algorithm design n ω s π n ATK (YTU) Signal Processing Fall 30 / 47
31 Parks-McClellan Method Matlab Implementation N = 5; wc = pi/2; dw = 0.1*pi; % Length of filter % Cut-off Frequency % Transition Band % Parks-McClellan Algorithm h = firpm(n*2, [0 (wc-dw)/pi (wc+dw)/pi 1], [ ]); %Plot the response [H,W] = freqz(h); plot(w/pi, abs(h)) ATK (YTU) Signal Processing Fall 31 / 47
32 Linear Programming Method Linear Programming Method ATK (YTU) Signal Processing Fall 32 / 47
33 Linear Programming Method What is Linear Programming? Linear programming is optimization Linear cost function Linear constraints General form of a linear program is minimize subject to c T x Ax b where x is the vector to be designed. ATK (YTU) Signal Processing Fall 33 / 47
34 Linear Programming Method The General Structure of a FIR Filter Let us consider a symmetric FIR filter such that h[ n] = h[n], for all n Suppose the coefficients of the filters are a k such that N h[n] = a k δ[n k] k= N Then, the frequency response of the filter is N H(e jω ) = a a k cos(kω) k=1 We would like H(e jω ) to approximate an desired filter H d (e jω ). ATK (YTU) Signal Processing Fall 34 / 47
35 Linear Programming Method Low-pass Filter Linear Programming Method: Low-pass Filter ATK (YTU) Signal Processing Fall 35 / 47
36 Linear Programming Method Low-pass Filter Filtering Specifications H(e jω ) 1 + δ 1 1 δ ω Transition band δ Ripple [ω s, π] Stop band [0, ω p] Pass band ω δ ω p ω c ω s π ω ATK (YTU) Signal Processing Fall 36 / 47
37 Linear Programming Method Low-pass Filter Constraints: Pass-band (1/2) In the pass-band: N a a k cos(kω) 1 δ k=1 for all 0 ω ω c. Equivalently, ( ) N a a k cos(kω) δ 1 k=1 ( ) N a a k cos(kω) δ 1 k=1 ATK (YTU) Signal Processing Fall 37 / 47
38 Linear Programming Method Low-pass Filter Constraints: Pass-band (2/2) If we re-write the conditions in vector form [ 1 2 cos(ω) 2 cos(2ω) 2 cos(nω) cos(ω) 2 cos(2ω) 2 cos(nω) 1 ] a 0 a 1. a N δ [ 1 1 ] for any 0 ω ω p. ATK (YTU) Signal Processing Fall 38 / 47
39 Linear Programming Method Low-pass Filter Constraints: Transition-band In the transition-band: No constraints ATK (YTU) Signal Processing Fall 39 / 47
40 Linear Programming Method Low-pass Filter Constraints: Stop-band (1/2) In the stop-band: a N a k cos(kω) δ k=1 Equivalently, ( ) N a a k cos(kω) δ 0 k=1 ( ) N a a k cos(kω) δ 0 k=1 ATK (YTU) Signal Processing Fall 40 / 47
41 Linear Programming Method Low-pass Filter Constraints: Stop-band (2/2) If we re-write the conditions in vector form [ 1 2 cos(ω) 2 cos(2ω) 2 cos(nω) cos(ω) 2 cos(2ω) 2 cos(nω) 1 ] a 0 a 1. a N δ [ 0 0 ] for any ω s ω π. ATK (YTU) Signal Processing Fall 41 / 47
42 Linear Programming Method Low-pass Filter Optimization Variables and Other Parameters Optimization variables to be designed are listed in the vector a 0 a 1 x =. a N δ Given parameters are Number of taps: N Cut-off frequency: ω c Transition bandwidth: ω From ω, pass-band freq: ω p = ω c ω/2 From ω, stop-band freq: ω s = ω c + ω/2 ATK (YTU) Signal Processing Fall 42 / 47
43 Linear Programming Method Low-pass Filter Cost Function and the Optimization Problem The cost function to be minimized is the amount of ripple δ and δ = [ ] x Then, the optimization problem is minimize subject to [ ] x a 0 a [ ] cos(ω) 2 cos(nω) cos(ω) 2 cos(nω) 1. a Nδ [ a 0 a ] cos(ω) 2 cos(nω) cos(ω) 2 cos(nω) 1. a Nδ [ 1 1 [ 0 0 ], ω [ 0, ω p ] ], ω [ω s, π] There are infinitely many constraints. Problem is not solvable. ATK (YTU) Signal Processing Fall 43 / 47
44 Linear Programming Method Low-pass Filter Discretization We can not evaluate at every ω point but instead in M samples among ω [0, π]. Let s choose M = 15N The samples are chosen as ω k = π k M, Then, the constraints become cos(ω 1 ) 2 cos(nω 1 ) cos(ω 1 ) 2 cos(nω 1 ) 1 k = 0, 1,, M. 1 2 cos(ω M ) 2 cos(nω M ) cos(ω M ) 2 cos(nω M ) 1. a N δ a 0 a ATK (YTU) Signal Processing Fall 44 / 47
45 Linear Programming Method Low-pass Filter Matlab Implementation (1/3) % % Set the given parameters L = 20; % Length of the filter N = 2L + 1 wc = pi/2; % Cut-off frequency dw = 0.1*pi; % Transition bandwidth wp = wc - dw/2; ws = wc + dw/2; % Pass-band and stop-band freqs. % % Discretize M = 15*L; wl = linspace(0, pi, M + 1); R = sum(wl <= wp) + sum(wl >=ws) - 2; % Number of constraints A = zeros(r, L+2); b = zeros(r,1); cnt = 1; ATK (YTU) Signal Processing Fall 45 / 47
46 Linear Programming Method Low-pass Filter Matlab Implementation (2/3) for k = 1:length(wl) w = wl(k); if w <= wp A(cnt,:) = [1, 2*cos([1:L]*w), -1]; A(cnt+1,:) = [-1, -2*cos([1:L]*w), -1]; b(cnt) = 1; b(cnt+1) = -1; cnt = cnt+2; elseif w >= ws A(cnt,:) = [1, 2*cos([1:L]*w), -1]; A(cnt+1,:) = [-1, -2*cos([1:L]*w), -1]; b(cnt) = 0; b(cnt+1) = 0; cnt = cnt+2; end end ATK (YTU) Signal Processing Fall 46 / 47
47 Linear Programming Method Low-pass Filter Matlab Implementation (3/3) c = [zeros(1,l+1), 1]; x = linprog(c, A, b); hhalf = x(1:l+1); h = [hhalf(end:-1:2); hhalf]; % Cost function to minimize c'*x % Minimize c'*x subject to Ax <= b % Construct the filter [H, W] = freqz(h); plot(w/pi, abs(h)) ATK (YTU) Signal Processing Fall 47 / 47
-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 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 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 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 informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 3 Signals & Systems Prof. ark Fowler Note Set #28 D-T Systems: DT Filters Ideal & Practical /4 Ideal D-T Filters Just as in the CT case we can specify filters. We looked at the ideal filter for the
More 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 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 information(Refer Slide Time: 01:28 03:51 min)
Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture 40 FIR Design by Windowing This is the 40 th lecture and our topic for
More 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 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 informationINF3440/INF4440. Design of digital filters
Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 INF3440/INF4440. Design of digital filters October 2004 Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 Last lectures:
More informationDSP. Chapter-3 : Filter Design. Marc Moonen. Dept. E.E./ESAT-STADIUS, KU Leuven
DSP Chapter-3 : Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ Filter Design Process Step-1 : Define filter specs Pass-band, stop-band,
More informationLABORATORY 3 FINITE IMPULSE RESPONSE FILTERS
LABORATORY 3 FINITE IMPULSE RESPONSE FILTERS 3.. Introduction A digital filter is a discrete system, used with the purpose of changing the amplitude and/or phase spectrum of a signal. The systems (filters)
More informationLecture 14: Windowing
Lecture 14: Windowing ECE 401: Signal and Image Analysis University of Illinois 3/29/2017 1 DTFT Review 2 Windowing 3 Practical Windows Outline 1 DTFT Review 2 Windowing 3 Practical Windows DTFT Review
More 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 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 informationFilter structures ELEC-E5410
Filter structures ELEC-E5410 Contents FIR filter basics Ideal impulse responses Polyphase decomposition Fractional delay by polyphase structure Nyquist filters Half-band filters Gibbs phenomenon Discrete-time
More 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 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 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 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 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 information1. FIR Filter Design
ELEN E4810: Digital Signal Processing Topic 9: Filter Design: FIR 1. Windowed Impulse Response 2. Window Shapes 3. Design by Iterative Optimization 1 1. FIR Filter Design! FIR filters! no poles (just zeros)!
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 informationExercises in Digital Signal Processing
Exercises in Digital Signal Processing Ivan W. Selesnick September, 5 Contents The Discrete Fourier Transform The Fast Fourier Transform 8 3 Filters and Review 4 Linear-Phase FIR Digital Filters 5 5 Windows
More 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 informationDSP-CIS. Chapter-4: FIR & IIR Filter Design. Marc Moonen
DSP-CIS Chapter-4: FIR & IIR Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ PART-II : Filter Design/Realization Step-1 : Define
More 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 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 informationContents. Digital Signal Processing, Part II: Power Spectrum Estimation
Contents Digital Signal Processing, Part II: Power Spectrum Estimation 5. Application of the FFT for 7. Parametric Spectrum Est. Filtering and Spectrum Estimation 7.1 ARMA-Models 5.1 Fast Convolution 7.2
More informationOptimal Design of Real and Complex Minimum Phase Digital FIR Filters
Optimal Design of Real and Complex Minimum Phase Digital FIR Filters Niranjan Damera-Venkata and Brian L. Evans Embedded Signal Processing Laboratory Dept. of Electrical and Computer Engineering The University
More informationProblem Set 9 Solutions
Problem Set 9 Solutions EE23: Digital Signal Processing. From Figure below, we see that the DTFT of the windowed sequence approaches the actual DTFT as the window size increases. Gibb s phenomenon is absent
More informationUNIVERSITY OF OSLO. Faculty of mathematics and natural sciences. Forslag til fasit, versjon-01: Problem 1 Signals and systems.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 1th, 016 Examination hours: 14:30 18.30 This problem set
More informationEFFICIENT REMEZ ALGORITHMS FOR THE DESIGN OF NONRECURSIVE FILTERS
EFFICIENT REMEZ ALGORITHMS FOR THE DESIGN OF NONRECURSIVE FILTERS Copyright 2003- Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 24, 2007 Frame # 1 Slide # 1 A. Antoniou EFFICIENT
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 informationYour solutions for time-domain waveforms should all be expressed as real-valued functions.
ECE-486 Test 2, Feb 23, 2017 2 Hours; Closed book; Allowed calculator models: (a) Casio fx-115 models (b) HP33s and HP 35s (c) TI-30X and TI-36X models. Calculators not included in this list are not permitted.
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 informationSchool of Information Technology and Electrical Engineering EXAMINATION. ELEC3004 Signals, Systems and Control
This exam paper must not be removed from the venue Venue Seat Number Student Number Family Name First Name School of Information Technology and Electrical Engineering EXAMINATION Semester One Final Examinations,
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 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 informationECSE 512 Digital Signal Processing I Fall 2010 FINAL EXAMINATION
FINAL EXAMINATION 9:00 am 12:00 pm, December 20, 2010 Duration: 180 minutes Examiner: Prof. M. Vu Assoc. Examiner: Prof. B. Champagne There are 6 questions for a total of 120 points. This is a closed book
More 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 informationVALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur- 603 203 DEPARTMENT OF INFORMATION TECHNOLOGY Academic Year 2016-2017 QUESTION BANK-ODD SEMESTER NAME OF THE SUBJECT SUBJECT CODE SEMESTER YEAR
More 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 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:
Digital Signal Processing: Mathematical and algorithmic manipulation of discretized and quantized or naturally digital signals in order to extract the most relevant and pertinent information that is carried
More informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set
More 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 informationQuestion Bank. UNIT 1 Part-A
FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY Senkottai Village, Madurai Sivagangai Main Road, Madurai -625 020 An ISO 9001:2008 Certified Institution Question Bank DEPARTMENT OF ELECTRONICS AND COMMUNICATION
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science Discrete-Time Signal Processing Fall 2005
1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.341 Discrete-Time Signal Processing Fall 2005 FINAL EXAM Friday, December 16, 2005 Walker (50-340) 1:30pm
More informationEEL3135: Homework #4
EEL335: Homework #4 Problem : For each of the systems below, determine whether or not the system is () linear, () time-invariant, and (3) causal: (a) (b) (c) xn [ ] cos( 04πn) (d) xn [ ] xn [ ] xn [ 5]
More informationSchool of Information Technology and Electrical Engineering EXAMINATION. ELEC3004 Signals, Systems & Control
This exam paper must not be removed from the venue Venue Seat Number Student Number Family Name First Name School of Information Technology and Electrical Engineering EXAMINATION Semester One Final Examinations,
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 informationVel Tech High Tech Dr.Ranagarajan Dr.Sakunthala Engineering College Department of ECE
Subject Code: EC6502 Course Code:C302 Course Name: PRINCIPLES OF DIGITAL SIGNAL PROCESSING L-3 : T-1 : P-0 : Credits 4 COURSE OBJECTIVES: 1. To learn discrete Fourier transform and its properties 2. To
More 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( ) 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 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 informationSchool of Information Technology and Electrical Engineering EXAMINATION. ELEC3004 Signals, Systems & Control
This exam paper must not be removed from the venue Venue Seat Number Student Number Family Name First Name School of Information Technology and Electrical Engineering EXAMINATION Semester One Final Examinations,
More informationEE 311 February 22, 2019 Lecture 13
EE 311 February, 019 Lecture 13 Minimum Phase FIR filters Geometric Interpretation of P/Z locations Minimum Phase System - A system is said to be a minimum phase system if the deviation from zero in the
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 18: Stability
Lecture 18: Stability ECE 401: Signal and Image Analysis University of Illinois 4/18/2017 1 Stability 2 Impulse Response 3 Z Transform Outline 1 Stability 2 Impulse Response 3 Z Transform BIBO Stability
More informationLecture 19 IIR Filters
Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class
More informationSignals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 2 Solutions
8-90 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 08 Midterm Solutions Name: Andrew ID: Problem Score Max 8 5 3 6 4 7 5 8 6 7 6 8 6 9 0 0 Total 00 Midterm Solutions. (8 points) Indicate whether
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 informationChap 2. Discrete-Time Signals and Systems
Digital Signal Processing Chap 2. Discrete-Time Signals and Systems Chang-Su Kim Discrete-Time Signals CT Signal DT Signal Representation 0 4 1 1 1 2 3 Functional representation 1, n 1,3 x[ n] 4, n 2 0,
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 informationREAL TIME DIGITAL SIGNAL PROCESSING
REAL TIME DIGITAL SIGNAL PROCESSING www.electron.frba.utn.edu.ar/dplab Digital Filters FIR and IIR. Design parameters. Implementation types. Constraints. Filters: General classification Filters: General
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 informationConvolution. Define a mathematical operation on discrete-time signals called convolution, represented by *. Given two discrete-time signals x 1, x 2,
Filters Filters So far: Sound signals, connection to Fourier Series, Introduction to Fourier Series and Transforms, Introduction to the FFT Today Filters Filters: Keep part of the signal we are interested
More informationIntroduction to DSP Time Domain Representation of Signals and Systems
Introduction to DSP Time Domain Representation of Signals and Systems Dr. Waleed Al-Hanafy waleed alhanafy@yahoo.com Faculty of Electronic Engineering, Menoufia Univ., Egypt Digital Signal Processing (ECE407)
More informationDiscrete Fourier Transform
Discrete Fourier Transform Virtually all practical signals have finite length (e.g., sensor data, audio records, digital images, stock values, etc). Rather than considering such signals to be zero-padded
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 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 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 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 informationECGR4124 Digital Signal Processing Final Spring 2009
ECGR4124 Digital Signal Processing Final Spring 2009 Name: LAST 4 NUMBERS of Student Number: Do NOT begin until told to do so Make sure that you have all pages before starting Open book, 2 sheet front/back
More informationTransform Analysis of Linear Time-Invariant Systems
Transform Analysis of Linear Time-Invariant Systems Discrete-Time Signal Processing Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-Sen University Kaohsiung, Taiwan ROC Transform
More informationEECS 123 Digital Signal Processing University of California, Berkeley: Fall 2007 Gastpar November 7, Exam 2
EECS 3 Digital Signal Processing University of California, Berkeley: Fall 7 Gastpar November 7, 7 Exam Last name First name SID You have hour and 45 minutes to complete this exam. he exam is closed-book
More informationDigital Signal Processing Lecture 4
Remote Sensing Laboratory Dept. of Information Engineering and Computer Science University of Trento Via Sommarive, 14, I-38123 Povo, Trento, Italy Digital Signal Processing Lecture 4 Begüm Demir E-mail:
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 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 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 informationDiscrete-Time Fourier Transform
C H A P T E R 7 Discrete-Time Fourier Transform In Chapter 3 and Appendix C, we showed that interesting continuous-time waveforms x(t) can be synthesized by summing sinusoids, or complex exponential signals,
More informationDigital Signal Processing Lecture 3 - Discrete-Time Systems
Digital Signal Processing - Discrete-Time Systems Electrical Engineering and Computer Science University of Tennessee, Knoxville August 25, 2015 Overview 1 2 3 4 5 6 7 8 Introduction Three components of
More informationProblem Value
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 30-Apr-04 COURSE: ECE-2025 NAME: GT #: LAST, FIRST Recitation Section: Circle the date & time when your Recitation
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 informationOptimum Design of Frequency-Response-Masking Filters Using Convex-Concave Procedure. Haiying Chen
Optimum Design of Frequency-Response-Masking Filters Using Convex-Concave Procedure by Haiying Chen B.Eng., University of Electrical Science and Technology of China, 2013 A Dissertation Submitted in partial
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,5 8,.7 Open access books available International authors and editors Downloads Our authors are
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 informationFINITE PRECISION EFFECTS 1. FLOATING POINT VERSUS FIXED POINT 3. TYPES OF FINITE PRECISION EFFECTS 4. FACTORS INFLUENCING FINITE PRECISION EFFECTS
FINITE PRECISION EFFECTS 1. FLOATING POINT VERSUS FIXED POINT 2. WHEN IS FIXED POINT NEEDED? 3. TYPES OF FINITE PRECISION EFFECTS 4. FACTORS INFLUENCING FINITE PRECISION EFFECTS 5. FINITE PRECISION EFFECTS:
More informationAnalog vs. discrete signals
Analog vs. discrete signals Continuous-time signals are also known as analog signals because their amplitude is analogous (i.e., proportional) to the physical quantity they represent. Discrete-time signals
More informationEE 313 Linear Signals & Systems (Fall 2018) Solution Set for Homework #7 on Infinite Impulse Response (IIR) Filters CORRECTED
EE 33 Linear Signals & Systems (Fall 208) Solution Set for Homework #7 on Infinite Impulse Response (IIR) Filters CORRECTED By: Mr. Houshang Salimian and Prof. Brian L. Evans Prolog for the Solution Set.
More informationMaximally Flat Lowpass Digital Differentiators
Maximally Flat Lowpass Digital Differentiators Ivan W. Selesnick August 3, 00 Electrical Engineering, Polytechnic University 6 Metrotech Center, Brooklyn, NY 0 selesi@taco.poly.edu tel: 78 60-36 fax: 78
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science. Fall Solutions for Problem Set 2
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Issued: Tuesday, September 5. 6.: Discrete-Time Signal Processing Fall 5 Solutions for Problem Set Problem.
More 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 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 informationProblem Value
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 30-Apr-04 COURSE: ECE-2025 NAME: GT #: LAST, FIRST Recitation Section: Circle the date & time when your Recitation
More informationChirp Transform for FFT
Chirp Transform for FFT Since the FFT is an implementation of the DFT, it provides a frequency resolution of 2π/N, where N is the length of the input sequence. If this resolution is not sufficient in a
More informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 2011) Final Examination December 19, 2011 Name: Kerberos Username: Please circle your section number: Section Time 2 11 am 1 pm 4 2 pm Grades will be determined by the correctness of your answers
More 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 information