Next Article in Journal
Relative Importance of Certain Factors Affecting Air Exchange in a High-Altitude Single-Heading Tunnels Based on the Numerical Simulation Method
Next Article in Special Issue
Kuramoto Model with Delay: The Role of the Frequency Distribution
Previous Article in Journal
Reducing the Reality Gap Using Hybrid Data for Real-Time Autonomous Operations
Previous Article in Special Issue
Improved Properties of Positive Solutions of Higher Order Differential Equations and Their Applications in Oscillation Theory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relaxation Oscillations in the Logistic Equation with Delay and Modified Nonlinearity

by
Alexandra Kashchenko
*,† and
Sergey Kashchenko
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(7), 1699; https://doi.org/10.3390/math11071699
Submission received: 28 February 2023 / Revised: 27 March 2023 / Accepted: 31 March 2023 / Published: 2 April 2023
(This article belongs to the Special Issue Advances in Delay Differential Equations)

Abstract

:
We consider the dynamics of a logistic equation with delays and modified nonlinearity, the role of which is to bound the values of solutions from above. First, the local dynamics in the neighborhood of the equilibrium state are studied using standard bifurcation methods. Most of the paper is devoted to the study of nonlocal dynamics for sufficiently large values of the ‘Malthusian’ coefficient. In this case, the initial equation is singularly perturbed. The research technique is based on the selection of special sets in the phase space and further study of the asymptotics of all solutions from these sets. We demonstrate that, for sufficiently large values of the Malthusian coefficient, a ‘stepping’ of periodic solutions is observed, and their asymptotics are constructed. In the case of two delays, it is established that there is attractor in the phase space of the initial equation, whose dynamics are described by special nonlinear finite-dimensional mapping.

1. Introduction

The logistic equation with delay is the name given to the equation
u ˙ = λ [ 1 u ( t T ) ] u .
Being a natural extension of the classical logistic equation, it arises in many applied problems, but most of all it is used in problems of mathematical ecology (see, for example [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15]). The value u describes the dynamics of changes in the number of specimens in the population or density of an isolated biological population living in a homogeneous environment. Therefore, it makes sense to consider only non-negative solutions of Equation (1). We note that the solution of Equation (1) with a non-negative initial function φ ( s ) C [ T , 0 ] specified at some value t 0 (i.e., u ( t 0 + s ) = φ ( s ) ) remains non-negative for all t > t 0 . The coefficient λ is called the Malthusian coefficient of linear growth and the parameter T > 0 is called the delay time. It is associated with the age of animal units capable of procreation in the corresponding population.
The dynamic properties of the solutions of Equation (1) are well understood. Under the condition λ T π 2 , the equilibrium state u 0 1 is asymptotically stable, while for λ T > π 2 it is unstable and the periodic solution U 0 ( t , λ ) is stable. Under the condition 0 < λ T π 2 1 , we obtain its asymptotics on the Andronov–Hopf bifurcation theory application basis (see, for example [16,17,18]). It is shown in [7] that for sufficiently large λ the periodic solution U 0 ( t , λ ) has a pronounced relaxation structure. There is a single wavelet of this function on the segment of the period T 0 ( λ ) length. The duration of the wavelet is close to the value of T and the period T 0 ( λ ) is equal to ( λ T ) 1 exp ( λ T ( 1 + o ( 1 ) ) ) . It is important to note the asymptotic equalities.
max t U 0 ( t , λ ) = exp ( λ T ( 1 + o ( 1 ) ) ) ,
min t U 0 ( t , λ ) = exp ( exp ( λ T ) ( 1 + o ( 1 ) ) ) .
We also note results of [19,20], in which important results on the dynamics of logistic equations with delay and diffusion are obtained, [21], where a randomized non-autonomous logistic equation is discussed, and Ref. [22], where interesting results were obtained for an equation in which the delay itself is a function of u.
Due to the applied significance, it is of interest to study equations with various types of nonlinearity. By virtue of biological meaning, solutions must remain positive. Similar to this restriction, it is natural to require solutions to be bounded not only from below, but also from above. In addition, unlimited asymptotic growth, as in the case of Equation (1), is not always convenient, for example, for describing the dynamics of changes in the number of species.
In this paper we study the dynamics of the new model—the logistic equation with delay and modified nonlinearity
u ˙ = λ [ 1 u ( t T ) ] u ( A u ) , A > 1 .
The main difference between the solutions of Equation (4) and the solutions of Equation (1) is that under the condition 0 φ ( s ) A the inequalities
0 u ( t , φ ) A
hold for the solution u ( t , φ ) of Equation (4) with the initial condition u ( t 0 + s ) = φ ( s ) as t > t 0 . Thus, the use of the factor ( A u ) on the right-hand side of Equation (4) leads to the restriction u ( t , φ ) A for the solutions of Equation (4). Therefore, the solutions of Equation (4) (in contrast to the solutions of Equation (1)) do not take asymptotically large values for large values of λ .
Model (4) is quite complex. At the moment, there are no analytical methods that would allow for absolutely any values of the parameters λ , T, and A to obtain results on the qualitative and quantitative behavior of all solutions of Equation (4) on the entire positive semiaxis t [ 0 , + ) . In addition, numerical methods are also not very effective in the study of such models, since, firstly, it is impossible to enumerate all possible initial conditions in order to draw a qualitative conclusion about the dynamics of the model, and secondly, for some values of the model parameters, the solutions are very close will approach zero (see Section 3 and Section 4), so a relatively small error in the calculations can lead to a completely wrong result.
Therefore, in the paper, the research is carried out by analytical (asymptotic) methods. At present, the most effective among them are the methods of the theory of bifurcations in the study in the neighbourhood of the equilibrium state and the methods of a large parameter (for λ 1 ) in the study of nonlocal processes. Large parameter methods are based on the use of singular perturbation theory, but even in cases where the λ parameter (Malthusian coefficient in model Equation (4)) is not large enough, nevertheless, conclusions can be drawn about the trends in fluctuations with an increase in this parameter. In addition, we point out that the relaxation nature of the oscillations obtained for model Equation (4) in the case of sufficiently large values of parameter λ corresponds to ideas about the dynamic behavior of population density dynamics in mathematical ecology [4,11].
The paper is organized as follows. In Section 2, the conditions for the parameter λ are found using the methods of bifurcation analysis. Under them, the Andronov–Hopf bifurcation occurs in Equation (4) and a cycle arises from the equilibrium state.
The main content of Section 3 and all subsequent sections is devoted to the study of the solutions of Equation (4) on the interval t [ 0 , + ) provided that the parameter λ is sufficiently large:
λ 1 .
In this case, the equation under consideration is singularly perturbed. After dividing by the left-hand and right-hand parts of Equation (4) by λ , we obtain an equation with a small parameter at the derivative:
λ 1 u ˙ = [ 1 u ( t T ) ] u ( A u ) .
However, the reduced equation
0 = [ 1 u ( t T ) ] u ( A u )
does not describe the behavior of the Equation (6) solutions (except the simplest equilibrium states) under the condition (5).
It is important to mention that from a computational point of view, this problem is difficult, since the relaxation solution approaches A and 0 very closely, so even a small error in the calculations will take us out of the class of solutions under consideration.
Therefore, a special analytical research method was developed. First, we briefly describe it for the simplest case, for studying slowly oscillating solutions. Recall that slowly oscillating solutions are those solutions for which the time distance between adjacent roots of the equation u ( t ) = 1 is greater than the delay time. Consider Equation (6). Denote by S C [ T , 0 ] the set of all such functions φ ( s ) S ( s [ T , 0 ] ) that satisfy the conditions (see Figure 1).
φ ( 0 ) = 1 ; 0 φ ( s ) exp λ 2 s .
Denote by u ( t , φ ) a solution of (6) with initial condition φ ( s ) : u ( s , φ ) = φ ( s ) ( s [ T , 0 ] ) . The method of steps will be used: for t [ 0 , T ] , taking into account the initial function φ ( s ) , to determine u ( t , φ ) we obtain an ordinary differential equation of the first order. It is easy to find the asymptotics of this solution as λ on the indicated interval. Then we consider u ( t , φ ) on the interval [ T , 2 T ] , on which we also obtain a first-order ordinary differential equation, and find the asymptotics of the solutions. After that, we will carry out the same actions on the segment [ 2 T , 3 T ] , then on the segment [ 3 T , 4 T ] , and so on.
In Section 3.1 we study the asymptotic behavior at λ of all solutions u ( t , φ ) with initial conditions from S. In particular, the asymptotics of the first two positive roots t 1 ( φ ) and t 2 ( φ ) of the equation u ( t , φ ) = 1 will be found. The main conclusion is that after the time interval t 2 ( φ ) the solution u ( t , φ ) will again fall on the set S: u ( s + t 2 ( φ ) , φ ) S . This gives reason to introduce the Π operator:
Π ( φ ( s ) ) = u ( t 2 ( φ ) + s , φ )
and justify the inclusion Π S S . From this and the well-known results of functional analysis [23], we conclude that S has a fixed point φ 0 ( s ) of the operator Π : Π ( φ 0 ( s ) ) = φ 0 ( s ) . Then the function u 0 ( t ) = u ( t , φ 0 ) is periodic with period t 2 ( φ 0 ) . In Section 3.1, asymptotic formulas for u 0 ( t ) will be given.
In Section 3.2 we investigate the asymptotic behavior of rapidly oscillating solutions. As a set of initial functions, we consider the set S ( τ 1 , τ 2 ) depending on two parameters τ 1 and τ 2 . It is shown schematically in Figure 3. In Section 3.2, a rigorous description of this set is given and the asymptotic behavior of all solutions (6) with initial conditions from S ( τ 1 , τ 2 ) is studied. Again, by t 1 ( φ ) and t 2 ( φ ) we denote the first and second positive roots of the equation u ( t , φ ) = 1 . We again introduce the operator Π : Π ( φ ( s ) ) = u ( t 2 ( φ ) + s , φ ) and it will be shown that up to o ( 1 ) the inclusion Π ( φ ( s ) ) S ( τ ¯ 1 , τ ¯ 2 ) , where τ ¯ 1 , 2 are explicitly expressed in terms of τ 1 , 2 . This gives grounds to take the function u ( s + t 2 ( φ ) , φ ) as the initial one and pass to the iteration process. It is possible to find a fixed point ( τ 1 0 , τ 2 0 ) and a function φ 0 ( s ) S ( τ 1 0 , τ 2 0 ) such that the solution u 0 ( t , φ 0 ) will be periodic with period t 2 ( φ 0 ) . It is important to note that the rapidly oscillating periodic solutions found in this way are unstable.
The results of Section 4 are more interesting. We consider equations with two delays
u ˙ = λ [ 1 α u ( t T ) ( 1 α ) u ( t h ) ] u ( A u ) ,
where α ( 0 , 1 ) and h < T .
It is commonly supposed that the use of two (or more) delays allows one to take into account the influence of the population age structure on the dynamics of population level changes [24,25,26,27].
One succeeded to clearly describe the nonlocal dynamic properties of Equation (8) under the condition (5). In contrast to the results of the Section 3, the relaxation step solutions for (8) can be stable.
We note that for the simpler logistic equation with two delays
u ˙ = λ [ 1 α u ( t T ) ( 1 α ) u ( t h ) ] u α ( 0 , 1 ) , h < T
the results were obtained in Ref. [7] as λ 1 . It follows from them that for sufficiently large λ in (9) there is an orbitally stable periodic solution which is similar in appearance to the solution of Equation (1) with the same initial conditions, and the asymptotic representation (2) holds for it.

