Fall 2011, EE123 Digital Signal Processing
|
|
- Erika Campbell
- 5 years ago
- Views:
Transcription
1 Lecture 6 Miki Lustig, UCB September 11, 2012 Miki Lustig, UCB
2 DFT and Sampling the DTFT X (e jω ) = e j4ω sin2 (5ω/2) sin 2 (ω/2) 5 x[n] 25 X(e jω ) n ω 5 reconstructed x[n] 25 X(e jω ) n Miki Lustig UCB Based on Course otes by JM Kahn ω
3 Circular Convolution as Matrix Operation Circular convolution: h[n] x[n] = = H c x h[0] h[ 1] h[1] h[1] h[0] h[2] h[ 1] h[ 2] h[0] x[0] x[1] x[] H c is a circulant matrix The columns of the DFT matrix are Eigen vectors of circulant matrices Eigen vectors are DFT coefficients How can you show? Miki Lustig UCB Based on Course otes by JM Kahn
4 Circular Convolution as Matrix Operation Diagonalize: W H c Wn 1 = Right-multiply by W W H c = Multiply both sides by x W H c x = H[0] H[1] 0 0 H[ 1] H[0] H[1] 0 0 H[ 1] H[0] H[1] 0 0 H[ 1] W W x Miki Lustig UCB Based on Course otes by JM Kahn
5 Miki Lustig UCB Based on Course otes by JM Kahn
6 Fast Fourier Transform Algorithms We are interested in efficient computing methods for the DFT and inverse DFT: X [k] = x[n] = 1 n=0 1 k=0 x[n]w kn, k = 0,, 1 X [k]w kn, n = 0,, 1 where W = e j( 2π ) Miki Lustig UCB Based on Course otes by JM Kahn
7 Recall that we can use the DFT to compute the inverse DFT: DFT 1 {X [k]} = 1 (DFT {X [k]}) Hence, we can just focus on efficient computation of the DFT Straightforward computation of an -point DFT (or inverse DFT) requires 2 complex multiplications Miki Lustig UCB Based on Course otes by JM Kahn
8 Fast Fourier transform algorithms enable computation of an -point DFT (or inverse DFT) with the order of just log 2 complex multiplications This can represent a huge reduction in computational load, especially for large 2 log 2 2 log , ,024 1,048,576 10, ,192 67,108, , * 6Mp image size Miki Lustig UCB Based on Course otes by JM Kahn
9 Most FFT algorithms exploit the following properties of W kn : Conjugate Symmetry W k( n) Periodicity in n and k: W kn = W kn = W k(n+) = (W kn ) = W (k+)n Power: W 2 = W /2 Miki Lustig UCB Based on Course otes by JM Kahn
10 Most FFT algorithms decompose the computation of a DFT into successively smaller DFT computations Decimation-in-time algorithms decompose x[n] into successively smaller subsequences Decimation-in-frequency algorithms decompose X [k] into successively smaller subsequences We mostly discuss decimation-in-time algorithms here Assume length of x[n] is power of 2 ( = 2 ν ) If smaller zero-pad to closest power Miki Lustig UCB Based on Course otes by JM Kahn
11 Decimation-in-Time Fast Fourier Transform We start with the DFT X [k] = 1 n=0 x[n]w kn, k = 0,, 1 Separate the sum into even and odd terms: X [k] = x[n]w kn + x[n]w kn n even n odd These are two DFT s, each with half of the samples Miki Lustig UCB Based on Course otes by JM Kahn
12 Decimation-in-Time Fast Fourier Transform Let n = 2r (n even) and n = 2r + 1 (n odd): X [k] = = (/2) 1 r=0 (/2) 1 r=0 x[2r]w 2rk x[2r]w 2rk (/2) 1 + r=0 + W k x[2r + 1]W (2r+1)k (/2) 1 r=0 x[2r + 1]W 2rk ote that: W 2rk ( ) = e j( 2π )(2rk) = e j 2π rk /2 = W rk /2 Remember this trick, it will turn up often Miki Lustig UCB Based on Course otes by JM Kahn
13 Decimation-in-Time Fast Fourier Transform Hence: X [k] = (/2) 1 r=0 (/2) 1 x[2r]w/2 rk + W k r=0 = G[k] + W k H[k], k = 0,, 1 where we have defined: x[2r + 1]W rk /2 G[k] H[k] = = (/2) 1 r=0 (/2) 1 r=0 x[2r]w rk /2 DFT of even idx x[2r + 1]W rk /2 DFT of odd idx Miki Lustig UCB Based on Course otes by JM Kahn
14 Decimation-in-Time Fast Fourier Transform An 8 sample DFT can then be diagrammed as Even Samples Odd Samples x[0] x[2] x[4] x[6] x[1] x[3] x[5] x[7] /2 - Point DFT /2 - Point DFT G[0] G[1] G[2] G[3] H[0] H[1] H[2] H[3] 0 W 1 W 2 W 3 W 4 W 5 W 6 W W 7 X[0] X[1] X[2] X[3] X[4] X[5] X[6] X[7] Miki Lustig UCB Based on Course otes by JM Kahn
15 Decimation-in-Time Fast Fourier Transform Both G[k] and H[k] are periodic, with period /2 For example so G[k + /2] = = = (/2) 1 r=0 (/2) 1 r=0 (/2) 1 r=0 = G[k] x[2r]w r(k+/2) /2 x[2r]w rk /2 W r(/2) /2 x[2r]w rk /2 G[k + (/2)] = G[k] H[k + (/2)] = H[k] Miki Lustig UCB Based on Course otes by JM Kahn
16 Decimation-in-Time Fast Fourier Transform The periodicity of G[k] and H[k] allows us to further simplify For the first /2 points we calculate G[k] and W k H[k], and then compute the sum X [k] = G[k] + W k H[k] {k : 0 k < 2 } How does periodicity help for 2 k <? Miki Lustig UCB Based on Course otes by JM Kahn
17 Decimation-in-Time Fast Fourier Transform X [k] = G[k] + W k H[k] {k : 0 k < 2 } for 2 k < : W k+(/2) =? X [k + (/2)] =? Miki Lustig UCB Based on Course otes by JM Kahn
18 Decimation-in-Time Fast Fourier Transform X [k + (/2)] = G[k] W k H[k] We previously calculated G[k] and W k H[k] ow we only have to compute their difference to obtain the second half of the spectrum o additional multiplies are required Miki Lustig UCB Based on Course otes by JM Kahn
19 Decimation-in-Time Fast Fourier Transform The -point DFT has been reduced two /2-point DFTs, plus /2 complex multiplications The 8 sample DFT is then: x[0] G[k] X[0] Even Samples Odd Samples x[2] x[4] x[6] x[1] x[3] x[5] x[7] /2 - Point DFT /2 - Point DFT 0 W 1 W 2 W 3 W H[k] W k X[1] X[2] X[3] X[4] X[5] X[6] X[7] Miki Lustig UCB Based on Course otes by JM Kahn
20 Decimation-in-Time Fast Fourier Transform ote that the inputs have been reordered so that the outputs come out in their proper sequence We can define a butterfly operation, eg, the computation of X [0] and X [4] from G[0] and H[0]: G[0] X[0] =G[0] + W 0 H[0] H[0] W 0 X[4] =G[0] - W 0 H[0] This is an important operation in DSP Miki Lustig UCB Based on Course otes by JM Kahn
21 Decimation-in-Time Fast Fourier Transform Still O( 2 ) operations What shall we do? x[0] G[k] X[0] Even Samples Odd Samples x[2] x[4] x[6] x[1] x[3] x[5] x[7] /2 - Point DFT /2 - Point DFT 0 W 1 W 2 W 3 W k H[k] W X[1] X[2] X[3] X[4] X[5] X[6] X[7] Miki Lustig UCB Based on Course otes by JM Kahn
22 Decimation-in-Time Fast Fourier Transform We can use the same approach for each of the /2 point DFT s For the = 8 case, the /2 DFTs look like x[0] x[4] x[2] x[6] /4 - Point DFT /4 - Point DFT W /2 0 W /2 1 G[0] G[1] G[2] G[3] *ote that the inputs have been reordered again Miki Lustig UCB Based on Course otes by JM Kahn
23 Decimation-in-Time Fast Fourier Transform At this point for the 8 sample DFT, we can replace the /4 = 2 sample DFT s with a single butterfly The coefficient is W /4 = W 8/4 = W 2 = e jπ = 1 The diagram of this stage is then x[0] x[0] + x[4] x[4] 1 x[0] - x[4] Miki Lustig UCB Based on Course otes by JM Kahn
24 Decimation-in-Time Fast Fourier Transform Combining all these stages, the diagram for the 8 sample DFT is: x[0] X[0] x[4] x[2] x[6] x[1] x[5] x[3] x[7] 0 W /2 1 W /2 W 0 /2 W 1 /2 W 0 W 1 W 2 W 3 X[1] X[2] X[3] X[4] X[5] X[6] X[7] This the decimation-in-time FFT algorithm Miki Lustig UCB Based on Course otes by JM Kahn
25 Decimation-in-Time Fast Fourier Transform In general, there are log 2 stages of decimation-in-time Each stage requires /2 complex multiplications, some of which are trivial The total number of complex multiplications is (/2) log 2 The order of the input to the decimation-in-time FFT algorithm must be permuted First stage: split into odd and even Zero low-order bit first ext stage repeats with next zero-lower bit first et effect is reversing the bit order of indexes Miki Lustig UCB Based on Course otes by JM Kahn
26 Decimation-in-Time Fast Fourier Transform This is illustrated in the following table for = 8 Decimal Binary Bit-Reversed Binary Bit-Reversed Decimal Miki Lustig UCB Based on Course otes by JM Kahn
27 Decimation-in-Frequency Fast Fourier Transform The DFT is X [k] = 1 n=0 x[n]w nk If we only look at the even samples of X [k], we can write k = 2r, X [2r] = 1 n=0 x[n]w n(2r) We split this into two sums, one over the first /2 samples, and the second of the last /2 samples X [2r] = (/2) 1 n=0 x[n]w 2rn (/2) 1 + n=0 x[n + /2]W 2r(n+/2) Miki Lustig UCB Based on Course otes by JM Kahn
28 Decimation-in-Frequency Fast Fourier Transform But W 2r(n+/2) = W 2rnW = W 2rn We can then write = W rn /2 X [2r] = = = (/2) 1 n=0 (/2) 1 n=0 (/2) 1 n=0 x[n]w 2rn x[n]w 2rn (/2) 1 + n=0 (/2) 1 + n=0 (x[n] + x[n + /2]) W rn /2 x[n + /2]W 2r(n+/2) x[n + /2]W 2rn This is the /2-length DFT of first and second half of x[n] summed Miki Lustig UCB Based on Course otes by JM Kahn
29 Decimation-in-Frequency Fast Fourier Transform X [2r] = DFT {(x[n] + x[n + /2])} 2 X [2r + 1] = DFT {(x[n] x[n + /2]) W n } 2 (By a similar argument that gives the odd samples) Continue the same approach is applied for the /2 DFTs, and the /4 DFT s until we reach simple butterflies Miki Lustig UCB Based on Course otes by JM Kahn
30 Decimation-in-Frequency Fast Fourier Transform The diagram for and 8-point decimation-in-frequency DFT is as follows x[0] X[0] x[1] x[2] x[3] x[4] x[5] x[6] x[7] W 0 1 W 2 W 3 W W 0 /2 W 1 /2 W 0 /2 W 1 /2 X[4] X[2] X[6] X[1] X[5] X[3] X[7] This is just the decimation-in-time algorithm reversed! The inputs are in normal order, and the outputs are bit reversed Miki Lustig UCB Based on Course otes by JM Kahn
31 on-power-of-2 FFT s A similar argument applies for any length DFT, where the length is a composite number For example, if = 6, a decimation-in-time FFT could compute three 2-point DFT s followed by two 3-point DFT s x[0] x[3] 2-Point DFT W Point DFT X[0] X[2] x[1] x[4] 2-Point DFT W 6 1 X[4] X[1] x[2] x[5] 2-Point DFT W Point DFT X[3] X[5] Miki Lustig UCB Based on Course otes by JM Kahn
32 on-power-of-2 FFT s Good component DFT s are available for lengths up to 20 or so Many of these exploit the structure for that specific length For example, a factor of W /4 = e j 2π (/4) = e j π 2 = j Why? just swaps the real and imaginary components of a complex number, and doesn t actually require any multiplies Hence a DFT of length 4 doesn t require any complex multiplies Half of the multiplies of an 8-point DFT also don t require multiplication Composite length FFT s can be very efficient for any length that factors into terms of this order Miki Lustig UCB Based on Course otes by JM Kahn
33 For example = 693 factors into = (7)(9)(11) each of which can be implemented efficiently We would perform 9 11 DFT s of length DFT s of length 9, and 7 9 DFT s of length 11 Miki Lustig UCB Based on Course otes by JM Kahn
34 Historically, the power-of-two FFTs were much faster (better written and implemented) For non-power-of-two length, it was faster to zero pad to power of two Recently this has changed The free FFTW package implements very efficient algorithms for almost any filter length Matlab has used FFTW since version 6 Miki Lustig UCB Based on Course otes by JM Kahn
35 001 run time [ms] Miki Lustig UCB Based on Course otes by JM Kahn
36 FFT as Matrix Operation X [0] X [k] X [ 1] = W 00 W 0n W 0( 1) W k0 W kn W k( 1) W ( 1)0 W ( 1)n W ( 1)( 1) x[0] x[n] x[ 1] W is fully populated 2 entries Miki Lustig UCB Based on Course otes by JM Kahn
37 FFT as Matrix Operation X [0] X [k] X [ 1] = W 00 W 0n W 0( 1) W k0 W kn W k( 1) W ( 1)0 W ( 1)n W ( 1)( 1) x[0] x[n] x[ 1] W is fully populated 2 entries FFT is a decomposition of W into a more sparse form: [ I/2 D F = /2 I /2 D /2 ] [ W/2 0 0 W /2 ] [ Even-Odd Perm Matrix ] I /2 is an identity matrix D /2 is a diagonal with entries 1, W,, W /2 1 Miki Lustig UCB Based on Course otes by JM Kahn
38 FFT as Matrix Operation Example: = 4 F 4 = W W Miki Lustig UCB Based on Course otes by JM Kahn
39 Miki Lustig UCB Based on Course otes by JM Kahn
EE123 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! 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 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 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 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 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 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 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 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 informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 5 based on slides by J.M. Kahn Info Last time Finished DTFT Ch. 2 z-transforms Ch. 3 Today: DFT Ch. 8 Reminders: HW Due tonight The effects of sampling What is going
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 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 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 informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:305:45 CBC C222 Lecture 8 Frequency Analysis 14/02/18 http://www.ee.unlv.edu/~b1morris/ee482/
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 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 informationECSE 512 Digital Signal Processing I Fall 2010 FINAL EXAMINATION
FINAL EXAMINATION 9:00 am 12:00 pm, December 20, 2010 Duration: 180 minutes Examiner: Prof. M. Vu Assoc. Examiner: Prof. B. Champagne There are 6 questions for a total of 120 points. This is a closed book
More informationINTRODUCTION TO THE DFS AND THE DFT
ITRODUCTIO TO THE DFS AD THE DFT otes: This brief handout contains in very brief outline form the lecture notes used for a video lecture in a previous year introducing the DFS and the DFT. This material
More informationChapter 6: Applications of Fourier Representation Houshou Chen
Chapter 6: Applications of Fourier Representation Houshou Chen Dept. of Electrical Engineering, National Chung Hsing University E-mail: houshou@ee.nchu.edu.tw H.S. Chen Chapter6: Applications of Fourier
More 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 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 informationEE123 Digital Signal Processing. M. Lustig, EECS UC Berkeley
EE123 Digital Signal Processing Today Last time: DTFT - Ch 2 Today: Continue DTFT Z-Transform Ch. 3 Properties of the DTFT cont. Time-Freq Shifting/modulation: M. Lustig, EE123 UCB M. Lustig, EE123 UCB
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 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 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 information! Review: Discrete Fourier Transform (DFT) ! DFT Properties. " Duality. " Circular Shift. ! Circular Convolution. ! Fast Convolution Methods
Toda ESE 531: Digital Signal Processing! Review: Discrete Fourier Transform (DFT)! DFT Properties Lec 20: April 11, 2017 Discrete Fourier Transform, Pt 2 " Dualit " Circular Shift! Circular Convolution!
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 1 Time-Dependent FT Announcements! Midterm: 2/22/216 Open everything... but cheat sheet recommended instead 1am-12pm How s the lab going? Frequency Analysis with
More information1. Calculation of the DFT
ELE E4810: Digital Signal Processing Topic 10: The Fast Fourier Transform 1. Calculation of the DFT. The Fast Fourier Transform algorithm 3. Short-Time Fourier Transform 1 1. Calculation of the DFT! Filter
More 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 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 informationEE123 Digital Signal Processing
Aoucemets HW solutios posted -- self gradig due HW2 due Friday EE2 Digital Sigal Processig ham radio licesig lectures Tue 6:-8pm Cory 2 Lecture 6 based o slides by J.M. Kah SDR give after GSI Wedesday
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 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 information1 1.27z z 2. 1 z H 2
E481 Digital Signal Processing Exam Date: Thursday -1-1 16:15 18:45 Final Exam - Solutions Dan Ellis 1. (a) In this direct-form II second-order-section filter, the first stage has
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 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 informationFall 2011, EE123 Digital Signal Processing
Lecture 5 Miki Lustig, UCB September 14, 211 Miki Lustig, UCB Motivatios for Discrete Fourier Trasform Sampled represetatio i time ad frequecy umerical Fourier aalysis requires a Fourier represetatio that
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 informationCSE 548: Analysis of Algorithms. Lecture 4 ( Divide-and-Conquer Algorithms: Polynomial Multiplication )
CSE 548: Analysis of Algorithms Lecture 4 ( Divide-and-Conquer Algorithms: Polynomial Multiplication ) Rezaul A. Chowdhury Department of Computer Science SUNY Stony Brook Spring 2015 Coefficient Representation
More informationESE 531: Digital Signal Processing
ESE 531: Digital Signal Processing Lec 22: April 10, 2018 Adaptive Filters Penn ESE 531 Spring 2018 Khanna Lecture Outline! Circular convolution as linear convolution with aliasing! Adaptive Filters Penn
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 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 informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 2B D. T. Fourier Transform M. Lustig, EECS UC Berkeley Something Fun gotenna http://www.gotenna.com/# Text messaging radio Bluetooth phone interface MURS VHF radio
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even
More informationDiscrete Time Fourier Transform (DTFT)
Discrete Time Fourier Transform (DTFT) 1 Discrete Time Fourier Transform (DTFT) The DTFT is the Fourier transform of choice for analyzing infinite-length signals and systems Useful for conceptual, pencil-and-paper
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Discrete Time Fourier Transform M. Lustig, EECS UC Berkeley A couple of things Read Ch 2 2.0-2.9 It s OK to use 2nd edition Class webcast in bcourses.berkeley.edu or linked
More informationLecture 7. Fourier Analysis
Lecture 7 Fourier Analysis Summary Lecture 6 Minima and maxima 1 dimension : Bracket minima : 3 values of f(x) : f(2) < f(1) and f(2)
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 information! Spectral Analysis with DFT. ! Windowing. ! Effect of zero-padding. ! Time-dependent Fourier transform. " Aka short-time Fourier transform
Lecture Outline ESE 531: Digital Signal Processing Spectral Analysis with DFT Windowing Lec 24: April 18, 2019 Spectral Analysis Effect of zero-padding Time-dependent Fourier transform " Aka short-time
More informationModule 3. Convolution. Aim
Module Convolution Digital Signal Processing. Slide 4. Aim How to perform convolution in real-time systems efficiently? Is convolution in time domain equivalent to multiplication of the transformed sequence?
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 informationChapter 8 The Discrete Fourier Transform
Chapter 8 The Discrete Fourier Transform Introduction Representation of periodic sequences: the discrete Fourier series Properties of the DFS The Fourier transform of periodic signals Sampling the Fourier
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 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 informationDigital Signal Processing Lecture 3 - Discrete-Time Systems
Digital Signal Processing - Discrete-Time Systems Electrical Engineering and Computer Science University of Tennessee, Knoxville August 25, 2015 Overview 1 2 3 4 5 6 7 8 Introduction Three components of
More informationContents. Digital Signal Processing, Part II: Power Spectrum Estimation
Contents Digital Signal Processing, Part II: Power Spectrum Estimation 5. Application of the FFT for 7. Parametric Spectrum Est. Filtering and Spectrum Estimation 7.1 ARMA-Models 5.1 Fast Convolution 7.2
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 informationE The Fast Fourier Transform
Fourier Transform Methods in Finance By Umberto Cherubini Giovanni Della Lunga Sabrina Mulinacci Pietro Rossi Copyright 2010 John Wiley & Sons Ltd E The Fast Fourier Transform E.1 DISCRETE FOURIER TRASFORM
More informationLecture 7 January 26, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 7 January 26, 26 Prof Emmanuel Candes Scribe: Carlos A Sing-Long, Edited by E Bates Outline Agenda:
More informationFourier Analysis of Signals Using the DFT
Fourier Analysis of Signals Using the DFT ECE 535 Lecture April 29, 23 Overview: Motivation Many applications require analyzing the frequency content of signals Speech processing study resonances of vocal
More informationDigital Signal Processing
DIGITAL SIGNAL PROCESSING SUBJECT CODE : NO. OF LECTURE HRS/WEEK : 04 TOTAL NO. OF LECTURE HRS. : 52 IA MARKS : 25 EXAM HOURS : 03 EXAMMARKS : 100 UNIT - 1 DISCRETE FOURIER TRANSFORMS (DFT): FREQUENCY
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 informationChapter 2: The Fourier Transform
EEE, EEE Part A : Digital Signal Processing Chapter Chapter : he Fourier ransform he Fourier ransform. Introduction he sampled Fourier transform of a periodic, discrete-time signal is nown as the discrete
More informationDFT-Based FIR Filtering. See Porat s Book: 4.7, 5.6
DFT-Based FIR Filtering See Porat s Book: 4.7, 5.6 1 Motivation: DTFT View of Filtering There are two views of filtering: * Time Domain * Frequency Domain x[ X f ( θ ) h[ H f ( θ ) Y y[ = h[ * x[ f ( θ
More informationBME 50500: Image and Signal Processing in Biomedicine. Lecture 2: Discrete Fourier Transform CCNY
1 Lucas Parra, CCNY BME 50500: Image and Signal Processing in Biomedicine Lecture 2: Discrete Fourier Transform Lucas C. Parra Biomedical Engineering Department CCNY http://bme.ccny.cuny.edu/faculty/parra/teaching/signal-and-image/
More informationChap 2. Discrete-Time Signals and Systems
Digital Signal Processing Chap 2. Discrete-Time Signals and Systems Chang-Su Kim Discrete-Time Signals CT Signal DT Signal Representation 0 4 1 1 1 2 3 Functional representation 1, n 1,3 x[ n] 4, n 2 0,
More informationDFT & Fast Fourier Transform PART-A. 7. Calculate the number of multiplications needed in the calculation of DFT and FFT with 64 point sequence.
SHRI ANGALAMMAN COLLEGE OF ENGINEERING & TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR,TRICHY-621105. DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING UNIT I DFT & Fast Fourier
More informationE : Lecture 1 Introduction
E85.2607: Lecture 1 Introduction 1 Administrivia 2 DSP review 3 Fun with Matlab E85.2607: Lecture 1 Introduction 2010-01-21 1 / 24 Course overview Advanced Digital Signal Theory Design, analysis, and implementation
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 informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 2A D.T Systems D. T. Fourier Transform A couple of things Read Ch 2 2.0-2.9 It s OK to use 2nd edition My office hours: posted on-line W 4-5pm Cory 506 ham radio
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 informationLTI H. the system H when xn [ ] is the input.
REVIEW OF 1D LTI SYSTEMS LTI xn [ ] y[ n] H Operator notation: y[ n] = H{ x[ n] } In English, this is read: yn [ ] is the output of the system H when xn [ ] is the input. 2 THE MOST IMPORTANT PROPERTIES
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 401 Digital Signal Processing Prof. Mark Fowler Note Set #21 Using the DFT to Implement FIR Filters Reading Assignment: Sect. 7.3 of Proakis & Manolakis Motivation: DTFT View of Filtering There are
More information( ) John A. Quinn Lecture. ESE 531: Digital Signal Processing. Lecture Outline. Frequency Response of LTI System. Example: Zero on Real Axis
John A. Quinn Lecture ESE 531: Digital Signal Processing Lec 15: March 21, 2017 Review, Generalized Linear Phase Systems Penn ESE 531 Spring 2017 Khanna Lecture Outline!!! 2 Frequency Response of LTI System
More informationLAB 2: DTFT, DFT, and DFT Spectral Analysis Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 2: DTFT, DFT, and DFT Spectral
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 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 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 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 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 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 informationFinite-length Discrete Transforms. Chapter 5, Sections
Fiite-legth Discrete Trasforms Chapter 5, Sectios 5.2-50 5.0 Dr. Iyad djafar Outlie The Discrete Fourier Trasform (DFT) Matrix Represetatio of DFT Fiite-legth Sequeces Circular Covolutio DFT Symmetry Properties
More informationWEEK10 Fast Fourier Transform (FFT) (Theory and Implementation)
WEEK10 Fast Fourier Transform (FFT) (Theory and Implementation) Learning Objectives DFT algorithm. Conversion of DFT to FFT algorithm. Implementation of the FFT algorithm. Chapter 19, Slide 2 DFT Algorithm
More informationDiscrete-Time Signals and Systems
ECE 46 Lec Viewgraph of 35 Discrete-Time Signals and Systems Sequences: x { x[ n] }, < n
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 informationL. Vandenberghe EE133A (Spring 2017) 3. Matrices. notation and terminology. matrix operations. linear and affine functions.
L Vandenberghe EE133A (Spring 2017) 3 Matrices notation and terminology matrix operations linear and affine functions complexity 3-1 Matrix a rectangular array of numbers, for example A = 0 1 23 01 13
More informationDigital Signal Processing Chapter 10. Fourier Analysis of Discrete- Time Signals and Systems CHI. CES Engineering. Prof. Yasser Mostafa Kadah
Digital Signal Processing Chapter 10 Fourier Analysis of Discrete- Time Signals and Systems Prof. Yasser Mostafa Kadah CHI CES Engineering Discrete-Time Fourier Transform Sampled time domain signal has
More information5 Irreducible representations
Physics 129b Lecture 8 Caltech, 01/1/19 5 Irreducible representations 5.5 Regular representation and its decomposition into irreps To see that the inequality is saturated, we need to consider the so-called
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 informationDSP Configurations. responded with: thus the system function for this filter would be
DSP Configurations In this lecture we discuss the different physical (or software) configurations that can be used to actually realize or implement DSP functions. Recall that the general form of a DSP
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 informationVery useful for designing and analyzing signal processing systems
z-transform z-transform The z-transform generalizes the Discrete-Time Fourier Transform (DTFT) for analyzing infinite-length signals and systems Very useful for designing and analyzing signal processing
More informationDigital Signal Processing. Lecture Notes and Exam Questions DRAFT
Digital Signal Processing Lecture Notes and Exam Questions Convolution Sum January 31, 2006 Convolution Sum of Two Finite Sequences Consider convolution of h(n) and g(n) (M>N); y(n) = h(n), n =0... M 1
More informationOLA and FBS Duality Review
MUS421/EE367B Lecture 10A Review of OverLap-Add (OLA) and Filter-Bank Summation (FBS) Interpretations of Short-Time Fourier Analysis, Modification, and Resynthesis Julius O. Smith III (jos@ccrma.stanford.edu)
More informationDigital Signal Processing. Midterm 1 Solution
EE 123 University of California, Berkeley Anant Sahai February 15, 27 Digital Signal Processing Instructions Midterm 1 Solution Total time allowed for the exam is 8 minutes Some useful formulas: Discrete
More informationDiscrete-time signals and systems
Discrete-time signals and systems 1 DISCRETE-TIME DYNAMICAL SYSTEMS x(t) G y(t) Linear system: Output y(n) is a linear function of the inputs sequence: y(n) = k= h(k)x(n k) h(k): impulse response of the
More informationComplex symmetry Signals and Systems Fall 2015
18-90 Signals and Systems Fall 015 Complex symmetry 1. Complex symmetry This section deals with the complex symmetry property. As an example I will use the DTFT for a aperiodic discrete-time signal. The
More informationFast Fourier Transform Discrete-time windowing Discrete Fourier Transform Relationship to DTFT Relationship to DTFS Zero padding
Fast Fourier Transform Discrete-time windowing Discrete Fourier Transform Relationship to DTFT Relationship to DTFS Zero padding Fourier Series & Transform Summary x[n] = X[k] = 1 N k= n= X[k]e jkω
More informationFOURIER - GENIUS OR TERRORISTE?
FOURIER - GEIUS OR TERRORISTE? Samir Chettri HUGHES-STX, Kanazawa University, GST, UMBC In this memo we discuss certain aspects of Fourier analysis This is by no means a complete document However, with
More informationIn this Lecture. Frequency domain analysis
In this Lecture Frequency domain analysis Introduction In most cases we want to know the frequency content of our signal Why? Most popular analysis in frequency domain is based on work of Joseph Fourier
More information