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
where
,
, and
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
of
N oscillators of the form (
1) with parameter
. Here, parameter
is positive and sufficiently large
, delay time
and coupling parameter
are positive, and feedback function
F is a compactly supported positive function:
where
p is some positive constant, and
is a piece-wise continuous and bounded function on the segment
and
for all
. In the segment
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
,
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
under condition
was studied [
17,
18]. It was proved that when
all oscillators were synchronized, and when
and
N was even two-cluster synchronization was observed. When
and
, 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
where
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
from (
4) must satisfy the equation with two delays
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
, where
under some additional requirements on the function
f.
Condition
lead us to
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
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
). 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
. 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
for any arbitrary values of parameters
,
,
, and
. 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
. 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
, where
T is the delay time of the model. Then we integrated the model using the method of steps [
24]. On the first step
function
is a known function, the initial condition
. That is why, on this segment, we considered our initial equation with delay as an ordinary differential equation with inhomogeneity depending on
. We found asymptotics of the solution to this equation at
. After that, we considered our model on the segment
as an ordinary differential equation with inhomogeneity depending on
, 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
there existed a moment
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
. That is why using the Schauder fixed-point theorem [
25] we concluded that there was a fixed point
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.
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 () and be the minimum of delays ().
Consider the following set of initial conditions:
Denote as
a solution to Equation (
5) with initial condition
. Let’s construct an asymptotic approximation of solution
.
1. Let
. For such
t inequality
holds, so
and Equation (
5) has the form of a linear differential equation with one delay (
8). Therefore, solution
satisfies the formula
From Lemma 8 we find that
for all
.
2. Let . For these values of t, the inequality holds; therefore, .
By Lemma 6 solution
on this segment becomes asymptotically large. It is of the order
at
. Let’s write out the exact solution formulas. First, consider Equation (
5) on the segment
.
In Equation (
5) is considered to be an ordinary differential equation, then its solution has the form
It follows from Formula (
16), Lemma 8, and conditions
and
that there exists a positive value
such that
and
. On the interval
, the solution
satisfies the formula
If
, then on segment
the following asymptotic formula for solution
holds:
Note, that the first term on the right side of (
18) is of the order
. Since
(see Lemma 1),
on the segment
, and
, then on the segment
solution
has the order
.
Figure 2a shows the solution of Equation (
5) on a segment
, and
Figure 2b shows an enlarged part in the neighbourhood of
.
3. Let
. Then Equation (
5) is the linear differential equation of the form (
8) while solution
. From Lemma 7, its solution is
Here
. This function tends to zero; therefore, for some
equality
holds. Moreover,
By Lemma 9 on the segment
, the solution satisfies the inequality
and on the segment
condition
holds; therefore,
. Moreover, the solution
for
satisfies the following inequality:
Consider a set
:
Note that
is a non-empty bounded closed convex subset of
S. In
Figure 3 the examples of
are shown.
Thus, we have proved that the Poincaré operator maps the set S to . Since is a subset of S, then maps to .
Since
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
such that
; therefore, solution
is a periodic solution.
Thus, we obtain the following statement.
Theorem 1. For all sufficiently large λ, Equation (5) has a periodic solution with initial conditions from the set with the period and amplitude . Asymptotics of this solution areon the segment ,on the segment ,on the segment ,on the segment , where and . Proof. We have proven that the compact Poincaré operator
maps the non-empty bounded closed convex set
to its pre-compact subset. That is why by the Schauder theorem there exists a fixed point of this operator
[
25]. From Formula (
21) we obtained asymptotics of this fixed point:
at
.
If we get
as an initial condition to Equation (
5), then we obtain a periodic solution of Equation (
5). Denote it as
. 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
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
that satisfies an exponential estimation
, where
. On the interval
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
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.
7. Discussion and Conclusions
We studied the periodic solutions to equation with two delays and compactly supported nonlinearity (
5) under condition
. 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
, and an amplitude of the order
. With additional requirements for the function
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
).
The results allowed for some generalizations.
1. The function
may take negative values on the segment
. The only reason we required positiveness of function
on the segment
while constructing asymptotics was to guarantee the positiveness of the solution on the interval
. Although for the positiveness of solution, it was sufficient to demand the positiveness of integral.
Therefore, under condition (
32) (instead of condition
on the
) there also existed a relaxation cycle of Equation (
5).
2. If we required that
for
and that
was Lipschitz continuous on the segment
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
.
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.