273 research outputs found
Linear theory for filtering nonlinear multiscale systems with model error
We study filtering of multiscale dynamical systems with model error arising
from unresolved smaller scale processes. The analysis assumes continuous-time
noisy observations of all components of the slow variables alone. For a linear
model with Gaussian noise, we prove existence of a unique choice of parameters
in a linear reduced model for the slow variables. The linear theory extends to
to a non-Gaussian, nonlinear test problem, where we assume we know the optimal
stochastic parameterization and the correct observation model. We show that
when the parameterization is inappropriate, parameters chosen for good filter
performance may give poor equilibrium statistical estimates and vice versa.
Given the correct parameterization, it is imperative to estimate the parameters
simultaneously and to account for the nonlinear feedback of the stochastic
parameters into the reduced filter estimates. In numerical experiments on the
two-layer Lorenz-96 model, we find that parameters estimated online, as part of
a filtering procedure, produce accurate filtering and equilibrium statistical
prediction. In contrast, a linear regression based offline method, which fits
the parameters to a given training data set independently from the filter,
yields filter estimates which are worse than the observations or even divergent
when the slow variables are not fully observed
Nonparametric Uncertainty Quantification for Stochastic Gradient Flows
This paper presents a nonparametric statistical modeling method for
quantifying uncertainty in stochastic gradient systems with isotropic
diffusion. The central idea is to apply the diffusion maps algorithm to a
training data set to produce a stochastic matrix whose generator is a discrete
approximation to the backward Kolmogorov operator of the underlying dynamics.
The eigenvectors of this stochastic matrix, which we will refer to as the
diffusion coordinates, are discrete approximations to the eigenfunctions of the
Kolmogorov operator and form an orthonormal basis for functions defined on the
data set. Using this basis, we consider the projection of three uncertainty
quantification (UQ) problems (prediction, filtering, and response) into the
diffusion coordinates. In these coordinates, the nonlinear prediction and
response problems reduce to solving systems of infinite-dimensional linear
ordinary differential equations. Similarly, the continuous-time nonlinear
filtering problem reduces to solving a system of infinite-dimensional linear
stochastic differential equations. Solving the UQ problems then reduces to
solving the corresponding truncated linear systems in finitely many diffusion
coordinates. By solving these systems we give a model-free algorithm for UQ on
gradient flow systems with isotropic diffusion. We numerically verify these
algorithms on a 1-dimensional linear gradient flow system where the analytic
solutions of the UQ problems are known. We also apply the algorithm to a
chaotically forced nonlinear gradient flow system which is known to be well
approximated as a stochastically forced gradient flow.Comment: Find the associated videos at: http://personal.psu.edu/thb11
- …