A Deterministic Sparse FFT for Functions with Structured Fourier Sparsity
share
In this paper a deterministic sparse Fourier transform algorithm is presented which breaks the quadratic-in-sparsity runtime bottleneck for a large class of periodic functions exhibiting structured frequency support. These functions include, e.g., the oft-considered set of block frequency sparse functions of the form f(x) = ∑^n_j=1∑^B-1_k=0 c_ω_j + k e^i(ω_j + k)x, {ω_1, ..., ω_n }⊂(-N/2, N/2]∩Z as a simple subclass. Theoretical error bounds in combination with numerical experiments demonstrate that the newly proposed algorithms are both fast and robust to noise. In particular, they outperform standard sparse Fourier transforms in the rapid recovery of block frequency sparse functions of the type above.
READ FULL TEXT