2. Andronov–Hopf Bifurcations in Equations with One Delay

We shall fix the parameters λ 0 , A 0 and T 0 so that the following equality is fulfilled:
λ 0 ( A 0 1 ) T 0 = π 2 .
Lemma 1. 
Suppose that in Equation (4)
λ ( A 1 ) T < π 2 .
Then the equilibrium state u 0 1 is asymptotically stable.
Proof. 
To prove this statement we linearize (4) on u 0 and consider the characteristic quasi-polynomial
μ = λ ( A 1 ) exp ( μ T ) .
Under the condition (11) all its roots have negative real parts. From this, and from known results (see Ref. [5]) on the stability in the first approximation, the proof of Lemma 1 follows. □
Next we introduce a small positive parameter ε :
0 < ε 1 ,
and assume that in Equation (4), for certain constant λ 1 , A 1 and T 1 , the following relations are fulfilled:
λ = λ 0 + ε λ 1 , A = A 0 + ε A 1 , T = T 0 + ε T 1 .
We set
b = 1 + π 2 4 1 π 2 + i ( λ 1 ( A 0 1 ) + λ 0 A 1 ) + λ 0 2 ( A 0 1 ) 2 T 1 1 i π 2 ,
d = λ 0 1 + π 2 4 1 π 2 + i + 3 ( A 0 2 ) 2 5 ( A 0 1 ) π 2 1 + i π 2 + 1 ,
ω = π ( 2 T 0 ) 1 ,
ξ 0 ( τ ) = ξ 0 exp ( i ψ τ ) , ψ = I m   b + ξ 0 2 I m   d , ξ 0 = π 2 ( λ 1 ( A 0 1 ) + λ 0 A 1 ) + λ 0 2 ( A 0 1 ) 2 T 1 λ 0 1 π 2 + 3 ( A 0 2 ) 2 5 ( A 0 1 ) π 2 1 1 1 2 .
Theorem 1. 
Let conditions (13)–(15), (17), and (18) hold.
1. Let R e   b < 0 . Then all solutions of Equation (4) with initial conditions from a sufficiently small (and independent of ε) neighbourhood of the equilibrium state u 1 tend to 1 as t .
2. Let
R e   b > 0 .
Then for all sufficiently small ε Equation (4) has in the neighbourhood of the unit equilibrium state a stable periodic solution u 0 ( t , ε ) for which the following asymptotic equality is fulfilled:
u 0 ( t , ε ) = 1 + ε 1 / 2 [ ξ 0 ( τ ) exp ( i ω t ) + ξ ¯ 0 ( τ ) exp ( i ω t ) ] + O ( ε ) ,
where τ = ε t .
Proof. 
Under the condition (10) the quasi-polynomial (12) has a pair of purely imaginary roots μ 1 , 2 = ± i ω , where the value of ω is given in (17), and all its other roots have negative real parts. Then, under the condition (14), and for all sufficiently small ε , Equation (4) in the neighbourhood of u 1 has (see Refs. [16,17,18,28]) a two-dimensional stable local invariant integral manifold, on which this equation can be written, to within higher-order terms, in the form of a complex scalar equation—a normal form:
d ξ d τ = B ξ ( τ ) + D ξ ( τ ) | ξ ( τ ) | 2 ,
where values B and D are to be determined.
In the case R e   B 0 and R e   D 0 normal form completely determines the behavior of all solutions in a neighbourhood of u 1 (see Refs. [16,17,18,28]). The solutions of (20) are connected with the solutions of (4) by the relation
u ( t , ε ) = 1 + ε 1 / 2 [ ξ ( τ ) exp ( i ω t ) + ξ ¯ ( τ ) exp ( i ω t ) ] + + ε u 2 ( t , τ ) + ε 3 / 2 u 3 ( t , τ ) +
Here, τ = ε t , the u j ( t , τ ) are periodic in t with period 2 π / ω . Substituting (21) into (4) and collecting coefficients of equal powers of ε , we successively find all the elements appearing there. At ε 1 / 2 we obtain true identity, and collecting coefficients of the first power of ε and taking into account (10), (12), and (17), we find equation for determining the function u 2 ( t , τ ) :
u 2 t = λ 0 ( A 0 1 ) u 2 ( t T 0 , τ ) + i λ 0 ( A 0 2 ) ξ 2 ( τ ) exp ( 2 i ω t ) i λ 0 ( A 0 2 ) ξ ¯ 2 ( τ ) exp ( 2 i ω t ) .
Thus, taking into account that u 2 ( t , τ ) is periodic in t with period 2 π / ω , we find that
u 2 ( t , τ ) = ( A 0 2 ) ( 2 i ) 5 ( A 0 1 ) ξ 2 ( τ ) exp ( 2 i ω t ) + ( A 0 2 ) ( 2 + i ) 5 ( A 0 1 ) ξ ¯ 2 ( τ ) exp ( 2 i ω t ) .
Collecting coefficients of ε 3 / 2 we obtain an equation for u 3 ( t , τ ) :
u 3 t = λ 0 ( A 0 1 ) u 3 ( t T 0 , τ ) + B 1 exp ( i ω t ) + B ¯ 1 exp ( i ω t ) + B 3 exp ( 3 i ω t ) + B ¯ 3 exp ( 3 i ω t ) ,
where B 1 and B 3 are complex values.
From the condition that Equation (22) be solvable in the class of functions periodic in t with period 2 π / ω (this condition is B 1 = 0 ) we arrive at Equation (20) for the unknown amplitude ξ ( τ ) , in which B = b and D = d (coefficient b is defined in (15) and the coefficient d is defined in (16)).
Multiplying both sides of Equation (20) by ξ ¯ ( τ ) and taking into account equalities B = b and D = d we obtain a scalar real ordinary differential equation
1 2 ρ ˙ ( τ ) = R e   b ρ ( τ ) + R e   d ρ 2 ( τ ) ,
where ρ ( τ ) = | ξ ( τ ) | 2 .
We note that R e   d < 0 .
That is why, if R e   b < 0 , all solutions of (23) tend to zero as t + . Thus, all solutions of (4) from the neighbourhood of u 1 tend to 1 as t + . This completes the proof of the first part of the Theorem.
Under the condition (19) Equation (20) with B = b and D = d has a stable periodic solution ξ 0 ( τ ) = ξ 0 exp ( i ψ τ ) , where the values ξ 0 and ψ are given in (18). Taking into account here the asymptotic equality (21) we have completed the proof of the Theorem. □
It is interesting to note that the bifurcation effect considered here can be realized for fixed λ = λ 0 and T = T 0 , i.e., for λ 1 = T 1 = 0 , and only upon variation of the parameter A.

