The Discrete Fourier Transform

Similar documents
Fall 2011, EE123 Digital Signal Processing

Finite-length Discrete Transforms. Chapter 5, Sections

Chapter 8. DFT : The Discrete Fourier Transform

Using Difference Equations to Generalize Results for Periodic Nested Radicals

Progression. CATsyllabus.com. CATsyllabus.com. Sequence & Series. Arithmetic Progression (A.P.) n th term of an A.P.

CHAPTER 5 : SERIES. 5.2 The Sum of a Series Sum of Power of n Positive Integers Sum of Series of Partial Fraction Difference Method

Practical Spectral Anaysis (continue) (from Boaz Porat s book) Frequency Measurement

Mapping Radius of Regular Function and Center of Convex Region. Duan Wenxi

Conditional Convergence of Infinite Products

Technical Report: Bessel Filter Analysis

A. Basics of Discrete Fourier Transform

Advanced Higher Formula List

On composite conformal mapping of an annulus to a plane with two holes

Greatest term (numerically) in the expansion of (1 + x) Method 1 Let T

EE123 Digital Signal Processing

Ch3 Discrete Time Fourier Transform

Block-by Block Convolution, FFT/IFFT, Digital Spectral Analysis

2D DSP Basics: 2D Systems

Signal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform

Finite q-identities related to well-known theorems of Euler and Gauss. Johann Cigler

M2.The Z-Transform and its Properties

By the end of this section you will be able to prove the Chinese Remainder Theorem apply this theorem to solve simultaneous linear congruences

THE ANALYTIC LARGE SIEVE

15 Solving the Laplace equation by Fourier method

Einstein Classes, Unit No. 102, 103, Vardhman Ring Road Plaza, Vikas Puri Extn., Outer Ring Road New Delhi , Ph. : ,

EE Midterm Test 1 - Solutions

ADVANCED DIGITAL SIGNAL PROCESSING

FIR Filter Design: Part II

Taylor Transformations into G 2

UNIT V: Z-TRANSFORMS AND DIFFERENCE EQUATIONS. Dr. V. Valliammal Department of Applied Mathematics Sri Venkateswara College of Engineering

MATH Midterm Solutions

2012 GCE A Level H2 Maths Solution Paper Let x,

Solutions. Number of Problems: 4. None. Use only the prepared sheets for your solutions. Additional paper is available from the supervisors.

BINOMIAL THEOREM An expression consisting of two terms, connected by + or sign is called a

BINOMIAL THEOREM NCERT An expression consisting of two terms, connected by + or sign is called a

The Discrete Fourier Transform

On ARMA(1,q) models with bounded and periodically correlated solutions

Math 166 Week-in-Review - S. Nite 11/10/2012 Page 1 of 5 WIR #9 = 1+ r eff. , where r. is the effective interest rate, r is the annual

Minimization of the quadratic test function

Recursion. Algorithm : Design & Analysis [3]

EE422G Homework #13 (12 points)

Chapter 9 Computation of the Discrete. Fourier Transform

IDENTITIES FOR THE NUMBER OF STANDARD YOUNG TABLEAUX IN SOME (k, l)-hooks

Introduction to Digital Signal Processing

Mathematical Description of Discrete-Time Signals. 9/10/16 M. J. Roberts - All Rights Reserved 1

Lesson 5. Chapter 7. Wiener Filters. Bengt Mandersson. r x k We assume uncorrelated noise v(n). LTH. September 2010

HOMEWORK #4 - MA 504

( ) ( ) ( ) ( ) Solved Examples. JEE Main/Boards = The total number of terms in the expansion are 8.

a) The average (mean) of the two fractions is halfway between them: b) The answer is yes. Assume without loss of generality that p < r.

The Pigeonhole Principle 3.4 Binomial Coefficients

KEY. Math 334 Midterm II Fall 2007 section 004 Instructor: Scott Glasgow

Lower Bounds for Cover-Free Families

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

