Next Article in Journal
Estimation of Mediation Effect on Zero-Inflated Microbiome Mediators
Previous Article in Journal
Sparse Support Tensor Machine with Scaled Kernel Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Travelling Waves in the Ring of Coupled Oscillators with Delayed Feedback

by
Alexandra Kashchenko
*,†,
Ilia Kashchenko
and
Sergey Kondratiev
Regional Scientific and Educational Mathematical Center «Centre of Integrable Systems», P. G. Demidov Yaroslavl State University, 150003 Yaroslavl, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(13), 2827; https://doi.org/10.3390/math11132827
Submission received: 17 May 2023 / Revised: 19 June 2023 / Accepted: 22 June 2023 / Published: 24 June 2023
(This article belongs to the Section Difference and Differential Equations)

Abstract

:
We studied travelling waves in N nonlinear differential equations with a delay and large parameter. This system is important because it can be regarded as a phenomenological model of N-coupled neuron-like oscillators with delay. The problem of the existence of travelling-wave-type solutions was reduced to the study of the dynamics of an auxiliary equation with two delays. Using a special asymptotic method for the large parameter we proved that this equation has a relaxation cycle, studied its properties (amplitude, period and asymptotics) and found the sufficient stability conditions. Based on this periodic solution the travelling waves of the initial model were constructed.

1. Introduction

Differential equations with delay arise in many applied problems in biology, physics, medicine and ecology [1,2,3,4,5,6,7].
A single equation of the form
x ˙ + ν x = λ F ( x ( t T 1 ) ) ,
where ν , λ , and T 1 are positive parameters and F is a nonlinear function plays an important role in radiophysical and biological applications. For example, in radiophysics this model simulates a generator with a first-order low-pass filter and delayed feedback [8,9]. Such generators are used in the manufacture of D-amplifiers, sonars and noise radars [8]. Moreover, model (1) simulates a biological process where the single state variable x decays at a present rate ν proportional to x, and is produced at a rate dependent on its the value at some time in the past [10]. Such processes arise in many areas in biology, such as in population biology, neurophysiology and metabolic regulation [10]. If the function F is compactly supported and λ is large enough equation (1) has a relaxation cycle [11]. Note, that despite the fact that only the first derivative is present in the equation (there is no second derivative), model (1) has oscillatory regimes for many nonlinear functions F. (See, for example, [2,8,9,10,11]).
In this paper we considered a unidirectionally coupled ring
x ˙ k + x k = λ F ( x k ( t T 1 ) ) + γ ( x k 1 x k ) , k = 1 , , N , x 0 x N ,
of N oscillators of the form (1) with parameter ν = 1 . Here, parameter λ is positive and sufficiently large ( λ 1 ) , delay time T 1 and coupling parameter γ are positive, and feedback function F is a compactly supported positive function:
F ( y ) = f ( y ) , if   y [ p , p ] , 0 , if   y < p   or   y > p ,
where p is some positive constant, and f ( y ) is a piece-wise continuous and bounded function on the segment y [ p , p ] and f ( y ) > 0 for all y ( 0 , p ] . In the segment y [ p , 0 ] it may change its sign, but it can’t be zero on any interval of the nonzero length.
The study of the dynamics of coupled oscillators is of great interest because they simulate a lot of processes in different areas of science; for example, they simulate the work of the heart or neurons in physiology and the rotation of generators in electrical grids in physics [12,13,14].
In the papers of [15,16] the nonlocal dynamics of model (2) when N = 2 , λ 1 and an asymptotically small γ were studied. It was shown that studying the existence and stability of the periodic solutions of this infinite-dimensional problem may be reduced to studying the dynamics of constructed finite-dimensional mappings, which were constructed for different orders of smallness of parameter γ on λ . The asymptotics of the periodic solutions of the initial model were found, and multistability was proven.
The nonlocal dynamics of the ring of diffusion-coupled oscillators
x ˙ k + x k = λ F ( x k ( t T 1 ) ) + γ ( x k 1 2 x k + x k + 1 ) , k = 1 , , N , x 0 x N , x N + 1 x 1 ,
under condition λ 1 was studied [17,18]. It was proved that when γ > 0 all oscillators were synchronized, and when 1 4 < γ < 0 and N was even two-cluster synchronization was observed. When 1 4 < γ < 0 and N = 3 , different periodic inhomogeneous regimes were found, and there were non-regular oscillations.
In this paper we studied the existence and properties of the travelling waves
x k ( t ) = x ( t ( k 1 ) T 2 ) , k = 1 , , N ,
where T 2 is some phase lag of model (2).
The existence of the travelling waves in the rings of coupled oscillators was studied in the papers of [19,20]. Note that for the case of unidirectionally coupled oscillators, solution (4) is often called “a rotating wave” (see [21,22]).
The function x ( t ) from (4) must satisfy the equation with two delays
x ˙ + x = λ F ( x ( t T 1 ) ) + γ ( x ( t T 2 ) x ) .
We showed that Equation (5) has a periodic solution with period P. Moreover, we proved that it is exponentially orbitally stable in the phase space C [ T , 0 ] , where T = max ( T 1 , T 2 ) under some additional requirements on the function f.
Condition x 0 x N lead us to
T 2 N = n P
for some integer n. Therefore, if equality (6) is true, then our initial model (2) has a travelling wave solution (4).
Furthermore, if we consider the ring of N oscillators with delayed coupling
x ˙ k + x k = λ F ( x k ( t T 1 ) ) + γ ( x k 1 ( t T 2 ) x k ) , k = 1 , , N , x 0 x N
instead of model (2), then its homogeneous regime is the solution to Equation (5).
In the paper, we proved the existence of several coexisting travelling waves of model (2) in the case of large numbers N (N of the order O ( ln λ ) ). These solutions are relaxation spike-like regimes. In addition, we proved that all solutions of system (2) with positive initial conditions were positive for all t [ 0 , + ) . That is why we can say that this system may be regarded as a phenomenological model of N-coupled neurons [23].
It is important to mention, that models (2), (5) and (7) are rather complicated. At present, there are no analytical methods to study their dynamics on the semiaxis t [ 0 , + ) for any arbitrary values of parameters λ , γ , T 1 , and T 2 . If we study the dynamics of these models numerically, we take several concrete functions F, but we cannot enumerate all infinite sets of compactly supported functions F and draw justified conclusions about the qualitative behavior of a model with an arbitrarily compact supported nonlinearity. That is why we analyzed the nonlocal dynamics of this model under the assumption that λ was a large parameter.
This assumption allowed us to use a special analytical method to study the dynamic properties of the solutions to models (2), (5) and (7) on the entire semiaxis t [ 0 , + ) . Let’s describe the essence of this method in the simplest case of a model with one delay. First, we selected a special set of initial conditions S from phase space C [ T , 0 ] , where T is the delay time of the model. Then we integrated the model using the method of steps [24]. On the first step t [ 0 , T ] function u ( t T ) is a known function, the initial condition φ ( t ) S . That is why, on this segment, we considered our initial equation with delay as an ordinary differential equation with inhomogeneity depending on φ ( t ) . We found asymptotics of the solution to this equation at λ + . After that, we considered our model on the segment t [ T , 2 T ] as an ordinary differential equation with inhomogeneity depending on u ( t T ) , which is known from the previous step, and constructed the asymptotics to this model at λ + . We did the same several times and then proved that for every initial function φ S there existed a moment t 0 such that the solution to our model with this initial condition returned to the set S. This meant that there existed a Poincaré operator Π of translation along the trajectories such that Π S S . That is why using the Schauder fixed-point theorem [25] we concluded that there was a fixed point φ * S of the operator Π : Π φ * = φ * . Therefore, if we took the function φ * as the initial condition to the equation with delay, then we derived a periodic solution to this equation.
Note, that this analytical method is rather general and may be applied to various systems of differential equations with delay, including mathematical models of radiophysical devices.
The paper is organized is as follows. In Section 2 we prove the positiveness of solutions to system (2) and Equation (5) with positive initial conditions. In Section 3 we state and prove some properties of the linear part of Equation (5). In Section 4 we construct asymptotics of the relaxation solutions to Equation (5) and prove that there exists a relaxation cycle of this equation. In Section 5 we give the sufficient conditions of the stability of this cycle. In Section 6 we discuss the conditions for the existence of the travelling wave solutions to model (2). In Section 7 we give some generalizations about the results and draw conclusions.

