Next Article in Journal
Unconditional Superconvergence Error Estimates of Semi-Implicit Low-Order Conforming Mixed Finite Element Method for Time-Dependent Navier–Stokes Equations
Previous Article in Journal
Some Properties of the χ-Curvature under Conformal Transformation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability and Neimark–Sacker Bifurcation of a Delay Difference Equation

1
Chinese-German Institute of Engineering, Zhejiang University of Science and Technology, Hangzhou 310023, China
2
School of Science, Zhejiang University of Science and Technology, Hangzhou 310023, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(8), 1942; https://doi.org/10.3390/math11081942
Submission received: 15 March 2023 / Revised: 6 April 2023 / Accepted: 17 April 2023 / Published: 20 April 2023

Abstract

:
In this paper, we revisit a delay differential equation. By using the semidiscretization method, we derive its discrete model. We mainly deeply dig out a Neimark–Sacker bifurcation of the discrete model. Namely, some results for the existence and stability of Neimark–Sacker bifurcation are derived by using the center manifold theorem and bifurcation theory. Some numerical simulations are also given to validate the existence of the Neimark–Sacker bifurcation derived.

1. Introduction

Delay differential equations are primarily used to describe dynamical systems that depend on current and past historical states and have important applications in physics, chemistry, engineering, economics, biology and other fields. Due to the extensive use of delay differential equations in real life, the study of their stability theory has received much attention [1,2,3,4].
Recently, Li et al. [5] studied the chaotic behavior of the following delay difference equation:
x ( n + 1 ) = α x ( n ) + β x ( n k ) ( 1 x 2 ( n k ) ) ,
where α and β are nonzero real parameters and k is a positive integer.
Equation (1) can be viewed as a discrete analogue of the following one-dimensional (or 1D) delay differential equation using the forward Euler scheme [6,7]
d w d t = a w ( t ) + b w ( t τ ) ( 1 w 2 ( t τ ) ) ,
where a > 0 , b is a real parameter and τ > 0 is a delay. Equation (4) is a special case of the following Mackey–Glass equation
d w d t = a w ( t ) + f ( w ( t τ ) ) ,
where a > 0 , τ > 0 is the delay and f is a 1D nonlinear function. Many applications of Equation (5) have been found in physics [8], population dynamics [9], physiology [1], medicine [2], neural control [3] and economics [4].
For the discrete equations of Equation (4), there should be many different forms. Correspondingly, there are also many problems to be considered. Although some good results about the chaotic behavior of discrete Equation (1) have been presented in [5], some other problems of the discrete version of Equation (4), such as bifurcation problems, have not been considered yet. In this paper, we will mainly study the bifurcation problems of its discrete version.
For a complicated ordinary differential equation, generally speaking, it is impossible to solve it accurately. Therefore, one often considers using a computer to derive its numerical solutions, which leads one to consider its corresponding discrete model. There are many different discrete methods and ways for an ordinary differential equation. We will utilize the method of semidiscretization to Equation (4) instead of the forward Euler scheme used in [5,6,7]. For related results for the method of discretization, refer to [8,9].
First, without loss of generality, we can assume τ = 1 in Equation (1). In fact, by taking s = t τ , and letting w ( t ) = w ( s τ ) = η ( s ) , Equation (1) is transformed into
d η d s = a τ η ( s ) + b τ η ( s 1 ) ( 1 η 2 ( s 1 ) ) .
By resetting a and b by a τ and b τ respectively, Equation (4) can be rewritten as
d η d s = a η ( s ) + b η ( s 1 ) ( 1 η 2 ( s 1 ) ) .
This is just (1) in the case of τ = 1 .
Suppose that [ s ] denotes the greatest integer not exceeding s. Consider the following change rate of (5) at the integer point
d η d s = a η ( [ s ] ) + b η ( [ s 1 ] ) ( 1 η 2 ( [ s 1 ] ) ) .
Obviously, the system (6) has piecewise constant arguments. For s [ 0 , + ) , a solution η ( s ) of (6) possesses the following features:
(1)
η ( s ) is continuous on [ 0 , + ) ;
(2)
d η ( s ) d s exists everywhere when s [ 0 , ) except for the points s { 0 , 1 , 2 , 3 , } .
For any s [ n , n + 1 ) with n = 0 , 1 , 2 , , integrating (6) from n to t, one obtains the following equation:
η ( s ) η ( n ) = a η ( n ) + b η ( n 1 ) ( 1 η 2 ( n 1 ) ) ( s n ) .
Letting s ( n + 1 ) in (7) leads to
η ( n + 1 ) = ( 1 a ) η ( n ) + b η ( n 1 ) ( 1 η 2 ( n 1 ) ) ,
which can be viewed as a discrete form of (5), a delay difference equation.
We are now in a position to change the delay difference Equation (8) into a discrete system. Using the transformation
x n = η ( n 1 ) , y n = η ( n ) ,
one has
x n + 1 = y n , y n + 1 = ( 1 a ) y n + b x n ( 1 x n 2 ) ,
which is a discrete version of system (1), and where the parameters
( a , b ) Ω = { ( a , b ) R 2 | a > 0 , b ( , ) } .
In this paper, our main aim is to consider the dynamics of the discrete system (10), namely, for its bifurcation problems except its stability. There have been some studies that consider Neimark–Sacker bifurcation in discrete mathematical models [8,9], whereas less is known about Neimark–Sacker bifurcation occurring in delay difference equations. In order to study the stability and local bifurcation of fixed points, we need a definition [8] and a key lemma [6]. For readers’ convenience, we list them in the Appendix A.
The rest of this paper is organized as follows. In Section 2, we analyze the existence and stability of fixed points of system (10). In Section 3, we discuss its Neimark–Sacker bifurcation. In Section 4, we present some numerical simulations to illustrate the corresponding theoretical analysis results. Finally, we discuss and draw some conclusions in Section 5.