Auchmuty High School Mathematics Department Sequences & Series Notes Teacher Version

Ch 3.4 Binomial Coefficients. Pascal's Identit y and Triangle. Chapter 3.2 & 3.4. South China University of Technology

Electron states in a periodic potential. Assume the electrons do not interact with each other. Solve the single electron Schrodinger equation: KJ =

ELEG3503 Introduction to Digital Signal Processing

is monotonically decreasing function of Ω, it is also called maximally flat at the

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course

( ) 1 Comparison Functions. α is strictly increasing since ( r) ( r ) α = for any positive real number c. = 0. It is said to belong to

physicsandmathstutor.com

Frequency Domain Filtering

Chapter 8 Complex Numbers

DANIEL YAQUBI, MADJID MIRZAVAZIRI AND YASIN SAEEDNEZHAD

Frequency Response of FIR Filters

Lecture 24: Observability and Constructibility

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

Absolute Specifications: A typical absolute specification of a lowpass filter is shown in figure 1 where:

Minimal order perfect functional observers for singular linear systems

Linear time invariant systems

Asymptotic Expansions of Legendre Wavelet

Filter banks. Separately, the lowpass and highpass filters are not invertible. removes the highest frequency 1/ 2and

SHIFTED HARMONIC SUMS OF ORDER TWO

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.

ON CERTAIN CLASS OF ANALYTIC FUNCTIONS

Generalizing the DTFT. The z Transform. Complex Exponential Excitation. The Transfer Function. Systems Described by Difference Equations

x[0] x[1] x[2] Figure 2.1 Graphical representation of a discrete-time signal.

8 The Discrete Fourier Transform (DFT)

ICS141: Discrete Mathematics for Computer Science I

Chapter 7: The z-transform. Chih-Wei Liu

The Z-Transform. (t-t 0 ) Figure 1: Simplified graph of an impulse function. For an impulse, it can be shown that (1)

14th European Signal Processing Conference (EUSIPCO 2006), Florence, Italy, September 4-8, 2006, copyright by EURASIP

SOME ARITHMETIC PROPERTIES OF OVERPARTITION K -TUPLES

Counting Functions and Subsets

DIGITAL SIGNAL PROCESSING LECTURE 3

SVD ( ) Linear Algebra for. A bit of repetition. Lecture: 8. Let s try the factorization. Is there a generalization? = Q2Λ2Q (spectral theorem!

Chapter 2 Systems and Signals

FIR Filter Design: Part I

Discrete-time Fourier transform (DTFT) of aperiodic and periodic signals

Spring 2014, EE123 Digital Signal Processing

Multivector Functions

Adiabatic evolution of the constants of motion in resonance (I)

from definition we note that for sequences which are zero for n < 0, X[z] involves only negative powers of z.

ADDITIONAL INTEGRAL TRANSFORMS

Spectral Analysis. This week in lab. Next classes: 3/26 and 3/28. Your next experiment Homework is to prepare

2D DSP Basics: Systems Stability, 2D Sampling

Signals & Systems Chapter3

The Discrete-Time Fourier Transform (DTFT)

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)

Review of Discrete-time Signals. ELEC 635 Prof. Siripong Potisuk

Transcription:

(7) The Discete Fouie Tasfom The Discete Fouie Tasfom hat is Discete Fouie Tasfom (DFT)? (ote: It s ot DTFT discete-time Fouie tasfom) A liea tasfomatio (mati) Samples of the Fouie tasfom (DTFT) of a apeiodic (with fiite duatio) sequece Etesio of Discete Fouie Seies (DFS) Review: FT, DTFT, FS, DFS Time sigal Tasfom Coeffs. (peiodic/apeiodic) Coeffs. (coti./discete) Aalog apeiodic FT Apeiodic Cotiuous Aalog peiodic FT FS Apeiodic Apeiodic Discete apeiodic DTFT Peiodic Cotiuous Discete peiodic DFS Peiodic Discete Discete fiite-duatio DFT Cotiuous (impulse) Discete CTU EE

