DFT-Based Filterbanks 1.1 2005/02/04 09:35:40 US/Central 2005/04/28 14:43:51 GMT-5 Douglas L. Jones dl-jones@uiuc.edu Douglas L. Jones dl-jones@uiuc.edu Charlet Reedstrom charlet@rice.edu Kyle Clarkson kclarks@rice.edu Multirate Signal Processing One common application of multirate processing arises in multirate, multi-channel filter banks ().

One application is separating frequency-division-multiplexed channels. If the filters are narrowband, the output channels can be decimated without significant aliasing. Such structures are especially attractive when they can be implemented efficiently. For example, if the filters are simply frequency modulated (by 2 k L n ) versions of each other, they can be efficiently implemented using FFTs! Furthermore, there are classes of filters called perfect reconstruction filters which are of finite length but from which the signal can be reconstructed exactly (using all M channels), even though the output of each channel experiences aliasing in the decimation step. These types of filterbanks have received a lot of research attention, culminating in wavelet theory and techniques.

Uniform DFT Filter Banks Suppose we wish to split a digital input signal into N frequency bands, uniformly spaced at center frequencies ω k 2 k N , for 0 k N 1 . Consider also a lowpass filter h n , H ω 1 ω N 0 . Bandpass filters can be constructed which have the frequency response H k ω H ω 2 k N from h k n h n 2 k n N The output of the kth bandpass filter is simply (assume h n are FIR) x n h k n m 0 M 1 x n m h m 2 k m N y k n This looks suspiciously like a DFT, except that M N , in general. However, if we fix M N , then we can compute all y k n outputs simultaneously using an FFT of x n m h m : The kth FFT frequency output y k n ! So the cost of computing all of these filter banks outputs is O N N , rather than N 2 , per a given n. This is very useful for efficient implementation of transmultiplexors (FDM to TDM). How would we implement this efficiently if we wanted to decimate the individual channels y k n by a factor of N, to their approximate Nyquist bandwidth? Simply step by N time samples between FFTs. Do you expect significant aliasing? If so, how do you propose to combat it? Efficiently? Aliasing should be expected. There are two ways to reduce it: Decimate by less ("oversample" the individual channels) such as decimating by a factor of N 2 . This is efficiently done by time-stepping by the appropriate factor. Design better (and thus longer) filters, say of length L N . These can be efficiently computed by producing only N (every Lth) FFT outputs using simplified FFTs. How might one convert from N input channels into an FDM signal efficiently? ()

Such systems are used throughout the telephone system, satellite communication links, etc. Use an FFT and an inverse FFT for the modulation (TDM to FDM) and demodulation (FDM to TDM), respectively.