Learning particle swarming models from data with Gaussian processes

Abstract

Interacting particle or agent systems that display a rich variety of swarming behaviours are ubiquitous in science and engineering. A fundamental and challenging goal is to understand the link between individual interaction rules and swarming. In this paper, we study the data-driven discovery of a second-order particle swarming model that describes the evolution of $N$ particles in ${\mathbb{R}}^d$ under radial interactions. We propose a learning approach that models the latent radial interaction function as Gaussian processes, which can simultaneously fulfill two inference goals: one is the nonparametric inference of {the} interaction function with pointwise uncertainty quantification, and the other one is the inference of unknown scalar parameters in the non-collective friction forces of the system. We formulate the learning problem as a statistical inverse problem and provide a detailed analysis of recoverability conditions, establishing that a coercivity condition is sufficient for recoverability. Given data collected from $M$ i.i.d trajectories with independent Gaussian observational noise, we provide a finite-sample analysis, showing that our posterior mean estimator converges in a Reproducing kernel Hilbert space norm, at an optimal rate in $M$ equal to the one in the classical 1-dimensional Kernel Ridge regression. As a byproduct, we show we can obtain a parametric learning rate in $M$ for the posterior marginal variance using $L^{\infty }$ norm, and the rate could also involve $N$ and $L$ (the number of observation time instances for each trajectory), depending on the condition number of the inverse problem. Numerical results on systems that exhibit different swarming behaviors demonstrate efficient learning of our approach from scarce noisy trajectory data.

arXiv comment44 pages; Appendix 5 pages

References in corpus (12)

• Interaction Ruling Animal Collective Behaviour Depends on Topological rather than Metric Distance: Evidence from a Field Study
• Machine Learning of Linear Differential Equations using Gaussian Processes
• State Transitions and the Continuum Limit for a 2D Interacting, Self-Propelled Particle System
• Gaussian Processes and Kernel Methods: A Review on Connections and Equivalences
• Nonparametric inference of interaction laws in systems of agents from trajectory data
• Learning interaction kernels in heterogeneous systems of agents from multiple trajectories
• Frequentist coverage and sup-norm convergence rate in Gaussian process regression
• Learning unknown ODE models with Gaussian processes
• A Laplacian Approach to Stubborn Agents and their Role in Opinion Formation on Influence Networks
• Solving and Learning Nonlinear PDEs with Gaussian Processes
• On the Estimation of Derivatives Using Plug-in Kernel Ridge Regression Estimators
• Learning interaction kernels in stochastic systems of interacting particles from multiple trajectories

Cited by in corpus (2)

• Learning Mean-Field Equations from Particle Data Using WSINDy
• Machine learning of continuous and discrete variational ODEs with convergence guarantee and uncertainty quantification