WEEK10 Fast Fourier Transform (FFT) (Theory and Implementation)
|
|
- Quentin Hodge
- 5 years ago
- Views:
Transcription
1 WEEK10 Fast Fourier Transform (FFT) (Theory and Implementation)
2 Learning Objectives DFT algorithm. Conversion of DFT to FFT algorithm. Implementation of the FFT algorithm. Chapter 19, Slide 2
3 DFT Algorithm The Fourier transform of an analogue signal x(t) is given by: The Discrete Fourier Transform (DFT) of a discrete-time signal x(nt) is given by: Where: Chapter 19, Slide 3
4 DFT Algorithm If we let: then: Chapter 19, Slide 4
5 DFT Algorithm x[n] = input X[k] = frequency bins W = twiddle factors X(0) X(1) X(k) X(N) = x[0]w N0 + x[1]w 0*1 N + + x[n]w 0*(N) N = x[0]w N0 + x[1]w 1*1 N + + x[n]w 1*(N) N : = x[0]w N0 + x[1]w k*1 N + + x[n]w k*(n) N : = x[0]w N0 + x[1]w (N)*1 N + + x[n]w (N)(N) N Note: For N samples of x we have N frequencies representing the signal. Chapter 19, Slide 5
6 Performance of the DFT Algorithm The DFT requires N 2 (NxN) complex multiplications: Each X(k) requires N complex multiplications. Therefore to evaluate all the values of the DFT ( X(0) to X(N) ) N 2 multiplications are required. The DFT also requires (N)*N complex additions: Each X(k) requires N additions. Therefore to evaluate all the values of the DFT (N)*N additions are required. Chapter 19, Slide 6
7 Performance of the DFT Algorithm Can the number of computations required be reduced? Chapter 19, Slide 7
8 DFT FFT A large amount of work has been devoted to reducing the computation time of a DFT. This has led to efficient algorithms which are known as the Fast Fourier Transform (FFT) algorithms. Chapter 19, Slide 8
9 DFT FFT [1] x[n] = x[0], x[1],, x[n] Lets divide the sequence x[n] into even and odd sequences: x[2n] = x[0], x[2],, x[n-2] x[2n+1] = x[1], x[3],, x[n] Chapter 19, Slide 9
10 DFT FFT Equation 1 can be rewritten as: [2] Since: Then: Chapter 19, Slide 10
11 DFT FFT The result is that an N-point DFT can be divided into two N/2 point DFT s: N-point DFT Where Y(k) and Z(k) are the two N/2 point DFTs operating on even and odd samples respectively: Two N/2- point DFTs Chapter 19, Slide 11
12 DFT FFT Periodicity and symmetry of W can be exploited to simplify the DFT further: [3] Or: : Symmetry And: Chapter 19, Slide 12 : Periodicity
13 DFT FFT Symmetry and periodicity: W 8 4 W 8 5 W 8 3 W 8 6 W 8 2 W 8 7 W 8 0 =W 8 8 W 8 1 = W 8 9 W k+n/2 N = - W k N W k+n/2 N/2 =W k N/2 W k+4 8 = - W k 8 W k+8 8 =W k 8 Chapter 19, Slide 13
14 DFT FFT Finally by exploiting the symmetry and periodicity, Equation 3 can be written as: [4] Chapter 19, Slide 14
15 DFT FFT Y(k) and W N k Z(k) only need to be calculated once and used for both equations. Note: the calculation is reduced from 0 to N to 0 to (N/2-1). Chapter 19, Slide 15
16 DFT FFT Y(k) and Z(k) can also be divided into N/4 point DFTs using the same process shown above: Chapter 19, Slide 16 The process continues until we reach 2 point DFTs.
17 DFT FFT x[0] x[2] x[4] N/2 point DFT y[0] y[1] y[2] X[0] = y[0]+ z[0] X[1] = y[1]+w 1 z[1] x[n-2] y[n-2] x[1] x[3] x[5] N/2 point DFT z[0] z[1] z[2] W 1 X[N/2] = y[0]- z[0] X[N/2+1] = y[1]-w 1 z[1] x[n] z[n/2] Illustration of the first decimation in time FFT. Chapter 19, Slide 17
18 FFT Implementation Calculation of the output of a butterfly : Y(k) = U r +ju i U =U r +ju i = X(k) W Nk Z(k) = (L r +jl i )(W r +jw i ) L =L r +jl i = X(k+N/2) Key: U = Upper r = real L = Lower i = imaginary Different methods are available for calculating the outputs U and L. The best method is the one with the least number of multiplications and additions. Chapter 19, Slide 18
19 FFT Implementation Calculation of the output of a butterfly : U r +ju i U =U r +ju i (L r +jl i )(W r +jw i ) L =L r +jl i Chapter 19, Slide 19
20 FFT Implementation Calculation of the output of a butterfly : U r +ju i U =U r +ju i = (L r W r -L i W i + U r )+j(l r W i + L i W r + U i ) (L r +jl i )(W r +jw i ) To further minimise the number of operations (* and +), the following are calculated only once: temp1 = L r W r temp2 = L i W i temp3 = L r W i temp4 = L i W r temp1_2 = temp1 - temp2 L =L r +jl i = (U r -L r W r + L i W i )+j(u i -L r W i -L i W r ) temp3_4 = temp3 + temp4 U r = temp1 - temp2 + U r = temp1_2 + U r U i = temp3 + temp4 + U i = temp3_4 + U i L r = U r - temp1 + temp2 = U r - temp1_2 L i = U i - temp3 - temp4 = U i - temp3_4 Chapter 19, Slide 20
21 FFT Implementation (Butterfly Calculation) Converting the butterfly calculation into C code: temp1 = (y[lower].real * WR); temp2 = (y[lower].imag * WI); temp3 = (y[lower].real * WI); temp4 = (y[lower].imag * WR); temp1_2 = temp1 - temp2; temp3_4 = temp 3 + temp4; y[upper].real = temp1_2 + y[upper].real; y[upper].imag = temp3_4 + y[upper].imag; y[lower].imag = y[upper].imag - temp3_4; y[lower].real = y[upper].real - temp1_2; Chapter 19, Slide 21
22 FFT Implementation To efficiently implement the FFT algorithm a few observations are made: Each stage has the same number of butterflies (number of butterflies = N/2, N is number of points). The number of DFT groups per stage is equal to (N/2 stage ). The difference between the upper and lower leg is equal to 2 stage. The number of butterflies in the group is equal to 2 stage. Chapter 19, Slide 22
23 FFT Implementation Example: 8 point FFT W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Chapter 19, Slide 23 Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
24 FFT Implementation Example: 8 point FFT (1) Number of stages: W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Chapter 19, Slide 24 Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
25 Stage 1 FFT Implementation Example: 8 point FFT (1) Number of stages: N stages = 1 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Chapter 19, Slide 25 Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
26 Stage 1 FFT Implementation Stage 2 Example: 8 point FFT (1) Number of stages: N stages = 2 W 1 W 3 Chapter 19, Slide 26 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
27 Stage 1 FFT Implementation Stage 2 Stage 3 Example: 8 point FFT (1) Number of stages: N stages = 3 W 1 W 3 Chapter 19, Slide 27 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
28 Stage 1 FFT Implementation Stage 2 Stage 3 Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: W 1 W 3 Chapter 19, Slide 28 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
29 Block 1 Stage 1 FFT Implementation Stage 2 Stage 3 Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 1 W 1 W 3 Chapter 19, Slide 29 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
30 Block 1 Block 2 Stage 1 FFT Implementation Stage 2 Stage 3 Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 2 W 1 W 3 Chapter 19, Slide 30 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
31 Block 1 Block 2 Block 3 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 3 W 3 Chapter 19, Slide 31 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
32 Block 1 Block 2 Block 3 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Block 4 W 3 Chapter 19, Slide 32 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
33 Block 1 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 1 W 3 Chapter 19, Slide 33 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
34 Block 1 Block 2 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 W 3 Chapter 19, Slide 34 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
35 Block 1 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 W 3 Chapter 19, Slide 35 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage
36 Chapter 19, Slide 36 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 (3) B flies/block: Stage 1:
37 Chapter 19, Slide 37 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 (3) B flies/block: Stage 1: N btf = 1
38 Chapter 19, Slide 38 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 (3) B flies/block: Stage 1: N btf = 1 Stage 2: N btf = 1
39 Chapter 19, Slide 39 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 (3) B flies/block: Stage 1: N btf = 1 Stage 2: N btf = 2
40 Chapter 19, Slide 40 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 (3) B flies/block: Stage 1: N btf = 1 Stage 2: N btf = 2 Stage 3: N btf = 1
41 Chapter 19, Slide 41 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 (3) B flies/block: Stage 1: N btf = 1 Stage 2: N btf = 2 Stage 3: N btf = 2
42 Chapter 19, Slide 42 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 (3) B flies/block: Stage 1: N btf = 1 Stage 2: N btf = 2 Stage 3: N btf = 3
43 Chapter 19, Slide 43 Stage 1 FFT Implementation Stage 2 Stage 3 W 1 W 3 Decimation in time FFT: Number of stages = log 2 N Number of blocks/stage = N/2 stage Number of butterflies/block = 2 stage Example: 8 point FFT (1) Number of stages: N stages = 3 (2) Blocks/stage: Stage 1: N blocks = 4 Stage 2: N blocks = 2 Stage 3: N blocks = 1 (3) B flies/block: Stage 1: N btf = 1 Stage 2: N btf = 2 Stage 3: N btf = 4
44 FFT Implementation Stage 1 Stage 2 Stage 3 W 1 W 3 Start Index Input Index Twiddle Factor Index N/2 = 4 Chapter 19, Slide 44
45 FFT Implementation Stage 1 Stage 2 Stage 3 W 1 W 3 Start Index Input Index Twiddle Factor Index N/2 = 4 4/2 = 2 Chapter 19, Slide 45
46 FFT Implementation Stage 1 Stage 2 Stage 3 W 1 W 3 Start Index Input Index Twiddle Factor Index N/2 = 4 4/2 = 2 2/2 = 1 Chapter 19, Slide 46
47 FFT Implementation Stage 1 Stage 2 Stage 3 W 1 W 3 Start Index 0 Input Index 1 Twiddle Factor Index N/2 = 4 4/2 = 2 2/2 = 1 Indicies Used W Chapter 19, Slide 47 W 3
48 FFT Implementation The most important aspect of converting the FFT diagram to C code is to calculate the upper and lower indicies of each butterfly: GS = N/4; /* Group step initial value */ step = 1; /* Initial value */ NB = N/2; /* NB is a constant */ for (k=n; k>1; k>>1) /* Repeat this loop for each stage */ { for (j=0; j<n; j+=gs) /* Repeat this loop for each block */ { for (n=j; n<(j+gs); n+=step) /* Repeat this loop for each butterfly */ { upperindex = n; lowerindex = n+step; } } /* Change the GS and step for the next stage */ } GS = GS << 1; step = step << 1; Chapter 19, Slide 48
49 FFT Implementation How to declare and access twiddle factors: (1) How to declare the twiddle factors: struct { short real; // * cos (2*pi*n) -> Q15 format short imag; // * sin (2*pi*n) -> Q15 format } w[] = { 32767,0, 32767,-201,... }; (2) How to access the twiddle factors: short temp_real, temp_imag; temp_real = w[i].real; temp_imag = w[i].imag; Chapter 19, Slide 49
50 Links: FFT Implementation Project location: \Code\Chapter 19 - Fast Fourier Transform\ Projects: FFT Decimation in Time: \FFT_C_Fixed_DIT\ FFT Decimation in Frequency: \FFT_C_Fixed_DIF\ Decimation in Time Code: FFTMain.c FFT code (fft1024.c) Laboratory sheet: FFTLabSheet.pdf Chapter 19, Slide 50
51 WEEK10 Fast Fourier Transform (FFT) (Theory and Implementation) - End -
EE482: 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 informationDigital Signal Processing ECE-307. Tarun Gulati Assoc. Prof., ECE Deptt. MMEC, MMU, Mullana
Digital Signal Processing ECE-307 Tarun Gulati Assoc. Prof., ECE Deptt. MMEC, MMU, Mullana The z-transform: Introduction Definition of the z-transform Definition:The z-transform of a discrete-time signal
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationTransforms and Orthogonal Bases
Orthogonal Bases Transforms and Orthogonal Bases We now turn back to linear algebra to understand transforms, which map signals between different domains Recall that signals can be interpreted as vectors
More informationThe 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 informationFinal Exam January 31, Solutions
Final Exam January 31, 014 Signals & Systems (151-0575-01) Prof. R. D Andrea & P. Reist Solutions Exam Duration: Number of Problems: Total Points: Permitted aids: Important: 150 minutes 7 problems 50 points
More informationi r-s THE MEMPHIS, TENN., SATURDAY. DEGfMBER
N k Q2 90 k ( < 5 q v k 3X3 0 2 3 Q :: Y? X k 3 : \ N 2 6 3 N > v N z( > > :}9 [ ( k v >63 < vq 9 > k k x k k v 6> v k XN Y k >> k < v Y X X X NN Y 2083 00 N > N Y Y N 0 \ 9>95 z {Q ]k3 Q k x k k z x X
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 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 information1 Pheasant Lane December 1990 Ithaca, NY (607) AN INTUITIVE APPROACH TO FFT ALGORITHMS
APPLICATION NOTE NO. Pheasant Lane December Ithaca, NY ()-- AN INTUITIVE APPROACH TO FFT ALGORITHMS Most readers are probably familiar with the general notion of an FFT (Fast Fourier Transform). An FFT
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 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 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 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 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 informationDFT and Matlab with some examples.
DFT and Matlab with some examples 1 www.icrf.nl Fourier transform Fourier transform is defined as: X + jωt = x( t e dt ( ω ) ) with ω= 2 π f Rad/s And x(t) a signal in the time domain. www.icrf.nl 2 Fourier
More informationMIS S BALLS, L.L.A.]
& N k k QY GN ( x - N N & N & k QY GN x 00 - XX N X ± - - - - ---------------- N N G G N N N Y NG 5 880N GN N X GN x ( G ) 8N ---- N 8 Y 8 - N N ( G () G ( ) (N) N? k [ x-k NNG G G k k N NY Y /( Q G (-)
More informationx (2) k even h n=(n) + (n+% x(6) X(3), x (5) 4 x(4)- x(), x (2), Decomposition of an N-point DFT into 2 N/2-point DFT's.
COMPUTATION OF DISCRETE FOURIER TRANSFORM - PART 2 1. Lecture 19-49 minutes k even n.o h n=(n) + (n+% x (2), x (2) x(4)- x(6) Decomposition of an N-point DFT into 2 N/2-point DFT's. X(3), x (5) 4 x(),
More informationL29: Fourier analysis
L29: Fourier analysis Introduction The discrete Fourier Transform (DFT) The DFT matrix The Fast Fourier Transform (FFT) The Short-time Fourier Transform (STFT) Fourier Descriptors CSCE 666 Pattern Analysis
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 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 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 informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture 19: Lattice Filters Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 2008 K. E. Barner (Univ. of Delaware) ELEG
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 informationAPPENDIX A. The Fourier integral theorem
APPENDIX A The Fourier integral theorem In equation (1.7) of Section 1.3 we gave a description of a signal defined on an infinite range in the form of a double integral, with no explanation as to how that
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 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 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 informationFourier Series Representation of
Fourier Series Representation of Periodic Signals Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline The response of LIT system
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 informationSummary of lecture 1. E x = E x =T. X T (e i!t ) which motivates us to define the energy spectrum Φ xx (!) = jx (i!)j 2 Z 1 Z =T. 2 d!
Summary of lecture I Continuous time: FS X FS [n] for periodic signals, FT X (i!) for non-periodic signals. II Discrete time: DTFT X T (e i!t ) III Poisson s summation formula: describes the relationship
More informationIntroduction to Digital Signal Processing
Introduction to Digital Signal Processing 1.1 What is DSP? DSP is a technique of performing the mathematical operations on the signals in digital domain. As real time signals are analog in nature we need
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 informationEDISP (NWL2) (English) Digital Signal Processing Transform, FT, DFT. March 11, 2015
EDISP (NWL2) (English) Digital Signal Processing Transform, FT, DFT March 11, 2015 Transform concept We want to analyze the signal represent it as built of some building blocks (well known signals), possibly
More 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 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 informationSistemas de Aquisição de Dados. Mestrado Integrado em Eng. Física Tecnológica 2015/16 Aula 6-26 de Outubro
Sistemas de Aquisição de Dados Mestrado Integrado em Eng. Física Tecnológica 2015/16 Aula 6-26 de Outubro Flash Decoder Thermometer code Wired NOR based decoder 2 Successive Approximation ADC (SAR) CONVERT
More informationDigital Signal Processing: Signal Transforms
Digital Signal Processing: Signal Transforms Aishy Amer, Mohammed Ghazal January 19, 1 Instructions: 1. This tutorial introduces frequency analysis in Matlab using the Fourier and z transforms.. More Matlab
More informationAppendix 3: FFT Computer Programs
onnexions module: m17397 1 Appendix 3: FFT omputer Programs. Sidney Burrus This work is produced by The onnexions Project and licensed under the reative ommons Attribution License Abstract Fortran programs
More informationMicrocontrollers. Fast - Fourier - Transformation. þ additional file AP EXE available
Microcontrollers Apote AP634 þ additional file AP634.EXE available ast - ourier - Transformation The ast ourier Transformation (T is an algorithm frequently used in various applications, like telecommunication,
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 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 information6.003 Homework #10 Solutions
6.3 Homework # Solutions Problems. DT Fourier Series Determine the Fourier Series coefficients for each of the following DT signals, which are periodic in N = 8. x [n] / n x [n] n x 3 [n] n x 4 [n] / n
More informationIT DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A
DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING IT6502 - DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A 1. What is a continuous and discrete time signal? Continuous
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 informationHEAGAN & CO., OPP. f>, L. & W. DEPOT, DOYER, N. J, OUR MOTTO! ould Iwv ia immediate vltlui. VEEY BEST NEW Creamery Butter 22c ib,
#4 NN N G N N % XX NY N Y FY N 2 88 N 28 k N k F P X Y N Y /» 2«X ««!!! 8 P 3 N 0»9! N k 25 F $ 60 $3 00 $3000 k k N 30 Y F00 6 )P 0» «{ N % X zz» «3 0««5 «N «XN» N N 00/ N 4 GN N Y 07 50 220 35 2 25 0
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 informationSolutions. Number of Problems: 10
Final Exam February 4th, 01 Signals & Systems (151-0575-01) Prof. R. D Andrea Solutions Exam Duration: 150 minutes Number of Problems: 10 Permitted aids: One double-sided A4 sheet. Questions can be answered
More information!Sketch f(t) over one period. Show that the Fourier Series for f(t) is as given below. What is θ 1?
Second Year Engineering Mathematics Laboratory Michaelmas Term 998 -M L G Oldfield 30 September, 999 Exercise : Fourier Series & Transforms Revision 4 Answer all parts of Section A and B which are marked
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
More informationMATH3283W LECTURE NOTES: WEEK 6 = 5 13, = 2 5, 1 13
MATH383W LECTURE NOTES: WEEK 6 //00 Recursive sequences (cont.) Examples: () a =, a n+ = 3 a n. The first few terms are,,, 5 = 5, 3 5 = 5 3, Since 5
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 informationBEE604 Digital Signal Processing
BEE64 Digital Signal Processing Copiled by, Mrs.S.Sherine Assistant Professor Departent of EEE BIHER. COTETS Sapling Discrete Tie Fourier Transfor Properties of DTFT Discrete Fourier Transfor Inverse Discrete
More informationNeed for transformation?
Z-TRANSFORM In today s class Z-transform Unilateral Z-transform Bilateral Z-transform Region of Convergence Inverse Z-transform Power Series method Partial Fraction method Solution of difference equations
More informationLOWELL JOURNAL CRESP0 REIGNS. MICE DESTROY #10,000 CASH. T O W E L L S T A T E B A N K,
N V 2 2 7 2 9 2 K K «KN < K KY KK K KN ( *
More information6.003 Signal Processing
6.003 Signal Processing Week 6, Lecture A: The Discrete Fourier Transform (DFT) Adam Hartz hz@mit.edu What is 6.003? What is a signal? Abstractly, a signal is a function that conveys information Signal
More informationImage Acquisition and Sampling Theory
Image Acquisition and Sampling Theory Electromagnetic Spectrum The wavelength required to see an object must be the same size of smaller than the object 2 Image Sensors 3 Sensor Strips 4 Digital Image
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 informationDIGITAL SIGNAL PROCESSING UNIT I
DIGITAL SIGNAL PROCESSING UNIT I CONTENTS 1.1 Introduction 1.2 Introduction about signals 1.3 Signals Processing 1.4 Classification of signals 1.5 Operations performed on signals 1.6 Properties of signals
More informationEE216B: VLSI Signal Processing. FFT Processors. Prof. Dejan Marković FFT: Background
4/30/0 EE6B: VLSI Signal Processing FFT Processors Prof. Dejan Marković ee6b@gmail.com FFT: Background A bit of history 805 - algorithm first described by Gauss 965 - algorithm rediscovered (not for the
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 informationMATH 412: Homework # 5 Tim Ahn July 22, 2016
MATH 412: Homework # 5 Tim Ahn July 22, 2016 Supplementary Exercises Supp #23 Find the Fourier matrices F for = 2 and = 4. set.seed(2) = 2 omega = 2*pi/ F = exp(omega*outer(0:(-1),0:(-1))*(0+1i)) F ##
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 informationQuiz. Good luck! Signals & Systems ( ) Number of Problems: 19. Number of Points: 19. Permitted aids:
Quiz November th, Signals & Systems (--) Prof. R. D Andrea Quiz Exam Duration: Min Number of Problems: Number of Points: Permitted aids: Important: None Questions must be answered on the provided answer
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 informationCanonic FFT flow graphs for real-valued even/odd symmetric inputs
Lao and Parhi EURASIP Journal on Advances in Signal Processing (017) 017:45 DOI 10.1186/s13634-017-0477-9 EURASIP Journal on Advances in Signal Processing RESEARCH Canonic FFT flow graphs for real-valued
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 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 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 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 informationLet H(z) = P(z)/Q(z) be the system function of a rational form. Let us represent both P(z) and Q(z) as polynomials of z (not z -1 )
Review: Poles and Zeros of Fractional Form Let H() = P()/Q() be the system function of a rational form. Let us represent both P() and Q() as polynomials of (not - ) Then Poles: the roots of Q()=0 Zeros:
More information6.003 Signal Processing
6.003 Signal Processing Week 6, Lecture A: The Discrete Fourier Transform (DFT) Adam Hartz hz@mit.edu What is 6.003? What is a signal? Abstractly, a signal is a function that conveys information Signal
More informationIntroduction to DFT. Deployment of Telecommunication Infrastructures. Azadeh Faridi DTIC UPF, Spring 2009
Introduction to DFT Deployment of Telecommunication Infrastructures Azadeh Faridi DTIC UPF, Spring 2009 1 Review of Fourier Transform Many signals can be represented by a fourier integral of the following
More informationLOWELL. MICHIGAN, THURSDAY, AUGUST 16, Specialty Co. Sells to Hudson Mfg. Company. Ada Farmer Dies When Boat Tips
K N» N - K V XXXV - N 22 N 5 V N N q N 0 " " - x Q- & V N -" Q" 24-?? 2530 35-6 9025 - - K ; ) ) ( _) ) N K : ; N - K K ) q ( K N ( x- 89830 z 7 K $2000 ( - N " " K : K 89335 30 4 V 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 informationDesign and Analysis of Algorithms
Design and Analysis of Algorithms CSE 5311 Lecture 5 Divide and Conquer: Fast Fourier Transform Junzhou Huang, Ph.D. Department of Computer Science and Engineering CSE5311 Design and Analysis of Algorithms
More informationMATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS
MATH 1372, SECTION 33, MIDTERM 3 REVIEW ANSWERS 1. We have one theorem whose conclusion says an alternating series converges. We have another theorem whose conclusion says an alternating series diverges.
More informationCopyright 2010 Pearson Education, Inc. Publishing as Prentice Hall.
.1 Limits of Sequences. CHAPTER.1.0. a) True. If converges, then there is an M > 0 such that M. Choose by Archimedes an N N such that N > M/ε. Then n N implies /n M/n M/N < ε. b) False. = n does not converge,
More informationComputational Methods CMSC/AMSC/MAPL 460
Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Balaji Vasan Srinivasan Dept of Computer Science Several slides from Prof Healy s course at UMD Last time: Fourier analysis F(t) = A 0 /2 + A
More information