Probability Measures for Numerical Solutions of Differential Equations

作者: Simo Särkkä , Konstantinos Zygalakis , Mark Girolami , Andrew Stuart , Patrick R. Conrad

DOI:

关键词: Numerical stabilityCollocation methodMathematical analysisDifferential algebraic equationNumerical partial differential equationsExplicit and implicit methodsMathematicsExponential integratorDelay differential equationStochastic partial differential equation

摘要: In this paper, we present a formal quantification of epistemic uncertainty induced by numerical solutions ordinary and partial differential equation models. Numerical equations contain inherent uncertainties due to the finite dimensional approximation an unknown implicitly defined function. When statistically analysing models based on describing physical, or other naturally occurring, phenomena, it is therefore important explicitly account for introduced method. This enables objective determination its importance relative uncertainties, such as those caused data contaminated with noise model error missing physical inadequate descriptors. To end show that wide variety existing solvers can be randomised, inducing probability measure over equations. These measures exhibit contraction Dirac around true solution, where rates convergence are consistent underlying deterministic Ordinary elliptic used illustrate approach quantifying in both statistical analysis forward inverse problems.

参考文章(15)
J. Andrés Christen, Marcos Capistrán, Sophie Donnet, Bayesian Analysis of ODE's: solver optimal accuracy and Bayes factors arXiv: Computation. ,(2013) , 10.1137/140976777
Oksana A. Chkrebtii, Ben Calderhead, Mark A. Girolami, David A. Campbell, Bayesian Solution Uncertainty Quantification for Differential Equations arXiv: Methodology. ,(2013)
Finn Lindgren, Håvard Rue, Johan Lindström, An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach Journal of The Royal Statistical Society Series B-statistical Methodology. ,vol. 73, pp. 423- 498 ,(2011) , 10.1111/J.1467-9868.2011.00777.X
E. Hairer, R. I. McLachlan, A. Razakarivony, ACHIEVING BROUWER'S LAW WITH IMPLICIT RUNGE-KUTTA METHODS Bit Numerical Mathematics. ,vol. 48, pp. 231- 243 ,(2008) , 10.1007/S10543-008-0170-3
Marc C. Kennedy, Anthony O'Hagan, Bayesian Calibration of computer models Journal of The Royal Statistical Society Series B-statistical Methodology. ,vol. 63, pp. 425- 464 ,(2001) , 10.1111/1467-9868.00294
Andrea Arnold, Daniela Calvetti, Erkki Somersalo, Linear multistep methods, particle filtering and sequential Monte Carlo Inverse Problems. ,vol. 29, pp. 085007- ,(2013) , 10.1088/0266-5611/29/8/085007
Jari Kaipio, Erkki Somersalo, Statistical inverse problems: discretization, model reduction and inverse crimes Journal of Computational and Applied Mathematics. ,vol. 198, pp. 493- 504 ,(2007) , 10.1016/J.CAM.2005.09.027
Gilbert Stengle, Error analysis of a randomized numerical method Numerische Mathematik. ,vol. 70, pp. 119- 128 ,(1995) , 10.1007/S002110050113
Ibrahim Coulibaly, Christian Lécot, A quasi-randomized Runge-Kutta method Mathematics of Computation. ,vol. 68, pp. 651- 659 ,(1999) , 10.1090/S0025-5718-99-01056-X