Algorithms of Scientific Computing
|
|
- Drusilla Hawkins
- 6 years ago
- Views:
Transcription
1 Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Michael Bader Technical University of Munich Summer 2018
2 Fast Fourier Transform Outline Discrete Fourier transform Fast Fourier transform Special Fourier transform: real-valued FFT sine/cosine transform Applications: Fast Poisson solver (FST) Computergraphics (FCT) Efficient Implementation Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
3 Discrete Fourier Transform (DFT) Definition: For a vector of complex numbers (f 0,..., f 1 ) T, the discrete Fourier transform (DFT) is given by the vector (F 0,..., F 1 ) T, where Interpretation: The DFT can be derived as F k = 1 1 f n e i2πnk/. n=0 trigonometric interpolation/approxmiation approximation of the coefficients of the Fourier series Applications: JPEG, MPEG, MP3, etc. signal processing in general, solvers in science and engineering Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
4 DFT as Interpolation (1) Interpolation problem: ansatz functions: g k (x) := e ikx in the interval [0, 2π], k = 0,..., 1 supporting points: x n := 2πn/, n = 0,..., 1 interpolation value f n, n = 0,..., 1 find weights F k such that at all supporting points 1 1 f n = F k g k (x n ) f n = F k e i2πnk/. trigonometric interpolation Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
5 DFT as Interpolation (2) Interpolation problem: ansatz functions: g k (z) := z k (complex unit polynomials), k = 0,..., 1 supporting points: z n := e i2πn/ = ω n, where ω := e i2π/ interpolation values f n, n = 0,..., 1, respectively. find the weights F k such that at all supporting points 1 1 f n = F k g k (z n ) f n = F k e i2πnk/. Polynomial interpolation at the complex unit roots ω n Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
6 Interpretation of the Interpolation Problem Starting from the first formulation, 1 f n = F k g k (x n ), g k (x n ) = e i2πnk/, we look for a representation of the signal f n or of a function f (x) of the form 1 f (x) = F k g k (x), g k (x) = e i2πkx. The ansatz functions are sine or cosine oscillations: e ikx = cos(kx) + i sin(kx) Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
7 Interpretation of the Interpolation Problem (2) Conclusions: we look for the representation of a periodic function as a sum of sines and cosines the F k are, thus, called Fourier coefficients: k represents the wave number the value of F k represents the amplitude of the corresponding frequency the Fourier transform leads to a frequency spectrum useful when a problem is easier to solve in the frequency domain than in the spatial domain. Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
8 Example: Spatial vs. Frequency Domain 3D Kurve des Signals s 1 3D Kurve des Signals s 2 3D Kurve des Signals s 3 Im(s 1 ) Im(s 2 ) Im(s 3 ) t Werte Re(s ) t Werte Re(s 1 2 ) t Werte Re(s ) 3 s 1 = e 3it, s 2 = 0.2 i e it, s 3 = e it e 12it Frequenzspektrum f Frequenzspektrum f 2 3 Frequenzspektrum f 1 Anteile Anteile Anteile Frequenzen Frequenzen Frequenzen Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
9 Solution of the Interpolation Problem Both interpolation problems lead to the identical linear systems of equations: 1 f n = F k ω nk, for all n = 0,..., 1; where ω := e i2π/, i.e. ω nk := ei2πnk/ ( unit roots : ω = 1) If we write the vectors of the f n and F k as f := (f 0,..., f 1 ) and F := (F 0,..., F 1 ), the linear system of equations can be formulated in matrix-vector notation WF = f, where the entries of the Fourier matrix W are given by W nk := ω nk. Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
10 Properties of the Fourier Matrix W W is symmetric: W = W T, and has the form W = ω 0 ω 0 ω 0... ω 0 ω 0 ω 1 ω 2... ω ( 1) ω 0 ω 2 ω 4... ω 2( 1). ω 0. ω ( 1) W ( W T ) = WW H = I, since [ WW H] 1 ( = ω kj ω lj kl.. ω 2( 1)... ω ( 1)( 1) ) 1 = { ω (k l)j if k = l = 0 if k l. Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
11 Properties of the Fourier Matrix W (details) important property #1: (ω k ) = e i2πk/ = e i2πk = 1 important property #2: (ω k ) = ( e i2πk/) = (cos(2πk/) + i sin(2πk/)) = cos(2πk/) i sin(2πk/) = cos( 2πk/) + i sin( 2πk/) = e i2πk/ = ω k W ( W T ) = WW H = I, since if k = l then 1 if k l then 1 [ WW H] 1 ( = ω kj ω lj kl ω (k l)j ω (k l)j 1 (1 ξ) = 1 ω 0 = = 1 1 ξ j = ) 1 = ξ j where ξ = ω (k l) 1 ξ j ω kj ω lj and then: ξ (j+1) = 1 ξ, thus Remember that (ω k ) = 1 for any k ( unit roots!), thus 1 1 = ω (k l)j 1 ξ j = 1 ξ 1 ξ ω (k l)j = 0 Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
12 Computation of the Fourier Coefficients F k Since WW H = I, the inverse of W is W 1 = 1 WH : ω 0 ω 0 ω 0... ω 0 W 1 = 1 ω 0 ω 1 ω 2... ω ( 1) ω 0 ω 2 ω 4... ω 2( 1).... ω 0 ω ( 1) ω 2( 1)... ω ( 1)( 1) the vector F of the Fourier coefficients can be computed easily as a matrix-vector product with computational effort O( 2 ): F = 1 WH f or F k = 1 1 f n ω nk. n=0 Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
13 Inverse Discrete Fourier Transform (IDFT) Definition: The inverse Discrete Fourier Transform (IDFT) of the vector (F 0,..., F 1 ) is given by the vector (f 0,..., f 1 ), where 1 f n = F k e i2πnk/. Observation: DFT and IDFT are inverse operations: F k = f n e i2πnk/, f n = F k e i2πnk/. n=0 F = DFT(IDFT(F)) or f = IDFT(DFT(f)). Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
14 The Pair DFT/IDFT as Matrix-Vector Product With the notation ω := e i2π/, i.e. ω nk := e i2πnk/, we formulate the DFT/IDFT as F k = f n ω nk f n = F k ω nk n=0 With the vectors f := (f 0,..., f 1 ) T and F := (F 0,..., F 1 ) T, we denote (and compute) the DFT und IDFT as matrix-vector products F = 1 WH f, f = WF, where the elements of the Fourier matrix W are W nk := ω nk. Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
15 Properties of the DFT DFT and IDFT are (as a matrix-vector product) linear: since ω nk DFT(αf + βg) = α DFT(f ) + β DFT(g) IDFT(αf + βg) = α IDFT(f ) + β IDFT(g) = ωn(k+) = ω (n+)k, the f n and the F k are periodic: f n+ = f n F k+ = F k for all k, n Z Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
16 Alternative Forms of the DFT Possible variants (in all imaginable combinations): Scaling factor 1 in the IDFT instead of the DFT; 1 alternatively a factor in DFT and IDFT. switched signs in the exponent of the exponential function in DFT and IDFT use j for the imaginary unit (electrical engineering) Shift of indices: periodic data: F k = F k+ aliasing of frequencies: e i2πnk/ = e i2πn(k±)/ Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
17 DFT with Shifted Indices Data and frequencies symmetric : In general: F k = f n e i2πnk/, f n = F k e i2πnk/ n= 2 +1 k= 2 +1 F k = 1 P+ n=p+1 f n e i2πnk/, f n = Q+ k=q+1 F k e i2πnk/ Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
18 DFT in Program Libraries Matlab, IMSL (Int. Math. and Stat. Library): F k+1 = 1 f n+1 e i2πnk/ k = 0,..., 1 n=0 f n+1 = 1 1 F k+1 e i2πnk/ n = 0,..., 1 Maple: 1 as factor for DFT and IDFT. Index shift by +1, since: Data/coefficients start at index 0 Arrays to store the numbers start at index 1 Michael Bader Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Summer
Fast 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 informationMA3232 Summary 5. d y1 dy1. MATLAB has a number of built-in functions for solving stiff systems of ODEs. There are ode15s, ode23s, ode23t, ode23tb.
umerical solutions of higher order ODE We can convert a high order ODE into a system of first order ODEs and then apply RK method to solve it. Stiff ODEs Stiffness is a special problem that can arise in
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 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 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 informationMARN 5898 Fourier Analysis.
MAR 5898 Fourier Analysis. Dmitriy Leykekhman Spring 2010 Goals Fourier Series. Discrete Fourier Transforms. D. Leykekhman - MAR 5898 Parameter estimation in marine sciences Linear Least Squares 1 Complex
More informationLecture 6: Discrete Fourier Transform
Lecture 6: Discrete Fourier Transform In the previous lecture we introduced the discrete Fourier transform as given either by summations or as a matrix vector product The discrete Fourier transform of
More information1 (2n)! (-1)n (θ) 2n
Complex Numbers and Algebra The real numbers are complete for the operations addition, subtraction, multiplication, and division, or more suggestively, for the operations of addition and multiplication
More informationChapter 2 Algorithms for Periodic Functions
Chapter 2 Algorithms for Periodic Functions In this chapter we show how to compute the Discrete Fourier Transform using a Fast Fourier Transform (FFT) algorithm, including not-so special case situations
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 informationFOURIER TRANSFORM METHODS David Sandwell, January, 2013
1 FOURIER TRANSFORM METHODS David Sandwell, January, 2013 1. Fourier Transforms Fourier analysis is a fundamental tool used in all areas of science and engineering. The fast fourier transform (FFT) algorithm
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 informationFourier Analysis. 19th October 2015
Fourier Analysis Hilary Weller 19th October 2015 This is brief introduction to Fourier analysis and how it is used in atmospheric and oceanic science, for: Analysing data (eg climate
More informationThe Discrete Fourier Transform
The Discrete Fourier Transform / 23 The Fourier Transform Any discrete time series exponentials: can be built up of a sum of complex = X m e i2πmn/, n = 0,, m=0 = [ + + + + ]. X 0 X e i2πn (/) X 2 e i2πn
More informationMath 365 Homework 5 Due April 6, 2018 Grady Wright
Math 365 Homework 5 Due April 6, 018 Grady Wright 1. (Polynomial interpolation, 10pts) Consider the following three samples of some function f(x): x f(x) -1/ -1 1 1 1/ (a) Determine the unique second degree
More informationDiscrete Fourier transform
Discrete Fourier transform Alejandro Ribeiro Dept. of Electrical and Systems Engineering University of Pennsylvania aribeiro@seas.upenn.edu http://www.seas.upenn.edu/users/~aribeiro/ January 2, 216 Signal
More informationThe moral of the story regarding discontinuities: They affect the rate of convergence of Fourier series
Lecture 7 Inner product spaces cont d The moral of the story regarding discontinuities: They affect the rate of convergence of Fourier series As suggested by the previous example, discontinuities of a
More information18.085: Summer 2016 Due: 3 August 2016 (in class) Problem Set 8
Problem Set 8 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 on relation between Fourier
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationReview Problems for Test 2
Review Problems for Test Math 0 009 These problems are meant to help you study. The presence of a problem on this sheet does not imply that there will be a similar problem on the test. And the absence
More informationWaves, the Wave Equation, and Phase Velocity. We ll start with optics. The one-dimensional wave equation. What is a wave? Optional optics texts: f(x)
We ll start with optics Optional optics texts: Waves, the Wave Equation, and Phase Velocity What is a wave? f(x) f(x-) f(x-) f(x-3) Eugene Hecht, Optics, 4th ed. J.F. James, A Student's Guide to Fourier
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 018 30 Lecture 30 April 7, 018 Fourier Series In this lecture,
More informationA Library of Functions
LibraryofFunctions.nb 1 A Library of Functions Any study of calculus must start with the study of functions. Functions are fundamental to mathematics. In its everyday use the word function conveys to us
More informationLES of Turbulent Flows: Lecture 3
LES of Turbulent Flows: Lecture 3 Dr. Jeremy A. Gibbs Department of Mechanical Engineering University of Utah Fall 2016 1 / 53 Overview 1 Website for those auditing 2 Turbulence Scales 3 Fourier transforms
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 informationThe Fourier spectral method (Amath Bretherton)
The Fourier spectral method (Amath 585 - Bretherton) 1 Introduction The Fourier spectral method (or more precisely, pseudospectral method) is a very accurate way to solve BVPs with smooth solutions on
More informationFOURIER TRANSFORM AND
Version: 11 THE DISCRETE TEX d: Oct. 23, 2013 FOURIER TRANSFORM AND THE FFT PREVIEW Classical numerical analysis techniques depend largely on polynomial approximation of functions for differentiation,
More information23. The Finite Fourier Transform and the Fast Fourier Transform Algorithm
23. The Finite Fourier Transform and the Fast Fourier Transform Algorithm 23.1 Introduction: Fourier Series Early in the Nineteenth Century, Fourier studied sound and oscillatory motion and conceived of
More informationFiltering in the Frequency Domain
Filtering in the Frequency Domain Dr. Praveen Sankaran Department of ECE NIT Calicut January 11, 2013 Outline 1 Preliminary Concepts 2 Signal A measurable phenomenon that changes over time or throughout
More informationData Processing and Analysis
Data Processing and Analysis Rick Aster and Brian Borchers September 27, 2013 Sampled Time Series umerical scientific data are commonly organized into series or matrices, i.e., sets of spatially or temporally
More informationSummary of lecture 1. E x = E x =T. X T (e i!t ) which motivates us to define the energy spectrum Φ xx (!) = jx (i!)j 2 Z 1 Z =T. 2 d!
Summary of lecture I Continuous time: FS X FS [n] for periodic signals, FT X (i!) for non-periodic signals. II Discrete time: DTFT X T (e i!t ) III Poisson s summation formula: describes the relationship
More informationImage Analysis. 3. Fourier Transform
Image Analysis Image Analysis 3. Fourier Transform Lars Schmidt-Thieme Information Systems and Machine Learning Lab (ISMLL) Institute for Business Economics and Information Systems & Institute for Computer
More information1D Wave PDE. Introduction to Partial Differential Equations part of EM, Scalar and Vector Fields module (PHY2064) Richard Sear.
Introduction to Partial Differential Equations part of EM, Scalar and Vector Fields module (PHY2064) November 12, 2018 Wave equation in one dimension This lecture Wave PDE in 1D Method of Separation of
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 informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
More informationFourier Transform with Rotations on Circles and Ellipses in Signal and Image Processing
Fourier Transform with Rotations on Circles and Ellipses in Signal and Image Processing Artyom M. Grigoryan Department of Electrical and Computer Engineering University of Texas at San Antonio One UTSA
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 informationTutorial 9 The Discrete Fourier Transform (DFT) SIPC , Spring 2017 Technion, CS Department
Tutorial 9 The Discrete Fourier Transform (DFT) SIPC 236327, Spring 2017 Technion, CS Department The DFT Matrix The DFT matrix of size M M is defined as DFT = 1 M W 0 0 W 0 W 0 W where W = e i2π M i =
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 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 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 information2. Interpolation. Closing the Gaps of Discretization...
2. Interpolation Closing the Gaps of Discretization... Numerisches Programmieren, Hans-Joachim Bungartz page 1 of 48 2.1. Preliminary Remark The Approximation Problem Problem: Approximate a (fully or partially
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 informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationNotes on FFT-based differentiation
Notes on FFT-based differentiation Steven G Johnson, MIT Applied Mathematics Created April, 2011, updated May 4, 2011 Abstract A common numerical technique is to differentiate some sampled function y(x)
More informationIntroduction to Signal Processing
to Signal Processing Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Intelligent Systems for Pattern Recognition Signals = Time series Definitions Motivations A sequence
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 informationThe Discrete Fourier Transform
In [ ]: cd matlab pwd The Discrete Fourier Transform Scope and Background Reading This session introduces the z-transform which is used in the analysis of discrete time systems. As for the Fourier and
More informationLecture 6: Differential Equations Describing Vibrations
Lecture 6: Differential Equations Describing Vibrations In Chapter 3 of the Benson textbook, we will look at how various types of musical instruments produce sound, focusing on issues like how the construction
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 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 informationMultimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2
Multimedia Signals and Systems - Audio and Video Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Kunio Takaya Electrical and Computer Engineering University of Saskatchewan December
More informationTutorial 2 - Learning about the Discrete Fourier Transform
Tutorial - Learning about the Discrete Fourier Transform This tutorial will be about the Discrete Fourier Transform basis, or the DFT basis in short. What is a basis? If we google define basis, we get:
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 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 informationMATH 167: APPLIED LINEAR ALGEBRA Chapter 3
MATH 167: APPLIED LINEAR ALGEBRA Chapter 3 Jesús De Loera, UC Davis February 18, 2012 Orthogonal Vectors and Subspaces (3.1). In real life vector spaces come with additional METRIC properties!! We have
More informationDigital Signal Processing Module 6 Discrete Fourier Transform (DFT)
Objective: Digital Signal Processing Module 6 Discrete Fourier Transform (DFT) 1. To analyze finite duration sequences using Discrete Fourier Transform. 2. To understand the characteristics and importance
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 informationSolution Set #1
05-738-0083 Solution Set # i. Evaluate δ [SINC []] = δ sin[π] π Recall relation for Dirac delta function wit functional argument (Fourier Metods, δ [g []] = δ [ 0] (6.38) dg =0 Te SINC function as zeros
More information18.085, PROBLEM SET 8 SOLUTIONS
18.085, PROBLEM SET 8 SOLUTIONS Question 1. (15 pts.) Question 4 of pset7, continued. Let f(x) be the 2π-periodic extension of the function g(x) = e x defined on the interval π x π. See the solutions to
More informationReview of Linear System Theory
Review of Linear System Theory The following is a (very) brief review of linear system theory and Fourier analysis. I work primarily with discrete signals. I assume the reader is familiar with linear algebra
More informationDiscrete Fourier Transform
Last lecture I introduced the idea that any function defined on x 0,..., N 1 could be written a sum of sines and cosines. There are two different reasons why this is useful. The first is a general one,
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 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 informationFourier transforms. c n e inπx. f (x) = Write same thing in an equivalent form, using n = 1, f (x) = l π
Fourier transforms We can imagine our periodic function having periodicity taken to the limits ± In this case, the function f (x) is not necessarily periodic, but we can still use Fourier transforms (related
More informationWorksheet #2. Use trigonometry and logarithms to model natural phenomena from a periodic behavior perspective.
math 112 Worksheet #2 Name Due Date: Objective: Use trigonometry and logarithms to model natural phenomena from a periodic behavior perspective. We revisit the varying of tree ring widths over time. We
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 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 informationDISCRETE FOURIER TRANSFORM
DISCRETE FOURIER TRANSFORM 1. Introduction The sampled discrete-time fourier transform (DTFT) of a finite length, discrete-time signal is known as the discrete Fourier transform (DFT). The DFT contains
More informationMaths 381 The First Midterm Solutions
Maths 38 The First Midterm Solutions As Given by Peter Denton October 8, 8 # From class, we have that If we make the substitution y = x 3/, we get # Once again, we use the fact that and substitute into
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 informationChapter 1. Discrete Fourier Transform. 1.1 Fourier series
Chapter Discrete Fourier Transform We usually think about processes around us as functions of time. However, it is often useful to think about them as functions of frequencies. We naturally do this without
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 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 informationCOMPRESSED SENSING IN PYTHON
COMPRESSED SENSING IN PYTHON Sercan Yıldız syildiz@samsi.info February 27, 2017 OUTLINE A BRIEF INTRODUCTION TO COMPRESSED SENSING A BRIEF INTRODUCTION TO CVXOPT EXAMPLES A Brief Introduction to Compressed
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 informationFourier Series and Fourier Transform
Fourier Series and Fourier Transform An Introduction Michael Figl Center for Medical Physics and Biomedical Engineering Medical University of Vienna 1 / 36 Introduction 2 / 36 Introduction We want to recall
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 informationIntroduction to HPC. Lecture 20
COSC Introduction to HPC Lecture Dept of Computer Science COSC Fast Fourier Transform COSC Image Processing In Electron Microscopy Medical Imaging Signal Processing In Radio Astronomy Modeling and Solution
More informationDigital Signal Processing: Signal Transforms
Digital Signal Processing: Signal Transforms Aishy Amer, Mohammed Ghazal January 19, 1 Instructions: 1. This tutorial introduces frequency analysis in Matlab using the Fourier and z transforms.. More Matlab
More 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 informationAppendix C: Recapitulation of Numerical schemes
Appendix C: Recapitulation of Numerical schemes August 31, 2009) SUMMARY: Certain numerical schemes of general use are regrouped here in order to facilitate implementations of simple models C1 The tridiagonal
More informationCircuits and Systems I
Circuits and Systems I LECTURE #2 Phasor Addition lions@epfl Prof. Dr. Volkan Cevher LIONS/Laboratory for Information and Inference Systems License Info for SPFirst Slides This work released under a Creative
More informationFOURIER TRANSFORMS. At, is sometimes taken as 0.5 or it may not have any specific value. Shifting at
Chapter 2 FOURIER TRANSFORMS 2.1 Introduction The Fourier series expresses any periodic function into a sum of sinusoids. The Fourier transform is the extension of this idea to non-periodic functions by
More informationReview of Linear Systems Theory
Review of Linear Systems Theory The following is a (very) brief review of linear systems theory, convolution, and Fourier analysis. I work primarily with discrete signals, but each result developed in
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 5. Complex Numbers and Euler s Formula
Lecture 5. Complex Numbers and Euler s Formula University of British Columbia, Vancouver Yue-Xian Li March 017 1 Main purpose: To introduce some basic knowledge of complex numbers to students so that they
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 informationC. Complex Numbers. 1. Complex arithmetic.
C. Complex Numbers. Complex arithmetic. Most people think that complex numbers arose from attempts to solve quadratic equations, but actually it was in connection with cubic equations they first appeared.
More information3. Lecture. Fourier Transformation Sampling
3. Lecture Fourier Transformation Sampling Some slides taken from Digital Image Processing: An Algorithmic Introduction using Java, Wilhelm Burger and Mark James Burge Separability ² The 2D DFT can be
More informationPhysics 142 Wave Optics 1 Page 1. Wave Optics 1. For every complex problem there is one solution that is simple, neat, and wrong. H.L.
Physics 142 Wave Optics 1 Page 1 Wave Optics 1 For every complex problem there is one solution that is simple, neat, and wrong. H.L. Mencken Interference and diffraction of waves The essential characteristic
More informationEEL3135: Homework #3 Solutions
EEL335: Homework #3 Solutions Problem : (a) Compute the CTFT for the following signal: xt () cos( πt) cos( 3t) + cos( 4πt). First, we use the trigonometric identity (easy to show by using the inverse Euler
More informationENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin
ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM Dr. Lim Chee Chin Outline Introduction Discrete Fourier Series Properties of Discrete Fourier Series Time domain aliasing due to frequency
More 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 informationSection 3: Complex numbers
Essentially: Section 3: Complex numbers C (set of complex numbers) up to different notation: the same as R 2 (euclidean plane), (i) Write the real 1 instead of the first elementary unit vector e 1 = (1,
More informationModified Fourier Spectral Method for Non-periodic CFD
Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu Chaoqun Liu Technical Report 2011-14 http://www.uta.edu/math/preprint/ Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu
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 informationLecture 4: Oscillators to Waves
Matthew Schwartz Lecture 4: Oscillators to Waves Review two masses Last time we studied how two coupled masses on springs move If we take κ=k for simplicity, the two normal modes correspond to ( ) ( )
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 informationNYS Algebra II and Trigonometry Suggested Sequence of Units (P.I's within each unit are NOT in any suggested order)
1 of 6 UNIT P.I. 1 - INTEGERS 1 A2.A.1 Solve absolute value equations and inequalities involving linear expressions in one variable 1 A2.A.4 * Solve quadratic inequalities in one and two variables, algebraically
More information