Multichannel blind deconvolution using a generalized Gaussian source model

In this paper, we present an algorithm for the problem of multi-channel blind deconvolution which can adapt to un-known sources with both sub-Gaussian and super-Gaussian probability density distributions using a generalized gaussian source model. We use a state space representation to model the mixer and demixer respectively, and show how the parameters of the demixer can be adapted using a gradient descent algorithm incorporating the natural gradient extension. We also present a learning method for the unknown parameters of the generalized Gaussian source model. The performance of the proposed generalized Gaussian source model on a typical example is compared with those of other algorithm, viz the switching nonlinearity algorithm proposed by Lee et al. [8]. © Association for Scientific Research

Of `Cocktail Parties' and Exoplanets

The characterisation of ever smaller and fainter extrasolar planets requires an intricate understanding of one's data and the analysis techniques used. Correcting the raw data at the 10^-4 level of accuracy in flux is one of the central challenges. This can be difficult for instruments that do not feature a calibration plan for such high precision measurements. Here, it is not always obvious how to de-correlate the data using auxiliary information of the instrument and it becomes paramount to know how well one can disentangle instrument systematics from one's data, given nothing but the data itself. We propose a non-parametric machine learning algorithm, based on the concept of independent component analysis, to de-convolve the systematic noise and all non-Gaussian signals from the desired astrophysical signal. Such a `blind' signal de-mixing is commonly known as the `Cocktail Party problem' in signal-processing. Given multiple simultaneous observations of the same exoplanetary eclipse, as in the case of spectrophotometry, we show that we can often disentangle systematic noise from the original light curve signal without the use of any complementary information of the instrument. In this paper, we explore these signal extraction techniques using simulated data and two data sets observed with the Hubble-NICMOS instrument. Another important application is the de-correlation of the exoplanetary signal from time-correlated stellar variability. Using data obtained by the Kepler mission we show that the desired signal can be de-convolved from the stellar noise using a single time series spanning several eclipse events. Such non-parametric techniques can provide important confirmations of the existent parametric corrections reported in the literature, and their associated results. Additionally they can substantially improve the precision exoplanetary light curve analysis in the future.Comment: ApJ accepte

Fetal electrocardiogram extraction by sequential source separation in the wavelet domain

This work addresses the problem of fetal electrocardiogram extraction using blind source separation (BSS) in the wavelet domain. A new approach is proposed, which is particularly advantageous when the mixing environment is noisy and time-varying, and that is shown, analytically and in simulation, to improve the convergence rate of the natural gradient algorithm. The distribution of the wavelet coefficients of the source signals is then modeled by a generalized Gaussian probability density, whereby in the time-scale domain the problem of selecting appropriate nonlinearities when separating mixtures of both sub- and super-Gaussian signals is mitigated, as shown by experimental results

Neural networks and separation of Cosmic Microwave Background and astrophysical signals in sky maps

The Independent Component Analysis (ICA) algorithm is implemented as a neural network for separating signals of different origin in astrophysical sky maps. Due to its self-organizing capability, it works without prior assumptions on the signals, neither on their frequency scaling, nor on the signal maps themselves; instead, it learns directly from the input data how to separate the physical components, making use of their statistical independence. To test the capabilities of this approach, we apply the ICA algorithm on sky patches, taken from simulations and observations, at the microwave frequencies, that are going to be deeply explored in a few years on the whole sky, by the Microwave Anisotropy Probe (MAP) and by the {\sc Planck} Surveyor Satellite. The maps are at the frequencies of the Low Frequency Instrument (LFI) aboard the {\sc Planck} satellite (30, 44, 70 and 100 GHz), and contain simulated astrophysical radio sources, Cosmic Microwave Background (CMB) radiation, and Galactic diffuse emissions from thermal dust and synchrotron. We show that the ICA algorithm is able to recover each signal, with precision going from 10% for the Galactic components to percent for CMB; radio sources are almost completely recovered down to a flux limit corresponding to $0.7\sigma_{CMB}$, where $\sigma_{CMB}$ is the rms level of CMB fluctuations. The signal recovering possesses equal quality on all the scales larger then the pixel size. In addition, we show that the frequency scalings of the input signals can be partially inferred from the ICA outputs, at the percent precision for the dominant components, radio sources and CMB.Comment: 15 pages; 6 jpg and 1 ps figures. Final version to be published in MNRA

An Introduction to Independent Component Analysis: InfoMax and FastICA algorithms

This paper presents an introduction to independent component analysis (ICA). Unlike principal component analysis, which is based on the assumptions of uncorrelatedness and normality, ICA is rooted in the assumption of statistical independence. Foundations and basic knowledge necessary to understand the technique are provided hereafter. Also included is a short tutorial illustrating the implementation of two ICA algorithms (FastICA and InfoMax) with the use of the Mathematica software