STJan 11, 2008
A new transform for solving the noisy complex exponentials approximation problemPiero Barone
The problem of estimating a complex measure made up by a linear combination of Dirac distributions centered on points of the complex plane from a finite number of its complex moments affected by additive i.i.d. Gaussian noise is considered. A random measure is defined whose expectation approximates the unknown measure under suitable conditions. An estimator of the approximating measure is then proposed as well as a new discrete transform of the noisy moments that allows to compute an estimate of the unknown measure. A small simulation study is also performed to experimentally check the goodness of the approximations.
STMay 24, 2011
A diffusion equation for the density of the ratio of two jointly distributed Gaussian variables and the numerical inversion of Laplace transformPiero Barone
It is shown that the density of the ratio of two random variables with the same variance and joint Gaussian density satisfies a non stationary diffusion equation. Implications of this result for kernel density estimation of the condensed density of the generalized eigenvalues of a random matrix pencil useful for the numerical inversion of the Laplace transform is discussed.
NAJan 28, 2008
Computational aspects and applications of a new transform for solving the complex exponentials approximation problemPiero Barone
Many real life problems can be reduced to the solution of a complex exponentials approximation problem which is usually ill posed. Recently a new transform for solving this problem, formulated as a specific moments problem in the plane, has been proposed in a theoretical framework. In this work some computational issues are addressed to make this new tool useful in practice. An algorithm is developed and used to solve a Nuclear Magnetic Resonance spectrometry problem, two time series interpolation and extrapolation problems and a shape from moments problem.
COMay 2, 2012
A black box method for solving the complex exponentials approximation problemPiero Barone
A common problem, arising in many different applied contexts, consists in estimating the number of exponentially damped sinusoids whose weighted sum best fits a finite set of noisy data and in estimating their parameters. Many different methods exist to this purpose. The best of them are based on approximate Maximum Likelihood estimators, assuming to know the number of damped sinusoids, which can then be estimated by an order selection procedure. As the problem can be severely ill posed, a stochastic perturbation method is proposed which provides better results than Maximum Likelihood based methods when the signal-to-noise ratio is low. The method depends on some hyperparameters which turn out to be essentially independent of the application. Therefore they can be fixed once and for all, giving rise to a black box method.