Fixed-point algorithms for frequency estimation and structured low rank approximation

Research output: Contribution to journalArticle

Abstract

We develop fixed-point algorithms for the approximation of structured matrices with rank penalties. In particular we use these fixed-point algorithms for making approximations by sums of exponentials, i.e., frequency estimation. For the basic formulation of the fixed-point algorithm we show that it converges to the solution of a related minimization problem, namely the one obtained by replacing the original objective function with its convex envelope and keeping the structured matrix constraint unchanged.It often happens that this solution agrees with the solution to the original minimization problem, and we provide a simple criterion for when this is true. We also provide more general fixed-point algorithms that can be used to treat the problems of making weighted approximations by sums of exponentials given equally or unequally spaced sampling. We apply the method to the case of missing data, although the above mentioned convergence results do not hold in this case. However, it turns out that the method often gives perfect reconstruction (up to machine precision) in such cases. We also discuss multidimensional extensions, and illustrate how the proposed algorithms can be used to recover sums of exponentials in several variables, but when samples are available only along a curve.

Details

Authors
Organisations
Research areas and keywords

Subject classification (UKÄ) – MANDATORY

  • Mathematics

Keywords

  • Convex envelope, Fenchel conjugate, Fixed-point algorithms, Frequency estimation, General domain Hankel matrices
Original languageEnglish
Pages (from-to)40-65
JournalApplied and Computational Harmonic Analysis
Volume46
Issue number1
Early online date2015 Apr 1
Publication statusPublished - 2019 Jan
Publication categoryResearch
Peer-reviewedYes