English

Deterministic Sparse Sublinear FFT with Improved Numerical Stability

Numerical Analysis 2021-03-09 v2 Numerical Analysis

Abstract

In this paper we extend the deterministic sublinear FFT algorithm in Plonka et al. (2018) for fast reconstruction of MM-sparse vectors x{\mathbf x} of length N=2JN= 2^J, where we assume that all components of the discrete Fourier transform x^=FNx\hat{\mathbf x}= {\mathbf F}_{N} {\mathbf x} are available. The sparsity of x{\mathbf x} needs not to be known a priori, but is determined by the algorithm. If the sparsity MM is larger than 2J/22^{J/2}, then the algorithm turns into a usual FFT algorithm with runtime O(NlogN){\mathcal O}(N \log N). For M2<NM^{2} < N, the runtime of the algorithm is O(M2logN){\mathcal O}(M^2 \, \log N). 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 M>20M>20 because of numerical instabilities, the modified algorithm is still numerically stable for M=200M=200.

Keywords

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