Fundamentals of the Discrete Fourier Transform
|
|
- Virgil Parrish
- 6 years ago
- Views:
Transcription
1 Seminar presentation at the Politecnico di Milano, Como, November 12, 2012 Fundamentals of the Discrete Fourier Transform Michael G. Sideris Department of Geomatics Engineering University of Calgary
2 Contents (1) From Fourier Series to the Continuous Fourier Transform Sinusoids Fourier Series The Continuous Fourier Transform (CFT) Important elementary functions and their CFTs The impulse function The CFT of sine and cosine The sampling function The rectangle and sinc functions Properties of the CFT Convolution and Correlation Properties of convolution
3 Contents (2) From the CFT to the DFT (Discrete Fourier Transform) Aliasing Leakage Periodicity The Discrete Fourier Transform Circular convolution Scale change Phase shifting Correlation, Covariance and Power Spectra Density functions The Fast Fourier Transform (FFT) The two-dimensional CFT and DFT
4 From Fourier Series to the Continuous Fourier Transform
5 Sinusoids s(t) = A 0 cos(ω 0 t +φ 0 ) : Sinusoid of frequency ω 0 A 0... amplitude ω 0... cyclic frequency t... time (or distance) φ 0... phase angle ω 0 = 2π T 0 = 2π f 0 T 0... period f 0...(linear) frequency Expansion s(t) = acosω 0 t + bsinω 0 t where a = A 0 cosφ 0, b = A 0 sinφ 0 A 0 = (a 2 + b 2 ) 1/2 φ 0 = arctan( b a ) Sinusoids in complex form s c (t) = acosω 0 t ± iasinω 0 t = ae ±iω 0 t Real sinusoids in complex form s(t) =A 0 cos(ω 0 t + φ 0 ) = A 0 e i(ω 0 t +φ 0 ) + e i(ω 0 t +φ 0 ) 2 = A 0 2 eiφ 0 e iω 0t + A 0 2 e iφ 0 e iω 0t
6 Fourier Series If g(t) = g(t +T ); T g(t)dt = g(t)dt, then 0 t 0 +T g(t) = (a n cos 2πn T t + b n sin 2πn T t ) n=0 a n = 2 T t 0 +T t 0 g(t)cosnt dt ; b n = 2 T t 0 t 0 +T t 0 g(t)sin nt dt Provided that: g(t) has a finite numbers of maxima and minima in a period and a finite number of finite discontinuities (Dirichlet s conditions) Complex form : g(t) = 1 T T /2 G n e iωnt, ω n = 2πn n= T G n = g(t)e iωnt dt, G n = 1 2 (a n ib n ), n = 0,±1, ± 2,... T /2 Call Δω = 2π T ) ω n = nδω + * 1 T = Δω, + 2π g(t) = G n n= 2π eiωnt Δω
7 The Continuous Fourier Transform (1) g(t) = 1 2π G(ω) = G(ω)e iωt dω Inverse CFT g(t) e iωt dt Direct CFT Since ω = 2πf g(t) = G( f )e i2π ft df = F -1 G( f ) { } G( f ) = g(t) e i2π ft dt = F g(t) { } G( f ) = G R ( f )+ ig I ( f ) = G( f ) e iθ ( f ) G( f ) is complex Amplitude: G( f ) = " # G 2 R ( f )+ G 2 I ( f ) $ % 1/2 Phase angle: θ( f ) = Arg{ G( f )} = arctan G ( f ) I G R ( f )
8 The Continuous Fourier Transform (2) Conditions for Existence : The integral of g(t) from to + exists (it is < ) g(t) has only finite discontinuities If g(t) is periodic or impulse, G( f ) does not exist
9 Important elementary functions and their CFTs
10 The Impulse Function (1) Definition : δ(t t 0 ) = 0, t t 0 ' ) ( δ(t t 0 )dt =1 ) * Definition as a distribution δ(t) = lim f (t,a) δ( t t ) φ( t )dt = a 0 0 φ ( t 0 )
11 The Impulse Function (2) Definition as a generlized limit: If lim f n (t)φ(t)dt = φ(0) δ(t) = n Properties: δ(t 0 )h(t) = h(t 0 )δ(t 0 ) δ(at) = a 1 δ(t) F{ Kδ(t) } = K sin at e.g. If we define it as δ(t) = lim, a πt then cos(2π ft) df = e i2π ft df = δ(t) lim f n (t) n Used, as a distribution, for the (otherwise nonexistent) CFT of periodic functions
12 The CTF of cosine and sine Acos(2πf o t) A 2 δ( f f o ) + A 2 δ( f + f o ) Asin(2πf o t) i A 2 δ( f + f o ) i A 2 δ( f f o )
13 The Sampling Function III(t) = δ(t nt ) F III(t) n= { } = 1 T n= δ( f n T ) III(t) f (t) = n= f (nt )δ(t nt ) Digitization Note that denser sampling in one domain transforms to sparser sampling in the other
14 The Rectangle and sinc Functions h(t) = " A, t = T 0 /2 $ # A /2, t = ±T 0 /2 $ % 0, t > T 0 /2 H( f ) = 2AT o sinc(2t o f ), sinc( f ) = sin(πf ) πf
15 Properties of the CFT
16 Properties of the CFT (1) Linearity ah(t) + bg(t) ah( f ) + bg( f ) Symmetry H(t ) h( f ) Time scaling h(at) 1 a H( f a ) Time shifting Differentiation Integration h(t t o ) H( f )e i 2πft o n h(t) t n (i 2πf ) n H( f ) t h(x )dx 1 i 2πf H( f ) H(0)δ( f )
17 Properties of the CFT (2) DC-value h(t)dt = H(0) Even function h E (t) H E ( f ) = R E ( f ) Odd function h O (t ) H O ( f ) = ii O ( f ) Real function h(t) = h R (t) H( f ) = R E ( f ) + ii O ( f ) Imaginary function h(t) = ih I (t) H( f ) = R O ( f ) + ii E ( f )
18 Convolution and Correlation
19 Convolution and Correlation (1) x(t) = g(t')h(t t')dt'= g(t) h(t) = h(t) g(t) = h(t')g(t t')dt' y(t) = g(t')h(t + t')dt'= g(t) h(t) h(t) g(t) Convolution theorem X( f ) = F{ g(t) h(t) } = F{ g(t) }F{ h(t) } = G( f )H ( f ) Convolution is the operation of filtering one of the functions by the other Correlation theorem Y ( f ) = F{ g(t) h(t) } = G( f )H * ( f )
20 Pictorial Representation of Convolution and Correlation
21 Convolution and Correlation (3) Properties: a) If either g(t) or h(t) is even, then g(t) h(t) = g(t) h(t) b) g(t) h(t) = h(t) g(t) convolution is commutative c) [g(t) h(t)] q(t) = h(t) [g(t) q(t)] convolution is associative d) k(t) [ag(t) bq(t)] = a[k(t) g(t)]+ b[k(t) q(t)] convolution is linear e) δ(t +τ ) h(t) = h(t +τ ), δ(t) h(t) = h(t) δ(t) is the identity operator in convolution f) x'(t) = (g(t) h(t))' = g'(t) h(t) = g(t) h'(t) the derivative of the g) F{ h(t)g(t) } = F{ h(t) } F{ g(t) } = H ( f ) G( f ) convolution
22 Convolution and Correlation (4) h) If h(t) and g(t) are time limited functions, i.e., non zero in the domain T 0 t T 0, then x(t) = h(x) g(t) is time limitedwith twice the support of h(t) or g(t), i.e., non zero in the domain 2T 0 t 2T 0 i) Parseval's thereom: h 2 (t)e 2πσ t dt = H ( f )H (σ f )df with σ = 0 and for h(t) real: h 2 (t)dt = H ( f ) 2 df
23 From the CFT to the DFT (Discrete Fourier Transform)
24 From the Continuous to the Discrete Fourier Transform aliasing leakage periodicity
25 The Discrete Fourier Transform
26 The Discrete Fourier Transform (1) H(mΔf ) = h(kδt) = N 1 h(kδt)e i2πkδtmδf Δt = h(kδt)e i2πkm / N Δt N 1 k= 0 N 1 k= 0 H(mΔf )e i2πkδtmδf Δf = H(mΔf )e i2πkm / N Δf m= 0 N 1 m= 0 h(kδt) H (mδf ) or h(t k ) H ( f m ) or h(k) H (m) T o = 1 Δf = NΔt, F o = 1 Δt = NΔf f N = F o 2 = 1 2Δt Nyquist frequency
27 The Discrete Fourier Transform (2) Illustration of Nyquist Frequency Highest frequencies existing in the original signal cannot be reproduced from the samples In this example, a sinusoid of lower frequency is recovered (the higher frequencies cannot be recovered; they are aliased due to the under-sampling)
28 Linear Convolution Discrete Convolution (1) For g(k) and h(k) defined in 0 k < N: x(k) = g(k) h(k) = g(l)h(k l)δt = g(l)h(k l)δt, 0 k 2N 1 l= M-point Circular Convolution k l=0 x c (k) = g(k) h(k) = M 1 l=0 g(l)h(k l) M Δt DFT && N F{g(k)}F{h(k)} Circularly shifted within M positions x c (k) = M 1 g(l)h(k l) M Δt = g(l)h(k l)δt + g(l)h(n + k l)δt l=0 k l=0 M 1 l=k+1 Linear convolution plus Aliasing
29 Discrete Convolution (2) Linear Convolution via FFT using Circular Convolution M 1 x c (k) = x(k) iff g(l)h(m + k l)δt = 0 l=k+1 which happens when M 2N 1 Zero-padding Procedure To compute the linear convolution of two functions given in 0 k N-1: 1. Append at least N-1 zeros to each function so that they are now 2N-1 samples long (preferably, 2 p 2N-1) each 2. Use FFT to compute their spectra 3. Use FFT to compute the Inverse DFT of the product of the two spectra. This, i.e., the resulting values for 0 k 2N-1, is the correct linear convolution of the two functions
30 Circular Convolution as Linear Convolution Plus Aliasing
31 The DFT in Computers Subroutines usually assume Δt = 1 and also ignore T 0 This requires rescaling as follows: H( f m ) = T o H c (m) = NΔtH c (m) x(t k ) = g(t k ) h(t k ) = T o x c (t k ) = T o F c 1 {G c (m)h c (m)} It also yields: H c (0) = h Subroutines also assume the origin at the left of the record This requires changing the phase of the spectrum by e i2π mδft o/2 = e iπ m = cos(mπ ) = ( 1) m h(t k T o / 2) ( 1) m H( f m ) End point of a period must be omitted (assumed due to periodicity)
32 Correlation, Covariance and Power Spectral Density functions
33 Definitions: CR, CV and PSD Functions R gh (t k ) = E{g(t l )h(t k + t l )} = lim 1 g(t N N l )h(t k + t l ) N 1 l=0 = lim To 1 T o g(t k ) h(t k ) C gh (t k ) = E{{g(t l ) g][h(t k + t l ) h ]} = lim 1 [g(t N N l ) g][h(t k + t l ) h ] 1 = lim g(t k ) h(t k ) gh = R gh (t k ) gh To T o 1 P gh ( f m ) = F{R gh (t k )} = lim G( f m )H ( f m ) To T o Computation by FFT: ˆP gh ( f m ) = 1 νt o ν λ=1 ˆR gh (t k ) = F 1 { ˆP gh ( f m )} G λ ( f m )H λ ( f m ) Ĉ gh (t k ) = F 1 { ˆP gh ( f m ) ghδ( f m )} N 1 l=0
34 The Fast Fourier Transform
35 The FFT - Flow Graph of Operations for N=4 H = Wh FFT achieves its speed by factorizing W and exploring its symmetries in a way that the required number of multiplications and additions is drastically reduced
36 Computational efficiency: FFT vs DFT 1D case DFT: O(N 2 ) FFT: O(NlogN) Bergland (1969)
37 The 2D Continuous and Discrete Fourier Transform
38 The Two-dimensional CFT H (u,v) = h(x, y) = h(x, y)e i2π (ux+vy) dx dy H (u,v)e i2π (ux+vy) du dv The Fourier Transform is a separable transform Very important for practical applications H (u,v) = h(x, y)e i2π (ux+vy) dx dy = h(x, y)e i2πux e i2πvy dx dy = ( h(x, y)e i2πux dx)e i2πvy dy = F x {h(x, y)}e i2πvy dy = F y {F x {h(x, y)}}
39 The Two-dimensional DFT M 1 N 1 H(u m,v n ) = h(x k, y l ) e i2π (mk / M +nl / N ) ΔxΔy k= 0 1= 0 M 1 N 1 h(x k, y l ) = H(u m,v n ) e i2π (mk / M +nl / N ) ΔuΔv k= 0 1= 0 SPACE DOMAIN FREQUENCY DOMAIN Δu = 1 T y = 1 MΔx, Δv = 1 T y = 1 NΔy Δx = 1 F u = 1 MΔu = 1 2u N, Δy = 1 F v = 1 NΔv = 1 2v N
40 Concluding Remarks
41 Concluding Remarks Fourier Transforms are important in geodetic applications for, e.g., Spectral analysis of data Can decide on the spectral content, appropriate sampling interval, and/or best computational method Digital signal/image processing Filtering of data Filter design Noise minimization Efficient evaluation of convolutions and of covariance functions Fourier Transforms are subject to errors due to Aliasing (sampling) Leakage (truncation of domain) Periodicity (digitization) Data noise which should be understood and minimized in practical applications
Computational 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 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 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 information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
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 informationA523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011
A523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011 Lecture 6 PDFs for Lecture 1-5 are on the web page Problem set 2 is on the web page Article on web page A Guided
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 informationFourier series. Complex Fourier series. Positive and negative frequencies. Fourier series demonstration. Notes. Notes. Notes.
Fourier series Fourier series of a periodic function f (t) with period T and corresponding angular frequency ω /T : f (t) a 0 + (a n cos(nωt) + b n sin(nωt)), n1 Fourier series is a linear sum of cosine
More informationSolutions to Problems in Chapter 4
Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave
More information2 Fourier Transforms and Sampling
2 Fourier ransforms and Sampling 2.1 he Fourier ransform he Fourier ransform is an integral operator that transforms a continuous function into a continuous function H(ω) =F t ω [h(t)] := h(t)e iωt dt
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 informationChapter 4 The Fourier Series and Fourier Transform
Chapter 4 The Fourier Series and Fourier Transform Representation of Signals in Terms of Frequency Components Consider the CT signal defined by N xt () = Acos( ω t+ θ ), t k = 1 k k k The frequencies `present
More information2 Frequency-Domain Analysis
2 requency-domain Analysis Electrical engineers live in the two worlds, so to speak, of time and frequency. requency-domain analysis is an extremely valuable tool to the communications engineer, more so
More informationBME 50500: Image and Signal Processing in Biomedicine. Lecture 2: Discrete Fourier Transform CCNY
1 Lucas Parra, CCNY BME 50500: Image and Signal Processing in Biomedicine Lecture 2: Discrete Fourier Transform Lucas C. Parra Biomedical Engineering Department CCNY http://bme.ccny.cuny.edu/faculty/parra/teaching/signal-and-image/
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 informationVer 3808 E1.10 Fourier Series and Transforms (2014) E1.10 Fourier Series and Transforms. Problem Sheet 1 (Lecture 1)
Ver 88 E. Fourier Series and Transforms 4 Key: [A] easy... [E]hard Questions from RBH textbook: 4., 4.8. E. Fourier Series and Transforms Problem Sheet Lecture. [B] Using the geometric progression formula,
More informationFigure 3.1 Effect on frequency spectrum of increasing period T 0. Consider the amplitude spectrum of a periodic waveform as shown in Figure 3.2.
3. Fourier ransorm From Fourier Series to Fourier ransorm [, 2] In communication systems, we oten deal with non-periodic signals. An extension o the time-requency relationship to a non-periodic signal
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 information2.1 Basic Concepts Basic operations on signals Classication of signals
Haberle³me Sistemlerine Giri³ (ELE 361) 9 Eylül 2017 TOBB Ekonomi ve Teknoloji Üniversitesi, Güz 2017-18 Dr. A. Melda Yüksel Turgut & Tolga Girici Lecture Notes Chapter 2 Signals and Linear Systems 2.1
More informationX(t)e 2πi nt t dt + 1 T
HOMEWORK 31 I) Use the Fourier-Euler formulae to show that, if X(t) is T -periodic function which admits a Fourier series decomposition X(t) = n= c n exp (πi n ) T t, then (1) if X(t) is even c n are all
More informationESS Dirac Comb and Flavors of Fourier Transforms
6. Dirac Comb and Flavors of Fourier ransforms Consider a periodic function that comprises pulses of amplitude A and duration τ spaced a time apart. We can define it over one period as y(t) = A, τ / 2
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 informationOSE801 Engineering System Identification. Lecture 05: Fourier Analysis
OSE81 Engineering System Identification Lecture 5: Fourier Analysis What we will study in this lecture: A short introduction of Fourier analysis Sampling the data Applications Example 1 Fourier Analysis
More informationPART 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.
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
More informationSystem Identification & Parameter Estimation
System Identification & Parameter Estimation Wb3: SIPE lecture Correlation functions in time & frequency domain Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE // Delft University
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 informationChapter 4 The Fourier Series and Fourier Transform
Chapter 4 The Fourier Series and Fourier Transform Fourier Series Representation of Periodic Signals Let x(t) be a CT periodic signal with period T, i.e., xt ( + T) = xt ( ), t R Example: the rectangular
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 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 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 informationa k cos kω 0 t + b k sin kω 0 t (1) k=1
MOAC worksheet Fourier series, Fourier transform, & Sampling Working through the following exercises you will glean a quick overview/review of a few essential ideas that you will need in the moac course.
More informationCHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L
CHAPTER 4 FOURIER SERIES 1 S A B A R I N A I S M A I L Outline Introduction of the Fourier series. The properties of the Fourier series. Symmetry consideration Application of the Fourier series to circuit
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 informationTHE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3
THE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3 Any periodic function f(t) can be written as a Fourier Series a 0 2 + a n cos( nωt) + b n sin n
More information3. Frequency-Domain Analysis of Continuous- Time Signals and Systems
3. Frequency-Domain Analysis of Continuous- ime Signals and Systems 3.. Definition of Continuous-ime Fourier Series (3.3-3.4) 3.2. Properties of Continuous-ime Fourier Series (3.5) 3.3. Definition of Continuous-ime
More informatione iωt dt and explained why δ(ω) = 0 for ω 0 but δ(0) =. A crucial property of the delta function, however, is that
Phys 531 Fourier Transforms In this handout, I will go through the derivations of some of the results I gave in class (Lecture 14, 1/11). I won t reintroduce the concepts though, so you ll want to refer
More informationSignal and systems. Linear Systems. Luigi Palopoli. Signal and systems p. 1/5
Signal and systems p. 1/5 Signal and systems Linear Systems Luigi Palopoli palopoli@dit.unitn.it Wrap-Up Signal and systems p. 2/5 Signal and systems p. 3/5 Fourier Series We have see that is a signal
More informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
More informationFunctions of a Complex Variable (S1) Lecture 11. VII. Integral Transforms. Integral transforms from application of complex calculus
Functions of a Complex Variable (S1) Lecture 11 VII. Integral Transforms An introduction to Fourier and Laplace transformations Integral transforms from application of complex calculus Properties of Fourier
More informationHow many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation?
How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation? (A) 0 (B) 1 (C) 2 (D) more than 2 (E) it depends or don t know How many of
More informationLINEAR SYSTEMS. J. Elder PSYC 6256 Principles of Neural Coding
LINEAR SYSTEMS Linear Systems 2 Neural coding and cognitive neuroscience in general concerns input-output relationships. Inputs Light intensity Pre-synaptic action potentials Number of items in display
More informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
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 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 informationSEISMIC 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 informationLecture 7 ELE 301: Signals and Systems
Lecture 7 ELE 30: Signals and Systems Prof. Paul Cuff Princeton University Fall 20-2 Cuff (Lecture 7) ELE 30: Signals and Systems Fall 20-2 / 22 Introduction to Fourier Transforms Fourier transform as
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 informationFourier Series. Fourier Transform
Math Methods I Lia Vas Fourier Series. Fourier ransform Fourier Series. Recall that a function differentiable any number of times at x = a can be represented as a power series n= a n (x a) n where the
More informationMIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain
10/31/05 wk9-a-1 The spatial frequency domain Recall: plane wave propagation x path delay increases linearly with x λ z=0 θ E 0 x exp i2π sinθ + λ z i2π cosθ λ z plane of observation 10/31/05 wk9-a-2 Spatial
More information23.6. The Complex Form. Introduction. Prerequisites. Learning Outcomes
he Complex Form 3.6 Introduction In this Section we show how a Fourier series can be expressed more concisely if we introduce the complex number i where i =. By utilising the Euler relation: e iθ cos θ
More informationEE303: Communication Systems
EE303: Communication Systems Professor A. Manikas Chair of Communications and Array Processing Imperial College London Introductory Concepts Prof. A. Manikas (Imperial College) EE303: Introductory Concepts
More informationFourier Analysis Linear transformations and lters. 3. Fourier Analysis. Alex Sheremet. April 11, 2007
Stochastic processes review 3. Data Analysis Techniques in Oceanography OCP668 April, 27 Stochastic processes review Denition Fixed ζ = ζ : Function X (t) = X (t, ζ). Fixed t = t: Random Variable X (ζ)
More informationContinuous-time Fourier Methods
ELEC 321-001 SIGNALS and SYSTEMS Continuous-time Fourier Methods Chapter 6 1 Representing a Signal The convolution method for finding the response of a system to an excitation takes advantage of the linearity
More informationUnstable Oscillations!
Unstable Oscillations X( t ) = [ A 0 + A( t ) ] sin( ω t + Φ 0 + Φ( t ) ) Amplitude modulation: A( t ) Phase modulation: Φ( t ) S(ω) S(ω) Special case: C(ω) Unstable oscillation has a broader periodogram
More information8: Correlation. E1.10 Fourier Series and Transforms ( ) Fourier Transform - Correlation: 8 1 / 11. 8: Correlation
E. Fourier Series and Transforms (25-5585) - Correlation: 8 / The cross-correlation between two signals u(t) and v(t) is w(t) = u(t) v(t) u (τ)v(τ +t)dτ = u (τ t)v(τ)dτ [sub: τ τ t] The complex conjugate,
More information! Circular Convolution. " Linear convolution with circular convolution. ! Discrete Fourier Transform. " Linear convolution through circular
Previously ESE 531: Digital Signal Processing Lec 22: April 18, 2017 Fast Fourier Transform (con t)! Circular Convolution " Linear convolution with circular convolution! Discrete Fourier Transform " Linear
More informationDiscrete Fourier Transform
Discrete Fourier Transform Virtually all practical signals have finite length (e.g., sensor data, audio records, digital images, stock values, etc). Rather than considering such signals to be zero-padded
More 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 informationNon-parametric identification
Non-parametric Non-parametric Transient Step-response using Spectral Transient Correlation Frequency function estimate Spectral System Identification, SSY230 Non-parametric 1 Non-parametric Transient Step-response
More informationProblem Sheet 1 Examples of Random Processes
RANDOM'PROCESSES'AND'TIME'SERIES'ANALYSIS.'PART'II:'RANDOM'PROCESSES' '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''Problem'Sheets' Problem Sheet 1 Examples of Random Processes 1. Give
More informationDISCRETE FOURIER TRANSFORM
DD2423 Image Processing and Computer Vision DISCRETE FOURIER TRANSFORM Mårten Björkman Computer Vision and Active Perception School of Computer Science and Communication November 1, 2012 1 Terminology:
More informationIntroduction to Fourier Transforms. Lecture 7 ELE 301: Signals and Systems. Fourier Series. Rect Example
Introduction to Fourier ransforms Lecture 7 ELE 3: Signals and Systems Fourier transform as a limit of the Fourier series Inverse Fourier transform: he Fourier integral theorem Prof. Paul Cuff Princeton
More informationFourier Analysis and Power Spectral Density
Chapter 4 Fourier Analysis and Power Spectral Density 4. Fourier Series and ransforms Recall Fourier series for periodic functions for x(t + ) = x(t), where x(t) = 2 a + a = 2 a n = 2 b n = 2 n= a n cos
More informationEE401: Advanced Communication Theory
EE401: Advanced Communication Theory Professor A. Manikas Chair of Communications and Array Processing Imperial College London Introductory Concepts Prof. A. Manikas (Imperial College) EE.401: Introductory
More informationSeries FOURIER SERIES. Graham S McDonald. A self-contained Tutorial Module for learning the technique of Fourier series analysis
Series FOURIER SERIES Graham S McDonald A self-contained Tutorial Module for learning the technique of Fourier series analysis Table of contents Begin Tutorial c 24 g.s.mcdonald@salford.ac.uk 1. Theory
More informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
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 informationChapter 6: Applications of Fourier Representation Houshou Chen
Chapter 6: Applications of Fourier Representation Houshou Chen Dept. of Electrical Engineering, National Chung Hsing University E-mail: houshou@ee.nchu.edu.tw H.S. Chen Chapter6: Applications of Fourier
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 information1.1 SPECIAL FUNCTIONS USED IN SIGNAL PROCESSING. δ(t) = for t = 0, = 0 for t 0. δ(t)dt = 1. (1.1)
SIGNAL THEORY AND ANALYSIS A signal, in general, refers to an electrical waveform whose amplitude varies with time. Signals can be fully described in either the time or frequency domain. This chapter discusses
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 informationSummary of Fourier Transform Properties
Summary of Fourier ransform Properties Frank R. Kschischang he Edward S. Rogers Sr. Department of Electrical and Computer Engineering University of oronto January 7, 207 Definition and Some echnicalities
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 informationDiscrete Fourier Transform
Discrete Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013 Mårten Björkman
More informationLecture 11: Spectral Analysis
Lecture 11: Spectral Analysis Methods For Estimating The Spectrum Walid Sharabati Purdue University Latest Update October 27, 2016 Professor Sharabati (Purdue University) Time Series Analysis October 27,
More informationToday s lecture. The Fourier transform. Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm
Today s lecture The Fourier transform What is it? What is it useful for? What are its properties? Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm Jean Baptiste Joseph Fourier
More informatione iωt dt and explained why δ(ω) = 0 for ω 0 but δ(0) =. A crucial property of the delta function, however, is that
Phys 53 Fourier Transforms In this handout, I will go through the derivations of some of the results I gave in class (Lecture 4, /). I won t reintroduce the concepts though, so if you haven t seen the
More informationLOPE3202: Communication Systems 10/18/2017 2
By Lecturer Ahmed Wael Academic Year 2017-2018 LOPE3202: Communication Systems 10/18/2017 We need tools to build any communication system. Mathematics is our premium tool to do work with signals and systems.
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 informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring
A653 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 15 http://www.astro.cornell.edu/~cordes/a653 Lecture 3 Power spectrum issues Frequentist approach Bayesian approach (some
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 informationWave Phenomena Physics 15c. Lecture 11 Dispersion
Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed
More informationImages have structure at various scales
Images have structure at various scales Frequency Frequency of a signal is how fast it changes Reflects scale of structure A combination of frequencies 0.1 X + 0.3 X + 0.5 X = Fourier transform Can we
More informationData Processing and Analysis
Data Processing and Analysis Rick Aster and Brian Borchers October 10, 2008 Sampled Time Series Commonly, scientific data consist of numerical series or matrices, i.e., sets of spatially or temporally
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 informationThe Fourier Transform
The Fourier Transorm Fourier Series Fourier Transorm The Basic Theorems and Applications Sampling Bracewell, R. The Fourier Transorm and Its Applications, 3rd ed. New York: McGraw-Hill, 2. Eric W. Weisstein.
More informationDigital Signal Processing
Digital Signal Processing Introduction Moslem Amiri, Václav Přenosil Embedded Systems Laboratory Faculty of Informatics, Masaryk University Brno, Czech Republic amiri@mail.muni.cz prenosil@fi.muni.cz February
More informationFourier analysis of discrete-time signals. (Lathi Chapt. 10 and these slides)
Fourier analysis of discrete-time signals (Lathi Chapt. 10 and these slides) Towards the discrete-time Fourier transform How we will get there? Periodic discrete-time signal representation by Discrete-time
More information( ) f (k) = FT (R(x)) = R(k)
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) = R(x)
More informationReview of Concepts from Fourier & Filtering Theory. Fourier theory for finite sequences. convolution/filtering of infinite sequences filter cascades
Review of Concepts from Fourier & Filtering Theory precise definition of DWT requires a few basic concepts from Fourier analysis and theory of linear filters will start with discussion/review of: basic
More informationFourier Series Example
Fourier Series Example Let us compute the Fourier series for the function on the interval [ π,π]. f(x) = x f is an odd function, so the a n are zero, and thus the Fourier series will be of the form f(x)
More information06EC44-Signals and System Chapter Fourier Representation for four Signal Classes
Chapter 5.1 Fourier Representation for four Signal Classes 5.1.1Mathematical Development of Fourier Transform If the period is stretched without limit, the periodic signal no longer remains periodic but
More informationVisual features: From Fourier to Gabor
Visual features: From Fourier to Gabor Deep Learning Summer School 2015, Montreal Hubel and Wiesel, 1959 from: Natural Image Statistics (Hyvarinen, Hurri, Hoyer; 2009) Alexnet ICA from: Natural Image Statistics
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 information1 Signals and systems
978--52-5688-4 - Introduction to Orthogonal Transforms: With Applications in Data Processing and Analysis Signals and systems In the first two chapters we will consider some basic concepts and ideas as
More informationDiscrete Systems & Z-Transforms. Week Date Lecture Title. 9-Mar Signals as Vectors & Systems as Maps 10-Mar [Signals] 3
http:elec34.org Discrete Systems & Z-Transforms 4 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: eek Date Lecture Title -Mar Introduction
More informationECE 301 Fall 2011 Division 1 Homework 10 Solutions. { 1, for 0.5 t 0.5 x(t) = 0, for 0.5 < t 1
ECE 3 Fall Division Homework Solutions Problem. Reconstruction of a continuous-time signal from its samples. Let x be a periodic continuous-time signal with period, such that {, for.5 t.5 x(t) =, for.5
More informationEE 438 Essential Definitions and Relations
May 2004 EE 438 Essential Definitions and Relations CT Metrics. Energy E x = x(t) 2 dt 2. Power P x = lim T 2T T / 2 T / 2 x(t) 2 dt 3. root mean squared value x rms = P x 4. Area A x = x(t) dt 5. Average
More informationThe formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) R(x)
More informationECE Digital Image Processing and Introduction to Computer Vision. Outline
ECE592-064 Digital mage Processing and ntroduction to Computer Vision Depart. of ECE, NC State University nstructor: Tianfu (Matt) Wu Spring 2017 1. Recap Outline 2. Thinking in the frequency domain Convolution
More information