English

Nearly Optimal Sparse Fourier Transform

Data Structures and Algorithms 2012-04-09 v2

Abstract

We consider the problem of computing the k-sparse approximation to the discrete Fourier transform of an n-dimensional signal. We show: * An O(k log n)-time randomized algorithm for the case where the input signal has at most k non-zero Fourier coefficients, and * An O(k log n log(n/k))-time randomized algorithm for general input signals. Both algorithms achieve o(n log n) time, and thus improve over the Fast Fourier Transform, for any k = o(n). They are the first known algorithms that satisfy this property. Also, if one assumes that the Fast Fourier Transform is optimal, the algorithm for the exactly k-sparse case is optimal for any k = n^{\Omega(1)}. We complement our algorithmic results by showing that any algorithm for computing the sparse Fourier transform of a general signal must use at least \Omega(k log(n/k)/ log log n) signal samples, even if it is allowed to perform adaptive sampling.

Keywords

Cite

@article{arxiv.1201.2501,
  title  = {Nearly Optimal Sparse Fourier Transform},
  author = {Haitham Hassanieh and Piotr Indyk and Dina Katabi and Eric Price},
  journal= {arXiv preprint arXiv:1201.2501},
  year   = {2012}
}

Comments

28 pages, appearing at STOC 2012

R2 v1 2026-06-21T20:03:34.770Z