2.161 Signal Processing: Continuous and Discrete
|
|
- Eustacia Green
- 6 years ago
- Views:
Transcription
1 MIT OpenCourseWare Signal Processing: Continuous and Discrete Fall 8 For information about citing these materials or our Terms of Use, visit:
2 MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING.6 Signal Processing - Continuous and Discrete Fall Term 8 Solution of Problem Set 7 Assigned: October 3, 8 Due: November 6, 8 Problem : There are a number of ways to prove this here are two: (a) Start with an acausal filter of the same length H a (z) = (N)/ n= (N)/ n h n z with N odd, and with real, odd symmetric coefficients h n = h n so that the impulse response is a real, odd function. Since DTFT {h n } H(e jω ), from the properties of the DTFT, namely that a real, odd function has a imaginary, odd DTFT, then H a (e jω ) is imaginary and odd. H a (z) may be made causal by adding a delay of (N )/, that is H(z) = z (N)/ H a (z) so that Since H a (e jω ) is imaginary H(e jω ) = e jω(n)/ H a (e jω ) H(e jω ) = je jω(n)/ = (N ) Ω ± π/ which is linear with frequency Ω = ωδt (apart from possible jump discontinuities where H a (e jω ) changes sign. Note that in practical filters, any jump in the phase response will occur outside the pass-band, and is generally ignored.) The equivalent delay is found from the slope, and represents a delay of (N )/ samples. (b) Assume the causal filter and let H(z) = N n h n z n= = h z (N) + h z h N z with N odd, and with odd symmetric coefficients about the mid-point, h n = h N n. Group the symmetric components together ( ) ( ) ( ) (N) (N ) (N ) (N)/ H(z) = h z z + h z z + h z z h (N)/ z
3 (N)/ ( ) n (N n) (N)/ = h n z z + h (N)/ z n= (N)/ ( ) (N)/ = z h n z (N)/ n ((N/) n) z n= since, as noted in the problem statement, for an odd-symmetric function about the mid-point h (N)/ =. The frequency response is H(e jω ) = H(z) z=e jω (N)/ ( ) H(e jω jω(n)/ ) = e h n e jω((n)/ n) e jω((n/) n) n= (N)/ = je jω(n)/ h n sin (Ω((N )/ + n)) n= = je jω(n)/ H a (e jω ), The function H a (e jω ) is purely real, therefore H(e jω ) = je jω(n)/ = (N ) Ω ± π/ which is linear with frequency Ω = ωδt (apart from possible jump discontinuities of π at sign changes in H a (e jω ), which as noted above occur outside the pass-band) and represents a delay of (N )/ samples. Problem : (a) The system has 7 zeros and 7 poles at the origin. It is therefore a non-recursive FIR system. (b) The zeros are at: ( ( ) ( )) ( ( ) ( )) π π 3 π π,,, cos ± j sin, cos ± j sin With the help of MATLAB: z = ; z = -4/3; z3 = -3/4; z4 = (4/3)*(cos(pi/3) + i* sin(pi/3)); z5 = (4/3)*(cos(pi/3) - i* sin(pi/3)); z6 = (3/4)*(cos(pi/3) + i* sin(pi/3)); z7 = (3/4)*(cos(pi/3) - i* sin(pi/3)); z = [z z z3 z4 z5 z6 z7]; p = [ ]; K=; H = tf(zpk(z, p, K, -)) The transfer function is z 7 z z 4.79z 3 + z z H(z) = z7 = z z +.79z z + z z
4 so that this is an odd symmetric filter with N = 8 and {h n } = {,,,.79,.79,,, } Then H(ze jω ) has a linear phase shift with frequency and the delay Δ is the slope, that is (c) The MATLAB command ( ) H(e jω ) = je jω(n)/) 7 = Ω + π/ d Δ = H(e jω ) = 7/ steps. dω freqz([ ], []) generates the following plot Magnitude (db) Phase (degrees) which shows a crude band-pass characteristic. The phase plot shows a change in phase of 7π/ over a range of π in frequency, corresponding to a delay of 7/ steps. Problem 3: The following MATLAB script was used for this problem: Problem Set 6 -- Problem 3 Enter the signal and contaminate it: t = :.:; signal = sin(4*pi*t); noise = *(rand(size(signal))-.5); 3
5 noisy_signal = signal + noise; Plot the waveforms figure() plot(t,signal, b ) hold on plot(t,noisy_signal, r ) title( Signal and Additive Noise ) hold off Design the filter w_c =.4; w_stop =.48; passband_ripple =.; stopband_ripple =. [N,Wn,BETA,FILTYPE] = kaiserord([w_c w_stop],[ ], [passband_ripple stopband_ripple],) B = fir(n,wn,filtype,kaiser(n+,beta)); A = ; Filter the noisy data and plot the output filtered_signal = filter(b,a,noisy_signal); figure() plot(t,filtered_signal, g ) title( Filtered Signal ) Plot the frequency response and pole-zero plots figure(3) freqz(b,a); figure(4) zplane(b,a) df = /.; f = -5:df:-5+*df; figure(5) plot(f,fftshift(abs(fft(noisy_signal)))) xlabel( Frequency (Hz) ) ylabel( Magnitude F ) title( Magnitude Spectrum of Noisy Signal ) figure(6) plot(f,fftshift(abs(fft(filtered_signal)))) xlabel( Frequency (Hz) ) ylabel( Magnitude F ) title( Magnitude Spectrum of Low-Pass Filtered Signal ) This script gave the filter parameters from kaiserord() as: N = 75, β = 3.395, ω n =.45. The following plots were produced: 4
6 Signal and Additive Noise Filtered Signal The following two spectra have been prepared using fftshift() to center them on Hz, 5 Magnitude Spectrum of Noisy Signal Magnitude F Frequency (Hz) 5
7 5 Magnitude Spectrum of Low Pass Filtered Signal Magnitude F Frequency (Hz) 5 Magnitude (db) Phase (degrees) Imaginary Part Real Part Note that the pole-zero plot is composed with zeros on the unit-circle at angles Ω corresponding to frequencies above the cut-off frequency thus forcing the magnitude response to zero at those 6
8 frequencies - this is therefore a low-pass filter. It is interesting to note that below the cut-off frequency, the zeros appear to be in pairs where their magnitudes are reciprocals - one inside and one outside the unit circle (in addition to be in complex conjugate pairs as well). Think about what this means (it also is evident in the pole-zero plot given in Problem, and we will see the same thing in the results of Problems 4 and 5). Problem 4: The MATLAB script from Problem 3 was modified slightly for this problem to design a band-pass filter: Problem Set 6 -- Problem 4 Enter the signal and contaminate it: t = :.:; signal = sin(4*pi*t); noise = *(rand(size(signal))-.5); noisy_signal = signal + noise; Plot the waveforms figure() plot(t,signal, b ) hold on plot(t,noisy_signal, r ) title( Signal and Additive Noise ) hold off Design the filter Frequency specs - normalized to the Nyquist frequency w_sl=.35; w_cl=.39; w_cu =.4; w_su =.45; Pass and stop band ripple p_rip =.; s_rip =.; [N,Wn,BETA,TYPE] = kaiserord([w_sl w_cl w_cu w_su],[ ], [s_rip p_rip s_rip],) B = fir(n,wn,type, kaiser(n+,beta)); A = ; Filter the noisy data and plot the output filtered_signal = filter(b,a,noisy_signal); figure() plot(t,filtered_signal, g ) title( Filtered Signal ) Plot the frequency response and pole-zero plots figure(3) freqz(b,a); figure(4) zplane(b,a) df = /.; 7
9 f = -5:df:-5+*df; figure(5) plot(f,fftshift(abs(fft(noisy_signal)))) xlabel( Frequency (Hz) ) ylabel( Magnitude F ) title( Magnitude Spectrum of Noisy Signal ) figure(6) plot(f,fftshift(abs(fft(filtered_signal)))) xlabel( Frequency (Hz) ) ylabel( Magnitude F ) title( Magnitude Spectrum of Band-Pass Filtered Signal ) This script gave the filter parameters fromkaiserord() as: N =, β = 3.395, ω n = [.37.43]. The following plots were produced: Signal and Additive Noise Filtered Signal
10 The following two spectra have been prepared using fftshift() to center them on Hz, 5 Magnitude Spectrum of Noisy Signal Magnitude F Frequency (Hz) 45 Magnitude Spectrum of Band Pass Filtered Signal Magnitude F Frequency (Hz) Magnitude (db) Phase (degrees)
11 .5 Imaginary Part Real Part Note that the pole-zero plot is composed with zeros on the unit-circle at angles Ω corresponding to frequencies above and below the cut-off frequency thus forcing the magnitude response to zero at those frequencies - this is therefore a band-pass filter. As in Problem 3, it is interesting to note that within the pass-band the four zeros appear to be in pairs where their magnitudes are reciprocals one inside and one outside the unit circle (in addition to be in complex conjugate pairs as well). Problem 5: The following MATLAB script was used to design and test a band-stop filter: Problem Set 6 -- Problem 5 Enter the signal and contaminate it: Dt = /3; t = :Dt:999*Dt; f_sample = /Dt; f_nyquist = f_sample/; signal =.5*sin(*pi*3*t) + *cos(*pi*9*t); noise = sin(*pi*6*t); noisy_signal = signal + noise; Plot the magnitude spectrum of the contaminated signal f = -5:/(*Dt):-5+999/(*Dt); figure() plot(f,fftshift(abs(fft(noisy_signal)))) title( Magnitude Spectrum of Signal with 6 Hz Contamination ) xlabel( Frequency (Hz) ) ylabel( Magnitude F_m ) f_norm = 6/f_nyquist; Design the filter w_sl = f_norm*.9; w_cl = f_norm*.7; w_cu = f_norm*.3; w_su = f_norm*.;
12 p_ripple =.; s_ripple =.; [N,Wn,BETA,TYPE] = kaiserord([w_cl w_sl w_su w_cu],[ ], [p_rip s_rip p_rip],) B = fir(n,wn,type, kaiser(n+,beta)); A = ; Filter the noisy data and plot the output magnitude spectrum figure() filtered_signal = filter(b,a,noisy_signal); plot(f, fftshift(abs(fft(filtered_signal)))) title( Magnitude Spectrum of the Filtered Signal ) xlabel( Frequency (Hz) ) ylabel( Magnitude F_m ) Plot the frequency response and pole-zero plots figure(3) freqz(b,a); figure(4) zplane(b,a) Extra plots - not asked for! Plot segments of the input and filtered signals First interpolate the signals ant time by a factor of eight to make the plots more readable. s_8 = resample(signal,8,); ns_8 = resample(noisy_signal,8,); fs_8 = resample(filtered_signal,8,); t_8 = resample(t,8,); Plot the original signal and the contaminated signal figure(5) plot(t_8(3:5), s_8(3:5), b--,t_8(3:5), ns_8(3:5), r ) xlabel( Time (sec) ) ylabel( Signal ) title( Segment of Signal and Contaminated Signal ) figure(6) plot(t_8(3:5), s_8(3:5), b--,t_8(3:5), fs_8(3:5), r ) xlabel( Time (sec) ) ylabel( Signal ) title( Segment of Uncontaminated Signal and Filter Output ) figure(7) plot(t_8(:5), s_8(:5), b--,t_8(:5), fs_8(:5), r ) xlabel( Time (sec) ) ylabel( Signal ) title( Initial Segment of Uncontaminated Signal and Filter Output ) This script gave the filter parameters from kaiserord() as: N = 56, β = 3.395, ω cl = 36 rad/s, ω cu = 48 rad/s. The following plots were produced
13 5 Magnitude Spectrum of Signal with 6 Hz Contamination Phase (degrees) Magnitude F Magnitude F m m Magnitude (db) Frequency (Hz) Magnitude Spectrum of the Filtered Signal Frequency (Hz)
14 .5 Imaginary Part Real Part The following plots were not asked for in the problem. To look at the waveforms it was necessary to interpolate the data records to show details of the waveforms, using resample(). See the script for details. 4 Segment of Signal and Contaminated Signal 3 Signal Time (sec) Desired signal Signal contaminated with 6 Hz component Segment of Uncontaminated Signal and Filter Output 3 Signal 3 Uncontaminated input signal Filter output Time (sec) 3
15 4 Initial Segment of Uncontaminated Signal and Filter Output 3 Signal Time (sec) Uncontaminated input signal Filter output Note the delay in the filter output. 4
Digital 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 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 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 informationDigital Signal Processing: Signal Transforms
Digital Signal Processing: Signal Transforms Aishy Amer, Mohammed Ghazal January 19, 1 Instructions: 1. This tutorial introduces frequency analysis in Matlab using the Fourier and z transforms.. More Matlab
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 informationTTT4120 Digital Signal Processing Suggested Solutions for Problem Set 2
Norwegian University of Science and Technology Department of Electronics and Telecommunications TTT42 Digital Signal Processing Suggested Solutions for Problem Set 2 Problem (a) The spectrum X(ω) can be
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 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 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 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. I Reading:
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 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 informationLab 4: Quantization, Oversampling, and Noise Shaping
Lab 4: Quantization, Oversampling, and Noise Shaping Due Friday 04/21/17 Overview: This assignment should be completed with your assigned lab partner(s). Each group must turn in a report composed using
More informationDepartment of Electrical and Computer Engineering Digital Speech Processing Homework No. 6 Solutions
Problem 1 Department of Electrical and Computer Engineering Digital Speech Processing Homework No. 6 Solutions The complex cepstrum, ˆx[n], of a sequence x[n] is the inverse Fourier transform of the complex
More information2.161 Signal Processing: Continuous and Discrete
MIT 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 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 information6.003 (Fall 2011) Quiz #3 November 16, 2011
6.003 (Fall 2011) Quiz #3 November 16, 2011 Name: Kerberos Username: Please circle your section number: Section Time 2 11 am 3 1 pm 4 2 pm Grades will be determined by the correctness of your answers (explanations
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 02 DSP Fundamentals 14/01/21 http://www.ee.unlv.edu/~b1morris/ee482/
More informationResponses of Digital Filters Chapter Intended Learning Outcomes:
Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing
More information2. Typical Discrete-Time Systems All-Pass Systems (5.5) 2.2. Minimum-Phase Systems (5.6) 2.3. Generalized Linear-Phase Systems (5.
. Typical Discrete-Time Systems.1. All-Pass Systems (5.5).. Minimum-Phase Systems (5.6).3. Generalized Linear-Phase Systems (5.7) .1. All-Pass Systems An all-pass system is defined as a system which has
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 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 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 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 informationx[n] = x a (nt ) x a (t)e jωt dt while the discrete time signal x[n] has the discrete-time Fourier transform x[n]e jωn
Sampling Let x a (t) be a continuous time signal. The signal is sampled by taking the signal value at intervals of time T to get The signal x(t) has a Fourier transform x[n] = x a (nt ) X a (Ω) = x a (t)e
More informationLab 3 & 4. LTI Systems, the Z-Transform, and an Introduction
E E 2 7 5 Lab January 28, 2015 Lab 3 & 4. LTI Systems, the Z-Transform, and an Introduction to Filtering Discrete Linear Time Invariant Systems - Short Review x[n] h[n] y[n] A system is: linear if y[n]
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 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 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 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 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
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 informationVery useful for designing and analyzing signal processing systems
z-transform z-transform The z-transform generalizes the Discrete-Time Fourier Transform (DTFT) for analyzing infinite-length signals and systems Very useful for designing and analyzing signal processing
More informationTransform Representation of Signals
C H A P T E R 3 Transform Representation of Signals and LTI Systems As you have seen in your prior studies of signals and systems, and as emphasized in the review in Chapter 2, transforms play a central
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 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 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 informationDigital Signal Processing Lecture 10 - Discrete Fourier Transform
Digital Signal Processing - Discrete Fourier Transform Electrical Engineering and Computer Science University of Tennessee, Knoxville November 12, 2015 Overview 1 2 3 4 Review - 1 Introduction Discrete-time
More informationFinal Exam of ECE301, Prof. Wang s section 1 3pm Tuesday, December 11, 2012, Lily 1105.
Final Exam of ECE301, Prof. Wang s section 1 3pm Tuesday, December 11, 2012, Lily 1105. 1. Please make sure that it is your name printed on the exam booklet. Enter your student ID number, e-mail address,
More informationThe Discrete-Time Fourier
Chapter 3 The Discrete-Time Fourier Transform 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 3-1-1 Continuous-Time Fourier Transform Definition The CTFT of
More informationComplex symmetry Signals and Systems Fall 2015
18-90 Signals and Systems Fall 015 Complex symmetry 1. Complex symmetry This section deals with the complex symmetry property. As an example I will use the DTFT for a aperiodic discrete-time signal. The
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 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 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 information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.161 Signal rocessing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts
More informationHomework 6 Solutions
8-290 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 208 Homework 6 Solutions. Part One. (2 points) Consider an LTI system with impulse response h(t) e αt u(t), (a) Compute the frequency response
More informationFinal Exam of ECE301, Prof. Wang s section 8 10am Tuesday, May 6, 2014, EE 129.
Final Exam of ECE301, Prof. Wang s section 8 10am Tuesday, May 6, 2014, EE 129. 1. Please make sure that it is your name printed on the exam booklet. Enter your student ID number, e-mail address, and signature
More informationThis homework will not be collected or graded. It is intended to help you practice for the final exam. Solutions will be posted.
6.003 Homework #14 This homework will not be collected or graded. It is intended to help you practice for the final exam. Solutions will be posted. Problems 1. Neural signals The following figure illustrates
More informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More 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 informationAssignment 4 Solutions Continuous-Time Fourier Transform
Assignment 4 Solutions Continuous-Time Fourier Transform ECE 3 Signals and Systems II Version 1.01 Spring 006 1. Properties of complex numbers. Let c 1 α 1 + jβ 1 and c α + jβ be two complex numbers. a.
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 informationChapter 13 Z Transform
Chapter 13 Z Transform 1. -transform 2. Inverse -transform 3. Properties of -transform 4. Solution to Difference Equation 5. Calculating output using -transform 6. DTFT and -transform 7. Stability Analysis
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 informationCMPT 889: Lecture 5 Filters
CMPT 889: Lecture 5 Filters Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University October 7, 2009 1 Digital Filters Any medium through which a signal passes may be regarded
More informationDigital Filters. Linearity and Time Invariance. Linear Time-Invariant (LTI) Filters: CMPT 889: Lecture 5 Filters
Digital Filters CMPT 889: Lecture 5 Filters Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University October 7, 29 Any medium through which a signal passes may be regarded as
More informationsummable Necessary and sufficient for BIBO stability of an LTI system. Also see poles.
EECS 206 DSP GLOSSARY c Andrew E. Yagle Fall 2005 absolutely impulse response: h[n] is finite. EX: n=0 ( 3 4 )n = 1 = 4 but 1 3 n=1 1 n. 4 summable Necessary and sufficient for BIBO stability of an LI
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 informationLAB 2: DTFT, DFT, and DFT Spectral Analysis 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 2: DTFT, DFT, and DFT Spectral
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 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 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 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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering Dynamics and Control II Fall 2007
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering.4 Dynamics and Control II Fall 7 Problem Set #9 Solution Posted: Sunday, Dec., 7. The.4 Tower system. The system parameters are
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 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 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 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 information2.161 Signal Processing: Continuous and Discrete
MIT OpenCourseWare http://ocw.mit.edu.6 Signal Processing: Continuous and Discrete Fall 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. M MASSACHUSETTS
More information7.17. Determine the z-transform and ROC for the following time signals: Sketch the ROC, poles, and zeros in the z-plane. X(z) = x[n]z n.
Solutions to Additional Problems 7.7. Determine the -transform and ROC for the following time signals: Sketch the ROC, poles, and eros in the -plane. (a) x[n] δ[n k], k > 0 X() x[n] n n k, 0 Im k multiple
More informationX. Chen More on Sampling
X. Chen More on Sampling 9 More on Sampling 9.1 Notations denotes the sampling time in second. Ω s = 2π/ and Ω s /2 are, respectively, the sampling frequency and Nyquist frequency in rad/sec. Ω and ω denote,
More information2.004 Dynamics and Control II Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 2.004 Dynamics and Control II Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8 I * * Massachusetts
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 informationGeneralizing the DTFT!
The Transform Generaliing the DTFT! The forward DTFT is defined by X e jω ( ) = x n e jωn in which n= Ω is discrete-time radian frequency, a real variable. The quantity e jωn is then a complex sinusoid
More informationFinal Exam of ECE301, Section 3 (CRN ) 8 10am, Wednesday, December 13, 2017, Hiler Thtr.
Final Exam of ECE301, Section 3 (CRN 17101-003) 8 10am, Wednesday, December 13, 2017, Hiler Thtr. 1. Please make sure that it is your name printed on the exam booklet. Enter your student ID number, and
More informationLecture 8 Finite Impulse Response Filters
Lecture 8 Finite Impulse Response Filters Outline 8. Finite Impulse Response Filters.......................... 8. oving Average Filter............................... 8.. Phase response...............................
More 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 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 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 information2.004 Dynamics and Control II Spring 2008
MT OpenCourseWare http://ocw.mit.edu 2.004 Dynamics and Control Spring 2008 or information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Reading: ise: Chapter 8 Massachusetts
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 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 informationLecture 2 OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE
OKAN UNIVERSITY FACULTY OF ENGINEERING AND ARCHITECTURE EEE 43 DIGITAL SIGNAL PROCESSING (DSP) 2 DIFFERENCE EQUATIONS AND THE Z- TRANSFORM FALL 22 Yrd. Doç. Dr. Didem Kivanc Tureli didemk@ieee.org didem.kivanc@okan.edu.tr
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 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 informationELEN E4810: Digital Signal Processing Topic 11: Continuous Signals. 1. Sampling and Reconstruction 2. Quantization
ELEN E4810: Digital Signal Processing Topic 11: Continuous Signals 1. Sampling and Reconstruction 2. Quantization 1 1. Sampling & Reconstruction DSP must interact with an analog world: A to D D to A x(t)
More 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 information! Introduction. ! Discrete Time Signals & Systems. ! Z-Transform. ! Inverse Z-Transform. ! Sampling of Continuous Time Signals
ESE 531: Digital Signal Processing Lec 25: April 24, 2018 Review Course Content! Introduction! Discrete Time Signals & Systems! Discrete Time Fourier Transform! Z-Transform! Inverse Z-Transform! Sampling
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 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 informationDigital Signal Processing. Midterm 1 Solution
EE 123 University of California, Berkeley Anant Sahai February 15, 27 Digital Signal Processing Instructions Midterm 1 Solution Total time allowed for the exam is 8 minutes Some useful formulas: Discrete
More 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 informationDSP Laboratory (EELE 4110) Lab#5 DTFS & DTFT
Islamic University of Gaza Faculty of Engineering Electrical Engineering Department EG.MOHAMMED ELASMER Spring-22 DSP Laboratory (EELE 4) Lab#5 DTFS & DTFT Discrete-Time Fourier Series (DTFS) The discrete-time
More information(i) Understanding the characteristics and properties of DTFT
Discrete-Time Fourier Transform (DTFT) Chapter Intended Learning Outcomes: (i) Understanding the characteristics and properties of DTFT (ii) Ability to perform discrete-time signal conversion between the
More information6.003: Signals and Systems
6.003: Signals and Systems Discrete-Time Frequency Representations November 8, 2011 1 Mid-term Examination #3 Wednesday, November 16, 7:30-9:30pm, No recitations on the day of the exam. Coverage: Lectures
More informationZ-TRANSFORMS. Solution: Using the definition (5.1.2), we find: for case (b). y(n)= h(n) x(n) Y(z)= H(z)X(z) (convolution) (5.1.
84 5. Z-TRANSFORMS 5 z-transforms Solution: Using the definition (5..2), we find: for case (a), and H(z) h 0 + h z + h 2 z 2 + h 3 z 3 2 + 3z + 5z 2 + 2z 3 H(z) h 0 + h z + h 2 z 2 + h 3 z 3 + h 4 z 4
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 informationECE-314 Fall 2012 Review Questions for Midterm Examination II
ECE-314 Fall 2012 Review Questions for Midterm Examination II First, make sure you study all the problems and their solutions from homework sets 4-7. Then work on the following additional problems. Problem
More information2.004 Dynamics and Control II Spring 2008
MT OpenCourseWare http://ocw.mit.edu.004 Dynamics and Control Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Massachusetts nstitute of Technology
More informationSinc Functions. Continuous-Time Rectangular Pulse
Sinc Functions The Cooper Union Department of Electrical Engineering ECE114 Digital Signal Processing Lecture Notes: Sinc Functions and Sampling Theory October 7, 2011 A rectangular pulse in time/frequency
More information