2. Positiveness of Solutions

In this section we prove the positiveness of all solutions to models (2) and (5) with positive initial conditions.
First, let’s prove the positiveness of all solutions to an equation with two delays (5) with positive initial conditions.
Lemma 1.
Let φ ( t ) C [ T , 0 ] and φ ( t ) > 0 for all t [ T , 0 ] , where T = max ( T 1 , T 2 ) . Then solution x ( t ) to the Equation (5) with initial conditions x ( t ) = φ ( t ) for t [ T , 0 ] satisfies inequality x ( t ) > 0 for all t > 0 .
Proof. 
Initial condition φ ( t ) is positive on the segment t [ T , 0 ] ; therefore, x ( 0 ) > 0 and x ( t T 2 ) > 0 for all t [ 0 , T 2 ] . Function F is non-negative for all values of its argument.
Suppose there exists a moment t = t 0 > 0 such that x ( t 0 ) = 0 and x ( t ) > 0 for all t [ 0 , t 0 ) . Then from Formula (5) we obtain
x ˙ ( t 0 ) = λ F ( x ( t 0 T 1 ) ) + γ x ( t 0 T 2 ) > 0 .
However, if for all t < t 0 the inequality x ( t ) > 0 = x ( t 0 ) holds, then x ˙ ( t 0 ) should be negative. This contradiction completes the proof. □
The next Lemma states that all solutions of system (2) with positive initial conditions are positive. This result is important because the research model may be regarded as a phenomenological model of coupled neurons only if its solutions with positive initial conditions stay positive on the whole semiaxis t [ 0 , + ) .
Lemma 2.
Let φ k ( t ) C [ T 1 , 0 ] and φ k ( t ) > 0 for all t [ T 1 , 0 ] and k = 1 , , N . Then solution x ( t ) = ( x 1 ( t ) , , x N ( t ) ) to the system (2) with initial conditions x k ( t ) = φ k ( t ) for t [ T 1 , 0 ] , k = 1 , , N satisfies inequalities x k ( t ) > 0 for all t > 0 .
Proof. 
The proof of this Lemma is very similar to the proof of Lemma 1. Let t 0 > 0 be the first moment when x m ( t 0 ) = 0 for some m. Thus, x k ( t ) > 0 for all t [ T 1 , t 0 ) and k = 1 , , N . From (2), we obtain
x m ( t 0 ) = x m ( 0 ) e ( 1 + γ ) t 0 + 0 t 0 e ( 1 + γ ) ( t s ) λ F ( x m ( s T 1 ) ) + γ x m 1 ( s ) d s .
The value x m ( t 0 ) is positive because x m ( 0 ) > 0 and inequalities F ( x ) 0 and x m 1 ( s ) > 0 hold for all t [ 0 , t 0 ) . This contradiction completes the proof. □

3. Some Properties of the Linear Equation