3. Step-like Solutions of the Logistic Equation with Delay and with a Restriction on the Nonlinear Function

All statements below refer to the case of λ + . For example, everywhere below o ( 1 ) means that the function tends to zero as λ + . Note that the phrase “for all sufficiently large λ ” means the existence of a value λ 0 such that for all λ > λ 0 the corresponding assertion is true.
We separately consider questions about slowly and rapidly oscillating solutions of Equation (4).

3.1. Asymptotic Behavior of Slowly Oscillating Solutions

We shall consider the asymptotic behavior of a slowly oscillating relaxation periodic solution of Equation (4).
In (4) it is convenient to replace the time t t T and to again denote the product λ T by λ . Then Equation (4) takes the form
u ˙ = λ [ 1 u ( t 1 ) ] u ( A u ) .
Let us introduce some notation.
By S we denote the set of all functions φ ( s ) C [ 1 , 0 ] that satisfy the conditions (7).
Let u ( t , φ ) be the solution of (24) with an initial condition φ ( s ) specified at t = 0 , i.e., u ( s , φ ) = φ ( s ) for s [ 1 , 0 ] . We shall construct the asymptotics for λ of all solutions u ( t , φ ) with φ ( s ) S .
Theorem 2. 
For all sufficiently large λ Equation (24) has an asymptotically orbitally stable periodic (with period T ( λ ) ) solution u 0 ( t , λ ) , for which u 0 ( 0 , λ ) = u 0 ( t 1 ( λ ) , λ ) = u 0 ( t 2 ( λ ) , λ ) = 1 , where t j ( λ ) are successive positive roots of the equation u 0 ( t , λ ) = 1 , and
t 1 ( λ ) = 1 + ( A 1 ) 1 + o ( 1 ) , T ( λ ) = t 2 ( λ ) = A + 1 + ( A 1 ) 1 + o ( 1 ) .
For every t from the intervals ( 0 , t 1 ( λ ) ) and ( t 1 ( λ ) , t 2 ( λ ) ) , respectively, the following equalities are fulfilled:
u 0 ( t , λ ) = A + o ( 1 ) , t ( 0 , t 1 ( λ ) ) , o ( 1 ) , t ( t 1 ( λ ) , t 2 ( λ ) ) .
Proof. 
We use the following formula for the solutions u ( t , φ ) of Equation (24):
u ( t , φ ) = A u ( τ , φ ) exp λ A τ t 1 u ( s 1 , φ ) d s · A u ( τ , φ ) + u ( τ , φ ) exp λ A τ t 1 u ( s 1 , φ ) d s 1 ,
where t and τ are arbitrary values such that 0 < τ < t .
We shall formulate some simple statements:
Lemma 2. 
Let t ( 0 , 1 ] . Then
u ( t , φ ) = A 1 + ( A 1 ) exp λ A 1 + o ( 1 ) t 1 .
The proof of (26) follows from (25) with τ = 0 .
Setting τ = 1 in (25), we immediately arrive at the following statement:
Lemma 3. 
Let t ( 1 , 1 + ( A 1 ) 1 + δ ] , where δ > 0 is some small but fixed value. Then
u ( t , φ ) = A 1 + ( A 1 ) exp λ A ( A 1 ) 1 + o ( 1 ) ( t 1 ) 1 1 .
Below, we denote by t 1 ( φ ) , t 2 ( φ ) , … the first, second, … positive roots of the equation u ( t , φ ) = 1 .
From (27) we arrive at the conclusion that the value t 1 ( φ ) exists and
t 1 ( φ ) = 1 + ( A 1 ) 1 + o ( 1 ) , u ( t , φ ) = A + o ( 1 ) for every t ( 0 , t 1 ( φ ) ) , u ( t , φ ) = o ( 1 ) for every t ( t 1 ( φ ) , t 1 ( φ ) + 1 ) .
The last equality is true, because formula (27) holds on the interval ( t 1 ( φ ) , t 1 ( φ ) + 1 ) .
Then we construct asymptotics of solution on the interval t > t 1 ( φ ) + 1 . We set τ = t 1 ( φ ) and take into account that u ( s 1 ) = A + o ( 1 ) on the interval s ( t 1 ( φ ) , t 1 ( φ ) + 1 ) . We use equality u ( s 1 ) = o ( 1 ) for t 1 ( φ ) + 1 < s < t . It follows from formula (25) that this equality holds for all t > t 1 ( φ ) + 1 while
t 1 ( φ ) t 1 u ( s 1 , φ ) d s < 0 .
We obtain the following result:
Lemma 4. 
For t [ 2 + ( A 1 ) 1 , A + 1 + ( A 1 ) 1 + δ ] , where δ > 0 is some small but fixed value
u ( t , φ ) = A exp λ A 1 A + t 2 + ( A 1 ) 1 + o ( 1 ) · A 1 + exp λ A 1 A + t 2 + ( A 1 ) 1 + o ( 1 ) 1 .
From (28) the existence of t 2 ( φ ) and the asymptotic equalities
t 2 ( φ ) = A + 1 + ( A 1 ) 1 + o ( 1 )
and
u ( t , φ ) = o ( 1 ) for t ( t 1 ( φ ) , t 2 ( φ ) )
follow.
We introduce the translation along the trajectories operator Π , which links the function u ( s + t 2 ( φ ) , φ ) , ( s [ 1 , 0 ] ) to a function φ ( s ) . From the formulas given above we arrive at the conclusion that
Π φ ( s ) = u s + t 2 ( φ ) , φ S ,
meaning that Π S S .
From this and from Ref. [23] it follows that this operator has in S a fixed point φ 0 ( s ) , i.e.,
Π φ 0 ( s ) = φ 0 ( s ) .
This means that the solution u 0 ( t , λ ) = u ( t , φ 0 ( s ) ) is periodic, with period T ( λ ) = t 2 ( φ 0 ( s ) ) .
The proof of the stability of u 0 ( t , λ ) is rather cumbersome. In the more complicated situation discussed in Refs. [7,29] a detailed proof was given of the stability of the solution constructed there, and so we shall not give the proof here. □
It is possible to obtain the asymptotic expansion of u 0 ( t , λ ) to any degree of precision. Here we note only that
max t u 0 ( t , λ ) = u 0 ( 1 , λ ) = A [ 1 + A 1 + o ( 1 ) exp ( λ A ) ) ] 1 , min t u 0 ( t , λ ) = u 0 t 1 ( φ 0 ) + 1 , λ = A ( A 1 ) 1 exp λ A ( 1 A ) 1 + o ( 1 ) .
It is interesting to compare the principal characteristics of u 0 ( t , λ ) with the corresponding characteristics of the stable periodic solution U 0 ( t , λ ) of Equation (1) for λ . To obtain formulas for max t U 0 ( t , λ ) and min t U 0 ( t , λ ) one should replace T by 1 in formulas (2) and (3) (see introduction of this paper and Ref. [7]). In Figure 2 we give graphs of periodic functions u 0 ( t , λ ) and U 0 ( t , λ ) .

