Almost-sure identifiability of multidimensional harmonic retrieval

Two-dimensional (2-D) and, more generally, multidimensional harmonic retrieval is of interest in a variety of applications, including transmitter localization and joint time and frequency offset estimation in wireless communications. The associated identifiability problem is key in understanding the fundamental limitations of parametric methods in terms of the number of harmonics that can be resolved for a given sample size. Consider a mixture of 2-D exponentials, each parameterized by amplitude, phase, and decay rate plus frequency in each dimension. Suppose that I equispaced samples are taken along one dimension and, likewise, J along the other dimension. We prove that if the number of exponentials is less than or equal to roughly IJ/4, then, assuming sampling at the Nyquist rate or above, the parameterization is almost surely identifiable. This is significant because the best previously known achievable bound was roughly (I+J)/2. For example, consider I=J=32; our result yields 256 versus 32 identifiable exponentials. We also generalize the result to N dimensions, proving that the number of exponentials that can be resolved is proportional to total sample size.

[1]  Arye Nehorai,et al.  Passive localization of near-field sources by path following , 1994, IEEE Trans. Signal Process..

[2]  Nikos D. Sidiropoulos,et al.  On communication diversity for blind identifiability and the uniqueness of low-rank decomposition of N-way arrays , 2000, 2000 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings (Cat. No.00CH37100).

[3]  Nikos D. Sidiropoulos,et al.  Identifiability results for blind beamforming in incoherent multipath with small delay spread , 2001, IEEE Trans. Signal Process..

[4]  V. Pisarenko The Retrieval of Harmonics from a Covariance Function , 1973 .

[5]  R. Lorentz,et al.  Multivariate Birkhoff Interpolation , 1992 .

[6]  Yingbo Hua,et al.  High resolution imaging of continuously moving object using stepped frequency radar , 1994, Signal Process..

[7]  Yingbo Hua,et al.  On rank of block Hankel matrix for 2-D frequency detection and estimation , 1996, IEEE Trans. Signal Process..

[8]  Nikos D. Sidiropoulos,et al.  Generalizing Carathéodory's uniqueness of harmonic parameterization to N dimensions , 2001, IEEE Trans. Inf. Theory.

[9]  J. Kruskal Three-way arrays: rank and uniqueness of trilinear decompositions, with application to arithmetic complexity and statistics , 1977 .

[10]  Petre Stoica,et al.  Parameter identifiability problem in signal processing , 1994 .

[11]  Emad S. Ebbini Region-adaptive motion tracking of speckle imagery , 2000, 2000 IEEE International Conference on Acoustics, Speech, and Signal Processing. Proceedings (Cat. No.00CH37100).

[12]  N. Sidiropoulos,et al.  On the uniqueness of multilinear decomposition of N‐way arrays , 2000 .

[13]  Petre Stoica,et al.  MUSIC, maximum likelihood and Cramer-Rao bound , 1988, ICASSP-88., International Conference on Acoustics, Speech, and Signal Processing.

[14]  Petre Stoica,et al.  Introduction to spectral analysis , 1997 .

[15]  Nikos D. Sidiropoulos,et al.  Parallel factor analysis in sensor array processing , 2000, IEEE Trans. Signal Process..

[16]  Josef A. Nossek,et al.  3-D unitary ESPRIT for joint 2-D angle and carrier estimation , 1997, 1997 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[17]  Michael D. Zoltowski,et al.  Real-time frequency and 2-D angle estimation with sub-Nyquist spatio-temporal sampling , 1994, IEEE Trans. Signal Process..

[18]  Josef A. Nossek,et al.  Simultaneous Schur decomposition of several nonsymmetric matrices to achieve automatic pairing in multidimensional harmonic retrieval problems , 1998, IEEE Trans. Signal Process..

[19]  Josef A. Nossek,et al.  Efficient high-resolution 3-D channel sounding , 1998, VTC '98. 48th IEEE Vehicular Technology Conference. Pathway to Global Wireless Revolution (Cat. No.98CH36151).

[20]  C. Carathéodory,et al.  Über den zusammenhang der extremen von harmonischen funktionen mit ihren koeffizienten und über den picard-landau’schen satz , 1911 .

[21]  T. Sauer,et al.  On multivariate Hermite interpolation , 1995 .

[22]  Arogyaswami Paulraj,et al.  Estimation of multipath parameters in wireless communications , 1998, IEEE Trans. Signal Process..

[23]  Nikos D. Sidiropoulos,et al.  Blind PARAFAC receivers for DS-CDMA systems , 2000, IEEE Trans. Signal Process..

[24]  Tomas Sauer,et al.  Polynomial interpolation in several variables , 2000, Adv. Comput. Math..

[25]  Jian Li,et al.  An efficient algorithm for two-dimensional frequency estimation , 1996, Multidimens. Syst. Signal Process..

[26]  Michael D. Zoltowski,et al.  Closed-form 2-D angle estimation with rectangular arrays in element space or beamspace via unitary ESPRIT , 1996, IEEE Trans. Signal Process..

[27]  Joe Brewer,et al.  Kronecker products and matrix calculus in system theory , 1978 .