Рассматривается модель полносвязной ассоциации нейронов с синаптической электрической связью, представляющая собой систему дифференциальных уравнений с запаздыванием. Специальной заменой эта система приводится к системе импульсных обыкновенных дифференциальных уравнений. Для соответствующей динамической системы в случае изучаются вопросы существования, устойчивости и асимптотического представления периодических решений на основании бифуркационного анализа двумерного отображения — оператора сдвига по траекториям решения специальной системы двух дифференциальных уравнений. Особое внимание уделяется числу сосуществующих устойчивых режимов. Исследуется задача нахождения параметров, для которых количество таких режимов максимально. Для поиска неподвижных точек полученного двумерного отображения используется численное исследование, основанное на следующей итерационной процедуре. Выбирается начальная точка, методом Рунге-Кутты с заданным шагом вычисляются значения решения на отрезке , в конечной точке этого отрезка значение решения сравнивается с начальным и, если отклонение превышает заданное значение, то значение в конечной точке принимается за начальное и цикл вычислений методом Рунге-Кутты повторяется. Вычисления заканчиваются, если достигнуто требуемое малое отклонение, т. е. найдена неподвижная точка оператора сдвига, и соответствующий устойчивый периодический режим, или если количество итераций достигает заданного большого числа, что свидетельствует об отсутствии неподвижной точки. В работе представлены результаты проведенного численного исследования, позволившего продемонстрировать основные перестройки, происходящие в фазовом пространстве двумерного отображения. Полученные неподвижные точки позволяют найти асимптотические устойчивые решения исходной задачи.
A model of a fully connected association of neurons with a synaptic electrical connection which is a system of differential equations with delay is considered. By a special substitution, this system is reduced to a system of impulsive ordinary differential equations. For the corresponding dynamical system in the case , we study the existence, stability, and asymptotic representation of periodic solutions on the basis of a bifurcation analysis of a two-dimensional mapping, a shift operator along trajectories of a solution to a special system of two differential equations. Particular attention is paid to the number of coexisting stable regimes. We study the problem of finding parameters for which the number of such modes is maximum. In order to search the fixed points of the resulting two-dimensional mapping, a numerical study is used based on the following iterative procedure. Selected the starting point, the Runge-Kutta method with a given step calculates the solution values on the segment . At the end point of this segment, the solution value is compared with the initial one and if the deviation exceeds the specified value, then the value at the end point is taken as the initial one and the calculation cycle by the Runge-Kutta method is repeated. The calculations terminate when the required small deviation is reached, i.e., a fixed point of the shift operator is found, and so is the corresponding stable periodic mode, or when the number of iterations reaches a given large number, and this indicates the absence of a fixed point. The paper presents the results of a numerical study that made it possible to demonstrate the main rearrangements occurring in the phase space of a two-dimensional mapping. The obtained fixed points allow us to find asymptotic stable solutions of the original problem.
The authors declare that there are no conflicts of interest present.