Table of Contents
Fetching ...

Closed-form Filtering for Non-linear Systems

Théophile Cantelobre, Carlo Ciliberto, Benjamin Guedj, Alessandro Rudi

TL;DR

This work proposes a new class of filters based on Gaussian PSD Models, which offer several advantages in terms of density approximation and computational efficiency, and shows that filtering can be efficiently performed in closed form when transitions and observations are Gaussian PSD Models.

Abstract

Sequential Bayesian Filtering aims to estimate the current state distribution of a Hidden Markov Model, given the past observations. The problem is well-known to be intractable for most application domains, except in notable cases such as the tabular setting or for linear dynamical systems with gaussian noise. In this work, we propose a new class of filters based on Gaussian PSD Models, which offer several advantages in terms of density approximation and computational efficiency. We show that filtering can be efficiently performed in closed form when transitions and observations are Gaussian PSD Models. When the transition and observations are approximated by Gaussian PSD Models, we show that our proposed estimator enjoys strong theoretical guarantees, with estimation error that depends on the quality of the approximation and is adaptive to the regularity of the transition probabilities. In particular, we identify regimes in which our proposed filter attains a TV $ε$-error with memory and computational complexity of $O(ε^{-1})$ and $O(ε^{-3/2})$ respectively, including the offline learning step, in contrast to the $O(ε^{-2})$ complexity of sampling methods such as particle filtering.

Closed-form Filtering for Non-linear Systems

TL;DR

This work proposes a new class of filters based on Gaussian PSD Models, which offer several advantages in terms of density approximation and computational efficiency, and shows that filtering can be efficiently performed in closed form when transitions and observations are Gaussian PSD Models.

Abstract

Sequential Bayesian Filtering aims to estimate the current state distribution of a Hidden Markov Model, given the past observations. The problem is well-known to be intractable for most application domains, except in notable cases such as the tabular setting or for linear dynamical systems with gaussian noise. In this work, we propose a new class of filters based on Gaussian PSD Models, which offer several advantages in terms of density approximation and computational efficiency. We show that filtering can be efficiently performed in closed form when transitions and observations are Gaussian PSD Models. When the transition and observations are approximated by Gaussian PSD Models, we show that our proposed estimator enjoys strong theoretical guarantees, with estimation error that depends on the quality of the approximation and is adaptive to the regularity of the transition probabilities. In particular, we identify regimes in which our proposed filter attains a TV -error with memory and computational complexity of and respectively, including the offline learning step, in contrast to the complexity of sampling methods such as particle filtering.
Paper Structure (47 sections, 26 theorems, 92 equations, 2 algorithms)

This paper contains 47 sections, 26 theorems, 92 equations, 2 algorithms.

Key Result

Proposition 1

Let $f(x, y; \theta_1)$ and $g(y, z; \theta_2)$ be two Gaussian PSD Models of order $M_1$ and $M_2$ respectively, as in def:psd-model. Then there exist some algorithms $\mathop{\mathrm{\textsc{Integral}}}\nolimits$, $\mathop{\mathrm{\textsc{PartialEval}}}\nolimits$, $\mathop{\mathrm{\textsc{Product}

Theorems & Definitions (61)

  • Definition 1
  • Example 1: Gaussian Mixture Model
  • Example 2: Squared linear Gaussian model
  • Proposition 1: Closed form operations for Gaussian PSD Models
  • Proposition 2: Constant order for Markov transition
  • Proposition 3: Generality of \ref{['assumption:target_function']}, Prop. 5 in ciliberto2021
  • Theorem 4
  • Corollary 5: $\hat{\pi}_n$ has constant order for any $n$
  • Theorem 6: PSD filter robustness and stability
  • Lemma 7
  • ...and 51 more