In this section, we establish several important properties of the linear part of Equation (5).
Consider the linear part of Equation (5)—a differential equation with one delay
x ˙ = ( 1 + γ ) x + γ x ( t T 2 ) ,
where γ and T 2 are positive parameters. The solutions and stability of Equation (8) are described by the roots of the characteristic equation
μ = ( 1 + γ ) + γ e μ T 2 .
Let’s denote as μ * a root of Equation (9) with a real part that is maximal from all roots of this equation.
Lemma 3.
Root μ * of Equation (9) is real and negative. Its multiplicity is equal to one.
Proof. 
Let’s represent μ in the form μ = Re   μ + i Im   μ . Then if we separate the real and imaginary parts of the Equation (9), we obtain the following system of equations
Re   μ + γ + 1 = γ exp ( Re   μ T 2 ) cos Im μ T 2 ,
Im   μ = γ exp ( Re   μ T 2 ) sin Im μ T 2 .
Consider the right and left sides of Equation (10) as two functions of the argument Re μ taking the value of Im μ as a parameter. Let h 1 ( Re   μ ) = Re   μ + γ + 1 and h 2 ( Re   μ ) = γ exp ( Re   μ T 2 ) cos Im μ T 2 . Then h 1 ( 0 ) = γ + 1 > γ h 2 ( 0 ) . If Re μ > 0 , then h 1 ( Re μ ) > h 1 ( 0 ) and h 2 ( Re μ ) < γ < h 1 ( 0 ) . Therefore, all roots Re μ of Equation (10) are less than zero; thus, Re μ * < 0 .
Let cos ( Im μ T 2 ) = 1 . Then Equation (10) has a root Re   μ = μ 0 < 0 , where μ 0 satisfies the equation μ 0 + γ + 1 = γ exp ( μ 0 T 2 ) . This equation has only one root because the left side an is increasing function and the right is a decreasing function. The multiplicity of this root is equal to one because the derivative on the left is positive whereas the derivative on the right is negative). If cos ( Im μ T 2 ) < 1 , then h 2 ( Re   μ ) < γ exp ( Re   μ T 2 ) ; therefore, all roots of Equation (10) with cos ( Im μ T 2 ) < 1 are less then μ 0 .
Let’s prove that μ 0 is the root of the system (10) and (11). If cos ( Im μ T 2 ) = 1 , then sin ( Im μ T 2 ) = 0 ; hence, Equation (11) in this case is equivalent to equation Im μ = 0 . Consequently, the real value μ 0 is the root of the system (10) and (11). Moreover, this root has the maximal real part from all roots of this system, there are no another complex roots with the same real part, and this root has a multiplicity equal to 1. □
Lemma 4.
Let φ ( t ) C [ t 0 T 2 , t 0 ] and φ ( t ) > 0 for t [ t 0 T 2 , t 0 ] . Then the solution to Equation (8) with initial condition φ ( t ) is positive for all t > t 0 .
Proof. 
This follows from Lemma 1 with F 0 . □
Next, we recall some important results from the monograph [24] on the form of the solutions to differential equations with delay.
Lemma 5
([24]). The solution to Equation (8) with initial condition x ( t ) = φ ( t ) for t [ t 0 T 2 , t 0 ] is
x ( t ) = L ( t t 0 ) φ ,
where L ( t t 0 ) is a linear compact operator such that for any μ * < μ 0 < 0 there exists c such that | L ( t t 0 ) φ | c ˜ e μ 0 ( t t 0 ) | | φ | | T 2 = c e μ 0 t | | φ | | T 2 , where | | φ | | T 2 = max s [ T 2 , 0 ] | φ ( t 0 + s ) | .
Lemma 6
(Variation-of-constants formula [24]). Let X ( t ) be the fundamental solution to Equation (8); i.e., solution with initial condition
X ( t ) = 0 , t < 0 ; 1 , t = 0 .
Then the solution to the linear inhomogeneous equation
x ˙ = ( 1 + γ ) x + γ x ( t T 2 ) + h ( t )
with initial conditions x ( t ) = φ ( t ) for t [ t 0 T 2 , t 0 ] is
x ( t ) = L ( t t 0 ) φ + t 0 t X ( t s ) h ( s ) d s .
Lemma 7.
Let φ ( t ) > 0 for all t [ t 0 T 2 , t 0 ] . Then there exists positive x 0 , c, c , and M such that for all t t 0 solution to Equation (8) with initial conditions x ( t ) = φ ( t ) ( t [ t 0 T 2 , t 0 ] ) is
x ( t ) = x 0 e μ * ( t t 0 ) ( 1 + w ( t ) ) , w ( t ) = Q ( t t 0 ) φ ,
where Q ( t t 0 ) is a compact linear operator such that
| w ( t ) | = | Q ( t t 0 ) φ | c ˜ e M ( t t 0 ) | | φ | | T 2 = c e M t | | φ | | T 2 and | w ˙ ( t ) | c e M t | | φ | | T 2 ,
where | | φ | | T 2 = max s [ T 2 , 0 ] | φ ( t 0 + s ) | .
Proof. 
The representation of x ( t ) and inequality | w ( t ) | c e M t | | φ | | T 2 follows from [24] (see Chapter 7). The value of x 0 is determined by the formula
x 0 = x ( t 0 ) + T 2 0 x ( s + t 0 ) e μ * ( s + T 2 ) d s > 0 .
To prove the inequality | w ˙ ( t ) | c e M t | | φ | | T 2 substitute the formula for x ( t ) into Equation (8), and after simplification we have
w ˙ ( t ) = ( 1 + γ + μ * ) w + γ e μ * T 2 w ( t T 2 ) ,
so
| w ˙ ( t ) | | 1 + γ + μ * | c e M t | | φ | | T 2 + γ e ( M μ * ) T 2 c e M t | | φ | | T 2 = c e M t | | φ | | T 2 .
Lemma 8.
Let 0 < φ ( t 0 ) p and 0 < φ ( t ) 1 + γ γ p for all t [ t 0 T 2 , t 0 ] . Then the solution x ( t ) to Equation (8) with initial condition x ( t ) = φ ( t ) ( t [ t 0 T 2 , t 0 ] ) is less than or equal to p for all t t 0 .
Proof. 
If Equation (8) is a linear inhomogeneous ODE, its solution has the form
x ( t ) = x ( t 0 ) e ( 1 + γ ) ( t t 0 ) + γ t 0 t e ( 1 + γ ) ( t τ ) x ( τ T 2 ) d τ .
The Second term in this formula can be estimated on the segment t [ t 0 , t 0 + T 2 ]
γ t 0 t e ( 1 + γ ) ( t τ ) x ( τ T 2 ) d τ γ 1 + γ γ p t 0 t e ( 1 + γ ) ( t τ ) d τ = p ( 1 e ( 1 + γ ) ( t t 0 ) ) .
Thus, on the segment t [ t 0 , t 0 + T 2 ]
x ( t ) x ( t 0 ) e ( 1 + γ ) ( t t 0 ) + p ( 1 e ( 1 + γ ) ( t t 0 ) ) p ,
and 0 < x ( t 0 + T 2 ) p . Then, if we replace t 0 by t 0 + T 2 in Formulas (13) and (14) we obtain the estimation of 0 < x ( t ) p on the segment t [ t 0 + T 2 , t 0 + 2 T 2 ] .
Acting like this, by induction we obtain estimation 0 < x ( t ) p for all t t 0 . □
Lemma 9.
Let λ be sufficiently large and x ( t ) = λ φ ( t ) for all t [ t 0 T 2 , t 0 ] , where φ ( t ) is a positive continuous function. Then there exists t 1 such that x ( t 1 ) = p , x ( t ) > p for all t [ t 0 , t 1 ) and for all t [ t 1 T 2 , t 1 ]
p 1 + γ γ x ( t ) p .
Proof. 
Using Lemma 7 we get x ( t ) = λ e μ * ( t t 0 ) x 0 ( 1 + w ( t ) ) , where x 0 is determined by (12); w ( t ) = Q ( t t 0 ) φ ; and | w ( t ) | c e M t | | φ | | T 2 for some positive c and M. Solving equation x ( t ) = p we find t 1 = t 0 + | μ * 1 | ln λ ( 1 + o ( 1 ) ) .
Thus, the values t [ t 1 T 2 , t 1 ] are large enough, so for such t
x ( t ) = p e μ * ( t t 1 ) ( 1 + o ( 1 ) ) .
Since μ * is a negative root of Equation (9), then e μ * T 2 = μ * γ 1 + 1 + γ γ < 1 + γ γ .
Let t [ t 1 T 2 , t 1 ] , then
x ( t ) = p e μ * ( t t 1 ) ( 1 + o ( 1 ) ) < p ( 1 + o ( 1 ) ) e μ * T 2 < p 1 + γ γ .
Thus, Lemma is proved. □

4. Relaxation Oscillations

