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Ω(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 Ω(k log (n/k) / log log n) signal samples, even if it is allowed to perform adaptive sampling.
Supplemental Material
- R. Agrawal, C. Faloutsos, and A. Swami. Efficient similarity search in sequence databases. Int. Conf. on Foundations of Data Organization and Algorithms, pages 69--84, 1993. Google ScholarDigital Library
- A. Akavia, S. Goldwasser, and S. Safra. Proving hard-core predicates using list decoding. Annual Symposium on Foundations of Computer Science, 44:146--159, 2003. Google ScholarDigital Library
- A. Akavia. Deterministic sparse Fourier approximation via fooling arithmetic progressions. COLT, pages 381--393, 2010.Google Scholar
- A. Chandrakasan, V. Gutnik, and T. Xanthopoulos. Data driven signal processing: An approach for energy efficient computing. International Symposium on Low Power Electronics and Design, 1996. Google ScholarDigital Library
- E. Candes, J. Romberg, and T. Tao. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. IEEE Transactions on Information Theory, 52:489--509, 2006. Google ScholarDigital Library
- Thomas Cover and Joy Thomas.Elements of Information Theory. Wiley Interscience, 1991. Google ScholarDigital Library
- D. Donoho. Compressed sensing. IEEE Transactions on Information Theory, 52(4):1289--1306, 2006. Google ScholarDigital Library
- I. Daubechies, O. Runborg, and J. Zou. A sparse spectral method for homogenization multiscale problems. Multiscale Model. Sim., 6(3):711--740, 2007.Google ScholarCross Ref
- A. Gilbert, S. Guha, P. Indyk, M. Muthukrishnan, and M. Strauss. Near-optimal sparse Fourier representations via sampling.STOC, 2002. Google ScholarDigital Library
- O. Goldreich and L. Levin. A hard-corepredicate for allone-way functions. STOC, pages 25--32, 1989. Google ScholarDigital Library
- Anna C. Gilbert, Yi Li, Ely Porat, and Martin J. Strauss. Approximate sparse recovery: optimizing time and measurements. In STOC, pages 475--484, 2010. Google ScholarDigital Library
- A. Gilbert, M. Muthukrishnan, and M. Strauss. Improved time bounds for near-optimal space Fourier representations.SPIE Conference, Wavelets, 2005.Google Scholar
- A.C. Gilbert, M.J. Strauss, and J. A. Tropp. A tutorial on fast Fourier sampling. Signal Processing Magazine, 2008.Google ScholarCross Ref
- H. Hassanieh, P. Indyk, D. Katabi, and E. Price. sFFT: Sparse Fast Fourier Transform. http://groups.csail.mit.edu/netmit/sFFT/, 2012.Google Scholar
- H. Hassanieh, P. Indyk, D. Katabi, and E. Price. Simple and practical algorithm for sparse Fourier transform. SODA, 2012. Google ScholarDigital Library
- Juha Heiskala and John Terry, Ph.D. OFDM Wireless LANs: A Theoretical and Practical Guide. Sams, Indianapolis, IN, USA, 2001. Google ScholarDigital Library
- P. Indyk, E. Price, and D. P. Woodruff. On the power of adaptivity in sparse recovery. FOCS, 2011. Google ScholarDigital Library
- M. A. Iwen. Combinatorial sublinear-time Fourier algorithms.Foundations of Computational Mathematics, 10:303--338, 2010. Google ScholarDigital Library
- J. Kahn, G. Kalai, and N. Linial. The influence of variables on boolean functions. FOCS, 1988. Google ScholarDigital Library
- E. Kushilevitz and Y. Mansour. Learning decision trees using the Fourier spectrum. STOC, 1991. Google ScholarDigital Library
- N. Linial, Y. Mansour, and N. Nisan. Constant depth circuits, Fourier transform, and learnability. Journal of the ACM (JACM), 1993. Google ScholarDigital Library
- Mengda Lin, A. P. Vinod, and Chong Meng Samson See. A new exible filter bank for low complexity spectrum sensing in cognitive radios. Journal of Signal Processing Systems, 62(2):205--215, 2011. Google ScholarDigital Library
- Y. Mansour. Randomized interpolation and approximation of sparse polynomials. ICALP, 1992. Google ScholarDigital Library
- G. Marsaglia. Evaluating the normal distribution.Journal of Statistical Software, 11(4):1--7, 2004.Google ScholarCross Ref
- A. Mueen, S. Nath, and J. Liu. Fast approximate correlation for massive time-series data. InProceedings of the 2010 international conference on Management of data, pages 171--182. ACM, 2010. Google ScholarDigital Library
- R. O'Donnell. Some topics in analysis of boolean functions (tutorial). STOC, 2008. Google ScholarDigital Library
- E. Price and D. P. Woodruff. (1+ε)-approximate sparse recovery. FOCS, 2011. Google ScholarDigital Library
Index Terms
- Nearly optimal sparse fourier transform
Recommendations
(Nearly) sample-optimal sparse Fourier transform
SODA '14: Proceedings of the twenty-fifth annual ACM-SIAM symposium on Discrete algorithmsWe consider the problem of computing a k-sparse approximation to the discrete Fourier transform of an n-dimensional signal. Our main result is a randomized algorithm that computes such an approximation using O(k log n(log log n)O(1)) signal samples in ...
Nearly Optimal Deterministic Algorithm for Sparse Walsh-Hadamard Transform
For every fixed constant α > 0, we design an algorithm for computing the k-sparse Walsh-Hadamard transform (i.e., Discrete Fourier Transform over the Boolean cube) of an N-dimensional vector x ∈ RN in time k1 + α(log N)O(1). Specifically, the algorithm ...
Sparse fourier transform in any constant dimension with nearly-optimal sample complexity in sublinear time
STOC '16: Proceedings of the forty-eighth annual ACM symposium on Theory of ComputingWe consider the problem of computing a k-sparse approximation to the Fourier transform of a length N signal. Our main result is a randomized algorithm for computing such an approximation (i.e. achieving ℓ2/ℓ2 sparse recovery guarantees using Fourier ...
Comments