14.08.2013 Views

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

SHOW MORE
SHOW LESS

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 />

-

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!