A multichannel generalization of the HAVOK method for the analysis of nonlinear dynamical systems
arXiv:2509.13424 · doi:10.1063/5.0303718
Abstract
By extending Takens' embedding theorem (1981), Deyle and Sugihara (2011) provided a theoretical justification for using parallel measurement time series to reconstruct a system's attractor. Building on Takens' framework, Brunton et al. (2017) introduced the Hankel alternative view of Koopman (HAVOK) algorithm, a data-driven approach capable of linearizing chaotic systems through delay embeddings. In this work, a modified version of the original algorithm is presented (mHAVOK), a practical realization of Deyle and Sugihara's generalized embedding theory. mHAVOK extends the original algorithm from one to multiple input time series and introduces a systematic approach to separating linear and nonlinear terms. An R2-informed quality score is introduced and shown to be a reliable guide for the selection of the reduced rank. The algorithm is tested on the familiar Lorenz system, as well as the more sophisticated Sprott system, which features different behaviors depending on the initial conditions. The quality of the reconstructions is assessed with the Chamfer distance, validating how mHAVOK allows for a more accurate reconstruction of the system dynamics. The new methodology generalizes HAVOK by allowing the analysis of multivariate time series, fundamental in real life data-driven applications.
References in corpus (7)
- A Data-Driven Approximation of the Koopman Operator: Extending Dynamic Mode Decomposition
- Linear predictors for nonlinear dynamical systems: Koopman operator meets model predictive control
- Applied Koopmanism
- Chaos as an Intermittently Forced Linear System
- Ergodic theory, Dynamic Mode Decomposition and Computation of Spectral Properties of the Koopman operator
- Extended dynamic mode decomposition with dictionary learning: a data-driven adaptive spectral decomposition of the Koopman operator
- Reduced-Order Models for Coupled Dynamical Systems: Data-driven Methods and the Koopman Operator