PART 1. Review of DSP. f (t)e iωt dt. F(ω) = f (t) = 1 2π. F(ω)e iωt dω. f (t) F (ω) The Fourier Transform. Fourier Transform.
|
|
- Madison Hopkins
- 5 years ago
- Views:
Transcription
1 PART 1 Review of DSP Mauricio Sacchi University of Alberta, Edmonton, AB, Canada The Fourier Transform F() = f (t) = 1 2π f (t)e it dt F()e it d Fourier Transform Inverse Transform f (t) F () Part 1 Review of DSP 2 1
2 Forward / Inverse Pairs FT -1 f (t) F() FT Part 1 Review of DSP 3 Amplitude and Phase Amplitude and Phase of a complex number X = R + i G α = tan 1 G R G A X = R + i G A = (R 2 + G 2 ) α R Part 1 Review of DSP 4 2
3 Amplitude and Phase Amplitude and Phase of the FT F() = R() + i G() α() = tan 1 G() R() G A X = R + i G A() = (R() 2 + G() 2 ) α R Part 1 Review of DSP 5 Symmetries of the FT If the signal is real, then R() = R( ) Even A() = A( ) G() = G( ) Odd α() = α() Part 1 Review of DSP 6 3
4 FT: a simple physical interpretation A signal can be represented as a superposition of elementary signals (complex exponentials) of frequency k scaled by a complex amplitude F( k ) f (t) = 1 2π F()e it d Δ 2π f (t) f k (t), f k (t) = Δ 2π F( k )ei k t k k F( k )e i k t Part 1 Review of DSP 7 Properties of the FT f (t) F () g(t) G() Linearity Scale Shifting Modulation in time Convolution Convolution in Frequency f (t) + g(t) F() + G() f (ta) 1 a F ( a ) f (t τ ) F ()e iτ f (t)e iϕ t F( ϕ) f (τ )g(t τ )dτ F()G() f (t)g(t) 1 2π F(ϕ)G( ϕ)dϕ Part 1 Review of DSP 8 4
5 Delta function f (t) F () δ (t) 1 f (t) = δ (t) F() = 1 t 1 represents the delta function which cannot be drawn Part 1 Review of DSP 9 Cosine f (t) = cos( 0 t), < t < f (t) R() t G() f (t) F() = R() + ig() cos( 0 t) πδ ( 0 ) + πδ ( + 0 ) Part 1 Review of DSP 10 5
6 Sine f (t) = sin( 0 t), < t < f (t) R() t f (t) F() = R() + ig() G() sin( 0 t) i πδ ( 0 ) + i πδ ( + 0 ) Part 1 Review of DSP 11 Truncation: zeros or??? Atmospheric pressure FCAG - UNLP - La Plata ( ) Part 1 Review of DSP 12 6
7 Boxcar FT of the truncation operator (Boxcar) t t Part 1 Review of DSP 13 The truncation problem Time Frequency Desire Aperture X t 0 * 2 Peaks Acquisition window t 0 = = 1 Peak Data t 0 Part 1 Review of DSP 14 7
8 The discrete world Analog signals (waveforms) are transformed into digital signals by acquisition systems How the FT of the true underlying continue signal/process relates to its discrete version?? This is answered by Nyquist theorem Acquisition System Time (s) Sample * Part 1 Review of DSP 15 The discrete world Analog signal Digital Signal sampled every Δt secs s(t) S() s(t n ) S d () t n = (n 1) Δt Part 1 Review of DSP 16 8
9 Nyquist theorem The FT of the discrete signal is a distorted version of the FT of the analog signal. The distortion is given by Poisson Formula: This formula can be found in any book on harmonic analiyis S d () = 1 Δt S( k 0 ), 0 = 2π Δt k= What you can measure S() is what you would have liked to measure Part 1 Review of DSP 17 Nyquist theorem Nyquist theorem or formula provides the sampling condition to compute the FT of the discrete signal is such a way that it is a perfect representation of the FT of the analog signal. The theorem is derived by simple inspection of Poisson formula. In a graphical manner: S d () = 1 S( k 0 ), 0 = 2π Δt Δt k= S() max S d () max S( ) S( + 0 ) S() S( 0 ) S( 2 0 ) max max Non-aliased Spectrum 2 max < 0 Part 1 Review of DSP 18 9
10 Nyquist theorem S d () = 1 S( k 0 ), 0 = 2π Δt Δt k= S() max S d () max S( ) S( + 0 ) S() S( 0 ) S( 2 0 ) max max Non-aliased Spectrum 2 max = 0 Part 1 Review of DSP 19 Nyquist theorem S d () = 1 S( k 0 ), 0 = 2π Δt Δt k= S() max S d () max S( ) S( + 0 ) S() S( 0 ) S( 2 0 ) Non-aliased Spectrum 2 max > 0 Part 1 Review of DSP 20 10
11 Nyquist theorem S d () = 1 S( k 0 ), 0 = 2π Δt Δt k= S() Alias, the true spectrum is distorted!! S( ) S( + 0 ) max S d () S() max S( 0 ) S( 2 0 ) Non-aliased Spectrum 2 max > 0 Part 1 Review of DSP 21 Nyquist theorem From the previous figures we have found the condition to avoid aliasing: 2 max < 0 0 = 2π /Δt max < π /Δt If we prefer to use frequency (Hz) rather than angular frequency (rad/sec): 1 2π f max < π /Δt Δt < 2 f max Famous Nyquist condition Part 1 Review of DSP 22 11
12 Nyquist theorem From now on, we consider signals where the sampling interval satisfies the Nyquist condition. It is clear that discrete signals must arise from the discretization of a band-limited analog signal. Electronic filters are often placed prior to discretization to guarantee that the signal to sample does contain not energy above a maximum frequency. Nyquist condition is easy to satisfied in the time domain (temporal sampling) Spatial sampling is often dictated by cost & logistics not by hardware!! Δx Multi-dimensional sampling in space is a problem of current research since prestack seismic data are often under sampled in one or more coordinates (4 spatial coordinates) x Part 1 Review of DSP 23 DFT When dealing with discrete time series or evenly sample data along the spatial domain we will use the Discrete Fourier Transform (DFT) S() = N 1 k= 0 s k e ik : angular frequency [rads, no dimensions] l = 2π l N S( l ) = N 1 k= 0, l = 0,...N 1 discrete angular frequency s k e i l k Example N=8 = π, k = 4 = π /2,k = 2 = 0,k = 0 = π /2,k = 6 Part 1 Review of DSP 24 12
13 IDFT We also need a transform to come back Inverse Discrete Fourier Transform (IDFT) s k = 1 N A note about frequency N 1 k= 0 S( l )e i l k l = 2π l N l = 2π l N Δt,, l = 0,...N 1 f l = l 2π = l N Δt, discrete angular frequency radians/secs Hertz Part 1 Review of DSP 25 Notes Wrong wording à The FFT Spectrum, You should say the DFT Spectrum because the FFT is just the tool that is used to compute the DFT in a fast way Remember that to apply the DFT is equivalent to multiply a Matrix times a Vector (N 2 operations) FFT is a simple matrix multiplication via a faster algorithm ( N log 2 N operations ) Part 1 Review of DSP 26 13
14 Linear systems Linear systems An easy way of describing physical phenomena A good approximation to some inverse problems in geophysics Given x 1 (t) y 1 (t) x 2 (t) y 2 (t) The system is linear if αx 1 (t) + β x 2 (t) α y 1 (t) + β y 2 (t) Part 1 Review of DSP 27 Linear systems x 1 (t) x 2 (t) LS LS y 1 (t) y 2 (t) αx 1 (t) + β x 2 (t) LS α y 1 (t) + β y 2 (t) LS: Linear System (the Earth, if you do not consider important phenomena!) Part 1 Review of DSP 28 14
15 Linear systems and Invariance Invariance [Linear Time Invariant System] Consider a system that is linear and also impose the condition of invariance: x (t) y(t) x(t τ ) y(t τ ) x(t) x(t τ ) τ LTIS LTIS τ y(t) y(t τ ) Example: deconvolution operator Part 1 Review of DSP 29 Linear systems and Invariance If the system is linear and time invariant, input and output are related by the following expression (it can be proven) y(t) = h(t τ )x(τ )dτ = h(t) * x(t) Convolution Symbol We are saying that if our process is represented by an LTIS then the I/O can be represented via a convolution integral The new signal h(t) is called the impulse response of the system Part 1 Review of DSP 30 15
16 Linear systems and Invariance Impulse response (hitting the system with an impulse) y(t) = h(t τ )x(τ )dτ = h(t) * x(t) x(t) = δ (t) x(t) h(t) h(t) h(t) y(t) Input δ (t) h(t) Input Part 1 Review of DSP 31 Linear systems and Invariance - Discrete case Convolution Sum y n = h k n x k = h n * x n k Signals are time series or vectors x = [1,0,0,0,0,0...] h h = [h 0,h 1,h 2,...] x = [ x 0, x 1, x 2,...] h y = [ y 0, y 1, y 2,...] Part 1 Review of DSP 32 16
17 Discrete convolution Formula Finite length signals y n = h n k x k = h n * x n k x k, k = 0,NX 1 y k, h k, k = 0,NY 1 k = 0,NH 1 How do we do the convolution with finite length signals? With paper and pencil Computer code Matrix times vector Polinomial multiplication DFT Part 1 Review of DSP 33 Discrete convolution % Initialize output y(1: NX + HH -1) = 0 % Do convolution sum for i = 1: NX for j = 1: NH y(i+ j -1) = y(i+ j -1) + x(i)h(j) end end Part 1 Review of DSP 34 17
18 Discrete convolution Example: x = [ x 0, x 1, x 2, x 3, x 4 ], NX = 5 h = [h 0,h 1,h 2 ], NH = 3 y 0 = x 0 h 0 y 1 = x 1 h 0 + x 0 h 1 y 2 = x 2 h 0 + x 1 h 1 + x 0 h 2 y 3 = x 3 h 0 + x 2 h 1 + x 1 h 2 y 4 = x 4 h 0 + x 3 h 1 + x 2 h 2 y 5 = x 4 h 1 + x 3 h 2 y 6 = x 4 h 2 " # y 0 y 1 y 2 y 3 y 4 y 5 y 6 y n = h k n x k = h n * x n k % " x % ' ' ' x 1 x 0 0 ' ' x 2 x 1 x 0 ' " ' ' ' = x 3 x 2 x 1 ' ' x 4 x 3 x 2 ' # ' ' 0 x 4 x 3 ' ' & # 0 0 x 4 & h 0 h 1 h 2 % ' ' & Part 1 Review of DSP 35 Transient-free Convolution Matrix " # y 0 y 1 y 2 y 3 y 4 y 5 y 6 % " x % ' ' ' x 1 x 0 0 ' ' x 2 x 1 x 0 ' " ' ' ' = x 3 x 2 x 1 ' ' x 4 x 3 x 2 ' # ' ' 0 x 4 x 3 ' ' & # 0 0 x 4 & Classical convolution h 0 h 1 h 2 % ' ' & " # y 2 y 3 y 4 % " x 2 x 1 x 0 % " ' ' = ' x 3 x 2 x 1 ' & # x 4 x 3 x 2 & # h 0 h 1 h 2 Transient-free convolution % ' ' & Part 1 Review of DSP 36 18
19 Discrete convolution and the z-transform Z-transform: a compact way of dealing with time series The z-transform of x = [ x 0, x 1, x 2, x 3, x 4 ], NX = 5 is given by X(z) = x 0 + x 1 z + x 2 z 2 + x 3 z 3 + x 4 z 4 Example: x = [2, 1,3] Casual x = [ 1,2, 1,3] Non-causal X(z) = 2 1z + 3z 2 X(z) = 1z z + 3z 2 Indicates sample n=0 Part 1 Review of DSP 37 What can we do with the z-transform? Convolve series x = [ x 0, x 1, x 2, x 3, x 4 ] h = [h 0,h 1,h 2 ] y = x * h X(z) = x 0 + x 1 z + x 2 z 2 + x 3 z 3 + x 4 z 4 H (z) = h 0 + h 1 z + h 2 z Y (z) = X(z).H (z) Design inverse filters (finally some seismology ) x Unknown Filter h y Let s see how one can use the z-transform to find Inverse Filters of simple signals Part 1 Review of DSP 38 19
20 Dipoles and inverse of a dipole Dipole: a signal made of two elements x = [1,a] h y = [1,0] Unknown Inverse filter that turns the dipole into a spike Find the inverse filter with the z-transform: X(z) = 1+ az, Y (z) = 1 Geometric Series y = x * h Y (z) = X (z).h (z) H (z) = Y (z) X (z) = 1 1+ az = 1 az + a 2 z 2 a 3 z 3 + a 4 z 4... h = [1, a,a 2, a 3,a 4,...] Part 1 Review of DSP 39 Inversion of a dipole using geometric series x = [1,a] h y a = 0.5 a = 0.5 Part 1 Review of DSP 40 20
21 Inversion of a dipole using geometric series: Truncation of the operator x = [1,a] h y a = 0.99 a = 0.99 Part 1 Review of DSP 41 Inversion of a dipole using geometric series: Deconvolution of a simple reflectivity series x = [1,a] h y Filter design a = 0.5 r s = x * r r = h * s Filter Application Part 1 Review of DSP 42 21
22 Inversion of a dipole using geometric series: Deconvolution of a simple reflectivity series Truncation in the operator introduces false reflections in the deconvolution output x = [1,a] h y Filter design a = 0.99 r s = x * r r = h * s Wrong estimate of the reflectivity Filter Application Part 1 Review of DSP 43 Minimum and Maximum Phase dipoles In simple terms Minimum Phase Maximum Phase x = [1,a], a < 1 x = [1,a], a > 1 Min Max Part 1 Review of DSP 44 22
23 Dipoles and Phase duality Take two dipoles x MIN = [1,a], a < 1 x MAX = [a,1] = a[1,1/a] = a[1,b], b > 1 You can show that X MAX () = X MIN () θ MAX () θ MIN () Same amplitude spectrum Different phase spectrum If only the amplitude spectrum is measured, one cannot uniquely determine the dipole (two dipoles produce the same amplitude) Part 1 Review of DSP 45 Dipoles and Phase duality x MIN = [1,a], a < 1 x MAX = [a,1] = a[1,1/a] = a[1,b], b > 1 Part 1 Review of DSP 46 23
24 Dipole filters - careful here Some signal processing schemes attempt o increase BW by convolution with dipole filters. The amplitude spectrum of the dipole filter can be Low Pass or High Pass according to the sign of a Examples: Low Pass x MIN = [1,a], a = 0.9 High Pass x MIN = [1,a], a = 0.9 Part 1 Review of DSP 47 Dipole filters - careful here Differentiator (Extreme High Pass dipole) Examples: Wavelet convolved n times with differentiator - Cosmetic freq. enhancement?? N=4 High Pass x = [1,a], a = 1 N=1 Wavelet Part 1 Review of DSP 48 24
25 Dipole filters - careful here N=2 (two differentiations) Reflectivity Data Data after 2 differentiations Part 1 Review of DSP 49 Dipole filters - careful here N=2 (two differentiations) Oops!! Reflectivity Data Data after 2 differentiations Part 1 Review of DSP 50 25
26 More about dipoles: Spectral Decomposition Some modern seismic interpretation methods are based on properties of dipoles filters Spectral Decomposition attempts to image thing layers by the spectral behaviour of signals similar to dipoles Impedance Reflectivity Trace Amplitude Spectrum Time Time This is the basis of spectral decomposition. Partyka, G., 2005, Spectral Decomposition: Recorder, 30 ( Time Frequency Spectral notch is proportional to layer thickness One, for instance, can map the amplitude at that particular frequency. This will provide an attribute for the x-y variability of layer thickness Interesting to point out that rather than whitening (flattening) the spectrum like in conventional decon, spectral decomposition attempts to track spectral features/attributes Part 1 Review of DSP 51 More about dipoles: Spectral Decomposition τ = 4.Δt Thin Layer r = [a,0,0,0,0,b,0,0,0...] R() = a + be iτ Spectrum R() 2 = a 2 + b a b cos(τ ) Min/Max condition Frequency at stationary point d R() 2 = 0 sin(τ ) = 0 τ = πk, k = 0,1,2,3,4 d f s = k /(2τ ) The second derivative can be used to determine if the stationary point is a min or max. Min or max depends on the signs of the reflection coefficients a and b. Part 1 Review of DSP 52 26
27 More about dipoles: Spectral Decomposition For trace #3: τ = 0.028s T (s) f (Hz) f s = 17.8, 35.7, 53.6, 71.4Hz Part 1 Review of DSP 53 More about dipoles: Spectral Decomposition For trace #3: τ = 0.028s T (s) f (Hz) f s = 17.8, 35.7, 53.6, 71.4Hz Part 1 Review of DSP 54 27
28 More about dipoles: Spectral Decomposition For trace #3: τ = 0.028s T (s) f (Hz) f s = 17.8, 35.7, 53.6, 71.4Hz Part 1 Review of DSP 55 28
SEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More 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 informationHomework 4. May An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt
Homework 4 May 2017 1. An LTI system has an input, x(t) and output y(t) related through the equation y(t) = t e (t t ) x(t 2)dt Determine the impulse response of the system. Rewriting as y(t) = t e (t
More informationESS Finite Impulse Response Filters and the Z-transform
9. Finite Impulse Response Filters and the Z-transform We are going to have two lectures on filters you can find much more material in Bob Crosson s notes. In the first lecture we will focus on some of
More informationTherefore the new Fourier coefficients are. Module 2 : Signals in Frequency Domain Problem Set 2. Problem 1
Module 2 : Signals in Frequency Domain Problem Set 2 Problem 1 Let be a periodic signal with fundamental period T and Fourier series coefficients. Derive the Fourier series coefficients of each of the
More informationTime and Spatial Series and Transforms
Time and Spatial Series and Transforms Z- and Fourier transforms Gibbs' phenomenon Transforms and linear algebra Wavelet transforms Reading: Sheriff and Geldart, Chapter 15 Z-Transform Consider a digitized
More informationInformation and Communications Security: Encryption and Information Hiding
Short Course on Information and Communications Security: Encryption and Information Hiding Tuesday, 10 March Friday, 13 March, 2015 Lecture 5: Signal Analysis Contents The complex exponential The complex
More informationFourier transform. Stefano Ferrari. Università degli Studi di Milano Methods for Image Processing. academic year
Fourier transform Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 27 28 Function transforms Sometimes, operating on a class of functions
More informationLecture 27 Frequency Response 2
Lecture 27 Frequency Response 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/12 1 Application of Ideal Filters Suppose we can generate a square wave with a fundamental period
More information06/12/ rws/jMc- modif SuFY10 (MPF) - Textbook Section IX 1
IV. Continuous-Time Signals & LTI Systems [p. 3] Analog signal definition [p. 4] Periodic signal [p. 5] One-sided signal [p. 6] Finite length signal [p. 7] Impulse function [p. 9] Sampling property [p.11]
More informationEE 224 Signals and Systems I Review 1/10
EE 224 Signals and Systems I Review 1/10 Class Contents Signals and Systems Continuous-Time and Discrete-Time Time-Domain and Frequency Domain (all these dimensions are tightly coupled) SIGNALS SYSTEMS
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 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 informationEC Signals and Systems
UNIT I CLASSIFICATION OF SIGNALS AND SYSTEMS Continuous time signals (CT signals), discrete time signals (DT signals) Step, Ramp, Pulse, Impulse, Exponential 1. Define Unit Impulse Signal [M/J 1], [M/J
More informationENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function
More informationThe Continuous Time Fourier Transform
COMM 401: Signals & Systems Theory Lecture 8 The Continuous Time Fourier Transform Fourier Transform Continuous time CT signals Discrete time DT signals Aperiodic signals nonperiodic periodic signals Aperiodic
More informationSummary of lecture 1. E x = E x =T. X T (e i!t ) which motivates us to define the energy spectrum Φ xx (!) = jx (i!)j 2 Z 1 Z =T. 2 d!
Summary of lecture I Continuous time: FS X FS [n] for periodic signals, FT X (i!) for non-periodic signals. II Discrete time: DTFT X T (e i!t ) III Poisson s summation formula: describes the relationship
More informationEE 3054: Signals, Systems, and Transforms Summer It is observed of some continuous-time LTI system that the input signal.
EE 34: Signals, Systems, and Transforms Summer 7 Test No notes, closed book. Show your work. Simplify your answers. 3. It is observed of some continuous-time LTI system that the input signal = 3 u(t) produces
More informationDIGITAL SIGNAL PROCESSING UNIT 1 SIGNALS AND SYSTEMS 1. What is a continuous and discrete time signal? Continuous time signal: A signal x(t) is said to be continuous if it is defined for all time t. Continuous
More informationNotes for Geophysics 426
Notes for Geophysics 426 Geophysical Signal Processing M. D. Sacchi Department of Physics University of Alberta 2 Contact: Dr M.D.Sacchi Department of Physics, University of Alberta, Edmonton, Canada,
More information4 The Continuous Time Fourier Transform
96 4 The Continuous Time ourier Transform ourier (or frequency domain) analysis turns out to be a tool of even greater usefulness Extension of ourier series representation to aperiodic signals oundation
More informationDEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010
[E2.5] IMPERIAL COLLEGE LONDON DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EXAMINATIONS 2010 EEE/ISE PART II MEng. BEng and ACGI SIGNALS AND LINEAR SYSTEMS Time allowed: 2:00 hours There are FOUR
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 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 informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Filtering in the Frequency Domain http://www.ee.unlv.edu/~b1morris/ecg782/ 2 Outline Background
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationThe (Fast) Fourier Transform
The (Fast) Fourier Transform The Fourier transform (FT) is the analog, for non-periodic functions, of the Fourier series for periodic functions can be considered as a Fourier series in the limit that the
More informationConvolution and Linear Systems
CS 450: Introduction to Digital Signal and Image Processing Bryan Morse BYU Computer Science Introduction Analyzing Systems Goal: analyze a device that turns one signal into another. Notation: f (t) g(t)
More informationComputational Methods for Astrophysics: Fourier Transforms
Computational Methods for Astrophysics: Fourier Transforms John T. Whelan (filling in for Joshua Faber) April 27, 2011 John T. Whelan April 27, 2011 Fourier Transforms 1/13 Fourier Analysis Outline: Fourier
More informationGabor Deconvolution. Gary Margrave and Michael Lamoureux
Gabor Deconvolution Gary Margrave and Michael Lamoureux = Outline The Gabor idea The continuous Gabor transform The discrete Gabor transform A nonstationary trace model Gabor deconvolution Examples = Gabor
More informationReview of Frequency Domain Fourier Series: Continuous periodic frequency components
Today we will review: Review of Frequency Domain Fourier series why we use it trig form & exponential form how to get coefficients for each form Eigenfunctions what they are how they relate to LTI systems
More informationLTI Systems (Continuous & Discrete) - Basics
LTI Systems (Continuous & Discrete) - Basics 1. A system with an input x(t) and output y(t) is described by the relation: y(t) = t. x(t). This system is (a) linear and time-invariant (b) linear and time-varying
More informationEECE 3620: Linear Time-Invariant Systems: Chapter 2
EECE 3620: Linear Time-Invariant Systems: Chapter 2 Prof. K. Chandra ECE, UMASS Lowell September 7, 2016 1 Continuous Time Systems In the context of this course, a system can represent a simple or complex
More informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationIntroduction to Digital Signal Processing
Introduction to Digital Signal Processing 1.1 What is DSP? DSP is a technique of performing the mathematical operations on the signals in digital domain. As real time signals are analog in nature we need
More informationQUESTION BANK SIGNALS AND SYSTEMS (4 th SEM ECE)
QUESTION BANK SIGNALS AND SYSTEMS (4 th SEM ECE) 1. For the signal shown in Fig. 1, find x(2t + 3). i. Fig. 1 2. What is the classification of the systems? 3. What are the Dirichlet s conditions of Fourier
More informationGabor multipliers for pseudodifferential operators and wavefield propagators
Gabor multipliers for pseudodifferential operators and wavefield propagators Michael P. Lamoureux, Gary F. Margrave and Peter C. Gibson ABSTRACT Gabor multipliers We summarize some applications of Gabor
More informationAspects of Continuous- and Discrete-Time Signals and Systems
Aspects of Continuous- and Discrete-Time Signals and Systems C.S. Ramalingam Department of Electrical Engineering IIT Madras C.S. Ramalingam (EE Dept., IIT Madras) Networks and Systems 1 / 45 Scaling the
More informationContinuous-Time Fourier Transform
Signals and Systems Continuous-Time Fourier Transform Chang-Su Kim continuous time discrete time periodic (series) CTFS DTFS aperiodic (transform) CTFT DTFT Lowpass Filtering Blurring or Smoothing Original
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 informationTutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling
2.39 utorial Sheet #2 discrete vs. continuous functions, periodicity, sampling We will encounter two classes of signals in this class, continuous-signals and discrete-signals. he distinct mathematical
More informationChapter 5. Section 5.1. Section ( x ) ( y ) 7. ( x ) ( y ) (0, 3 + 5) and (0, 3 5)
9 Chapter Section.. 0. ( x ) ( y ). ( x 7 ) + ( y+ ) = 9 7. ( x ) ( y ) 8 + + 0 = 8 + 8 = 9.. (0, + ) and (0, ). (.60786, 7.6887). (-.07,.8) 7. 9.87 miles Section. 70 0 -.. 00. 0 7. 9.. 8 9.. miles 7.
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 information3.2 Complex Sinusoids and Frequency Response of LTI Systems
3. Introduction. A signal can be represented as a weighted superposition of complex sinusoids. x(t) or x[n]. LTI system: LTI System Output = A weighted superposition of the system response to each complex
More informationEE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet
EE538 Final Exam Fall 005 3:0 pm -5:0 pm PHYS 3 Dec. 17, 005 Cover Sheet Test Duration: 10 minutes. Open Book but Closed Notes. Calculators ARE allowed!! This test contains five problems. Each of the five
More informationLine Spectra and their Applications
In [ ]: cd matlab pwd Line Spectra and their Applications Scope and Background Reading This session concludes our introduction to Fourier Series. Last time (http://nbviewer.jupyter.org/github/cpjobling/eg-47-
More information1 Understanding Sampling
1 Understanding Sampling Summary. In Part I, we consider the analysis of discrete-time signals. In Chapter 1, we consider how discretizing a signal affects the signal s Fourier transform. We derive the
More informationCorrelator I. Basics. Chapter Introduction. 8.2 Digitization Sampling. D. Anish Roshi
Chapter 8 Correlator I. Basics D. Anish Roshi 8.1 Introduction A radio interferometer measures the mutual coherence function of the electric field due to a given source brightness distribution in the sky.
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 informationLECTURE 12 Sections Introduction to the Fourier series of periodic signals
Signals and Systems I Wednesday, February 11, 29 LECURE 12 Sections 3.1-3.3 Introduction to the Fourier series of periodic signals Chapter 3: Fourier Series of periodic signals 3. Introduction 3.1 Historical
More informationSignal Processing Signal and System Classifications. Chapter 13
Chapter 3 Signal Processing 3.. Signal and System Classifications In general, electrical signals can represent either current or voltage, and may be classified into two main categories: energy signals
More informationENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin
ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM Dr. Lim Chee Chin Outline Introduction Discrete Fourier Series Properties of Discrete Fourier Series Time domain aliasing due to frequency
More informationChapter 8 The Discrete Fourier Transform
Chapter 8 The Discrete Fourier Transform Introduction Representation of periodic sequences: the discrete Fourier series Properties of the DFS The Fourier transform of periodic signals Sampling the Fourier
More informationFourier Transforms For additional information, see the classic book The Fourier Transform and its Applications by Ronald N. Bracewell (which is on the shelves of most radio astronomers) and the Wikipedia
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 informationRepresenting a Signal
The Fourier Series Representing a Signal The convolution method for finding the response of a system to an excitation takes advantage of the linearity and timeinvariance of the system and represents the
More informationENGIN 211, Engineering Math. Fourier Series and Transform
ENGIN 11, Engineering Math Fourier Series and ransform 1 Periodic Functions and Harmonics f(t) Period: a a+ t Frequency: f = 1 Angular velocity (or angular frequency): ω = ππ = π Such a periodic function
More informationModule 4. Related web links and videos. 1. FT and ZT
Module 4 Laplace transforms, ROC, rational systems, Z transform, properties of LT and ZT, rational functions, system properties from ROC, inverse transforms Related web links and videos Sl no Web link
More informationExploring Fourier Transform Techniques with Mathcad Document 1: Introduction to the Fourier Transform
Exploring Fourier Transform Techniques with Mathcad Document : Introduction to the Fourier Transform by Mark Iannone Department of Chemistry Millersville University Millersville, PA 7-3 miannone@maraudermillersvedu
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 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 informationFundamentals of the Discrete Fourier Transform
Seminar presentation at the Politecnico di Milano, Como, November 12, 2012 Fundamentals of the Discrete Fourier Transform Michael G. Sideris sideris@ucalgary.ca Department of Geomatics Engineering University
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 information2A1H Time-Frequency Analysis II
2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 209 For any corrections see the course page DW Murray at www.robots.ox.ac.uk/ dwm/courses/2tf. (a) A signal g(t) with period
More informationCommunications and Signal Processing Spring 2017 MSE Exam
Communications and Signal Processing Spring 2017 MSE Exam Please obtain your Test ID from the following table. You must write your Test ID and name on each of the pages of this exam. A page with missing
More informationCorrelation, discrete Fourier transforms and the power spectral density
Correlation, discrete Fourier transforms and the power spectral density visuals to accompany lectures, notes and m-files by Tak Igusa tigusa@jhu.edu Department of Civil Engineering Johns Hopkins University
More information9.4 Enhancing the SNR of Digitized Signals
9.4 Enhancing the SNR of Digitized Signals stepping and averaging compared to ensemble averaging creating and using Fourier transform digital filters removal of Johnson noise and signal distortion using
More informationBIOSIGNAL PROCESSING. Hee Chan Kim, Ph.D. Department of Biomedical Engineering College of Medicine Seoul National University
BIOSIGNAL PROCESSING Hee Chan Kim, Ph.D. Department of Biomedical Engineering College of Medicine Seoul National University INTRODUCTION Biosignals (biological signals) : space, time, or space-time records
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 informationTopic 7. Convolution, Filters, Correlation, Representation. Bryan Pardo, 2008, Northwestern University EECS 352: Machine Perception of Music and Audio
Topic 7 Convolution, Filters, Correlation, Representation Short time Fourier Transform Break signal into windows Calculate DFT of each window The Spectrogram spectrogram(y,1024,512,1024,fs,'yaxis'); A
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationLecture # 06. Image Processing in Frequency Domain
Digital Image Processing CP-7008 Lecture # 06 Image Processing in Frequency Domain Fall 2011 Outline Fourier Transform Relationship with Image Processing CP-7008: Digital Image Processing Lecture # 6 2
More informationMultidimensional digital signal processing
PSfrag replacements Two-dimensional discrete signals N 1 A 2-D discrete signal (also N called a sequence or array) is a function 2 defined over thex(n set 1 of, n 2 ordered ) pairs of integers: y(nx 1,
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 informationIT DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A
DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING IT6502 - DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A 1. What is a continuous and discrete time signal? Continuous
More informationUNIT 1. SIGNALS AND SYSTEM
Page no: 1 UNIT 1. SIGNALS AND SYSTEM INTRODUCTION A SIGNAL is defined as any physical quantity that changes with time, distance, speed, position, pressure, temperature or some other quantity. A SIGNAL
More informationHomework 9 Solutions
8-290 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 207 Homework 9 Solutions Part One. (6 points) Compute the convolution of the following continuous-time aperiodic signals. (Hint: Use the
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 informationThe Continuous-time Fourier
The Continuous-time Fourier Transform Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline Representation of Aperiodic signals:
More informationsinc function T=1 sec T=2 sec angle(f(w)) angle(f(w))
T=1 sec sinc function 3 angle(f(w)) T=2 sec angle(f(w)) 1 A quick script to plot mag & phase in MATLAB w=0:0.2:50; Real exponential func b=5; Fourier transform (filter) F=1.0./(b+j*w); subplot(211), plot(w,
More informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
More 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 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 informationFourier series for continuous and discrete time signals
8-9 Signals and Systems Fall 5 Fourier series for continuous and discrete time signals The road to Fourier : Two weeks ago you saw that if we give a complex exponential as an input to a system, the output
More informationLecture Schedule: Week Date Lecture Title
http://elec34.org Sampling and CONVOLUTION 24 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date Lecture Title 2-Mar Introduction 3-Mar
More informationDown-Sampling (4B) Young Won Lim 11/15/12
Down-Sampling (B) /5/ Copyright (c) 9,, Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version. or any later
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
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 informationReview of Fourier Transform
Review of Fourier Transform Fourier series works for periodic signals only. What s about aperiodic signals? This is very large & important class of signals Aperiodic signal can be considered as periodic
More informationContinuous Time Signal Analysis: the Fourier Transform. Lathi Chapter 4
Continuous Time Signal Analysis: the Fourier Transform Lathi Chapter 4 Topics Aperiodic signal representation by the Fourier integral (CTFT) Continuous-time Fourier transform Transforms of some useful
More informationChapter 2. Signals. Static and Dynamic Characteristics of Signals. Signals classified as
Chapter 2 Static and Dynamic Characteristics of Signals Signals Signals classified as. Analog continuous in time and takes on any magnitude in range of operations 2. Discrete Time measuring a continuous
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 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 informationModule 3. Convolution. Aim
Module Convolution Digital Signal Processing. Slide 4. Aim How to perform convolution in real-time systems efficiently? Is convolution in time domain equivalent to multiplication of the transformed sequence?
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
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. MASSACHUSETTS
More informationSquare Root Raised Cosine Filter
Wireless Information Transmission System Lab. Square Root Raised Cosine Filter Institute of Communications Engineering National Sun Yat-sen University Introduction We consider the problem of signal design
More informationContinuous and Discrete Time Signals and Systems
Continuous and Discrete Time Signals and Systems Mrinal Mandal University of Alberta, Edmonton, Canada and Amir Asif York University, Toronto, Canada CAMBRIDGE UNIVERSITY PRESS Contents Preface Parti Introduction
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationThe Discrete Fourier Transform. Signal Processing PSYCH 711/712 Lecture 3
The Discrete Fourier Transform Signal Processing PSYCH 711/712 Lecture 3 DFT Properties symmetry linearity shifting scaling Symmetry x(n) -1.0-0.5 0.0 0.5 1.0 X(m) -10-5 0 5 10 0 5 10 15 0 5 10 15 n m
More informationEA2.3 - Electronics 2 1
In the previous lecture, I talked about the idea of complex frequency s, where s = σ + jω. Using such concept of complex frequency allows us to analyse signals and systems with better generality. In this
More information