Abstract
The problem of extracting sinusoid signals from noisy observations made at equally spaced times is considered. Eigenanalysis methods, such as Pisarenko's method and the extended Prony method, find the eigenvector with minimum eigenvalue of a suitably chosen matrix and then obtain the complex sinusoids as the roots of the polynomial that has the components of the eigenvector as coefficients. For the sinusoids to be undamped, it is necessary that the roots lie on the unit circle and hence that the eigenvector be conjugate symmetric. It is shown how this symmetry constraint can be incorporated into eigenanalysis estimation methods in a routine way. The practical importance of the constraint is investigated by Cramer-Rao variance-bound calculations and by simulation. The following conclusions are made: (1) The symmetry constraint is straightforward to implement and reduces the amount of computation required. (2) The relative reduction in variance of the constrained over the unconstrained frequency estimators is arbitrarily large for frequencies close together or near a multiple of 4. There are also frequency values however, for which the symmetry constraint gives no reduction in variance at low noise-to-signal ratios. (3) The relative reduction in variance converges to 0 for large sample sizes. (4) The symmetry constraint increases the breakdown noise-to-signal ratios above which the various methods fail to give useful results.