Series Expansion Approximations of Brownian Motion for Non-Linear Kalman Filtering of Diffusion Processes

• Published in: IEEE transactions on signal processing 2014-03, Vol.62 (6), p.1514-1524
• Main Authors: Lyons, Simon M. J., Sarkka, Simo, Storkey, Amos J.
• Format: Article
• Language: English
• Subjects: Applied sciences; Approximation methods; Brownian motion; Decomposition; Density; Detection, estimation, filtering, equalization, prediction; Differential equations; Exact sciences and technology; Filtering; Filtration; Information, signal and communications theory; Kalman filters; Markov processes; Mathematical model; multidimensional signal processing; Noise; nonlinear filters; Nonlinearity; Signal and communications theory; Signal, noise; Statistics; Stochastic processes; Telecommunications and information theory; Transforms
• ISSN: 1053-587X; 1941-0476
• DOI: 10.1109/TSP.2014.2303430

In this paper, we describe a novel application of sigma-point methods to continuous-discrete filtering. The nonlinear continuous-discrete filtering problem is often computationally intractable to solve. Assumed density filtering methods attempt to match statistics of the filtering distribution to some set of more tractable probability distributions. Filters such as these are usually decompose the problem into two sub-problems. The first of these is a prediction step, in which one uses the known dynamics of the signal to predict its state at time tk+1 given observations up to time tk. In the second step, one updates the prediction upon arrival of the observation at time tk+1. The aim of this paper is to describe a novel method that improves the prediction step. We decompose the Brownian motion driving the signal in a generalised Fourier series, which is truncated after a number of terms. This approximation to Brownian motion can be described using a relatively small number of Fourier coefficients, and allows us to compute statistics of the filtering distribution with a single application of a sigma-point method. Assumed density filters that exist in the literature usually rely on discretisation of the signal dynamics followed by iterated application of a sigma point transform (or a limiting case thereof). Iterating the transform in this manner can lead to loss of information about the filtering distribution in highly non-linear settings. We demonstrate that our method is better equipped to cope with such problems.