In this section we prove that Equation (5) has a relaxation cycle, construct its asymptotics and study its main properties.
Let T be the maximum of delays ( T = max ( T 1 , T 2 ) ) and T m be the minimum of delays ( T m = min ( T 1 , T 2 ) ).
Consider the following set of initial conditions:
S = { φ ( t ) C [ T , 0 ] : φ ( 0 ) = p , φ ( t ) p t [ T , 0 ] , 1 + γ γ p φ ( t ) t [ T 2 , 0 ] } .
Examples of S are shown in Figure 1.
Denote as x φ ( t ) a solution to Equation (5) with initial condition φ S . Let’s construct an asymptotic approximation of solution x φ ( t ) .
1. Let t [ 0 , T 1 ] . For such t inequality x ( t T 1 ) > p holds, so F ( x ( t T 1 ) ) 0 and Equation (5) has the form of a linear differential equation with one delay (8). Therefore, solution x φ ( t ) satisfies the formula
x φ ( t ) = p e ( 1 + γ ) t + γ 0 t e ( 1 + γ ) ( t τ ) φ ( τ T 2 ) d τ .
From Lemma 8 we find that x φ ( t ) p for all t [ 0 , T 1 ] .
2. Let t [ T 1 , 2 T 1 ] . For these values of t, the inequality x φ ( t T 1 ) < p holds; therefore, F ( x φ ( t T 1 ) ) = f ( x φ ( t T 1 ) ) > 0 .
By Lemma 6 solution x φ ( t ) on this segment becomes asymptotically large. It is of the order O ( λ ) at λ + . Let’s write out the exact solution formulas. First, consider Equation (5) on the segment t [ T 1 , T 1 + T m ] .
In Equation (5) is considered to be an ordinary differential equation, then its solution has the form
x φ ( t ) = x ( T 1 ) e ( 1 + γ ) ( t T 1 ) + γ T 1 t e ( 1 + γ ) ( t τ ) x φ ( τ T 2 ) d τ + λ T 1 t e ( 1 + γ ) ( t τ ) F ( x φ ( τ T 1 ) ) d τ .
It follows from Formula (16), Lemma 8, and conditions f > 0 and λ 1 that there exists a positive value δ such that x ( T 1 + δ ) = p and δ = o ( 1 ) . On the interval t ( T 1 , T 1 + T m ] , the solution x φ ( t ) satisfies the formula
x φ ( t ) = λ T 1 t e ( 1 + γ ) ( t τ ) f ( x φ ( τ T 1 ) ) d τ ( 1 + o ( 1 ) ) .
If T m < T 1 , then on segment t [ T 1 + T m , 2 T 1 ] the following asymptotic formula for solution x φ ( t ) holds:
x φ ( t ) = γ T 1 + T m t e ( 1 + γ ) ( t τ ) x φ ( τ T 2 ) d τ + λ T 1 t e ( 1 + γ ) ( t τ ) f ( x φ ( τ T 1 ) ) d τ ( 1 + o ( 1 ) ) .
Note, that the first term on the right side of (18) is of the order O ( λ ) . Since x φ ( t T 2 ) > 0 (see Lemma 1), f ( x φ ( t T 1 ) ) > 0 on the segment t [ T 1 + T m , 2 T 1 ] , and 2 T 1 = O ( 1 ) , then on the segment t [ T 1 + T m , 2 T 1 + δ ] solution x φ ( t ) has the order O ( λ ) .
Figure 2a shows the solution of Equation (5) on a segment [ 0 , 2 T 1 ] , and Figure 2b shows an enlarged part in the neighbourhood of T 1 .
3. Let t 2 T 1 + δ . Then Equation (5) is the linear differential equation of the form (8) while solution x φ ( t T 1 ) > p . From Lemma 7, its solution is
x φ ( t ) = λ e μ * ( t 2 T 1 δ ) x 0 ( 1 + w ( t ) ) , | w ( t ) | c e M t max s [ T 2 , 0 ] | λ 1 x φ ( 2 T 1 + δ + s ) | .
Here x 0 , c , M > 0 . This function tends to zero; therefore, for some P > 2 T 1 + δ equality x φ ( P ) = p holds. Moreover,
P = | μ * 1 | ln λ + 2 T 1 + | μ * 1 | ( ln p ln x 0 ) | μ * 1 | ln ( 1 + w ( P ) ) + δ = | μ * 1 | ln λ ( 1 + o ( 1 ) ) .
By Lemma 9 on the segment t [ P T 2 , P ] , the solution satisfies the inequality p 1 + γ γ x φ ( t ) p and on the segment t [ P T , P ] condition x ( t ) p holds; therefore, x φ ( t + P ) S . Moreover, the solution x φ ( t + P ) for t [ T , 0 ] satisfies the following inequality:
x φ ( t + P ) = λ e μ * ( t + P 2 T 1 δ ) x 0 ( 1 + w ( t + P ) ) = e μ * t p ( 1 + w ( t + P ) ) ( 1 + w ( P ) ) 2 + 2 γ 2 + 2 γ + μ * p e μ * t .
Consider a set S 0 :
S 0 = { φ ( t ) C [ T , 0 ] : φ ( 0 ) = p , p φ ( t ) 2 + 2 γ 2 + 2 γ + μ * p e μ * t for t [ T , 0 ] } .
Note that S 0 is a non-empty bounded closed convex subset of S. In Figure 3 the examples of S 0 are shown.
Thus, we have proved that the Poincaré operator Π maps the set S to S 0 . Since S 0 is a subset of S, then Π maps S 0 to S 0 .
Since S 0 is a non-empty bounded closed convex set and Π is a compact operator (see [24]), we may use the Schauder fixed-point theorem [25]. As a result, we find that φ * S 0 S such that Π φ * = φ * ; therefore, solution x φ * ( t ) is a periodic solution.
Thus, we obtain the following statement.
Theorem 1.
For all sufficiently large λ, Equation (5) has a periodic solution x * ( t ) with initial conditions from the set S 0 with the period P = 2 T 1 + | μ * 1 | ln λ ( 1 + o ( 1 ) ) and amplitude O ( λ ) . Asymptotics of this solution are
x * ( t ) = p e μ * t ( 1 + o ( 1 ) )
on the segment t [ 0 , T 1 ] ,
x * ( t ) = λ T 1 t e ( 1 + γ ) ( t τ ) f ( p e μ * ( τ T 1 ) ) d τ ( 1 + o ( 1 ) )
on the segment t ( T 1 , T 1 + T m ] ,
x * ( t ) = λ T 1 t e ( 1 + γ ) ( t τ ) f ( p e μ * ( τ T 1 ) ) d τ ( 1 + o ( 1 ) ) + γ T 1 + T m t e ( 1 + γ ) ( t τ ) x * ( τ T 2 ) d τ
on the segment t ( T 1 + T m , 2 T 1 ] ,
x * ( t ) = x * ( 2 T 1 ) + T 2 0 x * ( 2 T 1 + s ) e μ * ( s + T 2 ) d s e μ * ( t 2 T 1 ) ( 1 + w ( t ) + o ( 1 ) )
on the segment t ( 2 T 1 , P ] , where | w ( t ) | c e M t and c > 0 , M > 0 .
Proof. 
We have proven that the compact Poincaré operator Π maps the non-empty bounded closed convex set S 0 to its pre-compact subset. That is why by the Schauder theorem there exists a fixed point of this operator φ * ( t ) S 0 [25]. From Formula (21) we obtained asymptotics of this fixed point: φ * ( t ) = p e μ * t ( 1 + o ( 1 ) ) at λ + .
If we get φ * ( t ) as an initial condition to Equation (5), then we obtain a periodic solution of Equation (5). Denote it as x * ( t ) . For all of the solutions of Equation (5) with initial conditions from the set S Formulas (15) and (17)–(19) hold. Hence, if we substitute φ * ( t ) into them and take into account Formula (12), then we obtain Formulas (22)–(25). □
Examples of periodic solutions to Equation (5) are shown in the Figure 4.
Note, that asymptotic Formula (25) contains a function w ( t ) that satisfies an exponential estimation | w ( t ) | c e M t , where c > 0 , M > 0 . On the interval t ( 2 T 1 , 2 T 1 + O ( 1 ) ) this function can make a significant contribution to the asymptotics of the solution. For example, the second (smaller) spikes on the period of the function x * ( t ) on the Figure 4a are described by this function.
Theorem 1 does not say anything about the stability of the constructed periodic relaxation solution. In the next section we show the sufficient conditions under which this cycle is exponentially orbitally stable.

5. The Stability of Relaxation Cycle

