Lecture 7 January 26, 2016
|
|
- Bethanie Bishop
- 5 years ago
- Views:
Transcription
1 MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 7 January 26, 26 Prof Emmanuel Candes Scribe: Carlos A Sing-Long, Edited by E Bates Outline Agenda: Discrete Fourier transform umerical accuracy of the trapezoidal rule 2 Fourier transform of finite signals 3 Introduction to Fast Fourier Transforms Last Time: We defined discrete convolutions for signals defined on Z Following a similar approach to that used for the continuous-time Fourier transform, we deduced that time-invariant maps on discrete signals are precisely the maps that can be expressed as convolutions by a fixed function We also showed that complex exponentials with again diagonalize these operators In turn, this motivated the definition of the Fourier series of a discrete signal, which is a periodic function The natural question is, is every periodic function a Fourier series? For square-integrable functions on the circle S, the answer is yes We showed that complex exponentials are a complete orthonormal set in this space, which immediately yields expressions for the inverse transform, as well as the Parseval-Plancherel theorem Furthermore, we saw why Shannon s sampling theorem can be interpreted as a restatement of the completeness of this orthonormal sequence 2 umerical accuracy of the trapezoidal rule In many applications one wants to estimate the integral of a function over a finite interval Under proper normalizations, we can state this problem as estimating the integral I f(t) dt This is typically done by using quadrature rules, that is, by using a finite sequence of nodes {t n } n [, ] and weights {w n} n such that I n w n f(t n ),
2 is a good approximation of I The numerical accuracy of the quadrature rule is characterized by the behavior of the error I I as increases One would like to choose nodes and weights so that the error decays as quickly as possible with The simplest method of quadrature is the trapezoidal rule, which uses uniformly spaced nodes with (almost) uniform weights: I ( ( ) ( ) 2 f() + f + + f + 2 ) f() Despite its simplicity, it is extraordinarily accurate when f is smooth and periodic (for instance, see Fig ) Let us see why this is the case Working on [, ], assume f is -periodic In particular, f() f() so that I n f ( n ) How close is I to I? To learn the answer, we express f as a linear combination of the basis functions we identified in Lecture 6: Let e k (t) e i2πkt and write f(t) k Z e k, f e k (t), which yields I e k, f e i2πk n n k Z k Z e k, f ( n e i2πk n ) ow, e i2πk/ is an th root of unity, equal to if and only if k is a multiple of We shall then use the following fact: If ξ is an th root of unity, then Indeed, ξ is a root of x, and so it is also a root of Consequently, We then have + ξ + ξ ξ () x x + x + x2 + + x { e i2πk n k is a multiple of otherwise n I ( ) k, f δ(k m) k Z e m Z e, f + I + m Z, m m Z, m 2 e m, f e m, f
3 If the Fourier coefficients of f can be computed, we could even have an explicit expression for the remainder I I e m, f m Z, m To determine how fast it decays as a function of, consider f that is p-times differentiable and -periodic In this case, integration by parts tells us e k, f f(t)e i2πkt dt i2πk f (t)e i2πkt dt (i2πk) p f (p) (t)e i2πkt dt The periodicity is necessary for the boundary terms to vanish during integration by parts We now see that if p is large and f (p) is small, then this Fourier coefficient will be negligible More precisely, f(t)e i2πkt dt f (p) 2πk p, In this case, I I m Z, m f (p) 2πm p 2 f (p) (2π) p m> m p 2 f (p) (2π) p ζ(p), where ζ(p) tends to as p ow, when f is smooth, this bound holds for every p, meaning the error decays with faster than p for any p Such behavior is called spectral accuracy We noted that if f is smooth but not periodic (ie f() f()), then the argument fails as f(t)e i2πkt f() f() dt + i2πk i2πk f (t)e i2πkt dt And the failure is seemingly quite dramatic, since now the argument from above would stop with p, giving O( ) error The reality is only slightly better: standard numerical analysis would show O( 2 ) error The failure in accuracy can be understood by considering the periodic domain S R/Z of period The Fourier series we are considering are functions on S If f() f(), then the periodic extension of f to all of R has a discontinuity at ( mod Z) So even though f is smooth on (, ), it is not even continuous on S, hence the loss of accuracy One workaround to this problem is to window the function f Again suppose f is smooth on (, ) but not periodic For any δ >, we can construct a smooth function χ such that χ(t) for t [δ, δ] and decays to at the boundaries t and t Then we can write I f(t) dt χ(t)f(t) dt + χ(t)f(t) dt + δ ( χ(t))f(t) dt ( χ(t))f(t) dt + ( χ(t))f(t) dt δ The function χf is smooth and periodic, and consequently the trapezoidal rule will be extremely accurate in estimating the first term For the second and third terms, we can use a more sophisticated quadrature rule Moreover, for small δ, the expense of having a fine mesh (ie large ) on the short intervals [, δ] and [ δ, δ] is comparatively small 3
4 f(t) e cos2 (t) Figure : Here on [, π] instead of [, ], the trapezoidal rule with only 8 points (shown in red) yields an estimate I of impressive accuracy: I I (ie 8 digits of accuracy) 2 An example Consider f(t) sin(πt) over [, ] A direct computation gives We thus have the bound m Z, m e k, f e m, f 4 π sin(πt)e i2πkt dt 2 π m> 4k 2 4(m) 2 π 2 m> m 2 π 6 2, showing the trapezoidal quadrature has O( 2 ) accuracy For instance, see Fig 2 to see that 4 gives a quadrature error on the order 2 evertheless, the same quadrature performs perfectly when f(t) sin(2πt) The difference is that when extended periodically, sin(πt) has a discontinuity in its first derivative while sin(2πt) is smooth 3 Fourier transform of finite signals Here we consider finite signals {f[n]} n, for which we use brackets instead of parentheses to indicate their argument Here is the length of the signal We define the discrete Fourier transform (DFT) of f as ˆf[k] n f[n]e i2πkn, k,, ote that the DFT of f is also a finite signal of the same length We can naturally represent f as a vector f C with f n f[n] and similarly for ˆf Then the DFT is just the matrix multiplication ˆf F f where F kn e i2πkn, k, n,, 4
5 (a) f(t) sin(πt) (b) f(t) sin(2πt) Figure 2: (a) Integrating sin(πt) on [, ] using the trapezoidal rule always underestimates the area With 4 points we obtain I whereas the true value is I π/ We achieve digit of accuracy (b) Integrating sin(2πt) on [, ] using the trapezoidal rule underestimates the area over [, 2 ] and overestimates it over [ 2, ] The two errors exactly cancel each other to give I 4 I We use roots of unity to simplify the notation Let z e i2π Then ˆf k n z kn f n, and, in matrix notation ˆf ˆf z z 2 z ˆf 2 z 2 z ( 2)( 2) z ( 2)() ˆf z z ()( 2) z ()() }{{}}{{} ˆf F f f f 2 f } {{ } f This shows that F is a Vandermonde matrix (a matrix in which each row is a geometric sequence) Furthermore, the DFT (when suitable rescaled) is unitary, meaning F F F F I We verify this directly: Since z z, we compute { m n, l F ml F ln where we have again used () l z ml zln This yields the discrete inversion formula l z (n m)l F F or f F ˆf otherwise, 5
6 In function notation, f[n] k ˆf[k]e i2πkn In what follows we will not make distinction between vectors and finite signals 3 DFT as a quadrature Consider a continuous-time signal f(t) on S The Fourier transform of f is given by ˆf(ω) f(t)e i2πωt dt In practice, this integral may be difficult to compute analytically Suppose we estimated the integral by a left-hand sum with intervals This would give ˆf(ω) n f ( n ) 2πωn i e When ω {,,, }, this estimate is exactly the DFT of the finite discretization f[n] f ( ) n That is, when the finite signal consists of uniformly spaced samples of an analog signal, the DFT is an estimate of the first (nonnegative) Fourier coefficients And we have seen that when f(t) is smooth and periodic, this estimate is astonishingly accurate Furthermore, when the signal is real-valued, the positive-frequency coefficients immediately give the negative-frequency ones: ˆf( k) ˆf(k) 32 Circular convolution For infinite, discrete signals, one has a well-defined convolution (f g)(n) m Z f(m)g(n m) We would like a notion of convolution for finite signals as well However, finite signals {f[n]} n and {g[n]} n are only defined on {,,, }, meaning the sum above would include terms that are not defined Consequently, we extend the signals in the natural way: f[n] : f[n mod ] and g[n] g[n mod ], n Z ow f[n] and g[n] can be regarded as periodic signals on Z We can now define circular convolution as (f g)[n] m f[m]g[n m] m f[m mod ]g[n m mod ] In analogy with the infinite case, the discrete Fourier transform diagonalizes circular convolutions (this is an indication that we have made a good definition of convolution) That is, h g f DFT ĥ ˆfĝ 6
7 Continuing the analogy, we have circular translations and the circular impulse f f τ [n] f[n τ mod ], δ[n] { n mod, otherwise A translation-invariant map is a linear map L satisfying Lf τ (Lf) τ By the same arguments as given in Lectures and 6, a map is translation-invariant if and only if it is a circular convolution In particular, the corresponding convolution kernel is the response to the impulse, h : Lδ Consequently, we have the diagram DFT f[n] ˆf[k] L ĥ[k] (f h)[n] DFT ˆf[k] ĥ[k] If we denote by L the matrix representing L, this result can be interpreted as a factorization of the form L F diag(ĥ)f, where ĥ[] ĥ[] diag(ĥ) ĥ[ 2] ĥ[ ] We thus have a simple method of applying L to a signal f More important, we have a fast method Multiplication by a fixed transfer function ĥ[k] is certainly fast, and we will see in the next section how the DFT and inverse DFT are carried out efficiently 4 Introduction to Fast Fourier Transforms The DFT of a finite signal x of size is given by ˆx[k] n x[n]e i2πkn How do we compute it? The direct approach is to assemble the DFT matrix, which is dense The computational cost of this method is O( 2 ) This has many drawbacks For instance: We can only take DFTs of relatively small signals The problem with this is that one typically wants a large number of samples from a continuous signal to obtain an accurate estimate of its frequency content When the accuracy is proportional to, we are performing a computation that is four times more expensive to be only twice as accurate 7
8 It is not possible to compute the DFT in (near) real-time The problem with this is the DFT allows us to extract information contained on a prescribed frequency band by filtering This is the case in wireless communications, where a delay in extracting the relevant information amounts, for instance, to a significant delay between receiving a signal on a cellphone and actually listening to the message For these reasons and others, having an efficient algorithm to compute the DFT is extremely important The de facto algorithm to compute the DFT is the Fast Fourier Transform (FFT) This algorithm was conceived by Gauss (Gauss, circa 85) and rediscovered and popularized by Cooley & Tukey (Cooley & Tukey, 965) It allows us to compute the DFT of a signal of length in O( log ) operations Here we will discuss the radix-2 algorithm which assumes is a power of two The main idea behind it is the following When is even, a direct computation shows that ˆx[k] n /2 n /2 n x[n]e i2π kn /2 x[2n]e i2π k(2n) + n x[2n]e i2π /2 kn + e i2π x[2n + ]e i2π k(2n+) /2 k n x[2n + ]e i2π /2 kn Consequently, we can decompose the DFT of a signal of size into two DFTs of signals of size /2 To formalize this idea, let z e i2π, be the th root of unity and define the even and odd parts of x as x even [k] x[2k] and x odd [k] x[2k + ], for k,, /2, that is, the even and odd parts are signals of half the length of the original With this notation, the expression above can be written as {ˆxeven [k] + z k ˆx[k] ˆx odd[k], k /2, ˆx even [k /2] z k /2 ˆx odd [k /2], /2 k, and the DFT of a signal of size can be decomposed into two DFTs of signals of size /2 Furthermore, if we knew ˆx even and ˆx odd, then the above shows how to reconstruct ˆx Suppose that T is the cost of computing a DFT of a signal of size Then the above shows T 2T /2 }{{} computing ˆx even and ˆx odd + /2 }{{} computing z k for k,, /2 + 3/2 }{{} forming the upper and lower part of ˆx For a historical account we refer to Heideman, M, Johnson, D H, Burrus, C S, Gauss and the history of the Fast Fourier Transform, IEEE ASSP Magazine, (4), 4-2 (984) 8
9 If 2 M for some M > then or T 2 M 2T 2 M + 2 M+ T 2 M 2 M T + 2M2 M, T T + 2 log 2, which asymptotically yields a complexity of O( log ) ote the FFT can be used to compute the inverse DFT efficiently too This is because x[n] k ˆx[k]e i2πkn k ˆx[ k]e i2π ( k)n e i2π n k ˆx[ k]e i2πkn, and consequently we only need to take the DFT of a reversed ˆx and then multiply by the factor above It is clear this has complexity O( log ) too In general, does not need to be a power of two, although FFTs are most efficient when it is Many techniques have been developed to reduce the DFT of a signal of size to DFTs of signals of smaller size, for instance, by attempting to use the prime factors of to achieve such a decomposition References [] Cooley, J W, Tukey, J W, An algorithm for the machine calculation of complex Fourier series, Mathematics of Computation, 9(2), 297-3, 965 9
Lecture 6 January 21, 2016
MATH 6/CME 37: Applied Fourier Analysis and Winter 06 Elements of Modern Signal Processing Lecture 6 January, 06 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda: Fourier
More informationLecture 1 January 5, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture January 5, 26 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long; Edited by E. Candes & E. Bates Outline
More informationLecture 9 February 2, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 9 February 2, 26 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda:
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 Virtually all practical signals have finite length (e.g., sensor data, audio records, digital images, stock values, etc). Rather than considering such signals to be zero-padded
More 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 informationDiscrete Fourier transform (DFT)
Discrete Fourier transform (DFT) Signal Processing 2008/9 LEA Instituto Superior Técnico Signal Processing LEA (IST) Discrete Fourier transform 1 / 34 Periodic signals Consider a periodic signal x[n] with
More informationDigital Signal Processing. Midterm 2 Solutions
EE 123 University of California, Berkeley Anant Sahai arch 15, 2007 Digital Signal Processing Instructions idterm 2 Solutions Total time allowed for the exam is 80 minutes Please write your name and SID
More informationThe discrete and fast Fourier transforms
The discrete and fast Fourier transforms Marcel Oliver Revised April 7, 1 1 Fourier series We begin by recalling the familiar definition of the Fourier series. For a periodic function u: [, π] C, we define
More informationChirp Transform for FFT
Chirp Transform for FFT Since the FFT is an implementation of the DFT, it provides a frequency resolution of 2π/N, where N is the length of the input sequence. If this resolution is not sufficient in a
More informationMA3232 Numerical Analysis Week 9. James Cooley (1926-)
MA umerical Analysis Week 9 James Cooley (96-) James Cooley is an American mathematician. His most significant contribution to the world of mathematics and digital signal processing is the Fast Fourier
More informationDSP Algorithm Original PowerPoint slides prepared by S. K. Mitra
Chapter 11 DSP Algorithm Implementations 清大電機系林嘉文 cwlin@ee.nthu.edu.tw Original PowerPoint slides prepared by S. K. Mitra 03-5731152 11-1 Matrix Representation of Digital Consider Filter Structures This
More 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 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 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 informationThe Fourier transform allows an arbitrary function to be represented in terms of simple sinusoids. The Fourier transform (FT) of a function f(t) is
1 Introduction Here is something I wrote many years ago while working on the design of anemometers for measuring shear stresses. Part of this work required modelling and compensating for the transfer function
More informationMultidimensional digital signal processing
PSfrag replacements Two-dimensional discrete signals N 1 A 2-D discrete signal (also N called a sequence or array) is a function 2 defined over thex(n set 1 of, n 2 ordered ) pairs of integers: y(nx 1,
More informationAn Introduction and Analysis of the Fast Fourier Transform. Thao Nguyen Mentor: Professor Ron Buckmire
An Introduction and Analysis of the Fast Fourier Transform Thao Nguyen Mentor: Professor Ron Buckmire Outline 1. Outline of the Talk 2. Continuous Fourier Transform (CFT) Visualization of the CFT Mathematical
More informationVII. Discrete Fourier Transform (DFT) Chapter-8. A. Modulo Arithmetic. (n) N is n modulo N, n is an integer variable.
1 VII. Discrete Fourier Transform (DFT) Chapter-8 A. Modulo Arithmetic (n) N is n modulo N, n is an integer variable. (n) N = n m N 0 n m N N-1, pick m Ex. (k) 4 W N = e -j2π/n 2 Note that W N k = 0 but
More informationFourier Analysis for Engineers
Fourier Analysis for Engineers M. Behrens January 12, 22 1 Motivation: Finite dimensional vector spaces The periodic theory of Fourier analysis may be likened to the decomposition of an ordinary vector
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 informationEE123 Digital Signal Processing
Announcements EE Digital Signal Processing otes posted HW due Friday SDR give away Today! Read Ch 9 $$$ give me your names Lecture based on slides by JM Kahn M Lustig, EECS UC Berkeley M Lustig, EECS UC
More 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 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 informationEE 225D LECTURE ON DIGITAL FILTERS. University of California Berkeley
University of California Berkeley College of Engineering Department of Electrical Engineering and Computer Sciences Professors : N.Morgan / B.Gold EE225D Digital Filters Spring,1999 Lecture 7 N.MORGAN
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 informationDigital Signal Processing Lecture 10 - Discrete Fourier Transform
Digital Signal Processing - Discrete Fourier Transform Electrical Engineering and Computer Science University of Tennessee, Knoxville November 12, 2015 Overview 1 2 3 4 Review - 1 Introduction Discrete-time
More informationLecture 11 FIR Filters
Lecture 11 FIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/4/12 1 The Unit Impulse Sequence Any sequence can be represented in this way. The equation is true if k ranges
More informationLecture 16 February 25, 2016
MTH 262/CME 372: pplied Fourier nalysis and Winter 2016 Elements of Modern Signal Processing Lecture 16 February 25, 2016 Prof. Emmanuel Candes Scribe: Carlos. Sing-Long, Edited by E. Bates 1 Outline genda:
More informationFrequency-domain representation of discrete-time signals
4 Frequency-domain representation of discrete-time signals So far we have been looing at signals as a function of time or an index in time. Just lie continuous-time signals, we can view a time signal as
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 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 informationLecture 10. Digital Signal Processing. Chapter 7. Discrete Fourier transform DFT. Mikael Swartling Nedelko Grbic Bengt Mandersson. rev.
Lecture 10 Digital Signal Processing Chapter 7 Discrete Fourier transform DFT Mikael Swartling Nedelko Grbic Bengt Mandersson rev. 016 Department of Electrical and Information Technology Lund University
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 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 informationHST.582J / 6.555J / J Biomedical Signal and Image Processing Spring 2007
MIT OpenCourseare http://ocw.mit.edu HST.58J / 6.555J / 16.56J Biomedical Signal and Image Processing Spring 7 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationCh.11 The Discrete-Time Fourier Transform (DTFT)
EE2S11 Signals and Systems, part 2 Ch.11 The Discrete-Time Fourier Transform (DTFT Contents definition of the DTFT relation to the -transform, region of convergence, stability frequency plots convolution
More 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 informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture 18: Applications of FFT Algorithms & Linear Filtering DFT Computation; Implementation of Discrete Time Systems Kenneth E. Barner Department of Electrical and
More information2 The Fourier Transform
2 The Fourier Transform The Fourier transform decomposes a signal defined on an infinite time interval into a λ-frequency component, where λ can be any real(or even complex) number. 2.1 Informal Development
More informationAn Introduction to Filterbank Frames
An Introduction to Filterbank Frames Brody Dylan Johnson St. Louis University October 19, 2010 Brody Dylan Johnson (St. Louis University) An Introduction to Filterbank Frames October 19, 2010 1 / 34 Overview
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 information6.003: Signal Processing
6.003: Signal Processing Discrete Fourier Transform Discrete Fourier Transform (DFT) Relations to Discrete-Time Fourier Transform (DTFT) Relations to Discrete-Time Fourier Series (DTFS) October 16, 2018
More informationThe Discrete Fourier transform
453.70 Linear Systems, S.M. Tan, The University of uckland 9- Chapter 9 The Discrete Fourier transform 9. DeÞnition When computing spectra on a computer it is not possible to carry out the integrals involved
More informationLecture 20: Discrete Fourier Transform and FFT
EE518 Digital Signal Processing University of Washington Autumn 2001 Dept of Electrical Engineering Lecture 20: Discrete Fourier Transform and FFT Dec 10, 2001 Prof: J Bilmes TA:
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationThe Fast Fourier Transform. Andreas Klappenecker
The Fast Fourier Transform Andreas Klappenecker Motivation There are few algorithms that had more impact on modern society than the fast Fourier transform and its relatives. The applications of the fast
More informationEE-210. Signals and Systems Homework 7 Solutions
EE-20. Signals and Systems Homework 7 Solutions Spring 200 Exercise Due Date th May. Problems Q Let H be the causal system described by the difference equation w[n] = 7 w[n ] 2 2 w[n 2] + x[n ] x[n 2]
More informationSection 6.6 Gaussian Quadrature
Section 6.6 Gaussian Quadrature Key Terms: Method of undetermined coefficients Nonlinear systems Gaussian quadrature Error Legendre polynomials Inner product Adapted from http://pathfinder.scar.utoronto.ca/~dyer/csca57/book_p/node44.html
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 informationHow to manipulate Frequencies in Discrete-time Domain? Two Main Approaches
How to manipulate Frequencies in Discrete-time Domain? Two Main Approaches Difference Equations (an LTI system) x[n]: input, y[n]: output That is, building a system that maes use of the current and previous
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 informationDigital Image Processing Lectures 13 & 14
Lectures 13 & 14, Professor Department of Electrical and Computer Engineering Colorado State University Spring 2013 Properties of KL Transform The KL transform has many desirable properties which makes
More informationDefinition. A signal is a sequence of numbers. sequence is also referred to as being in l 1 (Z), or just in l 1. A sequence {x(n)} satisfying
Signals and Systems. Definition. A signal is a sequence of numbers {x(n)} n Z satisfying n Z x(n)
More informationIntroduction to Signals and Systems Lecture #4 - Input-output Representation of LTI Systems Guillaume Drion Academic year
Introduction to Signals and Systems Lecture #4 - Input-output Representation of LTI Systems Guillaume Drion Academic year 2017-2018 1 Outline Systems modeling: input/output approach of LTI systems. Convolution
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione dei Segnali Multi-dimensionali e
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationThe Frequency Domain
Mathematical Tools for Data Science Spring 9 The Frequency Domain Motivation In data science we call any structured object of interest an image, an audio sequence, a video... a signal. Signals are usually
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 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 information! Introduction. ! Discrete Time Signals & Systems. ! Z-Transform. ! Inverse Z-Transform. ! Sampling of Continuous Time Signals
ESE 531: Digital Signal Processing Lec 25: April 24, 2018 Review Course Content! Introduction! Discrete Time Signals & Systems! Discrete Time Fourier Transform! Z-Transform! Inverse Z-Transform! Sampling
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 informationVALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur- 603 203 DEPARTMENT OF INFORMATION TECHNOLOGY Academic Year 2016-2017 QUESTION BANK-ODD SEMESTER NAME OF THE SUBJECT SUBJECT CODE SEMESTER YEAR
More informationComputational Methods for Astrophysics: Fourier Transforms
Computational Methods for Astrophysics: Fourier Transforms John T. Whelan (filling in for Joshua Faber) April 27, 2011 John T. Whelan April 27, 2011 Fourier Transforms 1/13 Fourier Analysis Outline: Fourier
More informationHow to Write Fast Numerical Code
How to Write Fast Numerical Code Lecture: Discrete Fourier transform, fast Fourier transforms Instructor: Markus Püschel TA: Georg Ofenbeck & Daniele Spampinato Rest of Semester Today Lecture Project meetings
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 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 information[ ], [ ] [ ] [ ] = [ ] [ ] [ ]{ [ 1] [ 2]
4. he discrete Fourier transform (DF). Application goal We study the discrete Fourier transform (DF) and its applications: spectral analysis and linear operations as convolution and correlation. We use
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 informationESS Finite Impulse Response Filters and the Z-transform
9. Finite Impulse Response Filters and the Z-transform We are going to have two lectures on filters you can find much more material in Bob Crosson s notes. In the first lecture we will focus on some of
More informationMultirate signal processing
Multirate signal processing Discrete-time systems with different sampling rates at various parts of the system are called multirate systems. The need for such systems arises in many applications, including
More informationModeling and Analysis of Systems Lecture #3 - Linear, Time-Invariant (LTI) Systems. Guillaume Drion Academic year
Modeling and Analysis of Systems Lecture #3 - Linear, Time-Invariant (LTI) Systems Guillaume Drion Academic year 2015-2016 1 Outline Systems modeling: input/output approach and LTI systems. Convolution
More informationLecture 19: Discrete Fourier Series
EE518 Digital Signal Processing University of Washington Autumn 2001 Dept. of Electrical Engineering Lecture 19: Discrete Fourier Series Dec 5, 2001 Prof: J. Bilmes TA: Mingzhou
More informationWavelets and applications
Chapter 3 Wavelets and applications 3. Multiresolution analysis 3.. The limits of Fourier analysis Although along this chapter the underlying Hilbert space will be L 2 (R), we start with a completely explicit
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 informationMAHALAKSHMI ENGINEERING COLLEGE-TRICHY
DIGITAL SIGNAL PROCESSING DEPT./SEM.: ECE&EEE /V DISCRETE FOURIER TRANFORM AND FFT PART-A 1. Define DFT of a discrete time sequence? AUC MAY 06 The DFT is used to convert a finite discrete time sequence
More informationDiscrete Time Signals and Systems Time-frequency Analysis. Gloria Menegaz
Discrete Time Signals and Systems Time-frequency Analysis Gloria Menegaz Time-frequency Analysis Fourier transform (1D and 2D) Reference textbook: Discrete time signal processing, A.W. Oppenheim and R.W.
More informationFourier Reconstruction from Non-Uniform Spectral Data
School of Electrical, Computer and Energy Engineering, Arizona State University aditya.v@asu.edu With Profs. Anne Gelb, Doug Cochran and Rosemary Renaut Research supported in part by National Science Foundation
More informationLinear Convolution Using FFT
Linear Convolution Using FFT Another useful property is that we can perform circular convolution and see how many points remain the same as those of linear convolution. When P < L and an L-point circular
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 informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
More informationMultiresolution schemes
Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Multi-dimensional signal processing Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione
More informationTopics in Fourier analysis - Lecture 2.
Topics in Fourier analysis - Lecture 2. Akos Magyar 1 Infinite Fourier series. In this section we develop the basic theory of Fourier series of periodic functions of one variable, but only to the extent
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationSNR Calculation and Spectral Estimation [S&T Appendix A]
SR Calculation and Spectral Estimation [S&T Appendix A] or, How not to make a mess of an FFT Make sure the input is located in an FFT bin 1 Window the data! A Hann window works well. Compute the FFT 3
More informationInterchange of Filtering and Downsampling/Upsampling
Interchange of Filtering and Downsampling/Upsampling Downsampling and upsampling are linear systems, but not LTI systems. They cannot be implemented by difference equations, and so we cannot apply z-transform
More information8 The Discrete Fourier Transform (DFT)
8 The Discrete Fourier Transform (DFT) ² Discrete-Time Fourier Transform and Z-transform are de ned over in niteduration sequence. Both transforms are functions of continuous variables (ω and z). For nite-duration
More informationGEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING
GEORGIA INSIUE OF ECHNOLOGY SCHOOL of ELECRICAL and COMPUER ENGINEERING ECE 6250 Spring 207 Problem Set # his assignment is due at the beginning of class on Wednesday, January 25 Assigned: 6-Jan-7 Due
More informationDHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A
DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A Classification of systems : Continuous and Discrete
More informationFast Polynomials Multiplication Using FFT
Li Chen lichen.xd at gmail.com Xidian University January 17, 2014 Outline 1 Discrete Fourier Transform (DFT) 2 Discrete Convolution 3 Fast Fourier Transform (FFT) 4 umber Theoretic Transform (TT) 5 More
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationDiscrete Fourier Transform
Discrete Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013 Mårten Björkman
More 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 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 informationDesigning Information Devices and Systems I Spring 2016 Official Lecture Notes Note 21
EECS 6A Designing Information Devices and Systems I Spring 26 Official Lecture Notes Note 2 Introduction In this lecture note, we will introduce the last topics of this semester, change of basis and diagonalization.
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 informationDISCRETE FOURIER TRANSFORM
DD2423 Image Processing and Computer Vision DISCRETE FOURIER TRANSFORM Mårten Björkman Computer Vision and Active Perception School of Computer Science and Communication November 1, 2012 1 Terminology:
More informationAnalog vs. discrete signals
Analog vs. discrete signals Continuous-time signals are also known as analog signals because their amplitude is analogous (i.e., proportional) to the physical quantity they represent. Discrete-time signals
More information