Nonlinear tracking in a diffusion process with a Bayesian filter and the finite element method

• Published in: Computational statistics & data analysis Vol. 55; no. 1; pp. 280 - 290
• Main Authors: Pedersen, M.W., Thygesen, U.H., Madsen, H.
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 2011; Elsevier
• Series: Computational Statistics & Data Analysis
• Subjects: Computer simulation; Density; Finite element method; Finite element method Hidden Markov model Nonlinear state estimation Point-mass filter Sequential Monte Carlo Stochastic differential equation; Hidden Markov model; Mathematical analysis; Mathematical models; Monte Carlo methods; Nonlinear state estimation; Nonlinearity; Point-mass filter; Sequential Monte Carlo; Sheet molding compounds; Stochastic differential equation; Finite element method; Hidden Markov model; Sequential Monte Carlo; Nonlinear state estimation; Point-mass filter; Stochastic differential equation
• ISSN: 0167-9473; 1872-7352
• DOI: 10.1016/j.csda.2010.04.018

A new approach to nonlinear state estimation and object tracking from indirect observations of a continuous time process is examined. Stochastic differential equations (SDEs) are employed to model the dynamics of the unobservable state. Tracking problems in the plane subject to boundaries on the state-space do not in general provide analytical solutions. A widely used numerical approach is the sequential Monte Carlo (SMC) method which relies on stochastic simulations to approximate state densities. For off-line analysis, however, accurate smoothed state density and parameter estimation can become complicated using SMC because Monte Carlo randomness is introduced. The finite element (FE) method solves the Kolmogorov equations of the SDE numerically on a triangular unstructured mesh for which boundary conditions to the state-space are simple to incorporate. The FE approach to nonlinear state estimation is suited for off-line data analysis because the computed smoothed state densities, maximum a posteriori parameter estimates and state sequence are deterministic conditional on the finite element mesh and the observations. The proposed method is conceptually similar to existing point-mass filtering methods, but is computationally more advanced and generally applicable. The performance of the FE estimators in relation to SMC and to the resolution of the spatial discretization is examined empirically through simulation. A real-data case study involving fish tracking is also analysed.