Vandermonde Decomposition of Multilevel Toeplitz Matrices with Application to Multidimensional Super-Resolution
arXiv:1505.02510
Abstract
The Vandermonde decomposition of Toeplitz matrices, discovered by Carathéodory and Fejér in the 1910s and rediscovered by Pisarenko in the 1970s, forms the basis of modern subspace methods for 1D frequency estimation. Many related numerical tools have also been developed for multidimensional (MD), especially 2D, frequency estimation; however, a fundamental question has remained unresolved as to whether an analog of the Vandermonde decomposition holds for multilevel Toeplitz matrices in the MD case. In this paper, an affirmative answer to this question and a constructive method for finding the decomposition are provided when the matrix rank is lower than the dimension of each Toeplitz block. A numerical method for searching for a decomposition is also proposed when the matrix rank is higher. The new results are applied to studying MD frequency estimation within the recent super-resolution framework. A precise formulation of the atomic norm is derived using the Vandermonde decomposition. Practical algorithms for frequency estimation are proposed based on relaxation techniques. Extensive numerical simulations are provided to demonstrate the effectiveness of these algorithms compared to the existing atomic norm and subspace methods.
17 pages, double column, 5 figures, to appear in IEEE Transactions on Information Theory
References in corpus (5)
- On Gridless Sparse Methods for Line Spectral Estimation From Complete and Incomplete Data
- Super-resolution on the Sphere using Convex Optimization
- Continuous Compressed Sensing With a Single or Multiple Measurement Vectors
- Exact Joint Sparse Frequency Recovery via Optimization Methods
- Super-Resolution Radar
Cited by in corpus (5)
- Gridless DOA Estimation and Root-MUSIC for Non-Uniform Arrays
- Joint Localization and Orientation Estimation in Millimeter-Wave MIMO OFDM Systems via Atomic Norm Minimization
- Super-Resolution Radar
- Off-the-grid Recovery of Time and Frequency Shifts with Multiple Measurement Vectors
- Exact Dimensionality Reduction for Partial Line Spectra Estimation Problems