Deterministic Sparse Sublinear FFT with Improved Numerical Stability
Abstract
In this paper we extend the deterministic sublinear FFT algorithm in Plonka et al. (2018) for fast reconstruction of -sparse vectors of length , where we assume that all components of the discrete Fourier transform are available. The sparsity of needs not to be known a priori, but is determined by the algorithm. If the sparsity is larger than , then the algorithm turns into a usual FFT algorithm with runtime . For , the runtime of the algorithm is . The proposed modifications of the approach in Plonka et al. (2018) lead to a significant improvement of the condition numbers of the Vandermonde matrices which are employed in the iterative reconstruction. Our numerical experiments show that our modification has a huge impact on the stability of the algorithm. While the algorithm in Plonka et al. (2018) starts to be unreliable for because of numerical instabilities, the modified algorithm is still numerically stable for .
Cite
@article{arxiv.2004.11097,
title = {Deterministic Sparse Sublinear FFT with Improved Numerical Stability},
author = {Gerlind Plonka and Therese von Wulffen},
journal= {arXiv preprint arXiv:2004.11097},
year = {2021}
}
Comments
15 pages, 3 figures