Physiological signaling is often oscillatory and shows nonlinearity due to complex interactions of underlying processes or signal propagation delays. This is particularly evident in case of brain activity which is subject to various feedback loop interactions between different brain structures, that coordinate their activity to support normal function. In order to understand such signaling in health and disease, methods are needed that can deal with such complex oscillatory phenomena. In this paper, a data-driven method for analyzing anharmonic oscillations is introduced. The KurSL model incorporates two well-studied components, which in the past have been used separately to analyze oscillatory behavior. The Sturm–Liouville equations describe a form of a general oscillation, and the Kuramoto coupling model represents a set of oscillators interacting in the phase domain. Integration of these components provides a flexible framework for capturing complex interactions of oscillatory processes of more general form than the most commonly used harmonic oscillators. The paper introduces a mathematical framework of the KurSL model and analyzes its behavior for a variety of parameter ranges. The significance of the model follows from its ability to provide information about coupled oscillators’ phase dynamics directly from the time series. KurSL offers a novel framework for analyzing a wide range of complex oscillatory behaviors, such as the ones encountered in physiological signals.