A Fourier State Space Model for Bayesian ODE Filters

07/17/2020
by   Hans Kersting, et al.
0

Gaussian ODE filtering is a probabilistic numerical method to solve ordinary differential equations (ODEs). It computes a Bayesian posterior over the solution from evaluations of the vector field defining the ODE. Its most popular version, which employs an integrated Brownian motion prior, uses Taylor expansions of the mean to extrapolate forward and has the same convergence rates as classical numerical methods. As the solution of many important ODEs are periodic functions (oscillators), we raise the question whether Fourier expansions can also be brought to bear within the framework of Gaussian ODE filtering. To this end, we construct a Fourier state space model for ODEs and a `hybrid' model that combines a Taylor (Brownian motion) and Fourier state space model. We show by experiments how the hybrid model might become useful in cheaply predicting until the end of the time domain.

READ FULL TEXT

page 1

page 2

page 3

page 4

10/08/2018

Probabilistic Solutions To Ordinary Differential Equations As Non-Linear Bayesian Filtering: A New Perspective

We formulate probabilistic numerical approximations to solutions of ordi...
07/25/2018

Convergence Rates of Gaussian ODE Filters

A recently-introduced class of probabilistic (uncertainty-aware) solvers...
11/24/2021

State-space deep Gaussian processes with applications

This thesis is mainly concerned with state-space approaches for solving ...
03/18/2021

A Probabilistic State Space Model for Joint Inference from Differential Equations and Data

Mechanistic models with differential equations are a key component of sc...
09/25/2017

Bayesian Filtering for ODEs with Bounded Derivatives

Recently there has been increasing interest in probabilistic solvers for...
05/25/2018

A Multi-Scan Labeled Random Finite Set Model for Multi-object State Estimation

State space models in which the system state is a finite set--called the...
09/02/2018

Stable approximation schemes for optimal filters

We explore a general truncation scheme for the approximation of (possibl...