Fast Fourier Transforms on Motorola's Digital Signal Processors
Fast Fourier Transforms on Motorola's Digital Signal Processors
Fast Fourier Transforms on Motorola's Digital Signal Processors
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
X=A+jB<br />
Y=C+jD<br />
X’=A+C+j(D+B)<br />
-<br />
Y’=A-C+j(D-B)<br />
X=A+jC<br />
W<br />
Y=B+jD<br />
6.1.2 Reordering<br />
X’=A+BWr+DWi<br />
+j(C+DWr-BWi)<br />
Y’=A-(BWr+DWi)<br />
-j[C+(DWr-BWi)]<br />
The output order of the Bergland algorithm is slightly<br />
different than the bit-reversed order, and the<br />
twiddle factor required in the calculati<strong>on</strong> is also in<br />
Bergland order. To get this special order, <strong>on</strong>e may<br />
use the following algorithm for doubling the length<br />
of each number sequence:<br />
1. Multiply the sec<strong>on</strong>d entry of the sequence by<br />
two, and make this product the sec<strong>on</strong>d entry of<br />
the new sequence<br />
2. Subtract each n<strong>on</strong>zero entry of the sequence<br />
from twice the product formed in step 1 (these<br />
differences form the rest of the even entries of<br />
the new sequence)<br />
3. Take the odd entries of the new sequence as<br />
the numbers of the original sequence<br />
6-6 MOTOROLA<br />
*<br />
* -<br />
* denotes c<strong>on</strong>jugate<br />
(a) (b)<br />
Figure 6-3 (a) Butterfly of Bergland Algorithm with W = 1<br />
(b) Butterfly of Bergland Algorithm with W ≠ 1<br />
-