A fast prediction-error detector for estimating sparse-spike sequences

Based on the Maximum-Likelihood principle, we develop a locally optimal method for detecting the location and estimating the amplitude of spikes in a sequence, which are considered the random input of a known ARMA model. A Bernoulli-Gaussian product model is adopted for the sparse-spike sequence, and the available data consist of a single, noisy, output record. By employing a Prediction-Error formulation our iterative algorithm guarantees the increase of a unique likelihood function used for the combined estimation/detection problem. Amplitude estimation is carried out with Kalman smoothing techniques, and event detection is performed in two ways, as an event adder and as an event remover. Synthetic examples verify that our algorithm is self-initialized, consistent, and fast.