cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: FFTW 3.3.5: The 1d Discrete Fourier Transform (DFT) cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: cannam@127:
cannam@127:

cannam@127: Next: , Previous: , Up: What FFTW Really Computes   [Contents][Index]

cannam@127:
cannam@127:
cannam@127: cannam@127:

4.8.1 The 1d Discrete Fourier Transform (DFT)

cannam@127: cannam@127: cannam@127: cannam@127:

The forward (FFTW_FORWARD) discrete Fourier transform (DFT) of a cannam@127: 1d complex array X of size n computes an array Y, cannam@127: where: cannam@127:

.
The backward (FFTW_BACKWARD) DFT computes: cannam@127:
.

cannam@127: cannam@127:

FFTW computes an unnormalized transform, in that there is no coefficient cannam@127: in front of the summation in the DFT. In other words, applying the cannam@127: forward and then the backward transform will multiply the input by cannam@127: n. cannam@127:

cannam@127: cannam@127:

From above, an FFTW_FORWARD transform corresponds to a sign of cannam@127: -1 in the exponent of the DFT. Note also that we use the cannam@127: standard “in-order” output ordering—the k-th output cannam@127: corresponds to the frequency k/n (or k/T, where T cannam@127: is your total sampling period). For those who like to think in terms of cannam@127: positive and negative frequencies, this means that the positive cannam@127: frequencies are stored in the first half of the output and the negative cannam@127: frequencies are stored in backwards order in the second half of the cannam@127: output. (The frequency -k/n is the same as the frequency cannam@127: (n-k)/n.) cannam@127:

cannam@127: cannam@127: cannam@127: cannam@127: cannam@127: