SSP'05 -- Proceedings

	IEEE/SP 13th workshop on Statistical Signal Processing July, 17-20, 2005 - Bordeaux - France

Information regarding the paper

Title

Bernoulli-Gaussian spectral analysis of unevenly spaced astrophysical data

Author(s)

Sébastien Bourguignon	Laboratoire d'Astrophysique de l'OMP, Université Paul Sabatier, Toulouse
Hervé Carfantan	Laboratoire d'Astrophysique de l'OMP, Université Paul Sabatier, Toulouse

Get the paper in PDF format

To obtain Acrobat Reader (version 5 minimum required) necessary to his read.

Abstract

We address the problem of line spectra detection and estimation from astrophysical data. As observations generally suffer sampling irregularities, false peaks may appear in the Fourier spectrum. We propose a linear spectral model with an arbitrarily large number of fixed frequencies and search for a sparse solution by modelling the spectrum as a Bernoulli-Gaussian process. The use of Markov Chain Monte Carlo methods to compute the posterior mean estimate is discussed in the unsupervised framework. The original work by Cheng et al. is modified to account for specificities of the spectral analysis problem. Simulations reveal the efficiency of the method and its relevance to the astrophysical frequency detection context is emphasized. Finally, an application to astrophysical data is presented.