DFT Matrix
DFT Matrix
DFT Matrix
In applied mathematics, a DFT matrix is an expression of a discrete Fourier transform (DFT) as a transformation matrix, which
can be applied to a signal through matrix multiplication.
Definition
An N-point DFT is expressed as the multiplication , where is the original input signal, is the N-by-N square DFT
matrix, and is the DFT of the signal.
where is a primitive Nth root of unity in which . We can avoid writing large exponents for using the fact
that for any exponent we have the identity This is the Vandermonde matrix for the roots of unity, up to the
normalization factor. Note that the normalization factor in front of the sum ( ) and the sign of the exponent in ω are merely
conventions, and differ in some treatments. All of the following discussion applies regardless of the convention, with at most minor
adjustments. The only important thing is that the forward and inverse transforms have opposite-sign exponents, and that the
product of their normalization factors be 1/N. However, the choice here makes the resulting DFT matrix unitary, which is
convenient in many circumstances.
Fast Fourier transform algorithms utilize the symmetries of the matrix to reduce the time of multiplying a vector by this matrix,
from the usual . Similar techniques can be applied for multiplications by matrices such as Hadamard matrix and the Walsh
matrix.
Examples
Two-point
The two-point DFT is a simple case, in which the first entry is the DC (sum) and the second entry is the AC (difference).
The first row performs the sum, and the second row performs the difference.
Four-point
Eight-point
The first non-trivial integer power of two case is for eight points:
where
(Note that .)
The following image depicts the DFT as a matrix multiplication, with elements of the matrix depicted by samples of complex
exponentials:
The real part (cosine wave) is denoted by a solid line, and the imaginary part (sine wave) by a dashed line.
The top row is all ones (scaled by for unitarity), so it "measures" the DC component in the input signal. The next row is
eight samples of negative one cycle of a complex exponential, i.e., a signal with a fractional frequency of −1/8, so it "measures"
how much "strength" there is at fractional frequency +1/8 in the signal. Recall that a matched filter compares the signal with a time
reversed version of whatever we're looking for, so when we're looking for fracfreq. 1/8 we compare with fracfreq. −1/8 so that is
why this row is a negative frequency. The next row is negative two cycles of a complex exponential, sampled in eight places, so it
has a fractional frequency of −1/4, and thus "measures" the extent to which the signal has a fractional frequency of +1/4.
The following summarizes how the 8-point DFT works, row by row, in terms of fractional frequency:
0 measures how much DC is in the signal
−1/8 measures how much of the signal has a fractional frequency of +1/8
−1/4 measures how much of the signal has a fractional frequency of +1/4
−3/8 measures how much of the signal has a fractional frequency of +3/8
−1/2 measures how much of the signal has a fractional frequency of +1/2
−5/8 measures how much of the signal has a fractional frequency of +5/8
−3/4 measures how much of the signal has a fractional frequency of +3/4
−7/8 measures how much of the signal has a fractional frequency of +7/8
Equivalently the last row can be said to have a fractional frequency of +1/8 and thus measure how much of the signal has a
fractional frequency of −1/8. In this way, it could be said that the top rows of the matrix "measure" positive frequency content in
the signal and the bottom rows measure negative frequency component in the signal.
Unitary transform
The DFT is (or can be, through appropriate selection of scaling) a unitary transform, i.e., one that preserves energy. The
appropriate choice of scaling to achieve unitarity is , so that the energy in the physical domain will be the same as the
energy in the Fourier domain, i.e., to satisfy Parseval's theorem. (Other, non-unitary, scalings, are also commonly used for
computational convenience; e.g., the convolution theorem takes on a slightly simpler form with the scaling shown in the discrete
Fourier transform article.)
Other properties
For other properties of the DFT matrix, including its eigenvalues, connection to convolutions, applications, and so on, see the
discrete Fourier transform article.
See also
Multidimensional transform
Clock and shift matrices
References
The Transform and Data Compression Handbook by P. C. Yip, K. Ramamohan Rao (https://www.amazon.com/g
p/reader/0849336929) – See chapter 2 for a treatment of the DFT based largely on the DFT matrix
External links
Fourier Operator and Decimation In Time (DIT) (http://wearcam.org/ece431/course_material/fourierop_and_dit.ht
m)