Skip to content

Commit 6c69955

Browse files
Unlambderleios
Unlambder
authored andcommitted
Fixed discrete fourier formula (#199)
added a missing "i" in some formulas
1 parent 0af1527 commit 6c69955

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

chapters/FFT/cooley_tukey.md

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -58,11 +58,11 @@ Truth be told, I didn't understand it fully until I discretized real and frequen
5858

5959
In principle, the Discrete Fourier Transform (DFT) is simply the Fourier transform with summations instead of integrals:
6060

61-
$$X_k = \sum_{n=0}^{N-1} x_n \cdot e^{-2 \pi k n / N}$$
61+
$$X_k = \sum_{n=0}^{N-1} x_n \cdot e^{-2 \pi i k n / N}$$
6262

6363
and
6464

65-
$$x_n = \frac{1}{N} \sum_{k=0}^{N-1} X_k \cdot e^{2 \pi k n / N}$$
65+
$$x_n = \frac{1}{N} \sum_{k=0}^{N-1} X_k \cdot e^{2 \pi i k n / N}$$
6666

6767
Where $$X_n$$ and $$x_n$$ are sequences of $$N$$ numbers in frequency and real space, respectively.
6868
In principle, this is no easier to understand than the previous case!
@@ -138,12 +138,12 @@ Butterfly Diagrams show where each element in the array goes before, during, and
138138
As mentioned, the FFT must perform a DFT.
139139
This means that even though we need to be careful about how we add elements together, we are still ultimately performing the following operation:
140140

141-
$$X_k = \sum_{n=0}^{N-1} x_n \cdot e^{-2 \pi k n / N}$$
141+
$$X_k = \sum_{n=0}^{N-1} x_n \cdot e^{-2 \pi i k n / N}$$
142142

143143
However, after shuffling the initial array (by bit reversing or recursive subdivision), we perform the matrix multiplication of the $$e^{-2 \pi k n / N}$$ terms in pieces.
144144
Basically, we split the array into a series of omega values:
145145

146-
$$\omega_N^k = e^{-2 \pi k / N}$$
146+
$$\omega_N^k = e^{-2 \pi i k / N}$$
147147

148148
And at each step, we use the appropriate term.
149149
For example, imagine we need to perform an FFT of an array of only 2 elements.

0 commit comments

Comments
 (0)