3.2. Rapidly Oscillating Solutions of Equation (24)

We shall consider rapidly oscillating solutions of Equation (24).
These are the set of those solutions u ( t , φ ) that become equal to 1 at t = 0 , and have one “step” on the segment [ 1 , 0 ] , as shown in Figure 3.
We shall describe the corresponding set of initial functions φ ( s ) . We arbitrarily fix values τ 1 and τ 2 for which
1 < τ 1 < τ 2 < 0 .
Below, we shall denote by δ > 0 an arbitrary sufficiently small quantity that does not depend on λ and whose precise value is unimportant. We introduce into the analysis three functions κ 1 , 2 , 3 ( s , λ ) :
κ 1 ( s , λ ) = exp λ δ ( s τ 1 ) , s [ 1 , τ 1 ] ; κ 2 ( s , λ ) = min A 1 + ( A 1 ) exp λ δ ( s τ 1 ) 1 , A 1 + ( A 1 ) exp ( λ δ ( s τ 2 ) 1 ) , s ( τ 1 , τ 2 ) ; κ 3 ( s , λ ) = max exp λ δ ( s τ 2 ) , exp ( λ δ s ) , s [ τ 2 , 0 ] .
Finally, by S ( τ 1 , τ 2 ) we denote the set of functions φ ( s ) from C [ 1 , 0 ] that satisfy the conditions
φ ( τ 1 ) = φ ( τ 2 ) = φ ( 0 ) = 1 ; 0 φ ( s ) κ 1 ( s , λ ) , for s [ 1 , τ 1 ] ; κ 2 ( s , λ ) φ ( s , λ ) A , for s ( τ 1 , τ 2 ) ; 0 φ ( s ) κ 3 ( s , λ ) for s [ τ 2 , 0 ] .
The first and the second positive roots of the equation u ( t , φ ) = 1 will be denoted by t 1 ( φ ) and t 2 ( φ ) . We set t 0 = A ( A 1 ) 1 ( τ 1 + 1 ) . Asymptotic analysis of u ( t , φ ) for t [ 0 , τ 2 + 1 ] leads to the following statement:
Lemma 5. 
Suppose that the condition
t 0 < τ 2 + 1
is fulfilled. Then for t 1 ( φ ) the following asymptotic (for λ ) equality holds:
t 1 ( φ ) = t 0 + o ( 1 ) .
Proof. 
As in Theorem 2, when constructing the asymptotics of the solution to Equation (24), we use the formula (25). It follows from (25), that on the segment t [ 0 , τ 1 + 1 ] solution has form (26). It follows from (26), that on the interval t ( 0 , τ 1 + 1 ] equality u ( t , φ ) = A + o ( 1 ) holds.
On the segment t [ τ 1 + 1 , τ 2 + 1 ] the function u ( t 1 , φ ) = A + o ( 1 ) . That is why on this time segment the solution to Equation (24) satisfies the formula
u ( t , φ ) = A 1 + ( A 1 ) exp λ A ( τ 1 + 1 ) + ( A 1 ) ( t ( τ 1 + 1 ) ) + o ( 1 ) 1 .
While τ 1 + 1 < t < t 0 , the right-hand side in formula (30) is asymptotically close to A at λ + , and for t > t 0 this expression is o ( 1 ) at λ + . Therefore, if t 0 < τ 2 + 1 , then the value t 0 is asymptotically close to t 1 ( φ ) at λ + . □
We note that for every t ( 0 , t 1 ( φ ) ) the relation
u ( t , φ ) = A + o ( 1 )
is fulfilled. Next, for t ( t 1 ( φ ) , τ 2 + 1 ) , we have u ( t , φ ) = o ( 1 ) , with
u ( τ 2 + 1 , φ ) = A ( A 1 ) 1 exp ( λ A ( A 1 ) ( τ 2 + 1 t 0 + o ( 1 ) ) ) .
We set t 0 = A ( τ 2 + 1 ) + ( 1 A ) t 0 . If
t 0 < 1 ,
then for t ( τ 2 + 1 , t 0 + δ 1 ) (where δ 1 is some sufficiently small quantity that does not depend on λ ) we find that
u ( t , φ ) = A / ( A 1 ) exp λ A t A ( τ 2 + 1 ) ( 1 A ) t 0 + o ( 1 ) + 1 .
From this follow both the existence of t 2 ( φ ) and the asymptotic formula t 2 ( φ ) = t 0 + o ( 1 ) .
It is easy to show, that if inequality t 0 > 1 holds, then t 2 ( φ ) exists, but t 2 ( φ ) t 1 ( φ ) > 1 . So, in the case t 0 > 1 we return to the case of the slowly oscillating solutions.
If inequality (31) holds, then as in the preceding section, we introduce into the analysis the operator Π , using the rule Π ( φ ( s ) ) = u ( s + t 2 ( φ ) , φ ) .
Note, that t 2 ( φ ) corresponds to 0, t 1 ( φ ) corresponds to τ 2 , 0 corresponds to τ 1 (see Figure 3). That is why τ ¯ 1 = 0 t 0 + o ( 1 ) τ ¯ 2 = t 0 t 0 + o ( 1 ) .
A consequence of the formulas given here is the following statement:
Theorem 3. 
Let the conditions (29) and (31) are fulfilled. Then for all sufficiently large λ for the Equation (24) the following inclusion holds:
Π ( φ ( s ) ) S ( τ ¯ 1 + o ( 1 ) , τ ¯ 2 + o ( 1 ) ) ,
with
τ ¯ 1 = A τ 1 A τ 2 , τ ¯ 2 = A 2 ( A 1 ) 1 τ 1 A τ 2 + A ( A 1 ) 1 .
The system of Equation (32) is a linear inhomogeneous difference system of equations. After τ ¯ 1 and τ ¯ 2 have been determined, the situation is repeated, i.e., while for current values of τ ¯ 1 and τ ¯ 2 inequalities (29) and (31) are fulfilled, using (32) we calculate τ = 1 and τ = 2 , and so on. This iteration process has a fixed point. By virtue of the fact that the determinant of the linear part in (32) is greater in modulus than 1 the iteration process is divergent, i.e., the equilibrium state in (32) is unstable.
We also remark that the number of periodic solutions of Equation (24) grows without limit as λ . We shall show this. It was established above that for sufficiently large λ Equation (24) has a slowly oscillating periodic solution u 0 ( t , λ ) with period T ( λ ) . This solution for each integer n satisfies the equation
u ˙ = λ [ 1 u ( t 1 n T ( λ ) ) ] u ( A u ) .
In (33) we replace the time t ( 1 + n T ( λ ) ) t . We then find that the function u n ( t , λ ) = u 0 ( ( 1 + n T ( λ ) ) t , λ ) is a solution of the equation
u ˙ = λ ( 1 + n T ( λ ) ) [ 1 u ( t 1 ) ] u ( A u ) .
The period of the function u n ( t , λ ) is equal to T ( λ ) [ 1 + n T ( λ ) ] 1 . In particular, the unstable periodic solution determinable from the fixed points τ 10 , τ 20 of the mapping (32) corresponds to the value n = 1 in (34), i.e., to the function u 1 ( t , λ ) .
We shall consider the more general case when on the segment [ 1 , 0 ] the initial functions φ ( s ) take the value 1 exactly ( 2 n + 1 ) times:
1 < τ 1 < τ 2 < < τ 2 n < 0 , φ ( τ j ) = φ ( 0 ) = 1 ( j = 1 , , 2 n ) .
Following the method proposed above, we arrive at the 2 n -dimensional map
τ ¯ 1 = τ 3 t 0 , t 0 = A τ 2 A τ 1 , τ ¯ 2 = τ 4 t 0 , t 0 = A ( A 1 ) 1 ( τ 1 + 1 ) , τ ¯ 2 n 1 = t 0 , τ ¯ 2 n = t 0 t 0 .
While we are in the class of initial conditions (conditions (35) are satisfied), we iterate this map. It can be shown that if the equilibrium state in (36) exists, then it is unstable. Since solution u n ( t , λ ( 1 + n T ( λ ) ) ) corresponds to equilibrium state in map (36), then this periodic solution of (24) is unstable.