Let’s discuss the stability of relaxation cycle from Theorem 1. To prove sufficient stability conditions, we need the following auxiliary result:
Theorem 2.
Let f ( x ) > 0 be Lipschitz continuous on the segment [ 0 , p ] :
| f ( x ) f ( y ) | c f | x y | , x , y [ 0 , p ] ,
Then the Poincaré operator Π for the Equation (5) is a contraction operator on the set S 0 .
Proof. 
Denote as | | x | | = max t [ T , 0 ] | x ( t ) | —norm in the space C [ T , 0 ] and | | x | | T 2 = max t [ T 2 , 0 ] | x ( t ) | —norm in the space C [ T 2 , 0 ] .
We know that Π ( S 0 ) S 0 ; therefore, to prove the Theorem it is enough to prove that there exists 0 < q < 1 such that for any φ , ψ S 0 the inequality
| | Π φ Π ψ | | q | | φ ψ | |
holds.
Let φ and ψ be two functions from set S 0 , so let’s construct solutions to Equation (5) x φ ( t ) and x ψ ( t ) with initial conditions φ and ψ respectively.
We need to repeat all constructions of this section.
1. Let t [ 0 , T 1 ] . For these t, Equation (5) takes form of linear differential Equation (8). From Lemma 5 we found that its solutions for each t [ 0 , T 1 ] satisfied the
| x φ ( t ) x ψ ( t ) | = | L ( t ) ( φ ψ ) | c 1 e μ 0 t | | φ ψ | | c 1 | | φ ψ | | ,
where c 1 is some positive constant, and μ * < μ 0 < 0 .
2. Consider x φ ( t ) for t [ T 1 , 2 T 1 + δ φ ] and x ψ ( t ) for t [ T 1 , 2 T 1 + δ ψ ] . As before, denote as δ φ and δ ψ the first positive roots of equations x φ ( T 1 + δ φ ) = p , and x ψ ( T 1 + δ ψ ) = p , respectively. From Formula (16) and inequality f ( p ) > 0 , the values of δ φ and δ ψ are O ( λ 1 ) at λ + . Without any loss of generality, we assumed that δ φ is less than δ ψ .
By Lemma 6 we may represent solutions as follows:
x φ ( t ) = L ( t T 1 ) x φ T 1 + λ T 1 t X ( t s ) f ( x φ ( s T 1 ) ) d s ,
x ψ ( t ) = L ( t T 1 ) x ψ T 1 + λ T 1 t X ( t s ) f ( x ψ ( s T 1 ) ) d s ,
where X ( t ) is the fundamental solution to (8). Here x φ τ and x ψ τ are functions x φ ( τ + t ) and x ψ ( τ + t ) for t [ T 2 , 0 ] . Note that L ( t T 1 ) x φ T 1 = L ( t ) φ and L ( t T 1 ) x ψ T 1 = L ( t ) ψ . From this formula, the difference between solutions on the segment [ T 1 , 2 T 1 ] satisfies the inequalities
| x φ ( t ) x ψ ( t ) | | L ( t ) ( φ ψ ) | + λ T 1 t | X ( t s ) ( f ( x φ ( s T 1 ) ) f ( x ψ ( s T 1 ) ) ) | d s c 1 | | φ ψ | | + λ c f c 2 T 1 t | x φ ( t T 1 ) x ψ ( t T 1 ) | d s ( c 1 + c 3 λ ( t T 1 ) ) | | φ ψ | | .
Here c 2 = max t [ 0 , 2 T 1 ] | X ( t ) | and c 3 = c f c 2 c 1 are positive constants. In particular,
| x φ ( T 1 + δ ψ ) x ψ ( T 1 + δ ψ ) | ( c 1 + λ c 3 δ ψ ) | | φ ψ | | c 4 | | φ ψ | | .
In the same way inequality (26) can be generalized to the segment t [ 2 T 1 , 2 T 1 + min ( δ φ , δ ψ ) ] .
Let’s estimate δ ψ δ φ . Using (16), we get
0 = p p = x φ ( T 1 + δ φ ) x ψ ( T 1 + δ ψ ) = x φ ( T 1 + δ φ ) x ψ ( T 1 + δ φ ) + x ψ ( T 1 + δ φ ) x ψ ( T 1 + δ ψ ) = x φ ( T 1 + δ φ ) x ψ ( T 1 + δ φ ) + x ˙ ψ ( T 1 + θ ) ( δ φ δ ψ ) .
Here θ is a value between δ ψ and δ φ . It is asymptotically small and positive, so x ψ ( θ ) = p + o ( 1 ) < p . Let’s estimate x ˙ ψ ( T 1 + θ ) .
x ˙ ψ ( T 1 + θ ) = ( 1 + γ ) x ψ ( T 1 + θ ) + γ x ψ ( T 1 + θ T 2 ) + λ f ( x ψ ( θ ) ) c 5 λ ,
where c 5 = 1 2 f ( p ) > 0 . Thus from (27),
| δ φ δ ψ | = | x ˙ ψ ( T 1 + θ ) | 1 | x φ ( T 1 + δ φ ) x ψ ( T 1 + δ φ ) | c δ λ 1 | | φ ψ | | , c δ = 2 c 4 c 5 1 > 0 .
Here we also note that for all t [ T , 0 ]
λ c m a x x φ ( 2 T 1 + δ φ + t ) , λ c m a x x ψ ( 2 T 1 + δ ψ + t ) ,
where
c m a x = 2 T 1 max t [ 0 , 2 T 1 ] | X ( t ) | max x [ 0 , p ] f ( x ) ,
and
x φ ( 2 T 1 + δ φ ) λ c 6 , x ψ ( 2 T 1 + δ ψ ) λ c 6 ,
where c 6 = 1 2 T 1 min t [ 0 , 2 T 1 ] | X ( t ) | min x [ 0 , p ] f ( x ) .
3. Consider each solution on the time interval t from 2 T 1 + δ φ and 2 T 1 + δ ψ till the moment when it returns to the set S 0 . By Lemma 7, both solutions can be represented in the form
x φ ( t ) = λ x 1 e μ * ( t 2 T 1 δ φ ) ( 1 + w φ ( t ) ) , x ψ ( t ) = λ x 2 e μ * ( t 2 T 1 δ ψ ) ( 1 + w ψ ( t ) ) ,
where w φ ( t ) = λ 1 Q ( t 2 T 1 δ φ ) x φ 2 T 1 + δ φ and w ψ ( t ) = λ 1 Q ( t 2 T 1 δ ψ ) x ψ 2 T 1 + δ ψ . For short, we denote Q ( t 2 T 1 δ φ ) as Q φ ( t ) and denote Q ( t 2 T 1 δ ψ ) as Q ψ ( t ) .
Values λ x 1 and λ x 2 are determined by formula (12); therefore
λ | x 1 x 2 |     | x φ ( 2 T 1 + δ φ ) x ψ ( 2 T 1 + δ φ ) | + | x ψ ( 2 T 1   +   δ φ ) x ψ ( 2 T 1 + δ ψ ) | + | T 2 0 e μ * ( s + T 2 ) ( x φ ( s + 2 T 1 + δ φ ) x ψ ( s + 2 T 1 + δ φ ) ) d s | + | T 2 0 e μ * ( s + T 2 ) ( x ψ ( s + 2 T 1 + δ φ ) x ψ ( s + 2 T 1 + δ ψ ) ) d s | λ c 7 | | φ ψ | | + c 8 λ | δ φ δ ψ | .
Here were used that for t [ T 1 , 2 T 1 + δ ψ ] from (28) and equality (5) it follows that
| x ˙ ψ ( t ) | ( 1 + γ ) | x ψ ( t ) | + γ | x ψ ( t T 2 ) |   +   λ f ( x ψ ( t T 1 ) ) λ c 9 ,
where c 9 = ( 1 + 2 γ ) c m a x + max x [ 0 , p ] f ( x ) .
Therefore, there exists c x > 0 such that
| x 1 x 2 |   c x | | φ ψ | | .
In addition, from Formula (12) it follows that x 1 , 2 c 6 > 0 .
Let P φ and P ψ be the first moments of time such that P φ > 2 T 1 , P ψ > 2 T 1 , x φ ( P φ ) = p , and x ψ ( P ψ ) = p . They are represented by the Formula (20). Let’s estimate the difference P φ P ψ :
| P φ P ψ |     | μ * 1 | | ln x 1 ln x 2 |   +   | μ * 1 | | ln ( 1 + w φ ( P φ ) ) ln ( 1 + w ψ ( P ψ ) ) | + | δ ψ δ φ |   ( | μ * 1 | c 6 1 c x + c δ λ 1 ) | | φ ψ | | + c 10 | w φ ( P φ ) w ψ ( P ψ ) | .
To estimate | w φ ( P φ ) w ψ ( P ψ ) | consider a more general inequality for t [ T , 0 ]
| w φ ( P φ + t ) w ψ ( P ψ + t ) |     | w φ ( P φ + t ) w ψ ( P φ + t ) |   +   | w ψ ( P φ + t ) w ψ ( P ψ + t ) |   λ 1 | Q φ ( P φ + t ) ( x φ 2 T 1 + δ φ x ψ 2 T 1 + δ ψ ) |   +   λ 1 | Q ψ ( P φ + t δ φ + δ ψ ) x ψ 2 T 1 + δ ψ Q ψ ( P φ + t ) x ψ 2 T 1 + δ ψ |   +   | w ˙ ψ ( θ ) | | P φ P ψ |   λ 1 | Q φ ( P φ + t ) ( x φ 2 T 1 + δ φ x ψ 2 T 1 + δ φ ) |   +   λ 1 | Q φ ( P φ + t ) ( x ψ 2 T 1 + δ φ x ψ 2 T 1 + δ ψ ) |   +   | w ψ ( P φ + t δ φ + δ ψ ) w ψ ( P φ + t ) |   + | w ˙ ψ ( θ ) | | P φ P ψ |   c e M ( P φ + t ) c 11 | | φ ψ | | + c λ 1 e M ( P φ + t ) | x ˙ ψ ( θ 1 ) | | δ φ δ ψ |   +   | w ˙ ψ ( θ 2 ) | | δ φ δ ψ |   +   | w ˙ ψ ( θ ) | | P φ P ψ | .
Here c 11 = 2 c 3 T 1 , θ 1 is a time moment between 2 T 1 + δ φ and 2 T 1 + δ ψ , θ 2 —between P φ + t and P φ + t δ φ + δ ψ , and θ is a time moment between P φ + t and P ψ + t .
From Formula (29) | x ˙ ψ ( θ 1 ) |   λ c 9 . By Lemma 7 and inequality (28) | w ˙ ψ ( t ) | λ 1 c e M t | | x ψ 2 T 1 + δ ψ | | c c m a x e M t and, since θ 2 and θ are large enough, | w ˙ ψ ( θ 2 ) | and | w ˙ ψ ( θ ) | are small, so 1 c 10 | w ˙ ψ ( θ ) | > 1 / 2 . Then we have
( 1 c 10 | w ˙ ψ ( θ ) | ) | P φ P ψ |   | μ * 1 | c 6 1 c x + c δ λ 1 + c c 10 ( c 11 + c 9 c δ λ 1 ) e M P φ + c 10 c c δ c m a x λ 1 e M θ 2 | | φ ψ | | ,
so there exists c P > 0 such that
| P φ P ψ |   c P | | φ ψ | | .
Moreover, from (30) it follows that there exists c w > 0 such that the next inequality is true:
| | w φ ( P φ ) w ψ ( P ψ ) | | e M P φ c w | | φ ψ | | .
Represent solutions x φ and x ψ in the form
x φ ( t + P φ ) = λ e μ * ( t + P φ 2 T 1 δ φ ) x 1 ( 1 + w φ ( t + P φ ) ) = x φ ( P φ ) e μ * t 1 + w φ ( t + P φ ) 1 + w φ ( P φ ) = p e μ * t 1 + w φ ( t + P φ ) 1 + w φ ( P φ )
and
x ψ ( t + P ψ ) = p e μ * t 1 + w ψ ( t + P ψ ) 1 + w ψ ( P ψ ) .
Let t [ T , 0 ] . For these t let’s estimate the difference
| x φ ( t + P φ ) x ψ ( t + P ψ ) | p e μ * t | 1 + w φ ( t + P φ ) 1 + w φ ( P φ ) 1 + w ψ ( t + P ψ ) 1 + w ψ ( P ψ ) | .
Here
| 1 + w φ ( t + P φ ) 1 + w φ ( P φ ) 1 + w ψ ( t + P ψ ) 1 + w ψ ( P ψ ) | = | ( w φ ( t + P φ ) w ψ ( t + P ψ ) ) ( 1 + w ψ ( P ψ ) ) + ( 1 + w ψ ( t + P ψ ) ) ( w ψ ( P ψ ) w φ ( P φ ) ) ( 1 + w φ ( P φ ) ) ( 1 + w ψ ( P ψ ) ) | 2 | w φ ( t + P φ ) w ψ ( t + P ψ ) | + 2 | w ψ ( P ψ ) w φ ( P φ ) | 4 | | w φ P φ w ψ P ψ | | 4 c w e M P φ | | φ ψ | | = 4 c w λ M / μ * ( 1 + o ( 1 ) ) | | φ ψ | | .
Here, as before, we denote as w φ τ and w ψ τ functions w φ ( τ + t ) and w ψ ( τ + t ) , respectively, for t [ T 2 , 0 ] .
So for all t [ T , 0 ]
| x φ ( t + P φ ) x ψ ( t + P ψ ) | 5 c w p e μ * T λ M / μ * | | φ ψ | | .
Let q = 5 c w p e μ * T λ M / μ * . Since M > 0 , μ * < 0 then for sufficiently large λ the value q is less than 1 and is small enough. Therefore, for all φ , ψ S 0
| | Π φ Π ψ | | q | | φ ψ | | .
Thus, Π : S 0 S 0 is a contraction operator. The theorem is proved. □
Theorem 3.
Let F satisfy conditions from Theorem 2. Then for all sufficiently large values of λ, Equation (5) has a unique periodic solution with initial conditions from S 0 . This solution is exponentially orbitally stable. Asymptotics of this solution are given in the Theorem 1.
Proof. 
Due to the contraction mapping principle, operator Π has a unique fixed point φ * S 0 that corresponds to the periodic solution x * ( t ) of (5). For any φ S 0 , the sequence Π n φ tends to the φ * exponentially fast, so solution x * ( t ) is exponentially orbitally stable. □

