19 research outputs found
Likelihood Analysis of Power Spectra and Generalized Moment Problems
We develop an approach to spectral estimation that has been advocated by
Ferrante, Masiero and Pavon and, in the context of the scalar-valued covariance
extension problem, by Enqvist and Karlsson. The aim is to determine the power
spectrum that is consistent with given moments and minimizes the relative
entropy between the probability law of the underlying Gaussian stochastic
process to that of a prior. The approach is analogous to the framework of
earlier work by Byrnes, Georgiou and Lindquist and can also be viewed as a
generalization of the classical work by Burg and Jaynes on the maximum entropy
method. In the present paper we present a new fast algorithm in the general
case (i.e., for general Gaussian priors) and show that for priors with a
specific structure the solution can be given in closed form.Comment: 17 pages, 4 figure
Frequency-Selective Vandermonde Decomposition of Toeplitz Matrices with Applications
The classical result of Vandermonde decomposition of positive semidefinite
Toeplitz matrices, which dates back to the early twentieth century, forms the
basis of modern subspace and recent atomic norm methods for frequency
estimation. In this paper, we study the Vandermonde decomposition in which the
frequencies are restricted to lie in a given interval, referred to as
frequency-selective Vandermonde decomposition. The existence and uniqueness of
the decomposition are studied under explicit conditions on the Toeplitz matrix.
The new result is connected by duality to the positive real lemma for
trigonometric polynomials nonnegative on the same frequency interval. Its
applications in the theory of moments and line spectral estimation are
illustrated. In particular, it provides a solution to the truncated
trigonometric -moment problem. It is used to derive a primal semidefinite
program formulation of the frequency-selective atomic norm in which the
frequencies are known {\em a priori} to lie in certain frequency bands.
Numerical examples are also provided.Comment: 23 pages, accepted by Signal Processin