4. Step-like Solutions of the Equation with Two Delays

We normalize the time t T t and rename λ T as λ and h T 1 as h ( h < 1 ) again in Equation (8) with two delays. As a result, we arrive at the equation
u ˙ = λ [ 1 α u ( t 1 ) ( 1 α ) u ( t h ) ] u ( A u ) .

4.1. Construction of Slowly Oscillating Solutions

We shall consider the construction of slowly oscillating periodic solutions. We introduce into the analysis a set of initial conditions consisting of all those functions φ ( s ) C [ 1 , 0 ] for which
φ ( 0 ) = 1 , 0 φ ( s ) 1 ; φ ( s ) exp ( λ δ s ) ,
where δ > 0 is some sufficiently small but fixed value. We investigate the asymptotic form of all solutions u ( t , φ ) with initial conditions from S. We denote by t 1 ( φ ) and t 2 ( φ ) the first and second positive roots of equation u ( t , φ ) = 1 . In the case when
t 2 ( φ ) 1 > t 1 ( φ ) ,
we can determine the translation along the trajectories operator Π : Π ( φ ( s ) ) = u ( s + t 2 ( φ ) , φ ) . From the inequality (38) we find that Π ( φ ( s ) ) S . If for all φ ( s ) S we have Π ( φ ( s ) ) S , i.e., Π S S , we can conclude that there exists an attractor with initial conditions from S and that there exists a fixed point φ 0 ( s ) S . The solution u 0 ( t , λ ) = u ( t , φ 0 ( s ) ) is thereby periodic, with period T ( λ ) = t 2 ( φ 0 ) . For the expressions t 1 , 2 ( φ ) it is possible to obtain asymptotic formulas, and this means that we can formulate conditions on the parameters of the problem that ensure fulfilment of the inequality (38). We confine ourselves here to giving only sufficient conditions for fulfilment of this inequality.
Lemma 6. 
Let
1 ( 1 α ) A > 0
and
1 α A < 0 .
Then for all sufficiently large λ for Equation (37) we have the following inclusion:
Π S S .
Proof. 
For each t ( 0 , h ] the asymptotic equality
u ( t , φ ) = A 1 + ( A 1 ) exp λ A t + o ( 1 ) 1
holds. But if t [ h , 1 ] , then
u ( t , φ ) = A 1 + ( A 1 ) exp λ A ( 1 ( 1 α ) A ) t + h ( 1 α ) A + o ( 1 ) 1 .
In particular,
u ( 1 , φ ) = A 1 + ( A 1 ) exp λ A 1 ( 1 α ) A ( 1 h ) + o ( 1 ) 1 .
It follows from inequality (39) that u ( t , φ ) = A + o ( 1 ) on the segment t [ h , 1 + δ 1 ] , where δ 1 is some positive and small (but independent on λ ) constant.
Therefore, for t > 1 while u ( t , φ ) = A + o ( 1 ) formula
u ( t , φ ) = A 1 + ( A 1 ) exp λ A 1 ( 1 α ) A ( 1 h ) + ( 1 A ) ( t 1 ) + o ( 1 ) 1
holds. From this formula we find that t 1 ( φ ) exists and that
t 1 ( φ ) = 1 + 1 ( 1 α ) A ( 1 h ) ( A 1 ) 1 + o ( 1 ) , ( t 1 ( φ ) > 1 ) .
Starting from τ = t 1 ( φ ) we obtain that on the interval t ( t 1 ( φ ) , t 1 ( φ ) + h ] equality u ( t , φ ) = o ( 1 ) holds and that on the segment t [ t 1 ( φ ) + h , t 1 ( φ ) + 1 ] formula
u ( t , φ ) = A 1 + ( A 1 ) exp λ A ( 1 A ) h + ( 1 α A ) ( t h t 1 ( φ ) ) + o ( 1 ) 1
is true. It follows from (40) that u ( t , φ ) = o ( 1 ) on the segment t [ t 1 ( φ ) + h , t 1 ( φ ) + 1 ] . Then for t > t 1 ( φ ) + 1 while u ( t , φ ) = o ( 1 ) we have the equality
u ( t , φ ) = A 1 + ( A 1 ) exp λ A ( 1 A ) h + ( 1 α A ) ( 1 h ) + ( t 1 t 1 ( φ ) ) + o ( 1 ) 1 .
From Formula (41) we obtain both the existence of t 2 ( φ ) and inequality t 2 ( φ ) t 1 ( φ ) > 1 .
It follows from this that Π ( φ ( s ) ) S and Π S S . The Lemma is proven. □
Theorem 4. 
Let inequalities (39) and (40) hold. Then for all sufficiently large λ Equation (37) has an asymptotically orbitally stable slowly oscillating periodic (with period T ( λ ) ) solution u 0 ( t , λ ) , for which u 0 ( 0 , λ ) = u 0 ( t 1 ( λ ) , λ ) = u 0 ( t 2 ( λ ) , λ ) = 1 , where t j ( λ ) are successive positive roots of the equation u 0 ( t , λ ) = 1 , and
t 1 ( λ ) = 1 + 1 ( 1 α ) A ( 1 h ) ( A 1 ) 1 + o ( 1 ) , T ( λ ) = t 2 ( λ ) = 2 + 1 ( 1 α ) A ( 1 h ) ( A 1 ) 1 + ( A 1 ) h ( 1 α A ) ( 1 h ) + o ( 1 ) .
For every t from the intervals ( 0 , t 1 ( λ ) ) and ( t 1 ( λ ) , t 2 ( λ ) ) , respectively, the following equalities are fulfilled:
u 0 ( t , λ ) = A + o ( 1 ) , t ( 0 , t 1 ( λ ) ) , o ( 1 ) , t ( t 1 ( λ ) , t 2 ( λ ) ) .
The proof of existence of the periodic solution in Theorem 4 follows from the proof of the Lemma 6, and proof of stability of this solution is rather cumbersome. In the more complicated situation discussed in Refs. [7,29] a detailed proof was given of the stability of the solution constructed there, and so we shall not give the proof here.