6. Travelling Waves

Let’s discuss the existence of travelling wave solutions in the system of N-coupled oscillators (2).
Theorem 4.
Let λ be large enough and N be a large integer of the order O ( ln λ ) , i.e., N = c ( λ ) ln λ , where c ( λ ) is a positive bounded function such that N is an integer. Let n > 0 be an integer and n c 1 ( λ ) < ln ( 1 + 1 γ ) . Then there exist positive T 2 , 1 , T 2 , 2 , …, T 2 , n such that system (2) has at least n travelling wave solutions in the form
x k , m = x * m ( t ( k 1 ) T 2 , m ) , m = 1 , , n , k = 1 , , N ,
where x * m ( t ) is a periodic relaxation solution to Equation (5) from Theorem 1 for corresponding T 2 , m .
Proof. 
From Theorem 1 for any T 2 we may construct a periodic solution x * ( t ) to Equation (5) with period P ( T 2 ) . We need T 2 such that condition (6) holds. Thus, we need to prove the existence of at least n roots of equation
T 2 N = m P ( T 2 )
for integers m = 1 , , n . We know that P ( T 2 ) = 1 μ * ln λ ( 1 + o ( 1 ) ) . To determine T 2 , therefore, we have the equation
μ * T 2 = m c ( λ ) ( 1 + o ( 1 ) ) .
The product | μ * T 2 | takes values in the interval ( 0 , ln ( 1 + 1 γ ) ) ; therefore, Equation (31) has n solutions T 2 = T 2 , m ( m = 1 , , n ) if and only if n c 1 ( λ ) < ln ( 1 + 1 γ ) . □
Figure 5 and Figure 6 show examples of stable travelling waves for some parameter values. In the Figure 5 there is solution to the system (2) for λ = 500 , γ = 0.38 , p = 5 , T 1 = 1 , f ( x ) 5 , N = 7 . Here, ln λ 6.21 , so, c ( λ ) 1.13 and c 1 ( λ ) 0.88 < 1.29 ln ( 1 + γ 1 ) . Thus, conditions of the Theorem 4 are fulfilled. In this case T 2 1.77 . The corresponding periodic solution x * ( t ) to Equation (5) is in Figure 4a. In Figure 6 there is solution of the system (2) for λ = 500 , γ = 0.15 , p = 5 , T 1 = 1 , f ( x ) 5 , N = 7 . As in Figure 5  c 1 ( λ ) 0.88 . It is less than 2.03 ln ( 1 + γ 1 ) . In this case T 2 1.33 . The corresponding periodic solution x * ( t ) to the Equation (5) is in Figure 4b.
Note, that system (2) has a homogeneous relaxation cycle. We may treat it as a travelling wave corresponding to T 2 = 0 . Asymptotics of this cycle coincided with the asymptotics of the homogeneous relaxation cycle of the system (3). For asymptotics of this cycle, see [18]. Note that this cycle does not depend on γ . An example of such a solution is in Figure 7.

