7,800 research outputs found
A new variational radial basis function approximation for inference in multivariate diffusions
In this paper we present a radial basis function based extension to a recently proposed variational algorithm for approximate inference for diffusion processes. Inference, for state and in particular (hyper-) parameters, in diffusion processes is a challenging and crucial task. We show that the new radial basis function approximation based algorithm converges to the original algorithm and has beneficial characteristics when estimating (hyper-)parameters. We validate our new approach on a nonlinear double well potential dynamical system
Variational approach for learning Markov processes from time series data
Inference, prediction and control of complex dynamical systems from time
series is important in many areas, including financial markets, power grid
management, climate and weather modeling, or molecular dynamics. The analysis
of such highly nonlinear dynamical systems is facilitated by the fact that we
can often find a (generally nonlinear) transformation of the system coordinates
to features in which the dynamics can be excellently approximated by a linear
Markovian model. Moreover, the large number of system variables often change
collectively on large time- and length-scales, facilitating a low-dimensional
analysis in feature space. In this paper, we introduce a variational approach
for Markov processes (VAMP) that allows us to find optimal feature mappings and
optimal Markovian models of the dynamics from given time series data. The key
insight is that the best linear model can be obtained from the top singular
components of the Koopman operator. This leads to the definition of a family of
score functions called VAMP-r which can be calculated from data, and can be
employed to optimize a Markovian model. In addition, based on the relationship
between the variational scores and approximation errors of Koopman operators,
we propose a new VAMP-E score, which can be applied to cross-validation for
hyper-parameter optimization and model selection in VAMP. VAMP is valid for
both reversible and nonreversible processes and for stationary and
non-stationary processes or realizations
Data-driven model reduction and transfer operator approximation
In this review paper, we will present different data-driven dimension
reduction techniques for dynamical systems that are based on transfer operator
theory as well as methods to approximate transfer operators and their
eigenvalues, eigenfunctions, and eigenmodes. The goal is to point out
similarities and differences between methods developed independently by the
dynamical systems, fluid dynamics, and molecular dynamics communities such as
time-lagged independent component analysis (TICA), dynamic mode decomposition
(DMD), and their respective generalizations. As a result, extensions and best
practices developed for one particular method can be carried over to other
related methods
A variational approach to modeling slow processes in stochastic dynamical systems
The slow processes of metastable stochastic dynamical systems are difficult
to access by direct numerical simulation due the sampling problem. Here, we
suggest an approach for modeling the slow parts of Markov processes by
approximating the dominant eigenfunctions and eigenvalues of the propagator. To
this end, a variational principle is derived that is based on the maximization
of a Rayleigh coefficient. It is shown that this Rayleigh coefficient can be
estimated from statistical observables that can be obtained from short
distributed simulations starting from different parts of state space. The
approach forms a basis for the development of adaptive and efficient
computational algorithms for simulating and analyzing metastable Markov
processes while avoiding the sampling problem. Since any stochastic process
with finite memory can be transformed into a Markov process, the approach is
applicable to a wide range of processes relevant for modeling complex
real-world phenomena
Representation Learning on Graphs: A Reinforcement Learning Application
In this work, we study value function approximation in reinforcement learning
(RL) problems with high dimensional state or action spaces via a generalized
version of representation policy iteration (RPI). We consider the limitations
of proto-value functions (PVFs) at accurately approximating the value function
in low dimensions and we highlight the importance of features learning for an
improved low-dimensional value function approximation. Then, we adopt different
representation learning algorithm on graphs to learn the basis functions that
best represent the value function. We empirically show that node2vec, an
algorithm for scalable feature learning in networks, and the Variational Graph
Auto-Encoder constantly outperform the commonly used smooth proto-value
functions in low-dimensional feature space
Stochastic collective dynamics of charged--particle beams in the stability regime
We introduce a description of the collective transverse dynamics of charged
(proton) beams in the stability regime by suitable classical stochastic
fluctuations. In this scheme, the collective beam dynamics is described by
time--reversal invariant diffusion processes deduced by stochastic variational
principles (Nelson processes). By general arguments, we show that the diffusion
coefficient, expressed in units of length, is given by ,
where is the number of particles in the beam and the Compton
wavelength of a single constituent. This diffusion coefficient represents an
effective unit of beam emittance. The hydrodynamic equations of the stochastic
dynamics can be easily recast in the form of a Schr\"odinger equation, with the
unit of emittance replacing the Planck action constant. This fact provides a
natural connection to the so--called ``quantum--like approaches'' to beam
dynamics. The transition probabilities associated to Nelson processes can be
exploited to model evolutions suitable to control the transverse beam dynamics.
In particular we show how to control, in the quadrupole approximation to the
beam--field interaction, both the focusing and the transverse oscillations of
the beam, either together or independently.Comment: 15 pages, 9 figure
Hydrodynamic Flows on Curved Surfaces: Spectral Numerical Methods for Radial Manifold Shapes
We formulate hydrodynamic equations and spectrally accurate numerical methods
for investigating the role of geometry in flows within two-dimensional fluid
interfaces. To achieve numerical approximations having high precision and level
of symmetry for radial manifold shapes, we develop spectral Galerkin methods
based on hyperinterpolation with Lebedev quadratures for -projection to
spherical harmonics. We demonstrate our methods by investigating hydrodynamic
responses as the surface geometry is varied. Relative to the case of a sphere,
we find significant changes can occur in the observed hydrodynamic flow
responses as exhibited by quantitative and topological transitions in the
structure of the flow. We present numerical results based on the
Rayleigh-Dissipation principle to gain further insights into these flow
responses. We investigate the roles played by the geometry especially
concerning the positive and negative Gaussian curvature of the interface. We
provide general approaches for taking geometric effects into account for
investigations of hydrodynamic phenomena within curved fluid interfaces.Comment: 14 figure
- …