22.09.2015 Views

of Microprocessors

Musical-Applications-of-Microprocessors-2ed-Chamberlin-H-1987

Musical-Applications-of-Microprocessors-2ed-Chamberlin-H-1987

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

448 MUSICAL ApPLICATIONS OF MICROPROCESSORS<br />

Fourier Multipliers Expressed as<br />

Arguments <strong>of</strong> the Function<br />

W(I) ~ COS(21TI/N), SIN(21TIIN)<br />

N = 16<br />

Time<br />

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15<br />

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0<br />

1 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15<br />

2 0 2 4 6 8 10 12 14 0 2 4 6 8 10 12 14<br />

3 0 3 6 9 12 15 2 5 8 11 14 1 4 7 10 13<br />

4 0 4 8 12 0 4 8 12 0 4 8 12 0 4 8 12<br />

F<br />

5 0 5 10 15 4 9 14 3 8 13 2 7 12 1 6 11 r<br />

6 0 6 12 2 8 14 4 10 0 6 12 2 8 14 4 10 e<br />

7 0 7 14 5 12 3 10 1 8 15 6 13 4 11 2 9 q<br />

8 0 8 0 8 0 8 0 8 0 8 0 8 0 8 0 8 u<br />

9 0 9 2 11 4 13 6 15 8 1 10 3 12 5 14 7 e<br />

10 0 10 4 14 8 2 12 6 0 10 4 14 8 2 12 6 n<br />

11 0 11 6 1 12 7 2 13 8 3 14 9 4 15 10 5 c<br />

12 0 12 8 4 0 12 8 4 0 12 8 4 0 12 8 4 Y<br />

13 0 13 10 7 4 1 14 11 8 5 2 15 12 9 6 3<br />

14 0 14 12 10 8 6 4 2 0 14 12 10 8 6 4 2<br />

15 0 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1<br />

Fig. 13-13. Fourier multipliers expressed as arguments <strong>of</strong> the function W(I) =<br />

cos(21TJIN),sin(21TIIN). N = 16.<br />

consideration which is constant throughout the FFT computation. If I ever<br />

exceeds the range <strong>of</strong> a to N -1, it is customary to use the principal value,<br />

which is simply I mod N. Finally, as if nothing is sacred, even the time<br />

samples are considered to be complex numbers with a real part and an<br />

imaginary part.<br />

As a result <strong>of</strong> going complex, the number <strong>of</strong> harmonics is equal to the<br />

number <strong>of</strong> complex samples. Thus, the forward FFT takes a set <strong>of</strong> N complex<br />

time samples and generates a set <strong>of</strong> N complex harmonics. As we will see<br />

later, all <strong>of</strong> the'intetmediate arithmetic in the FFT algorithm is also complex.<br />

In the real world, which includes audio, the time samples are, <strong>of</strong><br />

course, real, which means that the imaginary parts are zero. The corresponding<br />

spectrum has only (N/2)+ 1 unique harmonics; the remaining (N/2)-1<br />

are simply the complex conjugate (the signs <strong>of</strong> the sine components are reversed)<br />

<strong>of</strong> the others less de and Nyquist frequency components. Thus, in actual use<br />

with audio samples, one has N samples and (N/2)+ 1 harmonics, just as<br />

before. Although the preceding implies considerable waste in the computation,<br />

a method <strong>of</strong> eliminating it will be described later.

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

Saved successfully!

Ooh no, something went wrong!