2. Existence and Stability of Fixed Points

In this section, we study the existence and stability of fixed points of system (10).
For the existence of fixed points of system (10), one can easily derive the following results.
Theorem 1.
For the existence of fixed points of system (10), the following statements are valid.
1. 
If b < 0 , then system (10) has three fixed points E 0 = ( 0 , 0 ) , E = ( 1 a b , 1 a b ) and E + = ( 1 a b , 1 a b ) ;
2. 
If 0 b a , then system (10) has a unique fixed point E 0 = ( 0 , 0 ) ;
3. 
If b > a , then system (10) has three fixed points E 0 = ( 0 , 0 ) , E = ( 1 a b , 1 a b ) and E + = ( 1 a b , 1 a b ) .
The Jacobian matrix of system (10) at a fixed point E ( x , y ) is
J ( E ) = 0 1 b ( 1 3 x 2 ) 1 a .
The characteristic polynomial of the Jacobian matrix J ( E ) reads as
F ( λ ) = λ 2 p λ + q with p = 1 a , q = b ( 3 x 2 1 ) .
Because of the symmetry of E + and E , it follows from the characteristic polynomial (11) that their characteristic polynomials are the same. Thus, in the sequel, it suffices for one to only consider the properties of E 0 and E + . Now, we formulate some results for the stability of the fixed points E 0 and E + in the following theorems.
Theorem 2.
The following statements about the fixed point E 0 = ( 0 , 0 ) of system (10) are true.
1. 
If a > b , then,
(a) 
for a < 2 b ,
i. 
when b < 1 , E 0 is a sink;
ii. 
when b = 1 , E 0 is non-hyperbolic;
iii. 
when b > 1 , E 0 is a a source;
(b) 
for a = 2 b , E 0 is non-hyperbolic;
(c) 
for a > 2 b , E 0 is a saddle;
2. 
If a = b , then E 0 is non-hyperbolic;
3. 
If a < b , then,
(a) 
for a < 2 b , E 0 is a saddle;
(b) 
for a = 2 b , E 0 is non-hyperbolic;
(c) 
for a > 2 b , E 0 is a source.
Proof. 
The Jacobian matrix of system (10) at E 0 = ( 0 , 0 ) is
J ( E 0 ) = 0 1 b 1 a .
The characteristic polynomial of the Jacobian matrix J ( E 0 ) can be written as
F ( λ ) = λ 2 p λ + q with p = T r ( J ( E 0 ) ) = 1 a , q = D e t ( J ( E 0 ) ) = b .
Note that F ( 1 ) = a b and F ( 1 ) = 2 a b .
When a > b , F ( 1 ) > 0 . When a < 2 b , F ( 1 ) > 0 . Therefore, for b < 1 , q > 1 . It follows from Lemma A1 (i.1) that | λ 1 | < 1 and | λ 2 | < 1 , so E 0 is a sink. For b = 1 , q = 1 . Lemma A1 (i.5) reads | λ 1 | = | λ 2 | = 1 ; therefore, E 0 is non-hyperbolic. For b > 1 , q < 1 , which reads | λ 1 | > 1 and | λ 2 | > 1 by Lemma A1 (i.4), so E 0 is a source. When a = 2 b , F ( 1 ) = 0 , meaning 1 is one root of the characteristic polynomial; therefore, E 0 is non-hyperbolic. For a > 2 b , F ( 1 ) < 0 , Lemma A1 (i.3) tells us that | λ 1 | < 1 and | λ 2 | > 1 , so E 0 is a saddle.
When a = b , then F ( 1 ) = 0 , meaning 1 is a root of the characteristic equation; therefore, E 0 is non-hyperbolic.
When a < b , F ( 1 ) < 0 . If a < 2 b , then F ( 1 ) > 0 , meaning | λ 1 | > 1 and | λ 2 | < 1 by Lemma A1 (iii.2), so E 0 is a saddle; if a = 2 b , then F ( 1 ) = 0 , meaning 1 is a root of the characteristic equation; hence, E 0 is non-hyperbolic. If a > 2 b , then F ( 1 ) < 0 . In view of Lemma A1 (iii.1), | λ 1 | > 1 and | λ 2 | > 1 , so E 0 is a source. The proof is complete. □
Theorem 3.
For b ( , 0 ) ( a , + ) , the positive fixed point E + = ( 1 a b , 1 a b ) of system (10) occurs. Moreover, the following statements are valid about the positive fixed point E + .
1. 
If b < 0 , then,
(a) 
for a < b + 1 2 , E + is a saddle;
(b) 
for a = b + 1 2 , E + is non-hyperbolic;
(c) 
for a > b + 1 2 , E + is a source.
2. 
If b > a , then,
(a) 
for a < b + 1 2 ,
i. 
when a < 2 b 1 3 , E + is a sink;
ii. 
when a = 2 b 1 3 , E + is non-hyperbolic;
iii. 
when a > 2 b 1 3 , E + is a source;
(b) 
for a = b + 1 2 , E 1 is non-hyperbolic;
(c) 
for a > b + 1 2 , E 1 is a saddle.
Proof. 
The Jacobian matrix of system (10) at E + is given by
J ( E 1 ) = 0 1 3 a 2 b 1 a .
We can express the characteristic equation of J ( E 1 ) as
F ( λ ) = λ 2 p λ + q where p = 1 a , q = 3 a + 2 b .
Obviously, F ( 1 ) = 2 ( b a ) and F ( 1 ) = 2 ( b 2 a + 1 ) = 4 ( b + 1 2 a ) .
When b < 0 , F ( 1 ) < 0 . If a < b + 1 2 , then F ( 1 ) > 0 , meaning | λ 1 | > 1 and | λ 2 | < 1 by Lemma A1 (iii.2), so E + is a saddle; if a = b + 1 2 , then F ( 1 ) = 0 , meaning 1 is a root of the characteristic equation; hence, E + is non-hyperbolic. If a > b + 1 2 , then F ( 1 ) < 0 . In view of Lemma A1 (iii.1), | λ 1 | > 1 and | λ 2 | > 1 , so E + is a source.
When b > a , F ( 1 ) > 0 . For a < b + 1 2 , F ( 1 ) > 0 . Thus, for a < 2 b 1 3 , q > 1 . It follows from Lemma A1 (i.1) that | λ 1 | < 1 and | λ 2 | < 1 , so E + is a sink. For a = 2 b 1 3 , q = 1 . Lemma A1 (i.5) reads | λ 1 | = | λ 2 | = 1 ; therefore, E + is non-hyperbolic. For a > 2 b 1 3 , q < 1 , which reads | λ 1 | > 1 and | λ 2 | > 1 by Lemma A1 (i.4), so E + is a source.
For a = b + 1 2 , F ( 1 ) = 0 , meaning 1 is a root of the characteristic polynomial; therefore, E + is non-hyperbolic.
For a > b + 1 2 , F ( 1 ) < 0 , Lemma A1 (i.3) tells us that | λ 1 | < 1 and | λ 2 | > 1 , so E + is a saddle.
The proof is over. □

