|
|
- Julian Leonard
- 5 years ago
- Views:
Transcription
1 Low-Pass Filter _012 Triple Butterworth Filter to avoid numerical instability for high order filter. In [5]: import numpy as np import matplotlib.pyplot as plt %matplotlib notebook fsz = (7,5) # figure size import scipy.signal as ss Suppose you need to design a Butterworth low-pass filter with -3 db cutoff at 50 Hz and stopband attenuation of -40 db at 75 Hz on a signal with sampling rate F s = 8000 Hz. This requires a fairly high filter order because the transition bandwidth of 25 Hz is only a small fraction (25/4000 = ) of the maximum frequency that can be represented at the given F s. In [2]: Fs = 8000 tt = np.arange(fs)/float(fs)-1/2.0 # time axis sec uimp = np.zeros(tt.size) ix0 = np.argmin(np.abs(tt)) uimp[ix0] = Fs # unit impulse In [3]: # Trying a filter order of 15 fl = 50 f_ord = 15 b,a = ss.butter(f_ord,2*fl/float(fs)) ht = ss.lfilter(b,a,uimp) # DT Butterworth filter In [6]: fsz2 = (fsz[0],fsz[1]/2.0) plt.figure(1, figsize=fsz2) plt.plot(tt,ht,'-b') plt.xlabel('t [sec]') plt.ylabel('h(t)') strt1 = 'Butterworth h(t), $f_l=${} Hz, ord={}'.format(fl,f_ord) strt1 = strt1 + ', $F_s=${} Hz'.format(Fs) plt.title(strt1) 1 of 11 10/2/18, 6:11 PM
2 This is clearly not the right impulse response of a low-pass filter. It is in fact the impulse response of an unstable filter. Because of the scaling of the y-axis (1e220), we can assume that the filtering was done using double precision floating point numbers (float64). To find out what's wrong, let's look at the poles and zeros of the filter. In [7]: f_poles = np.roots(a) f_zeros = np.roots(b) u_circ = np.exp(1j*2*np.pi*np.arange(361)/360.0) # unit circle plt.figure(3,fsz) plt.plot(np.real(u_circ),np.imag(u_circ),':k') plt.plot(np.real(f_poles),np.imag(f_poles),'xr',label='poles') plt.plot(np.real(f_zeros),np.imag(f_zeros),'ob',label='zeros') plt.gca().set_aspect('equal', adjustable='box') plt.xlabel('real Part') plt.ylabel('imaginary Part') strt3 = 'Butterworth Poles/Zeros, $f_l=${} Hz, ord={}'.format(fl,f_ord) strt3 = strt3 + ', $F_s=${} Hz'.format(Fs) plt.title(strt3) plt.legend(loc=1) It is easy to see that there are poles outside the unit circle which cause the instability of the filter. Also since f L = 50 = 1/ a b, the poles seem to be too far spread out. The most likely reason for this is that the computation of the and coefficients of the filter would have required a higher precision that float64. Through trial and error we found that for the given parameters ( f L and F s ) the highest filter order that results in a stable filter is 9. So let's try that. 2 of 11 10/2/18, 6:11 PM
3 In [8]: # Trying a filter order of 9 fl = 50 f_ord = 9 b,a = ss.butter(f_ord,2*fl/float(fs)) ht = ss.lfilter(b,a,uimp) # DT Butterworth filter In [9]: plt.figure(5, figsize=fsz2) plt.plot(tt,ht,'-b') plt.xlabel('t [sec]') plt.ylabel('h(t)') strt5 = 'Butterworth h(t), $f_l=${} Hz, ord={}'.format(fl,f_ord) strt5 = strt5 + ', $F_s=${} Hz'.format(Fs) plt.title(strt5) Now we can take a look at the frequency response of the filter. 3 of 11 10/2/18, 6:11 PM
4 In [10]: def FTapprox(tt,xt,dlim): """ DFT/FFT approximation to Fourier transform of x(t) with time axis tt. >>>>> ff, Xf, absxf, argxf, Df = FTapprox(tt,xt,dlim) <<<<< where ff: Frequency axis absxf: X(f) argxf: angle(x(f)) in degrees Df: Frequency resolution tt: Time axis xt: Sampled CT waveform dlim = [f1,f2,llim]: display limits f1,f2: Lower/upper frequency limits llim = 0: display Xfabs linear and absolute llim > 0: same as llim=0 but phase is masked for Xfabs<llim llim < 0: display 20*log_{10}( X(f)/max( X(f))) in db with lower display limit llim db. phase is masked for display values < llim db """ # ***** prepare x(t), swap pos/neg parts along time axis ***** N = tt.size Fs = int(round((n-1)/(tt[-1]-tt[0]))) ixtp = np.where(tt>=0)[0] ixtn = np.where(tt<0)[0] xtp = np.hstack((xt[ixtp],xt[ixtn])) # ***** compute X(f), make frequency axis ***** Xf = 1/float(Fs)*np.fft.fft(xtp) ff = Fs/float(N)*np.arange(N) if dlim[0] < 0: ixfp = np.where(ff<fs/2.0)[0] ixfn = np.where(ff>=fs/2.0)[0] ff = np.hstack((ff[ixfn]-fs, ff[ixfp])) Xf = np.hstack((xf[ixfn],xf[ixfp])) # ***** compute X(f), arg[x(f)], trim to [f1,f2) ***** absxf = np.abs(xf) maxxf = np.amax(absxf) # max of X(f) ixf = np.where(np.logical_and(ff>=dlim[0],ff<dlim[1]))[0] ff = ff[ixf] Xf = Xf[ixf] absxf = absxf[ixf] argxf = np.angle(xf) # ***** Mask phase ***** if dlim[2] > 0: ix = np.where(absxf<dlim[2])[0] argxf[ix] = 0 elif dlim[2] < 0: absxf = absxf/float(maxxf) # normalized magnitude ix = np.where(absxf==0)[0] absxf[ix] = 1e-300 # avoid log of zero absxf = 20*np.log10(absXf) ix = np.where(absxf<dlim[2])[0] absxf[ix] = dlim[2] argxf[ix] = 0 # mask phase return ff, Xf, absxf, argxf, Fs/float(N) In [11]: f_disp, llim = 200, 1e-6 dlim = [-f_disp, f_disp, llim] ff,hf,abshf,arghf,df = FTapprox(tt,ht,dlim) N = tt.size 4 of 11 10/2/18, 6:11 PM
5 In [12]: plt.figure(7,figsize=fsz) plt.subplot(211) plt.plot(ff,abshf,'-b') plt.ylabel('$ H(f) $') strt7 = 'Butterworth $H(f)$, $f_l=${} Hz, ord={}'.format(fl, f_ord) strt7 = strt7 + ', $F_s=${} Hz'.format(Fs) plt.title(strt7) plt.subplot(212) plt.plot(ff,180/np.pi*arghf,'-r') plt.ylabel('arg[$h(f)$] [deg]') plt.xlabel('f [Hz]') One of the defining characteristics of Butterworth filters is that their magnitude is maximally flat in the passband. This is clearly not the case in the above magnitude frequency response. Let's try what happens if we reduce the filter order to 8. In [13]: # Trying a filter order of 8 fl = 50 f_ord = 8 b,a = ss.butter(f_ord,2*fl/float(fs)) ht = ss.lfilter(b,a,uimp) # DT Butterworth filter In [14]: ff,hf,abshf,arghf,df = FTapprox(tt,ht,dlim) 5 of 11 10/2/18, 6:11 PM
6 In [15]: plt.figure(9,figsize=fsz) plt.subplot(211) plt.plot(ff,abshf,'-b') plt.ylabel('$ H(f) $') strt9 = 'Butterworth $H(f)$, $f_l=${} Hz, ord={}'.format(fl, f_ord) strt9 = strt9 + ', $F_s=${} Hz'.format(Fs) plt.title(strt9) plt.subplot(212) plt.plot(ff,180/np.pi*arghf,'-r') plt.ylabel('arg[$h(f)$] [deg]') plt.xlabel('f [Hz]') By design, the -3 db frequency is 50 Hz. But what about the -40 db attenuation? At which frequency does it occur? We need to take a look at the magnitude of the frequency response in db. In [16]: f_dispdb, llimdb = 200, -80 dlimdb = [-f_dispdb, f_dispdb, llimdb] ffdb,hf,abshfdb,arghfdb,df = FTapprox(tt,ht,dlimdB) 6 of 11 10/2/18, 6:11 PM
7 In [17]: plt.figure(11,figsize=fsz) plt.subplot(211) plt.plot(ffdb,abshfdb,'-b') plt.ylabel('$ H(f) $ [db]') strt11 = 'Butterworth $H(f)$, $f_l=${} Hz, ord={}'.format(fl, f_ord) strt11 = strt11 + ', $F_s=${} Hz'.format(Fs) plt.title(strt11) plt.subplot(212) plt.plot(ffdb,180/np.pi*arghfdb,'-r') plt.ylabel('arg[$h(f)$] [deg]') plt.xlabel('f [Hz]') Zooming in, we find that the -40 db frequency is approximately 89 Hz. To make the transition from passband to stopband more steep we need to add more filtering. To avoid numerical instability, we need to use more than one filter to increase the filter order. Ideally, we could compute the poles and zeros from the and polynomials and then distribute them over the filters. But since the problem here seems to be the and polynomials themselves, we just have to specify -3 db frequencies and orders 8 for two or more filters in cascade. a b a b In [18]: # parameters for cascade of two Butterworth filters fl1, f_ord1 = 57, 8 fl2, f_ord2 = 51, 7 In [19]: b1,a1 = ss.butter(f_ord1,2*fl1/float(fs)) # DT Butterworth filter 1 h1t = ss.lfilter(b1,a1,uimp) b2,a2 = ss.butter(f_ord2,2*fl2/float(fs)) # DT Butterworth filter 2 h2t = ss.lfilter(b2,a2,h1t) 7 of 11 10/2/18, 6:11 PM
8 In [20]: plt.figure(13, figsize=fsz2) plt.plot(tt,h2t,'-b') plt.xlabel('t [sec]') plt.ylabel('h(t)') strt13 = 'Two Butterworth h(t), $f_{{l1}}=${} Hz, ord$_1$={}'.format(fl1,f_ ord1) strt13 = strt13 + ', $f_{{l2}}=${} Hz, ord$_2$={}'.format(fl2,f_ord2) strt13 = strt13 + ', $F_s=${} Hz'.format(Fs) plt.title(strt13) In [21]: ffdb,h2f,absh2fdb,argh2fdb,df = FTapprox(tt,h2t,dlimdB) 8 of 11 10/2/18, 6:11 PM
9 In [22]: plt.figure(15,figsize=fsz) plt.subplot(211) plt.plot(ffdb,absh2fdb,'-b') plt.ylabel('$ H(f) $ [db]') strt15 = 'Two Butterworth $H(f)$, $f_{{l1}}=${} Hz, ord$_1$={}'.format(fl1, f_ord1) strt15 = strt15 + ', $f_{{l2}}=${} Hz, ord$_2$={}'.format(fl2,f_ord2) strt15 = strt15 + ', $F_s=${} Hz'.format(Fs) plt.title(strt15) plt.subplot(212) plt.plot(ffdb,180/np.pi*argh2fdb,'-r') plt.ylabel('arg[$h(f)$] [deg]') plt.xlabel('f [Hz]') f L1 = 57 1 = 8 f L2 = 51 = 7 By trial and error, we found a cascade of two Butterworth filters, with Hz, ord, and Hz, ord2 with -3 db frequency of Hz and -40 db frequency Hz. The total filter order is thus To compensate for the filter delay, we measure the delay of the peak of the unit impulse response and then pad the filter input with zeros corresponding to that delay. After filtering the same amount of samples is removed from the beginning of the filter response. Small adjustments to the delay can then be made based on the phase of the frequency response. In [23]: t_dly = 28.5e-3 # Filter delay in seconds n_dly = int(np.round(t_dly*fs)) h1t_dly = ss.lfilter(b1,a1,np.hstack((uimp,np.zeros(n_dly)))) h2t_dly = ss.lfilter(b2,a2,h1t_dly) h2t_dly = h2t_dly[n_dly:] # delay compensated output 9 of 11 10/2/18, 6:11 PM
10 In [24]: plt.figure(17, figsize=fsz2) plt.plot(tt,h2t_dly,'-b') plt.xlabel('t [sec]') plt.ylabel('h(t)') strt17 = 'Two Butterworth h(t), $f_{{l1}}=${} Hz, ord$_1$={}'.format(fl1,f_ ord1) strt17 = strt17 + ', $f_{{l2}}=${} Hz, ord$_2$={}'.format(fl2,f_ord2) strt17 = strt17 + ', $F_s=${} Hz'.format(Fs) plt.title(strt17) In [25]: ffdb,h2f_dly,absh2fdb_dly,argh2fdb_dly,df = FTapprox(tt,h2t_dly,dlimdB) 10 of 11 10/2/18, 6:11 PM
11 In [26]: plt.figure(19,figsize=fsz) plt.subplot(211) plt.plot(ffdb,absh2fdb_dly,'-b') plt.ylabel('$ H(f) $ [db]') strt19 = 'Two Butterworth $H(f)$, $f_{{l1}}=${} Hz, ord$_1$={}'.format(fl1, f_ord1) strt19 = strt19 + ', $f_{{l2}}=${} Hz, ord$_2$={}'.format(fl2,f_ord2) strt19 = strt19 + ', $F_s=${} Hz'.format(Fs) plt.title(strt19) plt.subplot(212) plt.plot(ffdb,180/np.pi*argh2fdb_dly,'-r') plt.ylabel('arg[$h(f)$] [deg]') plt.xlabel('f [Hz]') After some fine tuning, the phase reponse near dc is essentially flat, indicating a proper delay compensation of the filter. In [ ]: 11 of 11 10/2/18, 6:11 PM
EE482: 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 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 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 informationw n = c k v n k (1.226) w n = c k v n k + d k w n k (1.227) Clearly non-recursive filters are a special case of recursive filters where M=0.
Random Data 79 1.13 Digital Filters There are two fundamental types of digital filters Non-recursive N w n = c k v n k (1.226) k= N and recursive N M w n = c k v n k + d k w n k (1.227) k= N k=1 Clearly
More 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 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 informationTime Series Analysis: 4. Digital Linear Filters. P. F. Góra
Time Series Analysis: 4. Digital Linear Filters P. F. Góra http://th-www.if.uj.edu.pl/zfs/gora/ 2018 Linear filters Filtering in Fourier domain is very easy: multiply the DFT of the input by a transfer
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 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 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 informationUNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter?
UNIT - III PART A. Mention the important features of the IIR filters? i) The physically realizable IIR filters does not have linear phase. ii) The IIR filter specification includes the desired characteristics
More informationDIGITAL 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 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 informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
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 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 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 informationUniversity of Illinois at Chicago Spring ECE 412 Introduction to Filter Synthesis Homework #4 Solutions
Problem 1 A Butterworth lowpass filter is to be designed having the loss specifications given below. The limits of the the design specifications are shown in the brick-wall characteristic shown in Figure
More informationRepresentation of Signals & Systems
Representation of Signals & Systems Reference: Chapter 2,Communication Systems, Simon Haykin. Hilbert Transform Fourier transform frequency content of a signal (frequency selectivity designing frequency-selective
More informationMITOCW watch?v=jtj3v Rx7E
MITOCW watch?v=jtj3v Rx7E The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free. To
More informationDigital Signal Processing
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 informationQuestion Paper Code : AEC11T02
Hall Ticket No Question Paper Code : AEC11T02 VARDHAMAN COLLEGE OF ENGINEERING (AUTONOMOUS) Affiliated to JNTUH, Hyderabad Four Year B. Tech III Semester Tutorial Question Bank 2013-14 (Regulations: VCE-R11)
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 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 informationOversampling Converters
Oversampling Converters David Johns and Ken Martin (johns@eecg.toronto.edu) (martin@eecg.toronto.edu) slide 1 of 56 Motivation Popular approach for medium-to-low speed A/D and D/A applications requiring
More 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 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 informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
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 informationHigher-Order Σ Modulators and the Σ Toolbox
ECE37 Advanced Analog Circuits Higher-Order Σ Modulators and the Σ Toolbox Richard Schreier richard.schreier@analog.com NLCOTD: Dynamic Flip-Flop Standard CMOS version D CK Q Q Can the circuit be simplified?
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 informationLectures: Lumped element filters (also applies to low frequency filters) Stub Filters Stepped Impedance Filters Coupled Line Filters
ECE 580/680 Microwave Filter Design Lectures: Lumped element filters (also applies to low frequency filters) Stub Filters Stepped Impedance Filters Coupled Line Filters Lumped Element Filters Text Section
More 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 informationOrdinary Differential Equations II: Runge-Kutta and Advanced Methods
Ordinary Differential Equations II: Runge-Kutta and Advanced Methods Sam Sinayoko Numerical Methods 3 Contents 1 Learning Outcomes 2 2 Introduction 2 2.1 Note................................ 4 2.2 Limitations
More informationImproved Interpolation
Improved Interpolation Interpolation plays an important role for motion compensation with improved fractional pixel accuracy. The more precise interpolated we get, the smaller our prediction residual becomes,
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 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 informationLab 5: Power Spectral Density, Noise, and Symbol Timing Information
ECEN 4652/5002 Communications Lab Spring 2017 2-20-17 P. Mathys Lab 5: Power Spectral Density, Noise, and Symbol Timing Information 1 Introduction The two concepts that are most fundamental to the realistic
More informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More 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 informationFourier Methods in Digital Signal Processing Final Exam ME 579, Spring 2015 NAME
Fourier Methods in Digital Signal Processing Final Exam ME 579, Instructions for this CLOSED BOOK EXAM 2 hours long. Monday, May 8th, 8-10am in ME1051 Answer FIVE Questions, at LEAST ONE from each section.
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 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 informationDriven, damped, pendulum
Driven, damped, pendulum Note: The notation and graphs in this notebook parallel those in Chaotic Dynamics by Baker and Gollub. (There's a copy in the department office.) For the driven, damped, pendulum,
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 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 informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 301 Signals & Systems Prof. Mark Fowler Note Set #15 C-T Systems: CT Filters & Frequency Response 1/14 Ideal Filters Often we have a scenario where part of the input signal s spectrum comprises what
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 informationGrades will be determined by the correctness of your answers (explanations are not required).
6.00 (Fall 20) Final Examination December 9, 20 Name: Kerberos Username: Please circle your section number: Section Time 2 am pm 4 2 pm Grades will be determined by the correctness of your answers (explanations
More informationExperiment 13 Poles and zeros in the z plane: IIR systems
Experiment 13 Poles and zeros in the z plane: IIR systems Achievements in this experiment You will be able to interpret the poles and zeros of the transfer function of discrete-time filters to visualize
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 informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationSystems & Signals 315
1 / 15 Systems & Signals 315 Lecture 13: Signals and Linear Systems Dr. Herman A. Engelbrecht Stellenbosch University Dept. E & E Engineering 2 Maart 2009 Outline 2 / 15 1 Signal Transmission through a
More informationEnhanced Steiglitz-McBride Procedure for. Minimax IIR Digital Filters
Enhanced Steiglitz-McBride Procedure for Minimax IIR Digital Filters Wu-Sheng Lu Takao Hinamoto University of Victoria Hiroshima University Victoria, Canada Higashi-Hiroshima, Japan May 30, 2018 1 Outline
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.011: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 2010 ANSWER BOOKLET
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 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 informationStability. X(s) Y(s) = (s + 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = +2 and s = -2. Y(s) 8X(s) G 1 G 2
Stability 8X(s) X(s) Y(s) = (s 2) 2 (s 2) System has 2 poles: points where Y(s) -> at s = 2 and s = -2 If all poles are in region where s < 0, system is stable in Fourier language s = jω G 0 - x3 x7 Y(s)
More informationChapter 2. Classical Control System Design. Dutch Institute of Systems and Control
Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral
More informationExercise 1 (A Non-minimum Phase System)
Prof. Dr. E. Frazzoli 5-59- Control Systems I (Autumn 27) Solution Exercise Set 2 Loop Shaping clruch@ethz.ch, 8th December 27 Exercise (A Non-minimum Phase System) To decrease the rise time of the system,
More informationFROM ANALOGUE TO DIGITAL
SIGNALS AND SYSTEMS: PAPER 3C1 HANDOUT 7. Dr David Corrigan 1. Electronic and Electrical Engineering Dept. corrigad@tcd.ie www.mee.tcd.ie/ corrigad FROM ANALOGUE TO DIGITAL To digitize signals it is necessary
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 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 informationSignal Conditioning Issues:
Michael David Bryant 1 11/1/07 Signal Conditioning Issues: Filtering Ideal frequency domain filter Ideal time domain filter Practical Filter Types Butterworth Bessel Chebyschev Cauer Elliptic Noise & Avoidance
More informationSpeaker: Arthur Williams Chief Scientist Telebyte Inc. Thursday November 20 th 2008 INTRODUCTION TO ACTIVE AND PASSIVE ANALOG
INTRODUCTION TO ACTIVE AND PASSIVE ANALOG FILTER DESIGN INCLUDING SOME INTERESTING AND UNIQUE CONFIGURATIONS Speaker: Arthur Williams Chief Scientist Telebyte Inc. Thursday November 20 th 2008 TOPICS Introduction
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 informationDistributed Real-Time Control Systems
Distributed Real-Time Control Systems Chapter 9 Discrete PID Control 1 Computer Control 2 Approximation of Continuous Time Controllers Design Strategy: Design a continuous time controller C c (s) and then
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 informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 21: Stability Margins and Closing the Loop Overview In this Lecture, you will learn: Closing the Loop Effect on Bode Plot Effect
More informationProblem Value
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 2-May-05 COURSE: ECE-2025 NAME: GT #: LAST, FIRST (ex: gtz123a) Recitation Section: Circle the date & time when
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 informationDHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A
DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A Classification of systems : Continuous and Discrete
More 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 informationActive Control? Contact : Website : Teaching
Active Control? Contact : bmokrani@ulb.ac.be Website : http://scmero.ulb.ac.be Teaching Active Control? Disturbances System Measurement Control Controler. Regulator.,,, Aims of an Active Control Disturbances
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 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 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 informationExercises for lectures 13 Design using frequency methods
Exercises for lectures 13 Design using frequency methods Michael Šebek Automatic control 2016 31-3-17 Setting of the closed loop bandwidth At the transition frequency in the open loop is (from definition)
More informationChapter 4: Filtering in the Frequency Domain. Fourier Analysis R. C. Gonzalez & R. E. Woods
Fourier Analysis 1992 2008 R. C. Gonzalez & R. E. Woods Properties of δ (t) and (x) δ : f t) δ ( t t ) dt = f ( ) f x) δ ( x x ) = f ( ) ( 0 t0 x= ( 0 x0 1992 2008 R. C. Gonzalez & R. E. Woods Sampling
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 informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationLECTURE 3 CMOS PHASE LOCKED LOOPS
Lecture 03 (8/9/18) Page 3-1 LECTURE 3 CMOS PHASE LOCKED LOOPS Topics The acquisition process unlocked state Noise in linear PLLs Organization: Systems Perspective Types of PLLs and PLL Measurements PLL
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 informationModule 3 : Sampling and Reconstruction Lecture 22 : Sampling and Reconstruction of Band-Limited Signals
Module 3 : Sampling and Reconstruction Lecture 22 : Sampling and Reconstruction of Band-Limited Signals Objectives Scope of this lecture: If a Continuous Time (C.T.) signal is to be uniquely represented
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 informationTopic 3: Fourier Series (FS)
ELEC264: Signals And Systems Topic 3: Fourier Series (FS) o o o o Introduction to frequency analysis of signals CT FS Fourier series of CT periodic signals Signal Symmetry and CT Fourier Series Properties
More informationThe Approximation Problem
EE 508 Lecture 3 The Approximation Problem Classical Approximating Functions - Thompson and Bessel Approximations Review from Last Time Elliptic Filters Can be thought of as an extension of the CC approach
More informationThe Approximation Problem
EE 508 Lecture The Approximation Problem Classical Approximating Functions - Elliptic Approximations - Thompson and Bessel Approximations Review from Last Time Chebyshev Approximations T Type II Chebyshev
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 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 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 informationExercise 1 (A Non-minimum Phase System)
Prof. Dr. E. Frazzoli 5-59- Control Systems I (HS 25) Solution Exercise Set Loop Shaping Noele Norris, 9th December 26 Exercise (A Non-minimum Phase System) To increase the rise time of the system, we
More informationFINAL EXAM. Problem Value Score
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 27-Apr-09 COURSE: ECE-2025 NAME: GT username: LAST, FIRST (ex: gpburdell3) 3 points 3 points 3 points Recitation
More informationS4 hardware injections
S4 hardware injections Gabriela González, Michael Landry, Vuk Mandic, Brian O Reilly September 16, 2005 Abstract We show that hardware injections in ETMX EXC produce a desired strain h 0 if EXC = h 0 A
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 informationEE422G Homework #9 Solution
EE4G Homework #9 Solution. (3 points) Sampling and Reconstruction a) Given the following discrete-time input x(n), sketch the outputs reconstructed continuous signal y(t) using Sample-And-Hold and Linear
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 information!Sketch f(t) over one period. Show that the Fourier Series for f(t) is as given below. What is θ 1?
Second Year Engineering Mathematics Laboratory Michaelmas Term 998 -M L G Oldfield 30 September, 999 Exercise : Fourier Series & Transforms Revision 4 Answer all parts of Section A and B which are marked
More informationLab 3: Poles, Zeros, and Time/Frequency Domain Response
ECEN 33 Linear Systems Spring 2 2-- P. Mathys Lab 3: Poles, Zeros, and Time/Frequency Domain Response of CT Systems Introduction Systems that are used for signal processing are very often characterized
More informationFINAL EXAM. Problem Value Score No/Wrong Rec 3
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL & COMPUTER ENGINEERING FINAL EXAM DATE: 1-May-08 COURSE: ECE-2025 NAME: GT username: LAST, FIRST (ex: gpburdell3) 3 points 3 points 3 points Recitation
More information