<< Chapter < Page Chapter >> Page >

Calculation of the fourier transform and fourier series using the fft

Most theoretical and mathematical analysis of signals and systems use the Fourier series, Fourier transform, Laplace transform, discrete-timeFourier transform (DTFT), or the z-transform, however, when we want to actually evaluate transforms, we calculate values at sample frequencies.In other words, we use the discrete Fourier transform (DFT) and, for efficiency, usually evaluate it with the FFT algorithm. An importantquestion is how can we calculate or approximately calculate these symbolic formula-based transforms with our practical finite numerical tool. Itwould certainly seem that if we wanted the Fourier transform of a signal or function, we could sample the function, take its DFT with the FFT, andhave some approximation to samples of the desired Fourier transform. We saw in the previous section that it is, in fact, possible provided somecare is taken.

Summary

For the signal that is a function of a continuous variable we have

FT: f ( t ) F ( ω ) DTFT: f ( T n ) 1 T F p ( ω ) = 1 T F ( ω + 2 π / T ) DFT: f p ( T n ) 1 T F p ( Δ k ) for Δ T N = 2 π

For the signal that is a function of a discrete variable we have

DTFT: h ( n ) H ( ω ) DFT: h p ( n ) H ( Δ k ) for Δ N = 2 π

For the periodic signal of a continuous variable we have

FS: g ˜ ( t ) C ( k ) DFT: g ˜ ( T n ) N C p ( k ) for T N = P

For the sampled bandlimited signal we have

Sinc: f ( t ) f ( T n ) f ( t ) = n f ( T n ) sinc ( 2 π t / T - π n ) if F ( ω ) = 0 for | ω | > 2 π / T

These formulas summarize much of the relations of the Fourier transforms of sampled signals and how they might be approximately calculate with theFFT. We next turn to the use of distributions and strings of delta functions as tool to study sampling.

Sampling functions — the shah function

Th preceding discussions used traditional Fourier techniques to develop sampling tools. If distributions or delta functions are allowed, theFourier transform will exist for a much larger class of signals. One should take care when using distributions as if they were functions but itis a very powerful extension.

There are several functions which have equally spaced sequences of impulses that can be used as tools in deriving a sampling formula. Theseare called “pitch fork" functions, picket fence functions, comb functions and shah functions. We start first with a finite length sequence to beused with the DFT. We define

⨿ M ( n ) = m = 0 L - 1 δ ( n - M m )

where N = L M .

D F T { ⨿ M ( n ) } = n = 0 N - 1 m = 0 L - 1 δ ( n - M m ) e - j 2 π n k / N
= m = 0 L - 1 n = 0 N - 1 δ ( n - M m ) e - j 2 π n k / N
= m = 0 L - 1 e - j 2 π M m k / N = m = 0 L - 1 e - j 2 π m k / L
= { L < k > L = 0 0 otherwise
= L l = 0 M - 1 δ ( k - L l ) = L ⨿ L ( k )

For the DTFT we have a similar derivation:

D T F T { ⨿ M ( n ) } = n = - m = 0 L - 1 δ ( n - M m ) e - j ω n
= m = 0 L - 1 n = - δ ( n - M m ) e - j ω n
= m = 0 L - 1 e - j ω M m
= { L ω = k 2 π / M 0 otherwise
= l = 0 M - 1 δ ( ω - 2 π l / M l ) = K ⨿ 2π/M ( ω )

where K is constant.

An alternate derivation for the DTFT uses the inverse DTFT.

I D T F T { ⨿ 2 π / M ( ω ) } = 1 2 π - π π ⨿ 2 π / M ( ω ) e j ω n d ω
= 1 2 π - π π l δ ( ω - 2 π l / M ) e j ω n d ω
= 1 2 π l - π π δ ( ω - 2 π l / M ) e j ω n d ω
= 1 2 π l = 0 M - 1 e 2 π l n / M = { M / 2 π n = M 0 otherwise
= ( M ) ⨿ 2π/M ( ω )

Therefore,

⨿ M ( n ) ( M ) ⨿ 2π/T ( ω )

For regular Fourier transform, we have a string of impulse functions in both the time and frequency. This we see from:

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Brief notes on signals and systems. OpenStax CNX. Sep 14, 2009 Download for free at http://cnx.org/content/col10565/1.7
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Brief notes on signals and systems' conversation and receive update notifications?

Ask