7. Discussion and Conclusions

We studied the periodic solutions to equation with two delays and compactly supported nonlinearity (5) under condition λ 1 . We chose a set of initial conditions and proved that translation along the trajectories operator Π mapped this set to itself. We proved that Π had a fixed point, and the solution to Equation (5) corresponding to this fixed point was periodic. This cycle had a large period of the order O ( ln λ ) , and an amplitude of the order O ( λ ) . With additional requirements for the function f ( x ) we proved that Π is a contraction operator and a found cycle is exponentially orbitally stable.
This cycle was used to construct travelling waves of the system of N-coupled neuron-like oscillators (2) when N was asymptotically large (on the order of O ( ln λ ) ).
The results allowed for some generalizations.
1. The function f ( x ) may take negative values on the segment x [ 0 , p ] . The only reason we required positiveness of function f ( x ) on the segment x [ 0 , p ] while constructing asymptotics was to guarantee the positiveness of the solution on the interval [ T 1 , 2 T 1 ] . Although for the positiveness of solution, it was sufficient to demand the positiveness of integral.
T 1 t e ( 1 + γ ) ( t τ ) f ( p e μ * ( τ T 1 ) ) d τ > 0 .
Therefore, under condition (32) (instead of condition f ( x ) > 0 on the x [ 0 , p ] ) there also existed a relaxation cycle of Equation (5).
2. If we required that f ( x ) < 0 for x [ p , 0 ) and that f ( x ) was Lipschitz continuous on the segment [ p , 0 ] and took negative initial conditions and repeated all the constructions and proofs in the paper then we obtained an exponentially orbitally stable negative relaxation cycle. Therefore, under these conditions we obtained multistability in Equation (5): positive and negative exponentially orbitally stable relaxation cycles coexisted. In Figure 8, the exponentially orbitally stable negative relaxation cycle of the model (5) is shown.
3. All previous results were given for sufficiently large values of λ . However, if λ were not large enough, the method still worked: the fixed points and cycles of the operator Π corresponded to the periodic solutions of (5).
In Figure 9 there is an example of a complex periodic solution of (5) for λ = 500 .
Note that this solution was not predicted by Theorem 1. This is not a contradiction because increasing λ leads to the disappearance of this cycle, which exists when λ is not large enough.

Author Contributions

Conceptualization, A.K. and I.K.; methodology, A.K. and I.K.; software, S.K.; validation, A.K., I.K. and S.K.; formal analysis, A.K., I.K. and S.K.; investigation, A.K., I.K. and S.K.; data curation, S.K.; writing—original draft preparation, A.K. and I.K.; writing—review and editing, A.K. and I.K.; visualization, S.K.; supervision, A.K.; project administration, A.K.; funding acquisition, A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Russian Science Foundation (project No. 22-11-00209), https://rscf.ru/en/project/22-11-00209/, accessed on 3 April 2023.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arino, J.; Van Den Driessche, P. Time delays in epidemic models. In Delay Differential Equations and Applications; Springer: Berlin/Heidelberg, Germany, 2006; pp. 539–578. [Google Scholar]
  2. Erneux, T. Applied Delay Differential Equations; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2009. [Google Scholar] [CrossRef]
  3. Smith, H. An Introduction to Delay Differential Equations with Applications to the Life Sciences; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar] [CrossRef]
  4. Brauer, F.; Castillo-Chavez, C. Mathematical Models in Population Biology and Epidemiology; Springer: New York, NY, USA, 2012. [Google Scholar] [CrossRef]
  5. Cushing, J.M. Integrodifferential Equations and Delay Models in Population Dynamics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  6. Ponomarenko, V.; Prokhorov, M.; Karavaev, A.; Kulminskiy, D. An experimental digital communication scheme based on chaotic time-delay system. Nonlinear Dyn. 2013, 74, 1013–1020. [Google Scholar] [CrossRef]
  7. Kaschenko, S. Models of Wave Memory; Springer: Cham, Switzerland, 2015. [Google Scholar] [CrossRef]
  8. Kilias, T.; Kelber, K.; Mogel, A.; Schwarz, W. Electronic chaos generators—Design and applications. Int. J. Electron. 1995, 79, 737–753. [Google Scholar] [CrossRef]
  9. Kilias, T.; Mogel, A.; Schwarz, W. Generation and application of broadband signals using chaotic electronic systems. In Nonlinear Dynamics: New Theoretical and Applied Results; Akademie: Berlin, Germany, 1995; pp. 92–111. [Google Scholar]
  10. an der Heiden, U.; Mackey, M.C. The dynamics of production and destruction: Analytic insight into complex behavior. J. Math. Biol. 1982, 16, 75–101. [Google Scholar] [CrossRef]
  11. Kaschenko, S.A. Asymptotic Behavior of Relaxation Oscillations in Systems of Difference-Differential Equations with a Compactly Supported Nonlinearity. I. Differ. Equ. 1995, 31, 1275–1285. [Google Scholar]
  12. Pikovsky, A.; Rosenblum, M.; Kurths, J. Synchronization: A Universal Concept in Nonlinear Science; Cambridge University Press: Cambridge, UK, 2002. [Google Scholar]
  13. Strogatz, S.H. Sync: How Order Emerges from Chaos in the Universe, Nature, and Daily Life; Hachette: London, UK, 2012. [Google Scholar]
  14. Rodrigues, F.A.; Peron, T.K.D.; Ji, P.; Kurths, J. The Kuramoto model in complex networks. Phys. Rep. 2016, 610, 1–98. [Google Scholar] [CrossRef] [Green Version]
  15. Kashchenko, A.A. Relaxation Cycles in a Model of Two Weakly Coupled Oscillators with Sign-Changing Delayed Feedback. Theor. Math. Phys. 2020, 202, 381–389. [Google Scholar] [CrossRef]
  16. Kashchenko, A.A. Relaxation modes of a system of diffusion coupled oscillators with delay. Commun. Nonlinear Sci. Numer. Simul. 2021, 93, 105488. [Google Scholar] [CrossRef]
  17. Kashchenko, A.A. Influence of coupling on the dynamics of three delayed oscillators. Izv. VUZ Appl. Nonlinear Dyn. 2021, 29, 869–891. [Google Scholar] [CrossRef]
  18. Kashchenko, A.A. Dependence of the dynamics of a model of coupled oscillators on the number of oscillators. Dokl. Math. 2021, 104, 355–359. [Google Scholar] [CrossRef]
  19. Glyzin, S.D.; Kolesov, A.Y.; Rozov, N.K. Periodic traveling-wave-type solutions in circular chains of unidirectionally coupled equations. Theor. Math. Phys. 2013, 175, 499–517. [Google Scholar] [CrossRef]
  20. Preobrazhenskaya, M.M. Discrete traveling waves in a relay system of Mackey–Glass equations with two delays. Theor. Math. Phys. 2021, 207, 827–840. [Google Scholar] [CrossRef]
  21. Marino, I.; Pérez-Muñuzuri, V.; Pérez-Villar, V.; Sánchez, E.; Matıas, M. Interaction of chaotic rotating waves in coupled rings of chaotic cells. Phys. D Nonlinear Phenom. 1999, 128, 224–235. [Google Scholar] [CrossRef]
  22. Klinshov, V.; Shchapin, D.; Yanchuk, S.; Wolfrum, M.; D’Huys, O.; Nekorkin, V. Embedding the dynamics of a single delay system into a feed-forward ring. Phys. Rev. E 2017, 96, 042217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Rabinovich, M.I.; Varona, P.; Selverston, A.I.; Abarbanel, H.D. Dynamical principles in neuroscience. Rev. Mod. Phys. 2006, 78, 1213. [Google Scholar] [CrossRef] [Green Version]
  24. Hale, J.; Sjoerd, M. Introdution to Functional Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  25. Lusternik, L.; Sobolev, V. Elements of Functional Analysis; Gordon and Breach Science Publishers: Philadelphia, PA, USA, 1968. [Google Scholar]
