FAST PREDICTION-ERROR DETECTOR FOR ESTIMATING SPARSE-SPIKE SEQUENCES.

Based on the maximum-likelihood principle, a locally optimal method for detecting the location and estimating the amplitude of spikes in a sequence is considered, based on a 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 using a prediction-error formulation the 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 this algorithm is self-initialized, consistent, and fast.