Sequential Fast Fourier Transform
|
|
- Florence O’Brien’
- 5 years ago
- Views:
Transcription
1 Sequential Fast Fourier Transform Departement Computerwetenschappen Room (Celestijnenlaan 200A) Includes material from slides by Prof. dr. Rob H. Bisseling
2 Applications of Fourier analysis Fourier analysis studies the decomposition of functions into their frequency components. Piano Concerto no. 9 by Mozart: enhance high frequencies. Chest picture by Computerised Tomography: reconstruct your interior without slicing you up. Star picture by pre-repair Hubble Space Telescope: remove blur. Sequential FFT 2 / 18
3 Fourier Series General idea: every signal is built out of harmonic parts. (For example, a linear combination of various sinusoids each with varying frequency and amplitude.)
4 Fourier Series Given a T -periodic function f : R C, with f piecewise continuous. Then the Fourier-coefficients {ˆf k } k Z are defined as: ˆf k = 1 T The Fourier series are given as: S n = T 0 f (t)e 2πik t T dt. n k= n ˆf k e 2πik t T. In accordance with the general idea, would like that lim S n = f. n
5 Fourier Inverse theorems f : R C is T -periodic and lim n S n = f, if
6 Fourier Inverse theorems f : R C is T -periodic and lim n S n = f, if f is piecewise differentiable, ( x R, 1 2 (lim y x f (y) + lim y x f (y)) = lim n S n )
7 Fourier Inverse theorems f : R C is T -periodic and lim n S n = f, if f is piecewise differentiable, ( x R, 1 2 (lim y x f (y) + lim y x f (y)) = lim n S n ) f is two-times differentiable (f C 2 ), ( x R, f (x) = lim n S n ; convergence in norm)
8 Fourier Inverse theorems f : R C is T -periodic and lim n S n = f, if f is piecewise differentiable, ( x R, 1 2 (lim y x f (y) + lim y x f (y)) = lim n S n ) f is two-times differentiable (f C 2 ), ( x R, f (x) = lim n S n ; convergence in norm) f is in L 1 T (i.e., T 0 f (t) dt converges) as well as in C 1, or (lim n f S n = 0)
9 Fourier Inverse theorems f : R C is T -periodic and lim n S n = f, if f is piecewise differentiable, ( x R, 1 2 (lim y x f (y) + lim y x f (y)) = lim n S n ) f is two-times differentiable (f C 2 ), ( x R, f (x) = lim n S n ; convergence in norm) f is in L 1 T (i.e., T 0 f (t) dt converges) as well as in C 1, or (lim n f S n = 0) f is in L 2 T (i.e., T 0 f (t) 2 dt converges). (lim n f S n 2 = 0) (Being piecewise continuous is not enough.)
10 Fourier Inverse theorems for non-periodic input If f lives on a bounded domain (f : [0, T ] C), as for example with a time-finite input signal, we can copy the signal to form an aritificial T -periodic function g: g : R C, with g(t) = f (t T /2) (for t [ T /2, T /2], and repeated elsewhere). The inverse is correct if g meets any of the previous conditions.
11 Fourier Inverse theorems for non-periodic input Let f, ĝ : R C (unbounded in time) bounded in 2-norm (f, ĝ L 2 ), (Ff ) s : k s (F g) s : t then s s f (t)e 2πitk dt, and s ĝ(t)e 2πitk dt. ˆf L 2 such that lim s ˆf (Ff ) s 2 = 0, and g L 2 such that lim s g (F ˆf ) s 2 = 0, and F Fg = g. (Depending on f, (F Ff ) s may converge uniformly or point-wise as s ).
12 It s a discrete world One second of audio on a compact disc contains 44,100 function values f (t j ) in regularly spaced sample points t j = jt n, 0 j < n. Sequential FFT 5 / 18
13 Approximation of Fourier coefficients Trapezoidal rule on interval [t j, t j+1 ] = tj+1 On the whole interval [0, T ]: c k = 1 T f (t) dt f (t j) + f (t j+1 ) t j 2 T 0 1 T T n [ ] jt n, (j+1)t n : T n. f (t)e 2πikt/T dt f (0) n f (t j )e 2πikt j /T + f (T ) 2 j=1 = 1 n 1 f (t j )e 2πijk/n (since f (0) = f (T ) = f (t 0 )). n Sequential FFT 6 / 18
14 Discrete Fourier Transform We define the discrete Fourier transform (DFT): Definition. DFT : C n C n with y = DFT (x) such that n 1 y k = x j e 2πijk/n, for each k [0, n 1] (x, y C n ). (Note this skips the factor 1/n from the repeated trapezoidal rule.)
15 Inverse discrete Fourier transform Definition. invdft : C n C n with y = invdft (x) such that x j = 1 n 1 y k e 2πijk/n, for each k [0, n 1] (x, y C n ). n k=0 (Note the factor 1/n reappears here.) We have that for all x C n, invdft (DFT (x)) = x. But how close is DFT (x) to the true Ff?
16 Matrix vector multiplication Define the n n Fourier matrix F n by Hence F n x = DFT (x): (F n ) jk = ω n jk, for 0 j, k < n. n 1 n 1 (F n x) j = (F n ) jk x k = x k ω jk n = (DFT (x)) j. k=0 k=0 Because ω 4 = e 2πi/4 = e πi/2 = i: ω 0 4 ω 0 4 ω 0 4 ω 0 4 F 4 = ω 0 4 ω 1 4 ω 2 4 ω 3 4 ω 0 4 ω 2 4 ω 4 4 ω 6 4 = ω 0 4 ω 3 4 ω 6 4 ω i 1 i i 1 i. Sequential FFT 10 / 18
17 Roots of unity ω 5 ω 6 = i ω 7 ω 4 ω 0 = 1 ω 3 ω 2 ω 1 ω 8 = e 2πi/8 = e πi/4 = i. ω n n = e 2πin/n = e 2πi = 1. ω n n/2 = e 2πi(n/2)/n = e πi = 1. ω n 2 = e 4πi/n = e 2πi/(n/2) = ω n/2. Sequential FFT 9 / 18
18 Cost of straightforward DFT Complex multiplication (a + bi)(c + di) = (ac bd) + (ad + bc)i requires 1 real addition, 1 real subtraction, 4 real multiplications, hence a total of 6 flops. Complex addition (a + bi) + (c + di) = (a + c) + (b + d)i requires 2 real additions. To compute y k, we need n complex multiplications and n 1 complex additions, so 6n + 2(n 1) = 8n 2 flops. To compute the n components of y, we need 8n 2 2n flops. Sequential FFT 11 / 18
19 Towards FFT: divide... Reducing the DFT flop-cost: y k = n 1 x jω jk n, split the sum in odd and even parts: y k = n/2 1 n/2 1 x 2j wn 2jk + x 2j+1 ω (2j+1)k n. Since ωn 2jk = ω jk n/2 (roots of unity: the 2c th root out of n roots is the c th root out of n/2 roots): y k = n/2 1 n/2 1 x 2j w jk n/2 + ωk n x 2j+1 ω jk n/2. These are DFTs of half length each!
20 Towards FFT: divide... Reducing the DFT flop-cost: y k = n 1 x jω jk n, split the sum in odd and even parts: y k = n/2 1 n/2 1 x 2j wn 2jk + x 2j+1 ω (2j+1)k n. Since ωn 2jk = ω jk n/2 (roots of unity: the 2c th root out of n roots is the c th root out of n/2 roots): y 0 k<n/2 = DFT (x even )+Ω n DFT (x odd ), with Ω n = diag(ω 0 n,..., ω n/2 1 n ). These are DFTs of half length each!
21 Towards FFT: divide..., divide... Split up y in two parts of length n/2 as well: y top = DFT (x even ) + Ω n DFT (x odd ), with Ω n = diag(ω 0 n,..., ω n/2 1 n ). Substitute k = k n/2 for the remaining half of y: n/2 1 y0 k bottom <n = n/2 1 x 2j w j(n/2+k ) n/2 + ωn k Note 1 that ω jn/2 n/2 = (ωn/2 n/2 )j = 1 and k = k + n/2, x 2j+1 ω j(n/2+k ) n/2. 1, roots of unity: the n th root out of n roots equals 1.
22 Towards FFT: divide..., divide... Split up y in two parts of length n/2 as well: y top = DFT (x even ) + Ω n DFT (x odd ), with Ω n = diag(ω 0 n,..., ω n/2 1 n ). Substitute k = k n/2 for the remaining half of y: n/2 1 y0 k bottom <n = n/2 1 x 2j w j(n/2+k ) n/2 + ωn k x 2j+1 ω j(n/2+k ) n/2. Note 1 that ω jn/2 n/2 = (ωn/2 n/2 )j = 1 and k = k + n/2, hence: n/2 1 y0 k bottom <n = n/2 1 x 2j w jk n/2 + +n/2 ωk n x 2j+1 ω jk n/2 ; 1, roots of unity: the n th root out of n roots equals 1.
23 Towards FFT: divide..., divide... Split up y in two parts of length n/2 as well: y top = DFT (x even ) + Ω n DFT (x odd ), with Ω n = diag(ω 0 n,..., ω n/2 1 n ). Substitute k = k n/2 for the remaining half of y: n/2 1 y0 k bottom <n = n/2 1 x 2j w jk n/2 + +n/2 ωk n Note 2 that ω k +n/2 n = ω n/2 n ω k n = ω k n, x 2j+1 ω jk n/2. 2, roots of unity: the n/2 th root out of n roots always ends up halfway the circle, at 1.
24 Towards FFT: divide..., divide... Split up y in two parts of length n/2 as well: y top = DFT (x even ) + Ω n DFT (x odd ), with Ω n = diag(ω 0 n,..., ω n/2 1 n ). Substitute k = k n/2 for the remaining half of y: n/2 1 y0 k bottom <n = n/2 1 x 2j w jk n/2 + +n/2 ωk n Note 2 that ω k +n/2 n = ω n/2 n ω k n = ω k n, hence: n/2 1 y0 k bottom <n = n/2 1 x 2j w jk n/2 ωk n x 2j+1 ω jk n/2. x 2j+1 ω jk n/2. 2, roots of unity: the n/2 th root out of n roots always ends up halfway the circle, at 1.
25 Towards FFT: divide..., divide..., combine... Splitting up y in two parts thus yields: y top = DFT (x even ) + Ω n DFT (x odd ) y bottom = DFT (x even ) Ω n DFT (x odd ), with Ω n = diag(ω 0 n,..., ω n/2 1 n ).
26 Towards FFT: divide..., divide..., combine... Splitting up y in two parts thus yields: y top = DFT (x even ) + Ω n DFT (x odd ) y bottom = DFT (x even ) Ω n DFT (x odd ), with Ω n = diag(ωn, 0..., ωn n/2 1 ). Thus a DFT of size n is split in two half-size DFTs plus n/2 multiplications and n additions (Ω-twiddling). a DFT of size n/2 takes 8n 2 /4 n flops, two half-size DFTs take 4n 2 2n flops, n/2 complex multiplications take 3n flops, n complex additionals take 2n flops, so the above divide & combine strategy costs 4n 2 + 3n flops! (Compared to the regular cost of 8n 2 n/2.)
27 FFT: divide, combine, and conquer! Algortihm FFT (x); input: x C n with n a power of two, output: y C n with y = DFT (x): if n = 1 then return y = x else ( y t Let y = y b ) be a vector of size n (x e, x o ) = odd-even-sort(x) τ = Ω n FFT (x o ) y t = FFT (x e ) + τ y b = FFT (x o ) τ return y end
28 How fast is fast? Cost of the recursive FFT: τ = Ω n FFT (x o ): n/2 complex multiplications. y t = FFT (x e ) + τ: n/2 complex additions, y b = FFT (x o ) τ: n/2 complex additions; totals to 5n flops per call to the FFT.
29 Sequential FFT Recursive computation of DFT The problem is split repeatedly, until the problem size is / 18
30 How fast is fast? Cost of the recursive FFT: 5n flops per call to the FFT. Two recursive FFT calls if n > 1: T (n) = 2T (n/2) + 5n. This is linear: for n > 8, T (n) =
31 How fast is fast? Cost of the recursive FFT: 5n flops per call to the FFT. Two recursive FFT calls if n > 1: T (n) = 2T (n/2) + 5n. This is linear: for n > 8, T (n) = 2T (n/2) + 5n =
32 How fast is fast? Cost of the recursive FFT: 5n flops per call to the FFT. Two recursive FFT calls if n > 1: T (n) = 2T (n/2) + 5n. This is linear: for n > 8, T (n) = 2T (n/2) + 5n = 4T (n/4) + 2 5n =
33 How fast is fast? Cost of the recursive FFT: 5n flops per call to the FFT. Two recursive FFT calls if n > 1: T (n) = 2T (n/2) + 5n. This is linear: for n > 8, T (n) = 2T (n/2) + 5n = 4T (n/4) + 2 5n = 8T (n/8) =... ;
34 How fast is fast? Cost of the recursive FFT: 5n flops per call to the FFT. Two recursive FFT calls if n > 1: T (n) = 2T (n/2) + 5n. This is linear: for n > 8, T (n) = 2T (n/2) + 5n = 4T (n/4) + 2 5n = 8T (n/8) =... ; that is, each level of the recursion tree contains a total of 5n flops. The tree height is exactly 1 + log 2 n,
35 How fast is fast? Cost of the recursive FFT: 5n flops per call to the FFT. Two recursive FFT calls if n > 1: T (n) = 2T (n/2) + 5n. This is linear: for n > 8, T (n) = 2T (n/2) + 5n = 4T (n/4) + 2 5n = 8T (n/8) =... ; that is, each level of the recursion tree contains a total of 5n flops. The tree height is exactly 1 + log 2 n, but T (1) = 0 (no flops, just copying of a single complex value).
36 How fast is fast? Cost of the recursive FFT: 5n flops per call to the FFT. Two recursive FFT calls if n > 1: T (n) = 2T (n/2) + 5n. This is linear: for n > 8, T (n) = 2T (n/2) + 5n = 4T (n/4) + 2 5n = 8T (n/8) =... ; that is, each level of the recursion tree contains a total of 5n flops. The tree height is exactly 1 + log 2 n, but T (1) = 0 (no flops, just copying of a single complex value). Hence the total cost of the recursive FFT is 5n log 2 n.
37 Summary The fast Fourier transform (FFT) idea was discovered by Gauss (1805), rediscovered by Danielson and Lanczos (1942), and is commonly attributed to Cooley and Tukey (1965), who rediscovered it in the digital era. The FFT is the computational workhorse in many applications, from weather forecasting to signal and image processing. Without the FFT, modern medicine would be impossible. The cost of an FFT of length n is 5n log 2 n flops. We have derived a recursive FFT algorithm, i.e., an algorithm that calls itself with a smaller problem size. Sequential FFT 18 / 18
Sequential Fast Fourier Transform (PSC ) Sequential FFT
Sequential Fast Fourier Transform (PSC 3.1 3.2) 1 / 18 Applications of Fourier analysis Fourier analysis studies the decomposition of functions into their frequency components. Piano Concerto no. 9 by
More informationComputational Methods CMSC/AMSC/MAPL 460
Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Balaji Vasan Srinivasan Dept of Computer Science Several slides from Prof Healy s course at UMD Last time: Fourier analysis F(t) = A 0 /2 + A
More informationDivide and Conquer algorithms
Divide and Conquer algorithms Another general method for constructing algorithms is given by the Divide and Conquer strategy. We assume that we have a problem with input that can be split into parts in
More informationThe Fast Fourier Transform. Andreas Klappenecker
The Fast Fourier Transform Andreas Klappenecker Motivation There are few algorithms that had more impact on modern society than the fast Fourier transform and its relatives. The applications of the fast
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More informationThe Fast Fourier Transform: A Brief Overview. with Applications. Petros Kondylis. Petros Kondylis. December 4, 2014
December 4, 2014 Timeline Researcher Date Length of Sequence Application CF Gauss 1805 Any Composite Integer Interpolation of orbits of celestial bodies F Carlini 1828 12 Harmonic Analysis of Barometric
More informationChapter 1 Divide and Conquer Polynomial Multiplication Algorithm Theory WS 2015/16 Fabian Kuhn
Chapter 1 Divide and Conquer Polynomial Multiplication Algorithm Theory WS 2015/16 Fabian Kuhn Formulation of the D&C principle Divide-and-conquer method for solving a problem instance of size n: 1. Divide
More information5.6 Convolution and FFT
5.6 Convolution and FFT Fast Fourier Transform: Applications Applications. Optics, acoustics, quantum physics, telecommunications, control systems, signal processing, speech recognition, data compression,
More informationE The Fast Fourier Transform
Fourier Transform Methods in Finance By Umberto Cherubini Giovanni Della Lunga Sabrina Mulinacci Pietro Rossi Copyright 2010 John Wiley & Sons Ltd E The Fast Fourier Transform E.1 DISCRETE FOURIER TRASFORM
More informationLecture 7 January 26, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 7 January 26, 26 Prof Emmanuel Candes Scribe: Carlos A Sing-Long, Edited by E Bates Outline Agenda:
More informationInterpolation on the unit circle
Lecture 2: The Fast Discrete Fourier Transform Interpolation on the unit circle So far, we have considered interpolation of functions at nodes on the real line. When we model periodic phenomena, it is
More informationFast Polynomial Multiplication
Fast Polynomial Multiplication Marc Moreno Maza CS 9652, October 4, 2017 Plan Primitive roots of unity The discrete Fourier transform Convolution of polynomials The fast Fourier transform Fast convolution
More informationDesign and Analysis of Algorithms
Design and Analysis of Algorithms CSE 5311 Lecture 5 Divide and Conquer: Fast Fourier Transform Junzhou Huang, Ph.D. Department of Computer Science and Engineering CSE5311 Design and Analysis of Algorithms
More informationPeriodic motions. Periodic motions are known since the beginning of mankind: Motion of planets around the Sun; Pendulum; And many more...
Periodic motions Periodic motions are known since the beginning of mankind: Motion of planets around the Sun; Pendulum; And many more... Periodic motions There are several quantities which describe a periodic
More informationCSE 548: Analysis of Algorithms. Lecture 4 ( Divide-and-Conquer Algorithms: Polynomial Multiplication )
CSE 548: Analysis of Algorithms Lecture 4 ( Divide-and-Conquer Algorithms: Polynomial Multiplication ) Rezaul A. Chowdhury Department of Computer Science SUNY Stony Brook Spring 2015 Coefficient Representation
More informationAlgorithm Design and Analysis
Algorithm Design and Analysis LECTURE 14 Divide and Conquer Fast Fourier Transform Sofya Raskhodnikova 10/7/2016 S. Raskhodnikova; based on slides by K. Wayne. 5.6 Convolution and FFT Fast Fourier Transform:
More informationFFT: Fast Polynomial Multiplications
FFT: Fast Polynomial Multiplications Jie Wang University of Massachusetts Lowell Department of Computer Science J. Wang (UMass Lowell) FFT: Fast Polynomial Multiplications 1 / 20 Overview So far we have
More informationChapter 1 Divide and Conquer Algorithm Theory WS 2016/17 Fabian Kuhn
Chapter 1 Divide and Conquer Algorithm Theory WS 2016/17 Fabian Kuhn Formulation of the D&C principle Divide-and-conquer method for solving a problem instance of size n: 1. Divide n c: Solve the problem
More informationFast Fourier Transform
Why Fourier Transform? Fast Fourier Transform Jordi Cortadella and Jordi Petit Department of Computer Science Polynomials: coefficient representation Divide & Conquer Dept. CS, UPC Polynomials: point-value
More informationThe Discrete Fourier Transform (DFT) Properties of the DFT DFT-Specic Properties Power spectrum estimate. Alex Sheremet.
4. April 2, 27 -order sequences Measurements produce sequences of numbers Measurement purpose: characterize a stochastic process. Example: Process: water surface elevation as a function of time Parameters:
More informationCS168: The Modern Algorithmic Toolbox Lecture #11: The Fourier Transform and Convolution
CS168: The Modern Algorithmic Toolbox Lecture #11: The Fourier Transform and Convolution Tim Roughgarden & Gregory Valiant May 8, 2015 1 Intro Thus far, we have seen a number of different approaches to
More informationFast Fourier Transform
Fast Fourier Transform December 8, 2016 FFT JPEG RGB Y C B C R (luma (brightness), chroma 2 (color)) chroma resolution is reduced image is split in blocks 8 8 pixels JPEG RGB Y C B C R (luma (brightness),
More informationLecture 19: Discrete Fourier Series
EE518 Digital Signal Processing University of Washington Autumn 2001 Dept. of Electrical Engineering Lecture 19: Discrete Fourier Series Dec 5, 2001 Prof: J. Bilmes TA: Mingzhou
More informationHow to Multiply. 5.5 Integer Multiplication. Complex Multiplication. Integer Arithmetic. Complex multiplication. (a + bi) (c + di) = x + yi.
How to ultiply Slides by Kevin Wayne. Copyright 5 Pearson-Addison Wesley. All rights reserved. integers, matrices, and polynomials Complex ultiplication Complex multiplication. a + bi) c + di) = x + yi.
More informationParallel Numerical Algorithms
Parallel Numerical Algorithms Chapter 13 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign CS 554 / CSE 512 Michael T. Heath Parallel Numerical Algorithms
More informationFig. 1: Fourier series Examples
FOURIER SERIES AND ITS SOME APPLICATIONS P. Sathyabama Assistant Professor, Department of Mathematics, Bharath collage of Science and Management, Thanjavur, Tamilnadu Abstract: The Fourier series, the
More informationKartsuba s Algorithm and Linear Time Selection
CS 374: Algorithms & Models of Computation, Fall 2015 Kartsuba s Algorithm and Linear Time Selection Lecture 09 September 22, 2015 Chandra & Manoj (UIUC) CS374 1 Fall 2015 1 / 32 Part I Fast Multiplication
More informationFrom Fourier to Wavelets in 60 Slides
From Fourier to Wavelets in 60 Slides Bernhard G. Bodmann Math Department, UH September 20, 2008 B. G. Bodmann (UH Math) From Fourier to Wavelets in 60 Slides September 20, 2008 1 / 62 Outline 1 From Fourier
More informationComputational Methods CMSC/AMSC/MAPL 460. Fourier transform
Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Ramani Duraiswami, Dept. of Computer Science Several slides from Prof. Healy s course at UMD Fourier Methods Fourier analysis ( harmonic analysis
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 informationDSP Algorithm Original PowerPoint slides prepared by S. K. Mitra
Chapter 11 DSP Algorithm Implementations 清大電機系林嘉文 cwlin@ee.nthu.edu.tw Original PowerPoint slides prepared by S. K. Mitra 03-5731152 11-1 Matrix Representation of Digital Consider Filter Structures This
More informationAn Adaptive Sublinear Time Block Sparse Fourier Transform
An Adaptive Sublinear Time Block Sparse Fourier Transform Volkan Cevher Michael Kapralov Jonathan Scarlett Amir Zandieh EPFL February 8th 217 Given x C N, compute the Discrete Fourier Transform (DFT) of
More informationThe Discrete Fourier transform
453.70 Linear Systems, S.M. Tan, The University of uckland 9- Chapter 9 The Discrete Fourier transform 9. DeÞnition When computing spectra on a computer it is not possible to carry out the integrals involved
More informationFundamentals of the DFT (fft) Algorithms
Fundamentals of the DFT (fft) Algorithms D. Sundararajan November 6, 9 Contents 1 The PM DIF DFT Algorithm 1.1 Half-wave symmetry of periodic waveforms.............. 1. The DFT definition and the half-wave
More informationNAME (1 pt): SID (1 pt): TA (1 pt): Name of Neighbor to your left (1 pt): Name of Neighbor to your right (1 pt):
CS 170 First Midterm 26 Feb 2010 NAME (1 pt): SID (1 pt): TA (1 pt): Name of Neighbor to your left (1 pt): Name of Neighbor to your right (1 pt): Instructions: This is a closed book, closed calculator,
More informationEDISP (NWL3) (English) Digital Signal Processing DFT Windowing, FFT. October 19, 2016
EDISP (NWL3) (English) Digital Signal Processing DFT Windowing, FFT October 19, 2016 DFT resolution 1 N-point DFT frequency sampled at θ k = 2πk N, so the resolution is f s/n If we want more, we use N
More informationChapter 4 Discrete Fourier Transform (DFT) And Signal Spectrum
Chapter 4 Discrete Fourier Transform (DFT) And Signal Spectrum CEN352, DR. Nassim Ammour, King Saud University 1 Fourier Transform History Born 21 March 1768 ( Auxerre ). Died 16 May 1830 ( Paris ) French
More informationThe O () notation. Definition: Let f(n), g(n) be functions of the natural (or real)
The O () notation When analyzing the runtime of an algorithm, we want to consider the time required for large n. We also want to ignore constant factors (which often stem from tricks and do not indicate
More informationThe Fastest Convolution in the West
The Fastest Convolution in the West Malcolm Roberts and John Bowman University of Alberta 2011-02-15 www.math.ualberta.ca/ mroberts 1 Outline Convolution Definition Applications Fast Convolutions Dealiasing
More informationThe Fourier Transform (and more )
The Fourier Transform (and more ) imrod Peleg ov. 5 Outline Introduce Fourier series and transforms Introduce Discrete Time Fourier Transforms, (DTFT) Introduce Discrete Fourier Transforms (DFT) Consider
More informationFall 2011, EE123 Digital Signal Processing
Lecture 6 Miki Lustig, UCB September 11, 2012 Miki Lustig, UCB DFT and Sampling the DTFT X (e jω ) = e j4ω sin2 (5ω/2) sin 2 (ω/2) 5 x[n] 25 X(e jω ) 4 20 3 15 2 1 0 10 5 1 0 5 10 15 n 0 0 2 4 6 ω 5 reconstructed
More information4.3 The Discrete Fourier Transform (DFT) and the Fast Fourier Transform (FFT)
CHAPTER. TIME-FREQUECY AALYSIS: FOURIER TRASFORMS AD WAVELETS.3 The Discrete Fourier Transform (DFT and the Fast Fourier Transform (FFT.3.1 Introduction In this section, we discuss some of the mathematics
More informationMA3232 Numerical Analysis Week 9. James Cooley (1926-)
MA umerical Analysis Week 9 James Cooley (96-) James Cooley is an American mathematician. His most significant contribution to the world of mathematics and digital signal processing is the Fast Fourier
More informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationTHE FAST FOURIER TRANSFORM III 10.3
10.3 The Fast Fourier Transform 439 THE FAST FOURIER TRANSFORM III 10.3 Many applications of linear algebra take time to develop. It is not easy to explain them in an hour. The teacher and the author must
More informationDIVIDE AND CONQUER II
DIVIDE AND CONQUER II master theorem integer multiplication matrix multiplication convolution and FFT Lecture slides by Kevin Wayne Copyright 2005 Pearson-Addison Wesley http://www.cs.princeton.edu/~wayne/kleinberg-tardos
More information12.1 Fourier Transform of Discretely Sampled Data
494 Chapter 12. Fast Fourier Transform now is: The PSD-per-unit-time converges to finite values at all frequencies except those where h(t) has a discrete sine-wave (or cosine-wave) component of finite
More informationTransforms and Orthogonal Bases
Orthogonal Bases Transforms and Orthogonal Bases We now turn back to linear algebra to understand transforms, which map signals between different domains Recall that signals can be interpreted as vectors
More informationSignal Processing COS 323
Signal Processing COS 323 Digital Signals D: functions of space or time e.g., sound 2D: often functions of 2 spatial dimensions e.g. images 3D: functions of 3 spatial dimensions CAT, MRI scans or 2 space,
More informationPhase recovery with PhaseCut and the wavelet transform case
Phase recovery with PhaseCut and the wavelet transform case Irène Waldspurger Joint work with Alexandre d Aspremont and Stéphane Mallat Introduction 2 / 35 Goal : Solve the non-linear inverse problem Reconstruct
More informationMath 56 Homework 5 Michael Downs
1. (a) Since f(x) = cos(6x) = ei6x 2 + e i6x 2, due to the orthogonality of each e inx, n Z, the only nonzero (complex) fourier coefficients are ˆf 6 and ˆf 6 and they re both 1 2 (which is also seen from
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 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 informationLecture 5. The Digital Fourier Transform. (Based, in part, on The Scientist and Engineer's Guide to Digital Signal Processing by Steven Smith)
Lecture 5 The Digital Fourier Transform (Based, in part, on The Scientist and Engineer's Guide to Digital Signal Processing by Steven Smith) 1 -. 8 -. 6 -. 4 -. 2-1 -. 8 -. 6 -. 4 -. 2 -. 2. 4. 6. 8 1
More informationFast Polynomials Multiplication Using FFT
Li Chen lichen.xd at gmail.com Xidian University January 17, 2014 Outline 1 Discrete Fourier Transform (DFT) 2 Discrete Convolution 3 Fast Fourier Transform (FFT) 4 umber Theoretic Transform (TT) 5 More
More informationCS483 Design and Analysis of Algorithms
CS483 Design and Analysis of Algorithms Lecture 6-8 Divide and Conquer Algorithms Instructor: Fei Li lifei@cs.gmu.edu with subject: CS483 Office hours: STII, Room 443, Friday 4:00pm - 6:00pm or by appointments
More informationCS711008Z Algorithm Design and Analysis
CS711008Z Algorithm Design and Analysis Lecture 5 FFT and Divide and Conquer Dongbo Bu Institute of Computing Technology Chinese Academy of Sciences, Beijing, China 1 / 56 Outline DFT: evaluate a polynomial
More informationDiscrete Fourier Transform
Discrete Fourier Transform Valentina Hubeika, Jan Černocký DCGM FIT BUT Brno, {ihubeika,cernocky}@fit.vutbr.cz Diskrete Fourier transform (DFT) We have just one problem with DFS that needs to be solved.
More informationEE123 Digital Signal Processing
Announcements EE Digital Signal Processing otes posted HW due Friday SDR give away Today! Read Ch 9 $$$ give me your names Lecture based on slides by JM Kahn M Lustig, EECS UC Berkeley M Lustig, EECS UC
More informationResponse to Periodic and Non-periodic Loadings. Giacomo Boffi.
Periodic and Non-periodic http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March 7, 218 Outline Undamped SDOF systems
More informationOutline. Introduction
Outline Periodic and Non-periodic Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano March 5, 4 A periodic loading is characterized by the identity p(t) = p(t + T ) where T is the period
More informationMATH 350: Introduction to Computational Mathematics
MATH 350: Introduction to Computational Mathematics Chapter VIII: The Fast Fourier Transform Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Spring 2008 Outline 1 The
More informationCS S Lecture 5 January 29, 2019
CS 6363.005.19S Lecture 5 January 29, 2019 Main topics are #divide-and-conquer with #fast_fourier_transforms. Prelude Homework 1 is due Tuesday, February 5th. I hope you ve at least looked at it by now!
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 informationMultiplying huge integers using Fourier transforms
Fourier transforms October 25, 2007 820348901038490238478324 1739423249728934932894??? integers occurs in many fields of Computational Science: Cryptography Number theory... Traditional approaches to
More informationFast Multipole Methods: Fundamentals & Applications. Ramani Duraiswami Nail A. Gumerov
Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov Week 1. Introduction. What are multipole methods and what is this course about. Problems from physics, mathematics,
More informationω 0 = 2π/T 0 is called the fundamental angular frequency and ω 2 = 2ω 0 is called the
he ime-frequency Concept []. Review of Fourier Series Consider the following set of time functions {3A sin t, A sin t}. We can represent these functions in different ways by plotting the amplitude versus
More informationIntroduction to the FFT
Introduction to the FFT 1 Introduction Assume that we have a signal fx where x denotes time. We would like to represent fx as a linear combination of functions e πiax or, equivalently, sinπax and cosπax
More informationEDISP (NWL2) (English) Digital Signal Processing Transform, FT, DFT. March 11, 2015
EDISP (NWL2) (English) Digital Signal Processing Transform, FT, DFT March 11, 2015 Transform concept We want to analyze the signal represent it as built of some building blocks (well known signals), possibly
More informationDivide-and-Conquer and Recurrence Relations: Notes on Chapter 2, Dasgupta et.al. Howard A. Blair
Divide-and-Conquer and Recurrence Relations: Notes on Chapter 2, Dasgupta et.al. Howard A. Blair 1 Multiplication Example 1.1: Carl Friedrich Gauss (1777-1855): The multiplication of two complex numbers
More informationResponse to Periodic and Non-periodic Loadings. Giacomo Boffi. March 25, 2014
Periodic and Non-periodic Dipartimento di Ingegneria Civile e Ambientale, Politecnico di Milano March 25, 2014 Outline Introduction Fourier Series Representation Fourier Series of the Response Introduction
More informationFast Convolution; Strassen s Method
Fast Convolution; Strassen s Method 1 Fast Convolution reduction to subquadratic time polynomial evaluation at complex roots of unity interpolation via evaluation at complex roots of unity 2 The Master
More information1. Calculation of the DFT
ELE E4810: Digital Signal Processing Topic 10: The Fast Fourier Transform 1. Calculation of the DFT. The Fast Fourier Transform algorithm 3. Short-Time Fourier Transform 1 1. Calculation of the DFT! Filter
More informationLecture 01 August 31, 2017
Sketching Algorithms for Big Data Fall 2017 Prof. Jelani Nelson Lecture 01 August 31, 2017 Scribe: Vinh-Kha Le 1 Overview In this lecture, we overviewed the six main topics covered in the course, reviewed
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationFourier Reconstruction from Non-Uniform Spectral Data
School of Electrical, Computer and Energy Engineering, Arizona State University aditya.v@asu.edu With Profs. Anne Gelb, Doug Cochran and Rosemary Renaut Research supported in part by National Science Foundation
More informationReview: Continuous Fourier Transform
Review: Continuous Fourier Transform Review: convolution x t h t = x τ h(t τ)dτ Convolution in time domain Derivation Convolution Property Interchange the order of integrals Let Convolution Property By
More informationSparse Fourier Transform (lecture 1)
1 / 73 Sparse Fourier Transform (lecture 1) Michael Kapralov 1 1 IBM Watson EPFL St. Petersburg CS Club November 2015 2 / 73 Given x C n, compute the Discrete Fourier Transform (DFT) of x: x i = 1 x n
More informationDown-Sampling (4B) Young Won Lim 10/25/12
Down-Sampling (4B) /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 informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
4 Interpolation 4.1 Polynomial interpolation Problem: LetP n (I), n ln, I := [a,b] lr, be the linear space of polynomials of degree n on I, P n (I) := { p n : I lr p n (x) = n i=0 a i x i, a i lr, 0 i
More informationLecture4 INF-MAT : 5. Fast Direct Solution of Large Linear Systems
Lecture4 INF-MAT 4350 2010: 5. Fast Direct Solution of Large Linear Systems Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo September 16, 2010 Test Matrix
More informationCh.11 The Discrete-Time Fourier Transform (DTFT)
EE2S11 Signals and Systems, part 2 Ch.11 The Discrete-Time Fourier Transform (DTFT Contents definition of the DTFT relation to the -transform, region of convergence, stability frequency plots convolution
More informationChapter 2: The Fourier Transform
EEE, EEE Part A : Digital Signal Processing Chapter Chapter : he Fourier ransform he Fourier ransform. Introduction he sampled Fourier transform of a periodic, discrete-time signal is nown as the discrete
More informationDiscrete Time Fourier Transform (DTFT)
Discrete Time Fourier Transform (DTFT) 1 Discrete Time Fourier Transform (DTFT) The DTFT is the Fourier transform of choice for analyzing infinite-length signals and systems Useful for conceptual, pencil-and-paper
More informationJim Lambers ENERGY 281 Spring Quarter Lecture 5 Notes
Jim ambers ENERGY 28 Spring Quarter 27-8 ecture 5 Notes These notes are based on Rosalind Archer s PE28 lecture notes, with some revisions by Jim ambers. Fourier Series Recall that in ecture 2, when we
More informationFourier Analysis and Spectral Representation of Signals
MIT 6.02 DRAFT Lecture Notes Last update: November 3, 2012 CHAPTER 13 Fourier Analysis and Spectral Representation of Signals We have seen in the previous chapter that the action of an LTI system on a
More informationDiscrete Fourier transform (DFT)
Discrete Fourier transform (DFT) Signal Processing 2008/9 LEA Instituto Superior Técnico Signal Processing LEA (IST) Discrete Fourier transform 1 / 34 Periodic signals Consider a periodic signal x[n] with
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 informationNONLINEAR COMPLEX MODULUS IN BITUMENS
NONLINEAR COMPLEX MODULUS IN BITUMENS J. Stastna, K. Jorshari and L. Zanzotto Bituminous Materials Chair University of Calgary 2500 University Drive NW Calgary, Alberta T2N 1N4 Canada Large Amplitude Oscillations
More informationReductions, Recursion and Divide and Conquer
Chapter 5 Reductions, Recursion and Divide and Conquer CS 473: Fundamental Algorithms, Fall 2011 September 13, 2011 5.1 Reductions and Recursion 5.1.0.1 Reduction Reducing problem A to problem B: (A) Algorithm
More informationRadar Systems Engineering Lecture 3 Review of Signals, Systems and Digital Signal Processing
Radar Systems Engineering Lecture Review of Signals, Systems and Digital Signal Processing Dr. Robert M. O Donnell Guest Lecturer Radar Systems Course Review Signals, Systems & DSP // Block Diagram of
More informationMathematics for Chemists 2 Lecture 14: Fourier analysis. Fourier series, Fourier transform, DFT/FFT
Mathematics for Chemists 2 Lecture 14: Fourier analysis Fourier series, Fourier transform, DFT/FFT Johannes Kepler University Summer semester 2012 Lecturer: David Sevilla Fourier analysis 1/25 Remembering
More informationAlgorithms and data structures
Algorithms and data structures Amin Coja-Oghlan LFCS Complex numbers Roots of polynomials A polynomial of degree d is a function of the form p(x) = d a i x i with a d 0. i=0 There are at most d numbers
More informationAlgorithms and Data Structures
Charles A. Wuethrich Bauhaus-University Weimar - CogVis/MMC June 1, 2017 1/44 Space reduction mechanisms Range searching Elementary Algorithms (2D) Raster Methods Shell sort Divide and conquer Quicksort
More information2 The Fourier Transform
2 The Fourier Transform The Fourier transform decomposes a signal defined on an infinite time interval into a λ-frequency component, where λ can be any real(or even complex) number. 2.1 Informal Development
More informationMATH 220 solution to homework 4
MATH 22 solution to homework 4 Problem. Define v(t, x) = u(t, x + bt), then v t (t, x) a(x + u bt) 2 (t, x) =, t >, x R, x2 v(, x) = f(x). It suffices to show that v(t, x) F = max y R f(y). We consider
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 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 informationlecture 7: Trigonometric Interpolation
lecture : Trigonometric Interpolation 9 Trigonometric interpolation for periodic functions Thus far all our interpolation schemes have been based on polynomials However, if the function f is periodic,
More informationChapter 6 THE SAMPLING PROCESS 6.1 Introduction 6.2 Fourier Transform Revisited
Chapter 6 THE SAMPLING PROCESS 6.1 Introduction 6.2 Fourier Transform Revisited Copyright c 2005 Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 14, 2018 Frame # 1 Slide # 1 A. Antoniou
More information