4.2. Construction of More Complicated Solutions

We shall consider the construction of rapidly oscillating periodic solutions. We arbitrarily fix τ , satisfying condition
τ max h , 1 h h ( 1 α A ) ( A 1 ) 1 , 1
and consider the set
S ( τ ) = { φ ( s ) C [ 1 , 0 ] : 0 φ ( s ) A ; φ ( τ ) = φ ( 0 ) = 1 , φ ( s ) max exp λ δ ( s + τ ) , exp ( λ δ s ) , s [ τ , 0 ] ; φ ( s ) min A 1 + ( A 1 ) exp λ δ ( s + 1 ) 1 , A 1 + ( A 1 ) exp ( λ δ ( s + τ ) 1 ) , s [ 1 , τ ] } .
We shall give the final asymptotic formulas for u ( t , φ ) ( φ S ( τ ) ) . In them we denote by g ( x ) the function
g ( x ) = x , if x > 0 ; 0 , if x 0 .
We shall assume that the inequalities
1 α A > 0 , 1 ( 1 α ) A < 0 ,
and
h < 1 h
are fulfilled. Let t ( 0 , h ] . Then u ( t , φ ) = A + o ( 1 ) , and
u ( h , φ ) = A 1 + ( A 1 ) exp λ A h α A min ( h , 1 τ ) + o ( 1 ) 1 .
It follows from condition (42), that for t [ h , max ( h , 1 τ ) ] we also have u ( t , φ ) = A + o ( 1 ) , and
u ( max ( h , 1 τ ) , φ ) = A [ 1 + ( A 1 ) exp ( λ A [ max ( h , 1 τ ) α A min ( h , 1 τ ) A g ( 1 τ h ) + o ( 1 ) ] ) ] 1 .
We set
t 1 0 = α A ( 1 τ + h ) A h 1 ( 1 α ) A .
From (42) and (43) we get that t 1 0 > max ( h , 1 τ ) . Then t 1 ( φ ) = t 1 0 + o ( 1 ) , and for every t ( 0 , t 1 ( φ ) ) we have the relation u ( t , φ ) = A + o ( 1 ) .
Suppose that the following inequality is fulfilled:
t 1 0 < 1 h .
Then on the interval ( t 1 ( φ ) , t 1 ( φ ) + h ] we have the equality u ( t , φ ) = o ( 1 ) , and
u t 1 ( φ ) + h , φ = A 1 + ( A 1 ) exp λ A ( 1 ( 1 α ) A ) h + o ( 1 ) 1 .
If condition
( 1 ( 1 α ) A ) h + 1 t 1 0 h < 0
is true, then on the interval ( t 1 ( φ ) + h , 1 ] equality u ( t , φ ) = o ( 1 ) holds, and
u 1 , φ = A 1 + ( A 1 ) exp λ A ( 1 ( 1 α ) A ) h + ( 1 t 1 0 h ) + o ( 1 ) 1 .
We set
t 2 0 = 1 + ( 1 α ) A h + t 1 0 1 [ 1 α A ] 1 .
Suppose
t 2 0 < t 1 0 + 1 ,
then t 2 ( φ ) = t 2 0 + o ( 1 ) , and for every t [ 1 , t 2 ( φ ) ) the equality u ( t , φ ) = o ( 1 ) is fulfilled. Finally, we introduce the quantity τ ¯ = t 2 ( φ ) t 1 ( φ ) . Then, to order o ( 1 ) (for λ ), we obtain the equality
τ ¯ = f ( τ ) ,
where f ( τ ) = t 2 0 t 1 0 .
After τ ¯ has been determined, the situation is repeated, i.e., while for current value of τ ¯ inequalities (42), (45)–(47) are fulfilled, using (48) we calculate τ = and so on.
The one-dimensional mapping (48) plays a central role in the study of the solutions from S ( τ ) . When the conditions (42)–(47) are fulfilled its attractors determine the structure of the periodic solutions from S ( τ ) . For example, to the equilibrium state τ 0 :
τ 0 = f ( τ 0 )
there corresponds, under the condition | f ( τ 0 ) | < 1 , a stable one-step periodic solution u 0 ( t , λ ) with the initial conditions in S ( τ 0 ) .
Lemma 7. 
Let the inequality
1 A + A 2 α < 0
hold. Then there exist values of h such that the map (48) has a stable fixed point τ = τ * , where
τ * = A ( A 1 ) ( α + ( α 1 ) h ) 1 A + A 2 α ,
such that inequalities (42)–(47) are fulfilled.
Proof. 
First, let us show, that if inequality (50) holds, then inequalities (43) are true.
It is easy to see, that inequalities (43) are equivalent to inequality α < min ( 1 / A , 1 1 / A ) . Inequality (50) is equivalent to α < 1 / A 1 / A 2 . This value is less than min ( 1 / A , 1 1 / A ) (parameter A > 1 ), therefore, if (50) is true, then conditions (43) hold.
Second, let us substitute the values of t 1 0 and t 2 0 into the map (48):
τ ¯ = A 2 α 2 ( 1 + α A ) ( 1 ( 1 α ) A ) τ + A ( A 1 ) ( h ( 1 α ) α ) ( 1 + α A ) ( 1 ( 1 α ) A ) .
This map has a fixed point τ * satisfying formula (51). We should check whether this value τ * satisfies conditions (42), (45)–(47), because map (52) corresponds to the initial equation only when all these inequalities hold, and we should study stability of the fixed point (51).
Let us begin with the study of stability. The fixed point τ * is stable if and only if
1 < A 2 α 2 ( 1 + α A ) ( 1 ( 1 α ) A ) < 1 .
From the condition (43) we get that this value is greater than zero, and this value is less than 1 if and only if inequality (50) holds. Therefore, this fixed point is stable.
Now let us prove that if τ = τ * satisfies (51), then conditions (42), (45)–(47) hold.
Below, we always assume that inequalities (43) and (50) hold.
Substituting τ * in inequality τ < 1 from (42) we obtain, that this condition is equivalent to
h < ( 1 α ) A 1 A ( A 1 ) ( 1 α ) = : γ 1 ,
inequality h < τ from (42) is equivalent to
h > α A ( A 1 ) ( A 1 ) 2 + α A = : β 1 ,
and condition τ > 1 h h ( 1 α A ) ( A 1 ) 1 from (42) is equivalent to
h > A 1 A 2 ( 1 α ) = : β 2 .
Condition (45) at τ = τ * can be rewritten in the form
h < ( A 1 ) ( 1 α A ) ( A ( 1 α A ) ( 1 ( 1 α ) A ) ) = : γ 2 .
Condition (46) holds, if τ = τ * and (54) is true, and inequality (47) holds, if τ ¯ = τ * < 1 (the last inequality is equivalent to (53)).
Thus, at τ = τ * and under condition (50) the system of inequalities (42), (43), (45)–(47) is equivalent to
max ( β 1 , β 2 ) < h < min ( γ 1 , γ 2 ) .
Parameter h must satisfy condition h < 1 / 2 (we obtain this inequality from condition (44)).
Let us prove that under condition (50) interval ( max ( β 1 , β 2 ) , min ( 1 / 2 , γ 1 , γ 2 ) ) is not empty.
β 1 β 2 = ( A 1 ) ( 1 ( 1 α ) A ) ( 1 A + A 2 α ) A 2 ( 1 + α ) ( ( A 1 ) 2 + α A ) < 0 ,
therefore max ( β 1 , β 2 ) = β 2 . Let us estimate the values of γ 1 β 2 , γ 2 β 2 , 1 / 2 β 2 .
γ 1 β 2 = 1 A + A 2 α ( A 1 ) A 2 ( 1 + α ) > 0 , γ 2 β 2 = ( A 1 ) ( 1 A + A α ) ( 1 A + A 2 α ) A 2 ( 1 α ) A ( 1 α A ) ( 1 ( 1 α ) A ) > 0 .
Inequality
1 / 2 β 2 = 2 2 A + A 2 A 2 α 2 A 2 ( 1 α ) > 0
is equivalent to inequality 2 2 A + A 2 A 2 α > 0 . Function q ( α ) = 2 2 A + A 2 A 2 α is decreasing on the interval α ( 0 , 1 / A 1 / A 2 ) . Therefore, q ( 1 / A 1 / A 2 ) q ( α ) for all α ( 0 , 1 / A 1 / A 2 ) .
q ( 1 / A 1 / A 2 ) = 2 2 A + A 2 A 2 ( 1 / A 1 / A 2 ) = A 2 3 A + 3 > 0 ,
thus inequality 1 / 2 β 2 > 0 holds for all α ( 0 , 1 / A 1 / A 2 ) . The Lemma is proven. □
From the constructions and the reasoning given above we obtain the following statement.
Theorem 5. 
Let the inequality (50) hold. Then for all sufficiently large λ there exist h and τ = τ * , satisfying (51), such that equation (37) has an asymptotically orbitally stable rapidly oscillating periodic (with period T ( λ ) ) solution u * ( t , λ ) , for which u * ( 0 , λ ) = u * ( t 1 ( λ ) , λ ) = u * ( t 2 ( λ ) , λ ) = 1 , where t j are successive positive roots of the equation u * ( t , λ ) = 1 , and
t 1 ( λ ) = A ( α + ( α 1 ) h ) 1 A + A 2 α + o ( 1 ) , T ( λ ) = t 2 ( λ ) = A 2 ( α + ( α 1 ) h ) 1 A + A 2 α + o ( 1 ) , t 2 ( λ ) t 1 ( λ ) = τ * + o ( 1 ) < 1 .
For every t from intervals ( 0 , t 1 ( λ ) ) and ( t 1 ( λ ) , t 2 ( λ ) ) , respectively, the following equalities are fulfilled:
u * ( t , λ ) = A + o ( 1 ) , t ( 0 , t 1 ( λ ) ) , o ( 1 ) , t ( t 1 ( λ ) , t 2 ( λ ) ) .
In a similar way, one can determine the existence conditions and find the asymptotics of multi-step periodic solutions on the segment [ 1 , 0 ] and obtain an n-dimensional (in terms of the number of steps on the segment [ 1 , 0 ] ) mapping for their description.

