Nyquist does not have a special data type corresponding to a sequence of FFT frames. This would be nice, but it would require creating a large set of operations suitable for processing frame sequences. Another approach, and perhaps the most "pure" would be to convert a single sound into a multichannel sound, with one channel per bin of the FFT.
Instead, Nyquist violates its "pure" functional model and resorts to objects
for FFT processing. A sequence of frames is represented by an XLISP object. Whenever you
send the selector
:next to the object, you get back either NIL, indicating the
end of the sequence, or you get an array of FFT coefficients.
The Nyquist function
snd-fft (mnemonic, isn't it?) returns one of the frame sequence
generating objects. You can pass any frame sequence generating object to another function,
snd-ifft, and turn the sequence back into audio.
snd-ifft, you can create all sorts of interesting processes. The main
idea is to create intermediate objects that both accept and generate sequences of frames.
These objects can operate on the frames to implement the desired spectral-domain
processes. Examples of this can be found in the file
which is part of the standard Nyquist release. The documentation for
snd-fft(sound, length, skip, window)[SAL]
(snd-fft sound length skip window)[LISP]
FLONUMs. The function modifies the sound, violating the normal rule that sounds are immutable in Nyquist, so it is advised that you copy the sound using
snd-copyif there are any other references to sound. The length of the FFT is specified by length, a
FIXNUM(integer) which must be a power of 2. After each FFT, the sound is advanced by skip samples, also of type
FIXNUM. Overlapping FFTs, where skip is less than length, are allowed. If window is not
NIL, it must be a sound. The first length samples of window are multiplied by length samples of sound before performing the FFT. When there are no more samples in sound to transform, this function returns
NIL. The coefficients in the returned array, in order, are the DC coefficient, the first real, the first imaginary, the second real, the second imaginary, etc. The last array element corresponds to the real coefficient at the Nyquist frequency.
snd-ifft(time, srate, iterator, skip, window)[SAL]
(snd-ifft time srate iterator skip window)[LISP]
(local-to-global 0). The sample rate is given by srate. Typically, this would be
*sound-srate*, but it might also depend upon the sample rate of the sound from which the spectral frames were derived. To obtain each frame, the function sends the message
:nextto the iterator object, using XLISP's primitives for objects and message passing. The object should return an array in the same format as obtained from
snd-fft, and the object should return
NILwhen the end of the sound is reached. After each frame is inverse transformed into the time domain, it is added to the resulting sound. Each successive frame is added with a sample offset specified by skip relative to the previous frame. This must be an integer greater than zero and less than the frame (FFT) size. If window is not
NIL, it must be a sound. This window signal is multiplied by the inverse transformed frame before the frame is added to the output sound. The length of each frame should be the same power of 2. The length is implied by the first array returned by iterator, so it does not appear as a parameter. This length is also the number of samples used from window. Extra samples are ignored, and window is padded with zeros if necessary, so be sure window is the right length. The resulting sound is computed on demand as with other Nyquist sounds, so
:nextmessages are sent to iterator only when new frames are needed. One should be careful not to reuse or modify iterator once it is passed to