MATH 350: Introduction to Computational Mathematics
|
|
- Rolf Montgomery
- 6 years ago
- Views:
Transcription
1 MATH 350: Introduction to Computational Mathematics Chapter VIII: The Fast Fourier Transform Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Spring 2008
2 Outline 1 The (Continuous) Fourier Transform 2 The Discrete Fourier Transform 3 The FFT in MATLAB 4 Examples
3 Outline The (Continuous) Fourier Transform 1 The (Continuous) Fourier Transform 2 The Discrete Fourier Transform 3 The FFT in MATLAB 4 Examples
4 Introduction The (Continuous) Fourier Transform From your differential equations class you know the Laplace transform F(s) = L{f (t)} = 0 e st f (t)dt and perhaps remember some of its nice properties,
5 Introduction The (Continuous) Fourier Transform From your differential equations class you know the Laplace transform F(s) = L{f (t)} = 0 e st f (t)dt and perhaps remember some of its nice properties, such as translation e at F(s a),
6 Introduction The (Continuous) Fourier Transform From your differential equations class you know the Laplace transform F(s) = L{f (t)} = 0 e st f (t)dt and perhaps remember some of its nice properties, such as translation e at F(s a), derivatives t n f (t) ( 1) n d n ds n F(s),
7 Introduction The (Continuous) Fourier Transform From your differential equations class you know the Laplace transform F(s) = L{f (t)} = 0 e st f (t)dt and perhaps remember some of its nice properties, such as translation e at F(s a), derivatives convolution (f g)(t) = t n f (t) ( 1) n d n t 0 ds n F(s), f (τ)g(t τ)dτ F(s)G(s). These properties allow us to convert systems of linear ODEs to systems of linear algebraic equations.
8 The (Continuous) Fourier Transform Another very important integral transform is given by the Fourier transform.
9 The (Continuous) Fourier Transform Another very important integral transform is given by the Fourier transform. The Fourier transform takes a function in the "time-domain" and transforms it into another function its Fourier transform in the "frequency domain".
10 The (Continuous) Fourier Transform Another very important integral transform is given by the Fourier transform. The Fourier transform takes a function in the "time-domain" and transforms it into another function its Fourier transform in the "frequency domain". This terminology is due to the fact that the Fourier transform plays a fundamental role in signal processing.
11 The (Continuous) Fourier Transform Definition The Fourier transform ŷ of a square-integrable function y is defined as: ŷ(ω) = e iωt y(t)dt, < ω <. (1)
12 The (Continuous) Fourier Transform Definition The Fourier transform ŷ of a square-integrable function y is defined as: ŷ(ω) = e iωt y(t)dt, < ω <. (1) On the other hand, the inverse Fourier transform lets us reconstruct y from its Fourier transform ŷ: y(t) = 1 2π e iωt ŷ(ω)dω, < t <. (2)
13 The (Continuous) Fourier Transform Definition The Fourier transform ŷ of a square-integrable function y is defined as: ŷ(ω) = e iωt y(t)dt, < ω <. (1) On the other hand, the inverse Fourier transform lets us reconstruct y from its Fourier transform ŷ: Remark y(t) = 1 2π e iωt ŷ(ω)dω, < t <. (2) There are other (equivalent) ways to define the Fourier transform: ŷ(ω) = 1 2π e iωt y(t)dt and y(t) = 1 2π eiωt ŷ(ω)dω ŷ(f ) = e 2πift y(t)dt and y(f ) = e2πift ŷ(f )df
14 The (Continuous) Fourier Transform The Fourier transform also satisfies many nice properties that make it particularly useful in the context of signal processing.
15 The (Continuous) Fourier Transform The Fourier transform also satisfies many nice properties that make it particularly useful in the context of signal processing. For example, convolution (x y)(t) = x(τ)y(t τ)dτ ˆx(ω)ŷ(ω),
16 The (Continuous) Fourier Transform The Fourier transform also satisfies many nice properties that make it particularly useful in the context of signal processing. For example, convolution (x y)(t) = x(τ)y(t τ)dτ ˆx(ω)ŷ(ω), modulation e iω 0t y(t) ŷ(ω ω 0 ),
17 The (Continuous) Fourier Transform The Fourier transform also satisfies many nice properties that make it particularly useful in the context of signal processing. For example, convolution (x y)(t) = x(τ)y(t τ)dτ ˆx(ω)ŷ(ω), modulation e iω 0t y(t) ŷ(ω ω 0 ), scaling y(at) 1 ( ω ). a ŷ a
18 The (Continuous) Fourier Transform Example Consider the function y(t) = and compute its Fourier transform. { 1, if 1 2 t 1 2 0, otherwise,
19 The (Continuous) Fourier Transform Example Consider the function y(t) = and compute its Fourier transform. Solution { 1, if 1 2 t 1 2 0, otherwise, The definition of the Fourier transform (1) and the definition of y imply ŷ(ω) = e iωt y(t)dt
20 The (Continuous) Fourier Transform Example Consider the function y(t) = and compute its Fourier transform. Solution { 1, if 1 2 t 1 2 0, otherwise, The definition of the Fourier transform (1) and the definition of y imply ŷ(ω) = = e iωt y(t)dt e iωt dt
21 The (Continuous) Fourier Transform Solution (cont.) Now, Euler s formula and the fact that cosine is even and sine is odd give ŷ(ω) = [cos(ωt) i sin(ωt)] dt
22 The (Continuous) Fourier Transform Solution (cont.) Now, Euler s formula and the fact that cosine is even and sine is odd give ŷ(ω) = = [cos(ωt) i sin(ωt)] dt cos(ωt)dt
23 The (Continuous) Fourier Transform Solution (cont.) Now, Euler s formula and the fact that cosine is even and sine is odd give ŷ(ω) = = = 2 sin(ωt) ω [cos(ωt) i sin(ωt)] dt cos(ωt)dt 1 2 0
24 The (Continuous) Fourier Transform Solution (cont.) Now, Euler s formula and the fact that cosine is even and sine is odd give ŷ(ω) = = = 2 sin(ωt) ω [cos(ωt) i sin(ωt)] dt cos(ωt)dt = sin ω 2 ω 2.
25 The (Continuous) Fourier Transform Solution (cont.) Now, Euler s formula and the fact that cosine is even and sine is odd give ŷ(ω) = = = 2 sin(ωt) ω [cos(ωt) i sin(ωt)] dt cos(ωt)dt = sin ω 2 ω 2 These functions play an important role in many applications (e.g., signal processing). The function y is known as a square pulse or characteristic function of the interval [ 1 2, 1 2 ]..
26 The (Continuous) Fourier Transform Solution (cont.) Now, Euler s formula and the fact that cosine is even and sine is odd give ŷ(ω) = = = 2 sin(ωt) ω [cos(ωt) i sin(ωt)] dt cos(ωt)dt = sin ω 2 ω 2 These functions play an important role in many applications (e.g., signal processing). The function y is known as a square pulse or characteristic function of the interval [ 1 2, 1 2 ]. Its Fourier transform ŷ is known as the sinc function..
27 The (Continuous) Fourier Transform
28 The (Continuous) Fourier Transform Other information about the continuous Fourier transform and its connection to Fourier series is provided in Chapter 8.9 of [NCM].
29 Outline The Discrete Fourier Transform 1 The (Continuous) Fourier Transform 2 The Discrete Fourier Transform 3 The FFT in MATLAB 4 Examples
30 The Discrete Fourier Transform In many applications the data is given in discrete form (for example, digital audio for MP3, CD or DVD players).
31 The Discrete Fourier Transform In many applications the data is given in discrete form (for example, digital audio for MP3, CD or DVD players). Therefore, we now look at the Fourier transform of the discrete and periodic data y = [y 0,..., y n 1 ] T with y j = y( 2jπ n ), j = 0,..., n 1. Here time is specified by the n-vector t = [0, 2π n, π n,..., 2(n 1)π n ] T with entries t j = 2jπ n.
32 The Discrete Fourier Transform In many applications the data is given in discrete form (for example, digital audio for MP3, CD or DVD players). Therefore, we now look at the Fourier transform of the discrete and periodic data y = [y 0,..., y n 1 ] T with y j = y( 2jπ n ), j = 0,..., n 1. Here time is specified by the n-vector t = [0, 2π n, π n,..., 2(n 1)π n ] T with entries t j = 2jπ n. We can think of the data y j being sampled from the continuous function y at the time instances t j.
33 The Discrete Fourier Transform Definition The discrete Fourier transform (DFT) is given by n 1 n 1 ŷ k = e ikt j y j = ωn jk y j, k = 0,..., n 1, (3) j=0 j=0 where the numbers ω n = e 2πi n are so-called n-th roots of unity which lie equally spaced on the unit circle.
34 The Discrete Fourier Transform Definition The discrete Fourier transform (DFT) is given by n 1 n 1 ŷ k = e ikt j y j = ωn jk y j, k = 0,..., n 1, (3) j=0 j=0 where the numbers ω n = e 2πi n are so-called n-th roots of unity which lie equally spaced on the unit circle. Remark The MATLAB program fftmatrix.m) from [NCM] not only plots the roots of unity, but I creates a graph connecting the elements of each column of the Fourier matrix F of order n. If n is prime, we get a complete graph, otherwise a subgraph. The sparsity of the graph serves as an indicator for the speed of the FFT.
35 The Discrete Fourier Transform Definition The formula for the inverse discrete Fourier transform (inverse DFT) is given by y j = 1 n 1 n 1 e ikt j ŷ n k = ω jkŷ n k, j = 1,..., n (4) k=0 k=0 with ω n = e 2πi n, the complex conjugate of ω n (also n-th roots of unity).
36 The Discrete Fourier Transform Definition The formula for the inverse discrete Fourier transform (inverse DFT) is given by y j = 1 n 1 n 1 e ikt j ŷ n k = ω jkŷ n k, j = 1,..., n (4) k=0 k=0 with ω n = e 2πi n, the complex conjugate of ω n (also n-th roots of unity). Remark Since the data are periodic the Fourier domain will be discrete as well (since only waves e ikt with integer wavenumber k have period 2π).
37 The Discrete Fourier Transform Since n 1 ŷ k = ωn jk y j j=0 we can express the DFT ŷ in matrix-vector form ŷ = Fy, where the Fourier matrix F has entries F kj = ωn jk.
38 The Discrete Fourier Transform Since n 1 ŷ k = ωn jk y j j=0 we can express the DFT ŷ in matrix-vector form ŷ = Fy, where the Fourier matrix F has entries F kj = ω jk n. The definition of the inverse DFT implies y = 1 n F H ŷ with the Hermitian conjugate F H of F, or where I is the n n identity matrix. F H F = ni,
39 Outline The FFT in MATLAB 1 The (Continuous) Fourier Transform 2 The Discrete Fourier Transform 3 The FFT in MATLAB 4 Examples
40 The FFT in MATLAB Since the Fourier matrix has entries F kj = ω jk n MATLAB via we can generate F in omega = exp(-2*pi*i/n); j = 0:n-1; k = j ; F = omega.^(k*j)
41 The FFT in MATLAB Since the Fourier matrix has entries F kj = ω jk n MATLAB via we can generate F in omega = exp(-2*pi*i/n); j = 0:n-1; k = j ; F = omega.^(k*j) Alternatively, we can just use F = fft(eye(n))
42 Outline Examples 1 The (Continuous) Fourier Transform 2 The Discrete Fourier Transform 3 The FFT in MATLAB 4 Examples
43 Examples touchtone.m
44 Examples fftgui.m displays plots of the elements of the four vectors real(y), imag(y), real(fft(y)), imag(fft(y)).
45 Examples fftgui.m displays plots of the elements of the four vectors real(y), imag(y), real(fft(y)), imag(fft(y)). Example Setting only y 0 0 should give a constant (real) Fourier transform since n 1 ŷ k = e ikt j y j j=0
46 Examples fftgui.m displays plots of the elements of the four vectors Example real(y), imag(y), real(fft(y)), imag(fft(y)). Setting only y 0 0 should give a constant (real) Fourier transform since n 1 ŷ k = e ikt j onlyy 0 0 y j = e 0 y 0. j=0
47 Examples Example Setting only y 1 0 should give ŷ k = n 1 e ikt j y j j=0
48 Examples Example Setting only y 1 0 should give ŷ k = n 1 e ikt j y j j=0 onlyy 0 0 = e 2πik/n y 1
49 Examples Example Setting only y 1 0 should give ŷ k = n 1 e ikt j y j j=0 onlyy 0 0 = e 2πik/n y 1 Euler = [cos(2kπ/n) i sin(2kπ/n)] y 0.
50 Examples SpectralDiffDemo.m dftfilter.m
51 References I Appendix References W. L. Briggs and V. E. Henson. The DFT: An Owners Manual for the Discrete Fourier Transform. SIAM, Philadelphia, C. Moler. Numerical Computing with MATLAB. SIAM, Philadelphia, Also M. Frigo and S. G. Johnson. FFTW: The fastest Fourier transform in the West (C library),
If s is a character that labels one of the buttons on the key pad, the corresponding row index k and column index j can be found with
Chapter 8 Fourier Analysis We all use Fourier analysis every day without even knowing it. Cell phones, disc drives, DVDs and JPEGs all involve fast finite Fourier transforms. This chapter discusses both
More informationFourier series. Complex Fourier series. Positive and negative frequencies. Fourier series demonstration. Notes. Notes. Notes.
Fourier series Fourier series of a periodic function f (t) with period T and corresponding angular frequency ω /T : f (t) a 0 + (a n cos(nωt) + b n sin(nωt)), n1 Fourier series is a linear sum of cosine
More informationFourier Series and Transforms. Revision Lecture
E. (5-6) : / 3 Periodic signals can be written as a sum of sine and cosine waves: u(t) u(t) = a + n= (a ncosπnft+b n sinπnft) T = + T/3 T/ T +.65sin(πFt) -.6sin(πFt) +.6sin(πFt) + -.3cos(πFt) + T/ Fundamental
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 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 informationChapter 2. Signals. Static and Dynamic Characteristics of Signals. Signals classified as
Chapter 2 Static and Dynamic Characteristics of Signals Signals Signals classified as. Analog continuous in time and takes on any magnitude in range of operations 2. Discrete Time measuring a continuous
More informationFrequency Response and Continuous-time Fourier Series
Frequency Response and Continuous-time Fourier Series Recall course objectives Main Course Objective: Fundamentals of systems/signals interaction (we d like to understand how systems transform or affect
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 34. Fourier Transforms
Lecture 34 Fourier Transforms In this section, we introduce the Fourier transform, a method of analyzing the frequency content of functions that are no longer τ-periodic, but which are defined over the
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 informationDESIGN OF CMOS ANALOG INTEGRATED CIRCUITS
DESIGN OF CMOS ANALOG INEGRAED CIRCUIS Franco Maloberti Integrated Microsistems Laboratory University of Pavia Discrete ime Signal Processing F. Maloberti: Design of CMOS Analog Integrated Circuits Discrete
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 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 informationMath 3313: Differential Equations Laplace transforms
Math 3313: Differential Equations Laplace transforms Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX Outline Introduction Inverse Laplace transform Solving ODEs with Laplace
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 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 informationX(t)e 2πi nt t dt + 1 T
HOMEWORK 31 I) Use the Fourier-Euler formulae to show that, if X(t) is T -periodic function which admits a Fourier series decomposition X(t) = n= c n exp (πi n ) T t, then (1) if X(t) is even c n are all
More 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 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 informationChapter 4 The Fourier Series and Fourier Transform
Chapter 4 The Fourier Series and Fourier Transform Fourier Series Representation of Periodic Signals Let x(t) be a CT periodic signal with period T, i.e., xt ( + T) = xt ( ), t R Example: the rectangular
More 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 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 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 information4. Sinusoidal solutions
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More information2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form
2. Oscillation So far, we have used differential equations to describe functions that grow or decay over time. The next most common behavior for a function is to oscillate, meaning that it increases and
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
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 information( ) f (k) = FT (R(x)) = R(k)
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) = R(x)
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 9: Conditionally Positive Definite Radial Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationSystem Identification & Parameter Estimation
System Identification & Parameter Estimation Wb3: SIPE lecture Correlation functions in time & frequency domain Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE // Delft University
More informationMusic 270a: Complex Exponentials and Spectrum Representation
Music 270a: Complex Exponentials and Spectrum Representation Tamara Smyth, trsmyth@ucsd.edu Department of Music, University of California, San Diego (UCSD) October 24, 2016 1 Exponentials The exponential
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 information10. Operators and the Exponential Response Formula
52 10. Operators and the Exponential Response Formula 10.1. Operators. Operators are to functions as functions are to numbers. An operator takes a function, does something to it, and returns this modified
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationMATH 350: Introduction to Computational Mathematics
MATH 350: Introduction to Computational Mathematics Chapter V: Least Squares Problems Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Spring 2011 fasshauer@iit.edu MATH
More informationIntroduction to Fourier Transforms. Lecture 7 ELE 301: Signals and Systems. Fourier Series. Rect Example
Introduction to Fourier ransforms Lecture 7 ELE 3: Signals and Systems Fourier transform as a limit of the Fourier series Inverse Fourier transform: he Fourier integral theorem Prof. Paul Cuff Princeton
More informationDetailed Solutions to Exercises
Detailed Solutions to Exercises Digital Signal Processing Mikael Swartling Nedelko Grbic rev. 205 Department of Electrical and Information Technology Lund University Detailed solution to problem E3.4 A
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 informationThe formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) R(x)
More informationCHEE 319 Tutorial 3 Solutions. 1. Using partial fraction expansions, find the causal function f whose Laplace transform. F (s) F (s) = C 1 s + C 2
CHEE 39 Tutorial 3 Solutions. Using partial fraction expansions, find the causal function f whose Laplace transform is given by: F (s) 0 f(t)e st dt (.) F (s) = s(s+) ; Solution: Note that the polynomial
More informationFourier Analysis and Power Spectral Density
Chapter 4 Fourier Analysis and Power Spectral Density 4. Fourier Series and ransforms Recall Fourier series for periodic functions for x(t + ) = x(t), where x(t) = 2 a + a = 2 a n = 2 b n = 2 n= a n cos
More informationFINAL EXAM MAY 20, 2004
18.034 FINAL EXAM MAY 20, 2004 Name: Problem 1: /10 Problem 2: /20 Problem 3: /25 Problem 4: /15 Problem 5: /20 Problem 6: /25 Problem 7: /10 Problem 8: /35 Problem 9: /40 Problem 10: /10 Extra credit
More informationEECS 20N: Midterm 2 Solutions
EECS 0N: Midterm Solutions (a) The LTI system is not causal because its impulse response isn t zero for all time less than zero. See Figure. Figure : The system s impulse response in (a). (b) Recall that
More informationDiscrete Systems & Z-Transforms. Week Date Lecture Title. 9-Mar Signals as Vectors & Systems as Maps 10-Mar [Signals] 3
http:elec34.org Discrete Systems & Z-Transforms 4 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: eek Date Lecture Title -Mar Introduction
More informationCMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation
CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is
More informationThe Laplace Transform
The Laplace Transform Syllabus ECE 316, Spring 2015 Final Grades Homework (6 problems per week): 25% Exams (midterm and final): 50% (25:25) Random Quiz: 25% Textbook M. Roberts, Signals and Systems, 2nd
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 informationFourier Sin and Cos Series and Least Squares Convergence
Fourier and east Squares Convergence James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University May 7, 28 Outline et s look at the original Fourier sin
More informationEE/ME/AE324: Dynamical Systems. Chapter 7: Transform Solutions of Linear Models
EE/ME/AE324: Dynamical Systems Chapter 7: Transform Solutions of Linear Models The Laplace Transform Converts systems or signals from the real time domain, e.g., functions of the real variable t, to the
More informationMATH 135: COMPLEX NUMBERS
MATH 135: COMPLEX NUMBERS (WINTER, 010) The complex numbers C are important in just about every branch of mathematics. These notes 1 present some basic facts about them. 1. The Complex Plane A complex
More informatione st f (t) dt = e st tf(t) dt = L {t f(t)} s
Additional operational properties How to find the Laplace transform of a function f (t) that is multiplied by a monomial t n, the transform of a special type of integral, and the transform of a periodic
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 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 informationCourse roadmap. ME451: Control Systems. Example of Laplace transform. Lecture 2 Laplace transform. Laplace transform
ME45: Control Systems Lecture 2 Prof. Jongeun Choi Department of Mechanical Engineering Michigan State University Modeling Transfer function Models for systems electrical mechanical electromechanical Block
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 informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
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 informationENSC327 Communications Systems 2: Fourier Representations. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 2: Fourier Representations Jie Liang School of Engineering Science Simon Fraser University 1 Outline Chap 2.1 2.5: Signal Classifications Fourier Transform Dirac Delta Function
More informationGeneralizing the DTFT!
The Transform Generaliing the DTFT! The forward DTFT is defined by X e jω ( ) = x n e jωn in which n= Ω is discrete-time radian frequency, a real variable. The quantity e jωn is then a complex sinusoid
More informationTime and Spatial Series and Transforms
Time and Spatial Series and Transforms Z- and Fourier transforms Gibbs' phenomenon Transforms and linear algebra Wavelet transforms Reading: Sheriff and Geldart, Chapter 15 Z-Transform Consider a digitized
More 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 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 informationOrdinary Differential Equations. Session 7
Ordinary Differential Equations. Session 7 Dr. Marco A Roque Sol 11/16/2018 Laplace Transform Among the tools that are very useful for solving linear differential equations are integral transforms. An
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 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 informationReview of Frequency Domain Fourier Series: Continuous periodic frequency components
Today we will review: Review of Frequency Domain Fourier series why we use it trig form & exponential form how to get coefficients for each form Eigenfunctions what they are how they relate to LTI systems
More informationIdentification Methods for Structural Systems
Prof. Dr. Eleni Chatzi System Stability Fundamentals Overview System Stability Assume given a dynamic system with input u(t) and output x(t). The stability property of a dynamic system can be defined from
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 informationThis is the number of cycles per unit time, and its units are, for example,
16 4. Sinusoidal solutions Many things in nature are periodic, even sinusoidal. We will begin by reviewing terms surrounding periodic functions. If an LTI system is fed a periodic input signal, we have
More informationnatural frequency of the spring/mass system is ω = k/m, and dividing the equation through by m gives
77 6. More on Fourier series 6.. Harmonic response. One of the main uses of Fourier series is to express periodic system responses to general periodic signals. For example, if we drive an undamped spring
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 informationCHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L
CHAPTER 4 FOURIER SERIES 1 S A B A R I N A I S M A I L Outline Introduction of the Fourier series. The properties of the Fourier series. Symmetry consideration Application of the Fourier series to circuit
More 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 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 informationENGIN 211, Engineering Math. Laplace Transforms
ENGIN 211, Engineering Math Laplace Transforms 1 Why Laplace Transform? Laplace transform converts a function in the time domain to its frequency domain. It is a powerful, systematic method in solving
More informationLaplace Transforms and use in Automatic Control
Laplace Transforms and use in Automatic Control P.S. Gandhi Mechanical Engineering IIT Bombay Acknowledgements: P.Santosh Krishna, SYSCON Recap Fourier series Fourier transform: aperiodic Convolution integral
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 informationSolutions to Problems in Chapter 4
Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave
More 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 informationSpecial Mathematics Laplace Transform
Special Mathematics Laplace Transform March 28 ii Nature laughs at the difficulties of integration. Pierre-Simon Laplace 4 Laplace Transform Motivation Properties of the Laplace transform the Laplace transform
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 informationGATE EE Topic wise Questions SIGNALS & SYSTEMS
www.gatehelp.com GATE EE Topic wise Questions YEAR 010 ONE MARK Question. 1 For the system /( s + 1), the approximate time taken for a step response to reach 98% of the final value is (A) 1 s (B) s (C)
More informationA523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011
A523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011 Lecture 6 PDFs for Lecture 1-5 are on the web page Problem set 2 is on the web page Article on web page A Guided
More informationD.S.G. POLLOCK: BRIEF NOTES
BIVARIATE SPECTRAL ANALYSIS Let x(t) and y(t) be two stationary stochastic processes with E{x(t)} = E{y(t)} =. These processes have the following spectral representations: (1) x(t) = y(t) = {cos(ωt)da
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 informationDefinitions of the Laplace Transform (1A) Young Won Lim 1/31/15
Definitions of the Laplace Transform (A) Copyright (c) 24 Young W. Lim. Permission is granted to copy, distriute and/or modify this document under the terms of the GNU Free Documentation License, Version.2
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationTHE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3
THE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3 Any periodic function f(t) can be written as a Fourier Series a 0 2 + a n cos( nωt) + b n sin n
More 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 informationBackground ODEs (2A) Young Won Lim 3/7/15
Background ODEs (2A) Copyright (c) 2014-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any
More informationSystems Analysis and Control
Systems Analysis and Control Matthew M. Peet Arizona State University Lecture 5: Calculating the Laplace Transform of a Signal Introduction In this Lecture, you will learn: Laplace Transform of Simple
More informationHarman Outline 1A CENG 5131
Harman Outline 1A CENG 5131 Numbers Real and Imaginary PDF In Chapter 2, concentrate on 2.2 (MATLAB Numbers), 2.3 (Complex Numbers). A. On R, the distance of any real number from the origin is the magnitude,
More informationA1 Time-Frequency Analysis
A 20 / A Time-Frequency Analysis David Murray david.murray@eng.ox.ac.uk www.robots.ox.ac.uk/ dwm/courses/2tf Hilary 20 A 20 2 / Content 8 Lectures: 6 Topics... From Signals to Complex Fourier Series 2
More informationa k cos kω 0 t + b k sin kω 0 t (1) k=1
MOAC worksheet Fourier series, Fourier transform, & Sampling Working through the following exercises you will glean a quick overview/review of a few essential ideas that you will need in the moac course.
More 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 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 informationENGIN 211, Engineering Math. Fourier Series and Transform
ENGIN 11, Engineering Math Fourier Series and ransform 1 Periodic Functions and Harmonics f(t) Period: a a+ t Frequency: f = 1 Angular velocity (or angular frequency): ω = ππ = π Such a periodic function
More 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 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 informationLecture 7 ELE 301: Signals and Systems
Lecture 7 ELE 30: Signals and Systems Prof. Paul Cuff Princeton University Fall 20-2 Cuff (Lecture 7) ELE 30: Signals and Systems Fall 20-2 / 22 Introduction to Fourier Transforms Fourier transform as
More information