Abstract
A particle filter method is presented for the discrete-time filtering problem with nonlinear Itô stochastic ordinary differential equations (SODE) with additive noise supposed to be analytically integrable as a function of the underlying vector-Wiener process and time. The Diffusion Kernel Filter is arrived at by a parametrization of small noise-driven state fluctuations within branches of prediction and a local use of this parametrization in the Bootstrap Filter. The method applies for small noise and short prediction steps. With explicit numerical integrators, the operations count in the Diffusion Kernel Filter is shown to be smaller than in the Bootstrap Filter whenever the initial state for the prediction step has sufficiently few moments. The established parametrization is a dual-formula for the analysis of sensitivity to gaussian-initial perturbations and the analysis of sensitivity to noise-perturbations, in deterministic models, showing in particular how the stability of a deterministic dynamics is modeled by noise on short times and how the diffusion matrix of an SODE should be modeled (i.e. defined) for a gaussian-initial deterministic problem to be cast into an SODE problem. From it, a novel definition of prediction may be proposed that coincides with the deterministic path within the branch of prediction whose information entropy at the end of the prediction step is closest to the average information entropy over all branches. Tests are made with the Lorenz-63 equations, showing good results both for the filter and the definition of prediction.
Similar content being viewed by others
References
Apte, A., Hairer, M., Stuart, A.M., Voss, J.: Sampling the posterior: an approach to non-Gaussian data assimilation. Physica D 230, 50–64 (2007)
Chorin, A.J., Hald, O.H.: Stochastic Tools in Mathematics and Science. Springer, Berlin (2005)
Chorin, A.J., Hald, O.H., Kupferman, R.: Optimal prediction with memory. Physica D 166, 239–257 (2002)
Crisan, D., Del Moral, P., Lyons, T.J.: Discrete filtering using branching and interacting particle systems. Markov Proc. Relat. Fields 5, 293–331 (1999)
Doucet, A., Godsill, S., Andrieu, C.: On sequential Monte Carlo sampling methods for Bayesian filtering. Stat. Comput. 10, 197–208 (2000)
Gordon, N.J., Salmond, D.J., Smith, A.F.M.: Novel approach to nonlinear/non-Gaussian Bayesian state estimation. IEE Proc. F 140(2), 107–113 (1993)
Higham, D.J., Mao, X., Stuart, A.M.: Strong convergence of Euler-type methods for nonlinear stochastic differential equations. SIAM J. Numer. Anal. 40(3), 1041–1063 (2002)
Jazwinski, A.H.: Stochastic Processes and Filtering Theory. Academic Press, New York (1970)
Kato, T.: Perturbation Theory for Linear Operators. Springer, Berlin (1995)
Kim, S., Eyink, G.L., Restrepo, J.M., Alexander, F.J., Johnson, G.: Ensemble filtering for nonlinear dynamics. Mon. Wea. Rev. 131, 2586–2594 (2003)
Kloeden, P.E., Platen, E.: Numerical Solution of Stochastic Differential Equations. Springer, Berlin (1999)
Miller, R.N., Carter, E.F., Blue, S.T.: Data assimilation into nonlinear stochastic models. Tellus A 51(2), 167–194 (1999)
Öksendal, B.K.: Stochastic Differential Equations: An Introduction with Applications, 6th edn. Springer, Berlin (2003)
Raupp, C.F.M., Dias, P.L.S.: Dynamics of resonantly interacting equatorial waves. Tellus A 58(2), 263–276 (2006)
Smith, A.F.M., Gelfand, A.E.: Bayesian statistics without tears: a sampling-resampling perspective. Am. Stat. 46, 84–88 (1992)
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Krause, P. The Diffusion Kernel Filter. J Stat Phys 134, 365–380 (2009). https://doi.org/10.1007/s10955-008-9673-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10955-008-9673-1