3. Bifurcation Analysis

In this section, we focus on the bifurcation problem of system (10), namely, for its Neimark–Sacker bifurcation by using the center manifold theorem and bifurcation theory in [1,2,3,4]. For the discrete bifurcation results, also refer to [10,11,12,13,14] and the references cited therein. It suffices for us to study system (10) at the positive fixed point E + . The conclusion for system (10) at the negative fixed point E is the same.
Suppose the parameters
( a , b ) Ω + = { ( a , b ) Ω | b ( , 0 ) ( a , ) } ,
ensuring the existence of the positive fixed point E + . Now let us consider ( a , b ) Ω + and a ( 0 , b ) . Let a 0 = 2 b 1 3 .
When the parameter a goes through the critical value a 0 , it follows from Theorem 3 1(b) that the dimensions of the unstable manifold and stable manifold of system (10) at fixed point E + change. Therefore, system (10) may undergo a bifurcation at the fixed point E + . Furthermore, at this time, system (10) has a pair of conjugate complex roots λ 1 and λ 2 satisfying | λ 1 | = | λ 2 | = 1 . Thus, a phenomenon of Neimark–Sacker bifurcation may occur. Now we analyze the process.
Take l n = x n 1 a b and m n = y n 1 a b . Then, system (10) reads as
l n + 1 = m n , m n + 1 = ( 3 a 2 b ) l n + ( 1 a ) m n 3 b 1 a b l n 2 b l n 3 .
Choose the parameter a as bifurcation parameter. Give a small perturbation a of the parameter a around a 0 , i.e., a = a a 0 , with 0 < | a | 1 . Then system (12) is perturbed into
l n + 1 m n + 1 = 0 1 1 + 3 a 4 2 b 3 a l n m n + 0 3 b ( b + 1 3 a ) l n 2 b l n 3 .
Denote the characteristic polynomial of the Jacobian matrix of the linearized equation associated with system (13) as F ( λ ) = λ 2 p ( a ) λ + q ( a ) , where
p ( a ) = 4 2 b 3 a , q ( a ) = 1 3 a .
Then the two roots of F ( λ ) = 0 are
λ 1 , 2 ( a ) = p ( a ) ± p 2 ( a ) 4 q ( a ) 2 .
Noticing the parameter vector ( a , b ) Ω + , it is not hard to derive 4 q ( a ) p 2 ( a ) > 0 when a = 0 and b < 5 ( a = 0 means 0 < a = a 0 = 2 b 1 3 . So, b > 1 2 ). Therefore,
λ 1 , 2 ( 0 ) = p ( 0 ) ± i 4 q ( 0 ) p 2 ( 0 ) 2 = 2 b ± i ( b + 1 ) ( 5 b ) 3 .
It is easy to obtain, for 0 < | a | 1 ,
| λ 1 ( a ) | = | λ 2 ( a ) | = q ( a ) ,
and hence,
d | λ 1 ( a ) | d a | a = 0 = d | λ 2 ( a ) | d a | a = 0 = 3 2 < 0 .
Moreover, it is obvious that λ 1 , 2 m ( 0 ) 1 for all m = 1 , 2 , 3 , 4 for b 2 . Thus, all of the conditions for Neimark–Sacker bifurcation to happen are satisfied.
Now we are in a position to look for the normal form of system (13) when a = 0 . Then system (13) can be regarded as
l n + 1 m n + 1 = 0 1 1 4 2 b 3 l n m n + 0 G ( l n , m n ) ,
where
G ( l n , m n ) = 3 b ( b + 1 ) l n 2 b l n 3 .
It is easy to derive the two eigenvalues of the matrix A = 0 1 1 4 2 b 3 to be λ 1 ( 0 ) and λ 2 ( 0 ) with corresponding eigenvectors ξ 1 = 1 2 b 3 and ξ 2 = 0 1 3 ( b + 1 ) ( 5 b ) . Let T = ( ξ 1 , ξ 2 ) , i.e.,
T = 1 0 2 b 3 1 3 ( b + 1 ) ( 5 b ) , then T 1 = 1 0 2 b ( b + 1 ) ( 5 b ) 3 ( b + 1 ) ( 5 b ) .
The transformation l n m n = T u n v n brings system (14) to
u n + 1 v n + 1 = 2 b 3 1 3 ( b + 1 ) ( 5 b ) 1 3 ( b + 1 ) ( 5 b ) 2 b 3 u n v n + f ¯ ( u n , v n ) g ¯ ( u n , v n ) ,
where
f ¯ ( u n , v n ) = 0 , g ¯ ( u n , v n ) = 3 ( b + 1 ) ( 5 b ) G ( u n , 2 b 3 u n 1 3 ( b + 1 ) ( 5 b ) v n ) = 3 b ( b + 1 ) ( 5 b ) 3 b + 1 3 b u n 2 u n 3 .
Furthermore,
f ¯ u u = f ¯ u v = f ¯ v v = f ¯ u u u = f ¯ u u v = f ¯ u v v = f ¯ v v v = 0 , g ¯ u v = g ¯ v v = g ¯ u u v = g ¯ u v v = g ¯ v v v = 0 , g ¯ u u = 18 b 3 b ( 5 b ) , g ¯ u u u = 18 b ( b + 1 ) ( 5 b ) .
Next, compute the quantity L that is used to judge the stability and direction of a Neimark–Sacker bifurcation (see [1,2,3,4]):
L = R e ( 1 2 λ 1 ) λ 2 2 1 λ 1 ξ 20 ξ 11 1 2 | ξ 11 | 2 | ξ 02 | 2 + R e ( λ 2 ξ 21 ) ,
where
ξ 20 = 1 8 [ f ¯ u u f ¯ v v + 2 g ¯ u v + i ( g ¯ u u g ¯ v v 2 f ¯ u v ) ] ,
ξ 11 = 1 4 [ f ¯ u u + f ¯ v v + i ( g ¯ u u + g ¯ v v ) ] ,
ξ 02 = 1 8 [ f ¯ u u f ¯ v v 2 g ¯ u v + i ( g ¯ u u g ¯ v v + 2 f ¯ u v ) ] ,
ξ 21 = 1 16 [ f ¯ u u u + f ¯ v v v + g ¯ u u v + g ¯ v v v + i ( g ¯ u u u + g ¯ u v v f ¯ u u v f ¯ v v v ) ] .
Some calculations display
ξ 20 = 9 b i 4 3 b ( 5 b ) , ξ 11 = 9 b i 2 3 b ( 5 b ) , ξ 02 = 9 b i 4 3 b ( 5 b ) , ξ 21 = 9 b i 8 ( b + 1 ) ( 5 b ) , L = 3 ( b 2 + 4 b + 5 ) 8 ( b + 1 ) < 0 .
Summarizing the above analysis, we have the following consequence.
Theorem 4.
Assume the parameters a , b Ω satisfy b ( 1 2 , 2 ) ( 2 , 5 ) and a ( 0 , b ) . Denote a 0 = 2 b 1 3 . Then system (10) undergoes a Neimark–Sacker bifurcation at the positive fixed point E + and the negative fixed point E , respectively, when the parameter a varies in a small neighborhood of a 0 . Moreover, an attracting invariant closed curve bifurcates from the positive fixed point E + and the negative fixed point E , respectively, for a > a 0 .

