The Fast Slepian Transform
For signal processing practitioners, this work makes the Slepian representation computationally feasible, addressing a key bottleneck that previously limited its use.
The paper introduces fast algorithms for approximate projection onto the Slepian basis, achieving complexity comparable to the FFT. It also establishes new nonasymptotic eigenvalue distribution results and demonstrates efficient solutions to least-squares problems in signal processing.
The discrete prolate spheroidal sequences (DPSS's) provide an efficient representation for discrete signals that are perfectly timelimited and nearly bandlimited. Due to the high computational complexity of projecting onto the DPSS basis - also known as the Slepian basis - this representation is often overlooked in favor of the fast Fourier transform (FFT). We show that there exist fast constructions for computing approximate projections onto the leading Slepian basis elements. The complexity of the resulting algorithms is comparable to the FFT, and scales favorably as the quality of the desired approximation is increased. In the process of bounding the complexity of these algorithms, we also establish new nonasymptotic results on the eigenvalue distribution of discrete time-frequency localization operators. We then demonstrate how these algorithms allow us to efficiently compute the solution to certain least-squares problems that arise in signal processing. We also provide simulations comparing these fast, approximate Slepian methods to exact Slepian methods as well as the traditional FFT based methods.