Operator-theoretic Analysis of Mutual Interactions in Synchronized Dynamics
Yuka Hashimoto, Masahiro Ikeda, Hiroya Nakao, Yoshinobu Kawahara
TL;DR
The paper tackles understanding mutual interactions in synchronized nonlinear oscillators by introducing KGME, a data-driven operator-theoretic framework based on the Koopman operator to estimate phase dynamics from time-series data. By recasting the phase-interaction estimation as a generalized multiparameter eigenvalue problem and leveraging multiple Koopman eigenvectors, KGME robustly recovers the phase function and a Fourier-based phase coupling function with guaranteed stability under data perturbations. The authors provide theoretical stability results showing the learned parameters are insensitive to the number of Fourier components M and harmonics j, and validate the approach on the van der Pol and FitzHugh-Nagumo models as well as real circadian data from cohabiting mice. The method advances robust, data-driven synchronization analysis and strengthens the connection between operator theory and nonlinear dynamics for practical applications in biology and engineering.
Abstract
Analyzing synchronized nonlinear oscillators is one of the most important and attractive topics in nonlinear science. By understanding the interactions between the oscillators, we can figure out the synchronization process. A promising approach to the analysis of interacting oscillators in nonlinear science is the application of the phase model. In this paper, we propose a data-driven approach to extract mutual interactions of synchronized oscillators based on the phase model. Recently, applying machine learning techniques to estimate models in physics has been actively investigated. We propose an operator-theoretic approach to estimate the phase model of interacting oscillators. We reduce the estimation problem to a multiparameter eigenvalue problem of the Koopman operator, a linear operator that describes a dynamical system. By reducing the problem to a linear algebraic problem, we can theoretically show that the proposed approach is stable with respect to perturbations in the given data.