5. Conclusions

The dynamics of the logistic equation with delays and with bounding modified nonlinearity are studied.
First, oscillations close to harmonic are studied. Using the methods of bifurcation analysis, we singled out the critical case in the problem of stability of stationary state and constructed the normal form. Its nonlocal dynamics determine the local behavior of the initial equation solutions in the neighborhood of the equilibrium state.
The solutions of the relaxation (step-like) type are studied for sufficiently large values of the parameter λ . We stress that, from a computational point of view, Equations (4) and (8) are difficult, since the relaxation solution approaches A and 0 very closely, so even a small error in the calculations will take us out of the class of solutions under consideration. Therefore, asymptotic methods play a special role. They not only allow one to find an approximation of the solution, but also reduce the problem of the dynamics of the original infinite-dimensional problem to the problem of the dynamics of the constructed finite-dimensional mapping (this object is much more simpler than the initial equation). Asymptotic formulas that couple the solutions of (4) and (8) with such mappings trajectories are obtained. The resulting formulas are suitable for engineering calculations. It is important to note that this method is applicable to equations with different types of nonlinearities [30,31] and to systems of two, three, and more singular perturbed equations with delay [30,32,33].
The use of the modified nonlinearity leads to the fact that the oscillations turn out to be ‘safer’ in comparison with Equation (1): their largest value A is significantly less than the value of exp ( λ T ) for Equation (1), and the smallest value of exp ( λ c o n s t ) is significantly greater than the value of exp ( exp λ T ) .
As a generalization of the obtained results, we note that they do not change if the function ( A u ) in (8) is replaced by a more general function F 0 ( u ) for which the conditions F 0 ( u ) > 0 as 0 < u < A , F 0 ( 0 ) = F 0 ( A ) = 0 and F 0 ( 0 ) 0 , F 0 ( A ) 0 are satisfied. We also note that the proposed method allows extension to the case of more than two delays in Equation (8).

Author Contributions