4. Numerical Simulations

In this section, to illustrate our theoretical results obtained and reveal some new dynamical behaviors in the system (10), we present the bifurcation diagrams, phase portraits and Lyapunov exponents for specific parameter values using Matlab software (R2023a). For similar numerical simulation work, we refer readers to [6,8,9,10,11,12,13,14].
For the fixed point E + , choose the parameters a ( 0.5 , 0.7 ) , b = 1.4 and the initial values ( x 0 , y 0 ) = ( 0.05 , 0.05 ) , ( 0.75 , 0.75 ) . From Figure 1a we can find that when a = a 0 = 0.6 the system (10) undergoes a Neimark–Sacker bifurcation. In order to bear out it, we take a near 0.6 and obtain Figure 1, Figure 2, Figure 3 and Figure 4a–d, which illustrate the existence of Neimark–Sacker bifurcation at the fixed point E + ( 4 / 7 , 4 / 7 ) ( 0.756 , 0.756 ) .  Figure 2a–d mean that the fixed point E + is a stable attractor when a 0 . 6 + . Moreover, Figure 2a,b and Figure 3a show the occurrence of Neimark–Sacker bifurcation when a 0 . 6 . Figure 3 and Figure 4a–d illustrate that increasing the value of a leads to a change of stability of the fixed point E + and the occurrence of an invariant closed curve around E + , which fits the results in Theorem 4. The spectrum of the maximum Lyapunov exponent with respect to the parameter a ( 0.5 , 0.7 ) when b = 1.4 is presented in Figure 1b.
Remark 1. 
Figure 3a–d show that the bifurcated closed orbit is stable outside whereas Figure 4a–d show that the bifurcated closed orbit is stable inside. Thus, the closed orbit bifurcated from E + is stable. This is in accord with our Theorem 4.

