Reconstructing signals from noisy data with unknown signal and noise covariance
arXiv:1107.2384 · doi:10.1103/PhysRevE.84.041118
Abstract
We derive a method to reconstruct Gaussian signals from linear measurements with Gaussian noise. This new algorithm is intended for applications in astrophysics and other sciences. The starting point of our considerations is the principle of minimum Gibbs free energy which was previously used to derive a signal reconstruction algorithm handling uncertainties in the signal covariance. We extend this algorithm to simultaneously uncertain noise and signal covariances using the same principles in the derivation. The resulting equations are general enough to be applied in many different contexts. We demonstrate the performance of the algorithm by applying it to specific example situations and compare it to algorithms not allowing for uncertainties in the noise covariance. The results show that the method we suggest performs very well under a variety of circumstances and is indeed qualitatively superior to the other methods in cases where uncertainty in the noise covariance is present.
12 pages, 6 figures; 1D example added; accepted for publication in Phys. Rev. E
References in corpus (2)
Cited by in corpus (18)
- An improved map of the Galactic Faraday sky
- Estimating extragalactic Faraday rotation
- BLOBCAT: Software to Catalogue Flood-Filled Blobs in Radio Images of Total Intensity and Linear Polarization
- Reconstruction of Gaussian and log-normal fields with spectral smoothness
- NIFTY - Numerical Information Field Theory - a versatile Python library for signal inference
- IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays
- RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
- Efficient Wiener filtering without preconditioning
- Faraday Tomography of the North Polar Spur: Constraints on the distance to the Spur and on the Magnetic Field of the Galaxy
- Information theory for fields
- Constraining Magnetic Fields in the Circumgalactic Medium
- Faraday synthesis: The synergy of aperture and rotation measure synthesis
- A compiled catalog of rotation measures of radio point sources
- Information field dynamics for simulation scheme construction
- Improving stochastic estimates with inference methods: calculating matrix diagonals
- The XENON100 exclusion limit without considering Leff as a nuisance parameter
- A path-integral approach to Bayesian inference for inverse problems using the semiclassical approximation
- Improving self-calibration