Conceptualization, A.K. and S.K.; Methodology, A.K. and S.K.; Validation, A.K. and S.K.; Formal analysis, A.K. and S.K.; Investigation, A.K. and S.K.; Writing—original draft, A.K. and S.K.; Writing—review and editing, A.K. and S.K.; Visualization, A.K. and S.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. 21-71-30011), https://rscf.ru/en/project/21-71-30011/, accessed on 28 February 2023.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Murray, J.D. Mathematical Biology II, 3rd ed.; Springer: New York, NY, USA, 2003. [Google Scholar] [CrossRef]
  2. Okuba, A. Dynamical aspects of animal grouping: Swarms, schools, flocks, and herds. Adv. Biophys. 1986, 22, 1–94. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, J. Theory and Applications of Partial Functional Differential Equations; Springer: New York, NY, USA, 1996. [Google Scholar] [CrossRef]
  4. Kuang, Y. Delay Differential Equations: With Applications in Population Dynamics; Academic Press: Boston, MA, USA, 1993; Available online: http://www.ams.org/mathscinet-getitem?mr=1218880 (accessed on 1 December 2022).
  5. Wright, E.M. A non-linear difference-differential equation. J. Reine Angew. Math. 1955, 194, 66–87. [Google Scholar] [CrossRef]
  6. Kashchenko, S.A.; Loginov, D.O. Estimation of the Region of Global Stability of the Equilibrium State of the Logistic Equation with Delay. Russ. Math. 2020, 64, 34–49. [Google Scholar] [CrossRef]
  7. Kashchenko, S.A. Asymptotics of the Solutions of the Generalized Hutchinson Equation. Autom. Control Comput. Sci. 2013, 47, 470–494. [Google Scholar] [CrossRef]
  8. May, R.M. Stability and Complexity in Model Ecosystems, 2nd ed.; Princeton University Press: Princeton, NJ, USA, 2001. [Google Scholar] [CrossRef]
  9. Dorini, F.A.; Cecconello, M.S.; Dorini, L.B. On the logistic equation subject to uncertainties in the environmental carrying capacity and initial population density. Commun. Nonlinear Sci. Numer. Simul. 2016, 33, 160–173. [Google Scholar] [CrossRef]
  10. Cushing, J.M. Integrodifferential Equations and Delay Models in Population Dynamics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  11. Brauer, F.; Castillo-Chavez, C. Mathematical Models in Population Biology and Epidemiology; Springer: New York, NY, USA, 2012. [Google Scholar] [CrossRef]
  12. Jin, Z.; Yuan, R. Hopf bifurcation in a reaction-diffusion-advection equation with nonlocal delay effect. J. Differ. Equ. 2021, 271, 533–562. [Google Scholar] [CrossRef]
  13. Goryunov, V.E. Dynamics of solutions of logistic equation with delay and diffusion in a planar domain. Theor. Math. Phys. 2022, 212, 1092–1110. [Google Scholar] [CrossRef]
  14. Sawada, Y.; Takeuchi, Y.; Dong, Y. Stability analysis of a single-species logistic model with time delay and constant inflow. Appl. Math. Lett. 2023, 138, 108514. [Google Scholar] [CrossRef]
  15. Sorokin, V.G.; Vyazmin, A.V. Nonlinear Reaction–Diffusion Equations with Delay: Partial Survey, Exact Solutions, Test Problems, and Numerical Integration. Mathematics 2022, 10, 1886. [Google Scholar] [CrossRef]
  16. Hale, J.K. Theory of Functional Differential Equations, 2nd ed.; Springer: New York, NY, USA, 1977. [Google Scholar] [CrossRef]
  17. Oster, G.; Guckenheimer, J. Chapter Bifurcation Phenomena in Population Models. In The Hopf Bifurcation and Its Applications; Springer: New York, NY, USA, 1976; pp. 327–353. [Google Scholar] [CrossRef]
  18. Marsden, J.E.; McCracken, M.F. The Hopf Bifurcation and Its Applications; Springer: New York, NY, USA, 1976. [Google Scholar] [CrossRef]
  19. Gourley, S.A.; Sou, J.W.H.; Wu, J.H. Nonlocality of Reaction-Diffusion Equations Induced by Delay: Biological Modeling and Nonlinear Dynamics. J. Math. Sci. 2004, 124, 5119–5153. [Google Scholar] [CrossRef]
  20. Kashchenko, S.A. Infinite-Dimensional Bifurcations in Spatially Distributed Delay Logistic Equation. Mathematics 2022, 10, 775. [Google Scholar] [CrossRef]
  21. Jiang, D.; Shi, N. A note on nonautonomous logistic equation with random perturbation. J. Math. Anal. Appl. 2005, 303, 164–172. [Google Scholar] [CrossRef] [Green Version]
  22. Golubenets, V.O. Relaxation oscillations in a logistic equation with nonconstant delay. Math. Notes 2020, 107, 920–932. [Google Scholar] [CrossRef]
  23. Edwards, R.E. Functional Analysis. Theory and Applications; Holt, Rinehart and Winston: New York, NY, USA, 1965. [Google Scholar]
  24. Chen, F.; Shi, J. Periodicity in a logistic type system with several delays. Comput. Math. Appl. 2004, 48, 35–44. [Google Scholar] [CrossRef] [Green Version]
  25. Smith, H. An Introduction to Delay Differential Equations with Applications to the Life Sciences; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar] [CrossRef]
  26. Kashchenko, I.; Kaschenko, S. Infinite process of forward and backward bifurcations in the logistic equation with two delays. Nonlinear Phenom. Complex Syst. 2019, 22, 407–412. [Google Scholar] [CrossRef]
  27. Yan, X.P.; Zhang, C.H. Bifurcation analysis in a diffusive Logistic population model with two delayed density-dependent feedback terms. Nonlinear Anal. Real World Appl. 2022, 63, 103394. [Google Scholar] [CrossRef]
  28. Hartman, P. Ordinary Differential Equations, 2nd ed.; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2002. [Google Scholar] [CrossRef]
  29. Kashchenko, S.A. Relaxation oscillations in a system with delays modeling the predator–prey problem. Autom. Control Comput. Sci. 2015, 49, 547–581. [Google Scholar] [CrossRef] [Green Version]
  30. Kashchenko, A.A. Relaxation modes of a system of diffusion coupled oscillators with delay. Commun. Nonlinear. Sci. Numer. Simulat. 2021, 93, 105488. [Google Scholar] [CrossRef]
  31. Kashchenko, A. Asymptotics of Solutions to a Differential Equation with Delay and Nonlinearity Having Simple Behaviour at Infinity. Mathematics 2022, 10, 3360. [Google Scholar] [CrossRef]
  32. 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]
  33. 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]
Figure 1. The set S.
Figure 1. The set S.
Mathematics 11 01699 g001
Figure 2. Form of the functions u 0 ( t , λ ) , U 0 ( t , λ ) .
Figure 2. Form of the functions u 0 ( t , λ ) , U 0 ( t , λ ) .
Mathematics 11 01699 g002
Figure 3. Graph of φ ( s ) S ( τ 1 , τ 2 ) for s [ 1 , 0 ] and u ( t , φ ) for t 0 .
Figure 3. Graph of φ ( s ) S ( τ 1 , τ 2 ) for s [ 1 , 0 ] and u ( t , φ ) for t 0 .
Mathematics 11 01699 g003
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, S. Relaxation Oscillations in the Logistic Equation with Delay and Modified Nonlinearity. Mathematics 2023, 11, 1699. https://doi.org/10.3390/math11071699

AMA Style

Kashchenko A, Kashchenko S. Relaxation Oscillations in the Logistic Equation with Delay and Modified Nonlinearity. Mathematics. 2023; 11(7):1699. https://doi.org/10.3390/math11071699

Chicago/Turabian Style

Kashchenko, Alexandra, and Sergey Kashchenko. 2023. "Relaxation Oscillations in the Logistic Equation with Delay and Modified Nonlinearity" Mathematics 11, no. 7: 1699. https://doi.org/10.3390/math11071699

APA Style

Kashchenko, A., & Kashchenko, S. (2023). Relaxation Oscillations in the Logistic Equation with Delay and Modified Nonlinearity. Mathematics, 11(7), 1699. https://doi.org/10.3390/math11071699

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