A derivation of the discrete Fourier transform

From Wikipedia, the free encyclopedia
Jump to: navigation, search

In mathematics, computer science, and electrical engineering, the discrete Fourier transform (DFT), occasionally called the finite Fourier transform, is a transform for Fourier analysis of finite-domain discrete-time signals. As with most Fourier analysis, it expresses an input function in terms of a sum of sinusoidal components by determining the amplitude and phase of each component. Unlike the Fourier transform, which operates upon continuous functions assumed to extend to infinity, the DFT operates upon discrete and finite sets of values: the input to the DFT is a finite sequence of real or complex numbers, which makes the DFT ideal for processing information stored in computers. In particular, the DFT is widely employed in signal processing and related fields to analyze the frequencies contained in a sampled signal, to solve partial differential equations, and to perform other operations such as convolutions.

The article discrete Fourier transform presents the definition of the transform, without derivation, as:

X[k] = \sum_{n=0}^{N-1} x[n] \cdot e^{-i \frac{2 \pi}{N} k n} \quad \quad k = 0, \dots, N-1

 

 

 

 

(Eq.1)

Here we take the view that the DFT is motivated by a desire to study continuous functions or waveforms and their continuous Fourier transforms using only a finite amount of data. When the sequence {x[n]} represents a subset of the samples of a waveform x(t), we can model the process that created {x[n]} as applying a window function to x(t), followed by sampling (or vice versa). It is instructive to envision how those operations affect our ability to observe the Fourier transform, X(ƒ). The window function widens every frequency component of X(ƒ) in a way that depends on the type of window used. That effect is called spectral leakage. We can think of it as causing X(ƒ) to blur... thus a loss of resolution. The sampling operation causes the Fourier transform to become periodic. More precisely, what happens is that {x[n]} has no Fourier transform. It is undefined. But using the Poisson summation formula a periodic function of continuous frequency can be constructed from the samples, and it comprises copies of the blurred X(ƒ) repeated at regular multiples of the sampling frequency (Fs = 1/T) and summed together where they overlap (called periodic summation):

X_{1/T}(f)\ \stackrel{\text{def}}{=}\ \sum_{k=-\infty}^{\infty} X(f-k/T)= \sum_{n=-\infty}^{\infty} \underbrace{T\cdot x(nT)}_{x[n]}\cdot e^{-i 2\pi f nT}    [note 1]

 

 

 

 

(Eq.2)

The copies are aliases of the original frequency components. In particular, due to the overlap, aliases can significantly distort the region containing the original X(ƒ) (if Fs is not sufficiently large enough to prevent it).  But if the windowing and sampling are done with sufficient care, the Poisson summation still contains a reasonable semblance of X(ƒ). It is therefore a common practice to compute an arbitrary number of samples (N) of one cycle of the periodic function X_{1/T} :

\underbrace{X_{1/T}\left(\frac{k}{NT}\right)}_{X[k]} = \sum_{n=-\infty}^{\infty} x[n]\cdot e^{-i 2\pi \frac{kn}{N}} \quad \quad k = 0, \dots, N-1

Since the kernel, e^{-i 2\pi \frac{kn}{N}}\,  is N-periodic, it can readily be shown that this is equivalent to the following DFT:

X[k] = \sum_{N} x_N[n]\cdot e^{-i 2\pi \frac{kn}{N}},

 

 

 

 

(Eq.3)

where \textstyle{\sum_{N}} is a summation over any interval of length N, and x_N\, is another periodic summation:

x_N[n]\ \stackrel{\text{def}}{=}\ \sum_{m=-\infty}^{\infty} x[n-mN].

Eq.1 (the standard DFT) is just a simplification of Eq.3 when the x[n] sequence is zero outside the interval [0, N-1].[note 2] But regardless of the duration of the x[n] sequence, the inverse DFT produces the periodic x_N\, sequence. That can be thought of as a consequence of substituting a discrete set of frequencies for the continuous X_{1/T}.

Notes[edit]

  1. ^ With this definition:  \textstyle{\omega = {2 \pi f \over F_s} = 2 \pi f T,}  the right-hand side of Eq.2 becomes \sum_{n=-\infty}^{\infty} x[n]\cdot e^{-i \omega n},  which is a normalized-frequency form of the discrete-time Fourier transform.
  2. ^ The discussion of longer sequences can be found at Sampling the DTFT.