Mathematical Model of Van der Pol-Airy Fractional Oscillator

Abstract

В работе предложена математическая модель нелинейного осциллятора Ван дер Поля-Эйри с учетом наследственности. Нелинейность осциллятора обусловлена наличием зависимости коэффициента трения от квадрата функции смещения, что характерно для осциллятора Ван дер Поля. Также собственная частота колебаний представляет собой функцию от времени, которая линейно возрастает при его возрастании. Последнее характерно для осциллятора Эйри. Эффекты наследственности вводятся в модельное уравнение посредством дробных производных в смысле Герасимова-Капуто. Они указывают на то, что колебательная система может обладать эффектами памяти, которые проявляются в зависимости текущего ее состояния от предыдущих. Для предложенной математической модели был разработан численный алгоритм, основанный на явной конечно-разностной схемы первого порядка. Численный алгоритм был реализован в компьютерной программе на языке Maple, с помощью которой была произведена визуализация результатов моделирования. Были построены осциллограммы и фазовые траектории при различных значениях параметров модели. Показано, что дробная математическая модель может обладать различными колебательными режимами: от автоколебательных, затухающих и хаотических. Дается интерпретация результатов моделирования. The paper proposes a mathematical model of the nonlinear Van der Pol-Airy oscillator taking into account heredity. The nonlinearity of the oscillator is due to the dependence of the friction coefficient on the square of the displacement function, which is typical for the Van der Pol oscillator. Also, the natural frequency of oscillations is a function of time, which increases linearly as it increases. The latter is typical for the Airy oscillator. Heredity effects are introduced into the model equation through fractional derivatives in the Gerasimov-Caputo sense. They indicate that the oscillatory system may have memory effects that manifest themselves depending on its current state from previous ones. For the proposed mathematical model, a numerical algorithm was developed based on an explicit first-order finite-difference scheme. The numerical algorithm was implemented in a computer program in the Maple language, with the help of which the simulation results were visualized. Oscillograms and phase trajectories were constructed for various values of the model parameters. It is shown that a fractional mathematical model can have various oscillatory modes: from self-oscillatory, damped and chaotic. An interpretation of the simulation results is given