Global dynamical structures from infinitesimal data
Scientists and engineers alike target modeling of complex, high dimensional, and nonlinear dynamical systems as a central goal. Machine learning breakthroughs alongside mounting computation and data advance the efficacy of learning from trajectory measurements. However scientifically interpreting data-driven models, e.g., localizing attracting sets and their basins, remains elusive. Such limitations particularly afflict identification of system-level regulatory mechanisms characteristic of living systems, e.g., stabilizing control for whole-body locomotion, where discontinuous, transient, and multiscale phenomena are common and prior models are rare. As a next step towards theory-grounded discovery of behavioral mechanisms in biology and beyond, we introduce VERT, a framework for discovering attracting sets from trajectories without recourse to any global model. Our infinitesimal-local-global (ILG) pipeline estimates the proximity of any sampled state to an attracting set, if one exists, with formal accuracy guarantees. We demonstrate our approach on phenomenological and physical oscillators with hierarchical and impulsive dynamics, finding sensitivity to both global and intermediate attractors composed in sequence and parallel. Application of VERT to human running kinematics data reveals insight into control modules that stabilize task-level dynamics, supporting a longstanding neuromechanical control hypothesis. The VERT framework promotes rigorous inference of underlying dynamical structure even for systems where learning a global dynamics model is impractical or impossible.