5. Conclusions

In this paper, we apply the semidiscretization method to a delay differential equation considered in [5] and derive a delay difference equation. By using the bifurcation theory, we mainly obtain some results for the Neimark–Sacker bifurcation of the discrete model.
Neimark–Sacker bifurcation is an important mechanic for one system to produce complicated dynamical behaviors. The occurrence of a Neimark–Sacker bifurcation often causes the system to jump from a stable window to chaotic states through periodic and quasi-periodic states and triggers a route to chaos. We find this mechanic in the delay difference equation. Whether or not there are other similar bifurcation problems, such as flip bifurcation, fold bifurcation, etc., in other delay difference equations, and how to find such bifurcations, will be an interesting and important topic and a direction worth future study.

Author Contributions

Methodology, X.L.; Software, S.J. All authors contributed equally and significantly in writing this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (grant: 61473340), Distinguished Professor Foundation of Qianjiang Scholar in Zhejiang Province (grant: F703108L02) and the Natural Science Foundation of Zhejiang University of Science and Technology (grant: F701108G14).

Data Availability Statement

Data is unavailable due to privacy or ethical restrictions.

Acknowledgments

This work was finished in discussion with Zhikang Pan.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Definition A1. 
Let E ( x , y ) be a fixed piont of system (10) with multipliers λ 1 and λ 2 .
(1) 
If | λ 1 | < 1 and | λ 2 | < 1 , a fixed point E ( x , y ) is called a sink, so a sink is locally asymptotically stable.
(2) 
If | λ 1 | > 1 and | λ 2 | > 1 , a fixed point E ( x , y ) is called a source, so a source is locally asymptotically unstable.
(3) 
If | λ 1 | < 1 and | λ 2 | > 1 (or | λ 1 | > 1 and | λ 2 | < 1 ), a fixed point E ( x , y ) is called a saddle.
(4) 
If either | λ 1 | = 1 or | λ 2 | = 1 , a fixed point E ( x , y ) is called non-hyperbolic.
Lemma A1. 
Let F ( λ ) = λ 2 + B λ + C , where B and C are two real constants. Suppose λ 1 and λ 2 are two roots of F ( λ ) = 0 . Then the following statements hold.
(i) 
If F ( 1 ) > 0 , then
(i.1) 
| λ 1 | < 1 and | λ 2 | < 1 if and only if F ( 1 ) > 0 and C < 1 ;
(i.2) 
λ 1 = 1 and λ 2 1 if and only if F ( 1 ) = 0 and B 2 ;
(i.3) 
| λ 1 | < 1 and | λ 2 | > 1 if and only if F ( 1 ) < 0 ;
(i.4) 
| λ 1 | > 1 and | λ 2 | > 1 if and only if F ( 1 ) > 0 and C > 1 ;
(i.5) 
λ 1 and λ 2 are a pair of conjugate complex roots and | λ 1 | = | λ 2 | = 1 if and only if 2 < B < 2 and C = 1 ;
(i.6) 
λ 1 = λ 2 = 1 if and only if F ( 1 ) = 0 and B = 2 .
(ii) 
If F ( 1 ) = 0 , namely, 1 is one root of F ( λ ) = 0 , then the other root λ satisfies | λ | = ( < , > ) 1 if and only if | C | = ( < , > ) 1 .
(iii) 
If F ( 1 ) < 0 , then F ( λ ) = 0 has one root lying in ( 1 , ) . Moreover,
(iii.1) 
the other root λ satisfies λ < ( = ) 1 if and only if F ( 1 ) < ( = ) 0 ;
(iii.2) 
the other root 1 < λ < 1 if and only if F ( 1 ) > 0 .

