Fast Fourier Transform
|
|
- Phillip Goodman
- 6 years ago
- Views:
Transcription
1 Fast Fourier Transform December 8, 2016 FFT
2 JPEG RGB Y C B C R (luma (brightness), chroma 2 (color)) chroma resolution is reduced image is split in blocks 8 8 pixels
3 JPEG RGB Y C B C R (luma (brightness), chroma 2 (color)) chroma resolution is reduced image is split in blocks 8 8 pixels DCT - Discrete Cosine Transform (1974) 2D spatial domain frequency domain Functions can be decomposed into sums of sines and cosines (i.e. waves or harmonics) of various frequencies.
4 JPEG RGB Y C B C R (luma (brightness), chroma 2 (color)) chroma resolution is reduced image is split in blocks 8 8 pixels DCT - Discrete Cosine Transform (1974) 2D spatial domain frequency domain Functions can be decomposed into sums of sines and cosines (i.e. waves or harmonics) of various frequencies. quantization (irreversible): discard high-frequency info (human eye is more sensitive to small variations in color or brightness over large areas than to the strength of high-frequency brightness variations) further compression with lossless algorithm
5 2D spatial domain frequency domain DCT - Discrete Cosine Transform, acts on the sample x k DFT - Discrete Fourier Transform, acts on the sample x k Fourier Transform, f (x) nonrepetitive signal, -period, f = Fourier series, f (x) repetitive signal, oscillation, 2π-period f (x) k=0
6 Spatial domain frequency domain Source: Wikipedia
7 History and Motivation 3rd century BC: empiric model of planetary motions, based on deferents and epicycles
8 History and Motivation 3rd century BC: empiric model of planetary motions, based on deferents and epicycles d Alembert: wave equation - string oscillations Bernoulli: general solution of wave equation as harmonic series Leonard Euler: relation between the coeffecient of the series and parameters of the string
9 History and Motivation 3rd century BC: empiric model of planetary motions, based on deferents and epicycles d Alembert: wave equation - string oscillations Bernoulli: general solution of wave equation as harmonic series Leonard Euler: relation between the coeffecient of the series and parameters of the string Jean-Baptiste Joseph Fourier: solutions of the heat equation in a metal plate as a superposition (linear combination) of simple sine and cosine function Tomson, Michelson: mechanical harmonic analysor for tides predictions and optics Dirichlet and Riemann: math precision
10 Fourier Series Any 2π-periodic function f (x) can be expanded in Fourier series f (x) a a n cos(nx) + b n sin(nx) n=1
11 Fourier Series Any 2π-periodic function f (x) can be expanded in Fourier series f (x) a a n cos(nx) + b n sin(nx) n=1 = a a 1 cos x + b 1 sin x + a 2 cos(2x) + b 2 sin(2x) a N cos(nx) + b N sin(nx) +...
12 Fourier Series Any 2π-periodic function f (x) can be expanded in Fourier series f (x) a a n cos(nx) + b n sin(nx) n=1 = a a 1 cos x + b 1 sin x + a 2 cos(2x) + b 2 sin(2x) a N cos(nx) + b N sin(nx) +... where a 0 = 1 π a n = 1 π π π π π f (x) dx, f (x) cos(nx) dx, b n = 1 π π π f (x) sin(nx) dx
13 Example f (x) = x 2 cos 2 (x) on [0, π]
14 Example N = 0: f (x)
15 Example N = 1: f (x) cos x
16 Example N = 2: f (x) cos x cos(2x)
17 Example N = 3: f (x) cos x cos(2x) cos(3x)
18 Example N = 4: f (x) cos(4x)
19 Example N = 5: f (x) cos(5x)
20 Example N = 6: f (x) cos(6x)
21 Example N = 7: f (x) cos(7x)
22 Example N = 8: f (x) cos(8x)
23 Example N = 9: f (x) cos(9x)
24 Generalization to complex-valued functions: where f (x) = n= c n e inx c n = 1 π f (x)e inx dx 2π π Fourier series converges at each point x for piecewise smooth functions!! FFT
25 Nonperiodic Functions: Fourier Transform: time (sec) > frequency (hertz) For f (x) integrable on (, ) f (ξ) = f (x)e 2πixξ dx Inverse Fourier Transform f (x) = f (ξ)e 2πixξ dξ FFT
26 Discrete Fourier Transform Sequence of N complex numbers x 0, x 1,..., x N 1 is transformed into N-periodic sequence of complex numbers X k = N 1 n=0 Thus the domain in k [0, N 1] x k e 2πikn/N, k Z FFT
27 Discrete Fourier Transform - Linear C n C n Euler s formula gives X k = N 1 n=0 x n (cos( 2πk n N ) + i sin( 2πk n N )), k Z Inverse Fourier Transform: x n = 1 N N 1 k=0 X k (cos(2πk n N ) + i sin(2πk n N )), n Z N = number of time samples we have n = current sample we re considering (0,..., N 1) x n = value of the signal at time n k = current frequency we re considering (0 Hertz up to N 1 Hertz) X k = amount of frequency k in the signal (Amplitude and Phase, a complex number) FFT
28 Direct computation - Fourier Matrix ω ω 2 ω 3 ω N 1 W N = 1 1 ω 2 ω 4 ω 6 ω 2(N 1) N 1 ω 3 ω 6 ω 9 ω 3(N 1), ω N 1 ω 2(N 1) ω 3(N 1) ω (N 1)(N 1) Notice W N = W T N, complex number ω = ei2π/n, so ω N = 1 (N-th root of unity) Then Fourier transform X = W N x FFT
29 Example for N = 4 ω 0 ω 0 ω 0 ω W 4 = ω 0 ω 1 ω 2 ω 3 ω 0 ω 2 ω 0 ω 2 = 1 i 1 i ω 0 ω 3 ω 2 ω 1 1 i 1 i where ω = e πi 2 = iω = e πi 2 = i. Property: 1 4 W H 4 W 4 = I (so called hermitian trasnform) FFT
30 Example Single impulse at time 0 is x = [ ] T DFT (x = X = W 4 x = i 1 i i 1 i = Meaning: A single impulse has all frequencies in equal amounts. Get back, inverse transform i 1 i i 1 i = = FFT
31 Problem Direct computation X = W N x requires N 2 multiplications and N(N 1) additions. That s A LOT!!! Hence: FFT - Fast Fourier Transform: J.W. Cooley (IBM) and J.W. Tukey (Princeton) in 1960s Idea may be traced back to Gauss. FFT
32 The discrete Fourier transform (DFT): X k = N 1 n=0 x n e 2πi N nk, where 0 k N 1. For N as power of 2, N = 2 l, calculate separately ( decimation in time ) (x 2m = x 0, x 2,..., x N 2 ) (x 2m+1 = x 1, x 3,..., x N 1 ) Recursively: runtime O(N log N) FFT
33 Split: X k = = N/2 1 m=0 N/2 1 x 2m e 2πi N (2m)k + N/2 1 m=0 x 2m e 2πi N/2 mk m=0 }{{} DFT of even indexed part of x m x 2m+1 e 2πi N (2m+1)k + e 2πi N k N/2 1 m=0 x 2m+1 e 2πi N/2 mk }{{} DFT of odd indexed part of x m = EVEN k + e 2πi N k ODD k. Periodicity: and E k+ N 2 O k+ N 2 = E k = O k. FFT
34 For 0 k N/2: size-2 DFT ( butterfly ) X k = E k + e 2πi N k O k 0, 1,..., N/2 1 X k+ N 2 = E k e 2πi N k O k N/2, N/2 + 1,..., N 1 Recursive procedure: split even and odd x[2n] and x[2n + 1] compute DFT split divide and conquer algorithm - multi-branched recursion FFT
35 FFT
36 Fourier unit pulse FFT
PHYS 502 Lecture 3: Fourier Series
PHYS 52 Lecture 3: Fourier Series Fourier Series Introduction In mathematics, a Fourier series decomposes periodic functions or periodic signals into the sum of a (possibly infinite) set of simple oscillating
More informationComputer Problems for Fourier Series and Transforms
Computer Problems for Fourier Series and Transforms 1. Square waves are frequently used in electronics and signal processing. An example is shown below. 1 π < x < 0 1 0 < x < π y(x) = 1 π < x < 2π... and
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 informationExamples of the Fourier Theorem (Sect. 10.3). The Fourier Theorem: Continuous case.
s of the Fourier Theorem (Sect. 1.3. The Fourier Theorem: Continuous case. : Using the Fourier Theorem. The Fourier Theorem: Piecewise continuous case. : Using the Fourier Theorem. The Fourier Theorem:
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 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 informationMore on Fourier Series
More on Fourier Series R. C. Trinity University Partial Differential Equations Lecture 6.1 New Fourier series from old Recall: Given a function f (x, we can dilate/translate its graph via multiplication/addition,
More informationPeriodic functions: simple harmonic oscillator
Periodic functions: simple harmonic oscillator Recall the simple harmonic oscillator (e.g. mass-spring system) d 2 y dt 2 + ω2 0y = 0 Solution can be written in various ways: y(t) = Ae iω 0t y(t) = A cos
More informationAnalysis II: Fourier Series
.... Analysis II: Fourier Series Kenichi Maruno Department of Mathematics, The University of Texas - Pan American May 3, 011 K.Maruno (UT-Pan American) Analysis II May 3, 011 1 / 16 Fourier series were
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 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 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 informationAlgorithms of Scientific Computing
Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Michael Bader Technical University of Munich Summer 2018 Fast Fourier Transform Outline Discrete Fourier transform Fast Fourier transform
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 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 informationWeighted SS(E) = w 2 1( y 1 y1) y n. w n. Wy W y 2 = Wy WAx 2. WAx = Wy. (WA) T WAx = (WA) T b
6.8 - Applications of Inner Product Spaces Weighted Least-Squares Sometimes you want to get a least-squares solution to a problem where some of the data points are less reliable than others. In this case,
More information8/19/16. Fourier Analysis. Fourier analysis: the dial tone phone. Fourier analysis: the dial tone phone
Patrice Koehl Department of Biological Sciences National University of Singapore http://www.cs.ucdavis.edu/~koehl/teaching/bl5229 koehl@cs.ucdavis.edu Fourier analysis: the dial tone phone We use Fourier
More informationMATH 124B: HOMEWORK 2
MATH 24B: HOMEWORK 2 Suggested due date: August 5th, 26 () Consider the geometric series ( ) n x 2n. (a) Does it converge pointwise in the interval < x
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 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 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 informationPower Series Solutions We use power series to solve second order differential equations
Objectives Power Series Solutions We use power series to solve second order differential equations We use power series expansions to find solutions to second order, linear, variable coefficient equations
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 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 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 informationf (t) K(t, u) d t. f (t) K 1 (t, u) d u. Integral Transform Inverse Fourier Transform
Integral Transforms Massoud Malek An integral transform maps an equation from its original domain into another domain, where it might be manipulated and solved much more easily than in the original domain.
More informationIntelligent Visual Prosthesis
Orientation sensor (IMU) Intelligent Visual Prosthesis Depth image-based obstacle detection Depth camera Wideangle RGB camera Simultaneous object recognition, localization, and hand tracking New projects:
More informationImage Filtering, Edges and Image Representation
Image Filtering, Edges and Image Representation Capturing what s important Req reading: Chapter 7, 9 F&P Adelson, Simoncelli and Freeman (handout online) Opt reading: Horn 7 & 8 FP 8 February 19, 8 A nice
More informationMath 3150 HW 3 Solutions
Math 315 HW 3 Solutions June 5, 18 3.8, 3.9, 3.1, 3.13, 3.16, 3.1 1. 3.8 Make graphs of the periodic extensions on the region x [ 3, 3] of the following functions f defined on x [, ]. Be sure to indicate
More information6.003: Signals and Systems. Sampling and Quantization
6.003: Signals and Systems Sampling and Quantization December 1, 2009 Last Time: Sampling and Reconstruction Uniform sampling (sampling interval T ): x[n] = x(nt ) t n Impulse reconstruction: x p (t) =
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 informationTime-Frequency Analysis
Time-Frequency Analysis Basics of Fourier Series Philippe B. aval KSU Fall 015 Philippe B. aval (KSU) Fourier Series Fall 015 1 / 0 Introduction We first review how to derive the Fourier series of a function.
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 informationESE 250: Digital Audio Basics. Week 4 February 5, The Frequency Domain. ESE Spring'13 DeHon, Kod, Kadric, Wilson-Shah
ESE 250: Digital Audio Basics Week 4 February 5, 2013 The Frequency Domain 1 Course Map 2 Musical Representation With this compact notation Could communicate a sound to pianist Much more compact than 44KHz
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 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 informationFOURIER SERIES, HAAR WAVELETS AND FAST FOURIER TRANSFORM
FOURIER SERIES, HAAR WAVELETS AD FAST FOURIER TRASFORM VESA KAARIOJA, JESSE RAILO AD SAMULI SILTAE Abstract. This handout is for the course Applications of matrix computations at the University of Helsinki
More informationEmily Jennings. Georgia Institute of Technology. Nebraska Conference for Undergraduate Women in Mathematics, 2012
δ 2 Transform and Fourier Series of Functions with Multiple Jumps Georgia Institute of Technology Nebraska Conference for Undergraduate Women in Mathematics, 2012 Work performed at Kansas State University
More informationFOURIER SERIES. Chapter Introduction
Chapter 1 FOURIER SERIES 1.1 Introduction Fourier series introduced by a French physicist Joseph Fourier (1768-1830), is a mathematical tool that converts some specific periodic signals into everlasting
More informationIntroduction to Fourier Analysis Part 2. CS 510 Lecture #7 January 31, 2018
Introduction to Fourier Analysis Part 2 CS 510 Lecture #7 January 31, 2018 OpenCV on CS Dept. Machines 2/4/18 CSU CS 510, Ross Beveridge & Bruce Draper 2 In the extreme, a square wave Graphic from http://www.mechatronics.colostate.edu/figures/4-4.jpg
More informationSequential 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 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 informationSYDE 575: Introduction to Image Processing. Image Compression Part 2: Variable-rate compression
SYDE 575: Introduction to Image Processing Image Compression Part 2: Variable-rate compression Variable-rate Compression: Transform-based compression As mentioned earlier, we wish to transform image data
More information3 rd class Mech. Eng. Dept. hamdiahmed.weebly.com Fourier Series
Definition 1 Fourier Series A function f is said to be piecewise continuous on [a, b] if there exists finitely many points a = x 1 < x 2
More informationELEG 305: Digital Signal Processing
ELEG 5: Digital Signal Processing Lecture 6: The Fast Fourier Transform; Radix Decimatation in Time Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 8 K.
More informationProperties of Fourier Cosine and Sine Transforms
International Journal of Sciences: Basic and Applied Research (IJSBAR) ISSN 237-4531 (Print & Online) http://gssrr.org/index.php?journal=journalofbasicandapplied ---------------------------------------------------------------------------------------------------------------------------
More informationChapter 7: Techniques of Integration
Chapter 7: Techniques of Integration MATH 206-01: Calculus II Department of Mathematics University of Louisville last corrected September 14, 2013 1 / 43 Chapter 7: Techniques of Integration 7.1. Integration
More informationCSCI 1290: Comp Photo
CSCI 1290: Comp Photo Fall 2018 @ Brown University James Tompkin Many slides thanks to James Hays old CS 129 course, along with all of its acknowledgements. Capture Frequency - Rolling `Shutter James
More informationFOURIER ANALYSIS. (a) Fourier Series
(a) Fourier Series FOURIER ANAYSIS (b) Fourier Transforms Useful books: 1. Advanced Mathematics for Engineers and Scientists, Schaum s Outline Series, M. R. Spiegel - The course text. We follow their notation
More informationSequential Fast Fourier Transform
Sequential Fast Fourier Transform Departement Computerwetenschappen Room 03.018 (Celestijnenlaan 200A) albert-jan.yzelman@cs.kuleuven.be Includes material from slides by Prof. dr. Rob H. Bisseling Applications
More informationFourier Series and the Discrete Fourier Expansion
2 2.5.5 Fourier Series and the Discrete Fourier Expansion Matthew Lincoln Adrienne Carter sillyajc@yahoo.com December 5, 2 Abstract This article is intended to introduce the Fourier series and the Discrete
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 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 informationMa 221 Eigenvalues and Fourier Series
Ma Eigenvalues and Fourier Series Eigenvalue and Eigenfunction Examples Example Find the eigenvalues and eigenfunctions for y y 47 y y y5 Solution: The characteristic equation is r r 47 so r 44 447 6 Thus
More informationVU Signal and Image Processing. Torsten Möller + Hrvoje Bogunović + Raphael Sahann
052600 VU Signal and Image Processing Torsten Möller + Hrvoje Bogunović + Raphael Sahann torsten.moeller@univie.ac.at hrvoje.bogunovic@meduniwien.ac.at raphael.sahann@univie.ac.at vda.cs.univie.ac.at/teaching/sip/17s/
More informationThe wave equation in one dimension
The wave equation in one dimension Prof. Joyner 1 The theory of the vibrating string touches on musical theory and the theory of oscillating waves, so has likely been a concern of scholars since ancient
More informationLab Fourier Analysis Do prelab before lab starts. PHSX 262 Spring 2011 Lecture 5 Page 1. Based with permission on lectures by John Getty
Today /5/ Lecture 5 Fourier Series Time-Frequency Decomposition/Superposition Fourier Components (Ex. Square wave) Filtering Spectrum Analysis Windowing Fast Fourier Transform Sweep Frequency Analyzer
More informationLECTURE Fourier Transform theory
LECTURE 3 3. Fourier Transform theory In section (2.2) we introduced the Discrete-time Fourier transform and developed its most important properties. Recall that the DTFT maps a finite-energy sequence
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 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 informationMathematics for Engineers II. lectures. Power series, Fourier series
Power series, Fourier series Power series, Taylor series It is a well-known fact, that 1 + x + x 2 + x 3 + = n=0 x n = 1 1 x if 1 < x < 1. On the left hand side of the equation there is sum containing
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 information3150 Review Problems for Final Exam. (1) Find the Fourier series of the 2π-periodic function whose values are given on [0, 2π) by cos(x) 0 x π f(x) =
350 Review Problems for Final Eam () Find the Fourier series of the 2π-periodic function whose values are given on [0, 2π) by cos() 0 π f() = 0 π < < 2π (2) Let F and G be arbitrary differentiable functions
More informationThe Gibbs Phenomenon
The Gibbs Phenomenon Eli Dean & Britton Girard Math 572, Fall 2013 December 4, 2013 1. Introduction 1 S N f(x) := f(n)e inx := 1 2π n N n N denotes the partial Fourier sum for f at the point x f(x)e inx
More informationThe Frequency Domain : Computational Photography Alexei Efros, CMU, Fall Many slides borrowed from Steve Seitz
The Frequency Domain 15-463: Computational Photography Alexei Efros, CMU, Fall 2008 Somewhere in Cinque Terre, May 2005 Many slides borrowed from Steve Seitz Salvador Dali Gala Contemplating the Mediterranean
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 informationProblem Set 8 - Solution
Problem Set 8 - Solution Jonasz Słomka Unless otherwise specified, you may use MATLAB to assist with computations. provide a print-out of the code used and its output with your assignment. Please 1. More
More informationMath 489AB A Very Brief Intro to Fourier Series Fall 2008
Math 489AB A Very Brief Intro to Fourier Series Fall 8 Contents Fourier Series. The coefficients........................................ Convergence......................................... 4.3 Convergence
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 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 informationUniversity of Connecticut Lecture Notes for ME5507 Fall 2014 Engineering Analysis I Part III: Fourier Analysis
University of Connecticut Lecture Notes for ME557 Fall 24 Engineering Analysis I Part III: Fourier Analysis Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
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 informationFourier Series and Recent Developments in Analysis
Fourier Series and Recent Developments in Analysis Karlstad, June 2003 Javier Soria (U. Barcelona) 1 Jean Baptiste Joseph Fourier (1768-1830) It was around 1804 that Fourier did his important mathematical
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 informationPARTIAL DIFFERENTIAL EQUATIONS (MATH417) SOLUTIONS FOR THE FINAL EXAM
PARTIAL DIFFERENTIAL EQUATIONS (MATH417) SOLUTIONS FOR THE FINAL EXAM Problem 1 (1 pts.) Classify each equation as linear homogeneous, linear inhomogeneous, or nonlinear: a) u = u x + x 3 b) u + u c) u
More informationPREMED COURSE, 14/08/2015 OSCILLATIONS
PREMED COURSE, 14/08/2015 OSCILLATIONS PERIODIC MOTIONS Mechanical Metronom Laser Optical Bunjee jumping Electrical Astronomical Pulsar Biological ECG AC 50 Hz Another biological exampe PERIODIC MOTIONS
More informationencoding without prediction) (Server) Quantization: Initial Data 0, 1, 2, Quantized Data 0, 1, 2, 3, 4, 8, 16, 32, 64, 128, 256
General Models for Compression / Decompression -they apply to symbols data, text, and to image but not video 1. Simplest model (Lossless ( encoding without prediction) (server) Signal Encode Transmit (client)
More informationChapter 6: Fast Fourier Transform and Applications
Chapter 6: Fast Fourier Transform and Applications Michael Hanke Mathematical Models, Analysis and Simulation, Part I Read: Strang, Ch. 4. Fourier Sine Series In the following, every function f : [,π]
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 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 information5 Trigonometric Functions
5 Trigonometric Functions 5.1 The Unit Circle Definition 5.1 The unit circle is the circle of radius 1 centered at the origin in the xyplane: x + y = 1 Example: The point P Terminal Points (, 6 ) is on
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 information/ (2π) X(e jω ) dω. 4. An 8 point sequence is given by x(n) = {2,2,2,2,1,1,1,1}. Compute 8 point DFT of x(n) by
Code No: RR320402 Set No. 1 III B.Tech II Semester Regular Examinations, Apr/May 2006 DIGITAL SIGNAL PROCESSING ( Common to Electronics & Communication Engineering, Electronics & Instrumentation Engineering,
More 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 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 informationBernoulli Numbers and their Applications
Bernoulli Numbers and their Applications James B Silva Abstract The Bernoulli numbers are a set of numbers that were discovered by Jacob Bernoulli (654-75). This set of numbers holds a deep relationship
More informationFourier Series User Guide
Fourier Series User Guide K.N. Toosi University Of Technology Electrical And Computer Engineering Faculty Dr. Hadi Aliakbarian Autumn 2017 Contents 1 Requirements 1 2 Brief Description Of Fourier Series
More informationFourier and Partial Differential Equations
Chapter 5 Fourier and Partial Differential Equations 5.1 Fourier MATH 294 SPRING 1982 FINAL # 5 5.1.1 Consider the function 2x, 0 x 1. a) Sketch the odd extension of this function on 1 x 1. b) Expand the
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 informationFall Math 3410 Name (Print): Solution KEY Practice Exam 2 - November 4 Time Limit: 50 Minutes
Fall 206 - Math 340 Name (Print): Solution KEY Practice Exam 2 - November 4 Time Limit: 50 Minutes This exam contains pages (including this cover page) and 5 problems. Check to see if any pages are missing.
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 informationFourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012)
Fourier Transform in Image Processing CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) Basis Decomposition Write a function as a weighted sum of basis functions f ( x) wibi(
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 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 information7.2 Trigonometric Integrals
7.2 1 7.2 Trigonometric Integrals Products of Powers of Sines and Cosines We wish to evaluate integrals of the form: sin m x cos n xdx where m and n are nonnegative integers. Recall the double angle formulas
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 informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 43 Outline 1 Introduction Motivation History
More information2. the basis functions have different symmetries. 1 k = 0. x( t) 1 t 0 x(t) 0 t 1
In the next few lectures, we will look at a few examples of orthobasis expansions that are used in modern signal processing. Cosine transforms The cosine-i transform is an alternative to Fourier series;
More informationMath 5587 Midterm II Solutions
Math 5587 Midterm II Solutions Prof. Jeff Calder November 3, 2016 Name: Instructions: 1. I recommend looking over the problems first and starting with those you feel most comfortable with. 2. Unless otherwise
More information