(7) The Discete Fouie Tasfom The Discete Fouie Seies Popeties of e jπ, thus e π j -- is peiodic with peiod. (It is essetially cos ad si) : ± ± L -- l, if l m, if l m l m (Pf) (i) If l m, l l (ii) If l m, l l l l -- l Y [] l δ [ l m] m DFS fo peiodic sequeces [ ] [] [ ], peiod + Its DFS epesetatio is defied as follows: Sythesis equatio: j [] [] e [] o Aalysis equatio: [] [] π ote: The tilde i idicates a peiodic sigal. [ ] is peiodic of peiod. CTU EE

(7) The Discete Fouie Tasfom CTU EE 3 Pf) [] [] Pic a ( < ) [] [] [] [] [] [] [] [ ] [] ) ( ) ( + + + + L L That is, [] []. QED Eample: Peiodic Rectagula Pulse Tai [] 4 4 5 si si j K e π π π

(7) The Discete Fouie Tasfom Samplig the Fouie Tasfom Compae two cases: () Peiodic sequece [] [ ] () Fiite duatio sequece [ ] oe peiod of [ ] A apeiodic sequece: jω [] FT ( e ) b? samples jω [] IDFS [] ( e ) π ω Compae: ( t) FT ( jω) samples b? jω [] DTFT ( e ) Eample:, 4, 4 + [] [], othewiase, 5 + 9 + itege CTU EE 4

The Discete Fouie Tasfom CTU EE 5 [] [] ( ) [ ] [ ] [ ] [] [ ] [ ] + + m m m m j m m j j K m m e m e e δ π π ω ω π ω ω ) (Itechage (FT) (Samplig) (IDFS) If [] has fiite legth ad we tae a sufficiet umbe of equally spaced samples of its Fouie Tasfom ( a umbe geate tha o equal to the legth of []), the [] is ecoveable fom [].

The Discete Fouie Tasfom Two ways (equivaletly) to defie DFT: () samples of the DTFT of a fiite duatio sequece [ ] () Mae the peiodic eplica of [ ] [ ] Tae the DFS of [] Pic up oe segmet of [ ] [] DFT peiodic [] DFS [] [ ] oe segmet Popeties of the Discete Fouie Seies -- Simila to those of FT ad z-tasfom Lieaity [] [] [] [] a [] + b [] a [] + b [] Shift m [] [] > [ m] [ ] l [] [ l] Duality Def: [] [] [] [] [] [] [] [ ] () (# ) CTU EE 6

The Discete Fouie Tasfom CTU EE 7 Symmety [] [] [] { } [] [] [ ] ( ) [] { } [] [] [ ] ( ) + j o e Im Re [] [] [ ] ( ) [] { } [] [] [ ] ( ) [] { } j o e Im Re + [] [] [ ] is eal, If. [ ] [ ] [] [ ] [] { } [ ] { } [] { } [ ] { } Im Im Re Re Peiodic Covolutio [] [], ae peiodic sequeces with peiod [ ] [ ] [ ] [ ] [] [] [] [] [ ] 3 l m l l m m

The Discete Fouie Tasfom Discete Fouie Tasfom Defiitio [] : legth, Maig the peiodic eplica: [] [ + ] [] [] [( modulo )] [(( ) ) ] Keep oe segmet (fiite duatio) [], [], othewise This fiite duatio sequece [] Aalysis eq : Sythesis eq : That is, [ ] [(( )) ] is the discete Fouie tasfom (DFT) of [] [] [], [] [], Rema: DFT fomula is the same as DFS fomula. Ideed, may popeties of DFT ae deived fom those of DFS. CTU EE 8

The Discete Fouie Tasfom Popeties of Discete Fouie Tasfom Lieaity [] [] [] [] a legth ma [] + b [] a [] + b [] [, ] Cicula Shift [] [] > m [(( m) ) ] [ ] l [] [(( l) ) ] (Pf) Fom the ight side of the d eq. Km b DFT j m j m [] e [] e [] [(( m) ) ] [ (( m) ) ] [ m] π π IDFS QED Rema: This is cicula shift ot liea shift. (Liea shift of a peiodic sequece cicula shift of a fiite sequece.) CTU EE 9

The Discete Fouie Tasfom Duality [] [] [] [ (( ) ) ], - Symmety Popeties ep [] peiodic cojugate - symmetic e[] { [ (( ) ) ] + [(( ) ) ]}, op [] If ep Re { [] + [ ] }, { [] }, peiodic { [] [ ] }, - - cojugate - atisymmetic - Im{ [] }, [] Re{ [] } [ ] j Im [ ] [] eal, [] [(( )) ], [] [(( ) ) ] { [] } [(( ) ) ] op Re Im { } { [] } Re{ [(( ) ) ]} { [] } Im{ [(( ) ) ]} Re Im { } { [] } [] [(( ) ) ] + [(( ) ) ] ep { } { [] } [] [(( ) ) ] [(( ) ) ] op CTU EE

The Discete Fouie Tasfom Cicula Covolutio Θ 3 [] [] [] m [ m] [ (( m) ) ] [] Θ [] [] [] -poit cicula covolutio Eample: -poit cicula covolutio of two costat sequeces of legth L-poit cicula covolutio of two costat sequeces of legth L CTU EE

The Discete Fouie Tasfom Liea Covolutio Usig DFT hy usig DFT? Thee ae fast DFT algoithms (FFT) How to do it? () Compute the -poit DFT of [ ] ad [ ] [] ad [] () Compute the poduct [] [ ] [ ] 3 (3) Compute the -poit IDFT of 3 [ ] [ ] Poblems: (a) Aliasig (b) Vey log sequece sepaately 3 CTU EE

The Discete Fouie Tasfom Aliasig [], legth L (ozeo values) [], legth P I ode to avoid aliasig, L + P (hat do we mea avoid aliasig? The pecedig pocedue is cicula covolutio but we wat liea covolutio. That is, [ ] 3 equals to the liea covolutio of [] ad [ ] ) CTU EE 3

The Discete Fouie Tasfom [] [] pad with zeos legth pad with zeos legth Itepetatio: (hy call it aliasig?) 3 has a (time domai) badwidth of size L + P [] (That is, the ozeo values of [ ] Theefoe, [ ] 3 ca be at most L + P ) 3 should have at least L + P samples. If the samplig ate is isufficiet, aliasig occus o 3 [ ]. CTU EE 4

The Discete Fouie Tasfom Vey log sequece (FIR filteig) Bloc covolutio Method ovelap ad add Patitio the log sequece ito sectios of shote legth. Fo eample, the filte impulse espose [ ] ealy ifiite. [ + ] [] L, L- Let [ L] whee [], othewise The system (filte) output is a liea covolutio: y h has fiite legth P ad the iput data [ ] [] [] h[] y [ L] whee y [] [] h[] is Rema: The covolutio legth is L + P. That is, the L + P poit DFT is used. y [] has L + P data poits; amog them, (P-) poits should be added to the et sectio. This is called ovelap-add method. (Key: The iput data ae patitioed ito oovelappig sectios the sectio outputs ae ovelapped ad added togethe.) CTU EE 5

The Discete Fouie Tasfom Method ovelap ad save Patitio the log sequece ito ovelappig sectios. Afte computig DFT ad IDFT, thow away some (icoect) outputs. Fo each sectio (legth L, which is also the DFT size), we wat to etai the coect data of legth ( L ( P )) poits Let h [], legth P [], legth L (L>P) The, y [] cotais (P-) icoect poits at the begiig. Theefoe, we divide ito sectios of legth L but each sectio ovelaps the pecedig sectio by (P-) poits. + L P + P, L - [] [ ( ) ( )] This is called ovelap-save method. CTU EE 6