References

  1. Carr, J. Application of Center Manifold Theorem; Springer: New York, NY, USA, 1981. [Google Scholar]
  2. Kuzenetsov, Y.A. Elements of Applied Bifurcation Theory, 2nd ed.; Springer: New York, NY, USA, 1998. [Google Scholar]
  3. Robinson, C. Dynamical Systems: Stability, Symbolic Dynamics, and Chaos, 2nd ed.; Amazon: Boca Raton, FL, USA; London, UK; New York, NY, USA, 1999. [Google Scholar]
  4. Winggins, S. Introduction to Applied Nonlinear Dynamical Systems and Chaos; Springer: New York, NY, USA, 2003. [Google Scholar]
  5. Li, Z.C.; Zhao, Q.L.; Liang, D. Chaotic behavior in a class of delay difference equations. Adv. Differ. Equ. 2013, 2013, 99. [Google Scholar] [CrossRef]
  6. Wang, C.; Li, X.Y. Further investigations into the stability and bifurcation of a discrete predator-prey model. J. Math. Anal. Appl. 2015, 422, 920–939. [Google Scholar] [CrossRef]
  7. Din, Q. Compexity and chaos control in a discrete-time prey-predator model. Commun. Nonlinear Sci. Numer. Simul. 2017, 49, 113–134. [Google Scholar] [CrossRef]
  8. Yao, W.B.; Li, X.Y. Complicate bifurcation behaviors of a discrete predator–prey model with group defense and nonlinear harvesting in prey. Appl. Anal. 2022. [Google Scholar] [CrossRef]
  9. Pan, Z.K.; Li, X.Y. Stability and Neimark–Sacker bifurcation for a discrete Nicholson’s blowflies model with proportional delay. J. Differ. Equ. Appl. 2021, 27, 250–260. [Google Scholar] [CrossRef]
  10. Rana, S.S. Bifurcations and chaos control in a discrete-time predator-prey system of Leslie type. J. Appl. Anal. Comput. 2019, 9, 31–44. [Google Scholar]
  11. Atabaigi, A. Bifurcation and chaos in a discrete time predator-prey system of Leslie type with generalized Holling type III functional response. J. Appl. Anal. Comput. 2017, 7, 411–426. [Google Scholar]
  12. Singh, A.; Deolia, P. Dynamical analysis and chaos control in discrete-time prey-predator model. Commun. Nonlinear Sci. Numer. Simul. 2020, 90, 105313. [Google Scholar] [CrossRef]
  13. Naik, P.A.; Eskandari, Z.; Avazzadeh, Z.; Zu, J. Multiple Bifurcations of a discrete-time prey-predator model with mixed functional response. Int. J. Bifurc. Chaos 2022, 32, 2250050. [Google Scholar] [CrossRef]
  14. Naik, P.A.; Eskandari, Z.; Yavuz, M.; Zu, J. Complex dynamics of a discrete-time Bazykin-Berezovskaya prey-predator model with a strong Allee effect. J. Comput. Appl. Math. 2022, 413, 114401. [Google Scholar] [CrossRef]
Figure 1. Bifurcation of system (10) in ( a , x ) plane and maximal Lyapunov exponent for b = 1.4 .
Figure 1. Bifurcation of system (10) in ( a , x ) plane and maximal Lyapunov exponent for b = 1.4 .
Mathematics 11 01942 g001
Figure 2. Phase portraits for system (10) with b = 1.4 and different a when the initial value ( x 0 , y 0 ) = ( 0.05 , 0.05 ) .
Figure 2. Phase portraits for system (10) with b = 1.4 and different a when the initial value ( x 0 , y 0 ) = ( 0.05 , 0.05 ) .
Mathematics 11 01942 g002
Figure 3. Phase portraits for system (10) with b = 1.4 and different a when the initial value ( x 0 , y 0 ) = ( 0.05 , 0.05 ) is outside the closed orbit.
Figure 3. Phase portraits for system (10) with b = 1.4 and different a when the initial value ( x 0 , y 0 ) = ( 0.05 , 0.05 ) is outside the closed orbit.
Mathematics 11 01942 g003aMathematics 11 01942 g003b
Figure 4. Phase portraits for system (10) with b = 1.4 and different a when the initial value ( x 0 , y 0 ) = ( 0.75 , 0.75 ) is inside the closed orbit.
Figure 4. Phase portraits for system (10) with b = 1.4 and different a when the initial value ( x 0 , y 0 ) = ( 0.75 , 0.75 ) is inside the closed orbit.
Mathematics 11 01942 g004
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

Jin, S.; Li, X. Stability and Neimark–Sacker Bifurcation of a Delay Difference Equation. Mathematics 2023, 11, 1942. https://doi.org/10.3390/math11081942

AMA Style

Jin S, Li X. Stability and Neimark–Sacker Bifurcation of a Delay Difference Equation. Mathematics. 2023; 11(8):1942. https://doi.org/10.3390/math11081942

Chicago/Turabian Style

Jin, Shaoxia, and Xianyi Li. 2023. "Stability and Neimark–Sacker Bifurcation of a Delay Difference Equation" Mathematics 11, no. 8: 1942. https://doi.org/10.3390/math11081942

APA Style

Jin, S., & Li, X. (2023). Stability and Neimark–Sacker Bifurcation of a Delay Difference Equation. Mathematics, 11(8), 1942. https://doi.org/10.3390/math11081942

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