Richard A. Norton

Richard A. Norton, G. R. W. Quispel

In this paper we consider discrete gradient methods for approximating the solution and preserving a first integral (also called a constant of motion) of autonomous ordinary differential equations. We prove under mild conditions for a large class of discrete gradient methods that the numerical solution exists and is locally unique, and that for arbitrary $p \in \mathbb{N}$ we may construct a method that is of order $p$. In the proofs of these results we also show that the constants in the time step constraint and the error bounds may be chosen independently from the distance to critical points of the first integral. In the case when the first integral is quadratic, for arbitrary $p \in \mathbb{N}$, we have devised a new method that is linearly implicit at each time step and of order $p$. This new method has significant advantages in terms of efficiency. We illustrate our theory with a numerical example.

Richard A. Norton, Colin Fox, Malcolm E. Morrison

We develop a finite-dimensional approximation of the Frobenius-Perron operator using the finite volume method applied to the continuity equation for the evolution of probability. A Courant-Friedrichs-Lewy condition ensures that the approximation satisfies the Markov property, while existing convergence theory for the finite volume method guarantees convergence of the discrete operator to the continuous operator as mesh size tends to zero. Properties of the approximation are demonstrated in a computed example of sequential inference for the state of a low-dimensional mechanical system when observations give rise to multi-modal distributions.