Figure 1. Examples of set S for (a) T 1 < T 2 ; (b) T 2 < T 1 .
Figure 1. Examples of set S for (a) T 1 < T 2 ; (b) T 2 < T 1 .
Mathematics 11 02827 g001
Figure 2. Solution of the Equation (5) for (a) t [ 0 , 2 T 1 ] ; (b) t in neighbourhood of T 1 .
Figure 2. Solution of the Equation (5) for (a) t [ 0 , 2 T 1 ] ; (b) t in neighbourhood of T 1 .
Mathematics 11 02827 g002
Figure 3. Example of set S 0 for (a) T 1 < T 2 ; (b) T 2 < T 1 .
Figure 3. Example of set S 0 for (a) T 1 < T 2 ; (b) T 2 < T 1 .
Mathematics 11 02827 g003
Figure 4. Solution of the equation with two delays (5). Parameters: (a) λ = 500 , γ = 0.38 , p = 5 , T 1 = 1 , T 2 = 1.77 , f ( x ) = 5 ; (b) λ = 500 , γ = 0.15 , p = 5 , T 1 = 1 , T 2 = 1.33 , f ( x ) = 5 .
Figure 4. Solution of the equation with two delays (5). Parameters: (a) λ = 500 , γ = 0.38 , p = 5 , T 1 = 1 , T 2 = 1.77 , f ( x ) = 5 ; (b) λ = 500 , γ = 0.15 , p = 5 , T 1 = 1 , T 2 = 1.33 , f ( x ) = 5 .
Mathematics 11 02827 g004
Figure 5. Example of travelling wave solution to (2). Parameters: N = 7 ; λ = 500 ; γ = 0.38 ; p = 5 ; T 1 = 1 ; f ( x ) = 5 .
Figure 5. Example of travelling wave solution to (2). Parameters: N = 7 ; λ = 500 ; γ = 0.38 ; p = 5 ; T 1 = 1 ; f ( x ) = 5 .
Mathematics 11 02827 g005
Figure 6. Example of travelling wave solution to (2). Parameters: N = 7 ; λ = 500 ; γ = 0.15 ; p = 5 ; T 1 = 1 ; f ( x ) = 5 .
Figure 6. Example of travelling wave solution to (2). Parameters: N = 7 ; λ = 500 ; γ = 0.15 ; p = 5 ; T 1 = 1 ; f ( x ) = 5 .
Mathematics 11 02827 g006
Figure 7. Example of travelling wave solution (homogeneous cycle) of (2). Parameters: N = 7 ; λ = 500 ; p = 5 ; T 1 = 1 ; f ( x ) = 5 .
Figure 7. Example of travelling wave solution (homogeneous cycle) of (2). Parameters: N = 7 ; λ = 500 ; p = 5 ; T 1 = 1 ; f ( x ) = 5 .
Mathematics 11 02827 g007
Figure 8. Solution of the equation with two delays (5). λ = 550 ; γ = 0.38 ; p = 5 ; T 1 = 1 ; T 2 = 1.77 ; f ( x ) = 5 .
Figure 8. Solution of the equation with two delays (5). λ = 550 ; γ = 0.38 ; p = 5 ; T 1 = 1 ; T 2 = 1.77 ; f ( x ) = 5 .
Mathematics 11 02827 g008
Figure 9. Solution of the equation with two delays (5). λ = 500 ; γ = 0.3 ; p = 5 ; T 1 = 1 ; T 2 = 3 ; f ( x ) = 5 .
Figure 9. Solution of the equation with two delays (5). λ = 500 ; γ = 0.3 ; p = 5 ; T 1 = 1 ; T 2 = 3 ; f ( x ) = 5 .
Mathematics 11 02827 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kashchenko, A.; Kashchenko, I.; Kondratiev, S. Travelling Waves in the Ring of Coupled Oscillators with Delayed Feedback. Mathematics 2023, 11, 2827. https://doi.org/10.3390/math11132827

AMA Style

Kashchenko A, Kashchenko I, Kondratiev S. Travelling Waves in the Ring of Coupled Oscillators with Delayed Feedback. Mathematics. 2023; 11(13):2827. https://doi.org/10.3390/math11132827

Chicago/Turabian Style

Kashchenko, Alexandra, Ilia Kashchenko, and Sergey Kondratiev. 2023. "Travelling Waves in the Ring of Coupled Oscillators with Delayed Feedback" Mathematics 11, no. 13: 2827. https://doi.org/10.3390/math11132827

APA Style

Kashchenko, A., Kashchenko, I., & Kondratiev, S. (2023). Travelling Waves in the Ring of Coupled Oscillators with Delayed Feedback. Mathematics, 11(13), 2827. https://doi.org/10.3390/math11132827

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop