Dynamical Behavior of Discrete Ratio-Dependent Predator-Prey System ()
1. Introduction and Preliminaries
The interaction between predators and their prey has long been the studied topics of ecology and mathematical ecology. In recent years, discrete dynamic models have been widely studied. Compared with continuous dynamic systems, discrete dynamic systems have the following three advantages: Firstly, when populations have non-overlapping generations or the population number is small, discrete dynamic systems are more appropriate than continuous dynamic systems. Secondly, the numerical simulation results obtained from the discrete dynamic system are more accurate than the continuous dynamic system. Thirdly, discrete models exhibit more complex dynamic characteristics than corresponding continuous models [1]. For example, the single-species discrete Logistic model has more complex dynamic behaviors such as bifurcation and chaos.
There are two ways to get discrete system. On the one hand, we start from the continuous model and transform the continuous model into a discrete model by some method, such as Euler iteration [2] and semi-discrete methods [3]. On the other hand, we start from the discrete model directly. Elettreby [4] [5] studied the discrete prey-predator model with mixed functional response from the following discrete model.
(1.1)
where
and
represent population densities of prey and predator at a discrete time step n, respectively; d is the natural death rate of the predator;
and
represent the predator-prey interaction,
is the predator functional responses,
describes how predator transforms the consumed prey into the growth of predator. Because the function
only depends on the density of prey, it is often called a prey-dependent response function. However, more and more evidence shows that in many situations, especially when predators have to search for food, a more appropriate functional response is a function which is dependent on both prey and predator densities.
Differing from the prey-dependent predator-prey systems, the ratio-dependent predator-prey models have two main predictions: a) the equilibrium abundance is positively correlated with the enrichment gradient (see Arditi [6] ) and b) the “paradox of enrichment” (see Rosenzweig [7] ) either completely disappears or enrichment is linked to stability in a more complex way. In this paper, we will study some mathematical characteristic rather than discuss the ecological significance of the model. In 2001, Xiao [8] [9] [10] further studied the ratio-dependent predator-prey model based on Kuang [11]. Xiao made a global qualitative analysis of the model depending on all parameters, and gave the conditions for the existence and nonexistence of limit cycles in the model. Later, the discrete-time dynamical models described by difference equations are widely investigated. Chen [3] showed that positive equilibrium is globally asymptotically stable and established a new sufficient condition. Sohel Rana [12] discretized the continuous model using the forward Euler scheme. It showed that under some parametric conditions, the system undergoes flip and Neimark-Sacker bifurcation.
Arditi [6] proposed the ratio-dependent function in the following form
Then we get the following discrete predator-prey model with the ratio-dependent:
(1.2)
where
is the carrying capacity of the prey, d is the natural death rate of the predator, a is prey intrinsic growth rate, c is capturing rate, m is half saturation constant and f is conversion rate. So a, b, c, d, m and f are positive constants.
The paper is organized as follows. In Section 1, as preliminaries, an important lemma is given. In Section 2, we devote to the existence and local stability of the fixed points of system (1.2). In Sections 3 - 4, we analyze the bifurcation of system (1.2) at the boundary equilibrium point by using the bifurcation theory. Moreover, we proved that the system (1.2) undergoes fold and flip bifurcations for some parameters, respectively. In Section 5, numerical simulations are performed to illustrate the obtained theoretical results and reveal some new dynamical properties of the system (1.2). In Section 6, some conclusions close the paper.
Before analyzing the fixed points of the system (1.2), we introduce an important lemma [2], which will be useful later.
Lemma 1.1 Let
, where B and C are two real constants. Suppose
and
are two roots of
. Then the following statements hold. If
, then
1)
and
if and only if
and
;
2)
and
(or
and
) if and only if
;
3)
and
if and only if
and
;
4)
and
if and only if
and
;
5)
and
if and only if
and
;
6)
and
are a pair of conjugate complex roots and
if and only if
and
.
2. The Existence and Stability of Fixed Points
In this section, we give the conditions for existence of fixed points of system (1.2) in
and discuss the property of these fixed points.
It is clear that the fixed points of the system (1.2) satisfy the following equations:
(2.1)
By a simple computation, we have
Lemma 2.1 System (1.2) has two fixed points:
(i) the boundary equilibrium point
exists if
;
(ii) the unique positive fixed point
, where
and
. The fixed point exists if and only if any one of the following conditions holds:
(ii.1)
and
;
(ii.2)
and
.
Next, we study the stability of the system (1.2) at fixed points. The Jacobian matrix of the system (1.2) evaluated at
is as follows:
(2.2)
The characteristic equation of the matrix can be written as
where
Proposition 2.2 When
, system (1.2) has a boundary equilibrium point
and
(i) it is a sink if
and
;
(ii) it is a source if one of the following conditions holds:
(ii.1)
and
;
(ii.2)
and
;
(ii.3)
and
;
(ii.4)
and
;
(iii) it is non-hyperbolic if
or
or
;
(iv) it is a saddle except if conditions (1)-(3) are satisfied.
Proof. The Jacobian matrix of the fixed point
is
(2.3)
The two eigenvalues of the characteristic equation are:
Then, the fixed point
is a sink if
, where
.
So, the fixed point
is stable if the following condition are satisfied:
The other three items of the propositions are proved by using the above method.
Condition (iii) can be written as
(2.4)
(2.5)
(2.6)
For simple calculation, we denote
and
.
Proposition 2.3. When
and
, system (1.2) has a unique positive fixed point
and
(i) it is a sink if one of the following conditions holds:
(i.1)
,
and
;
(i.2)
,
,
and
;
(i.3)
,
and
;
(i.4)
,
and
;
(i.5)
,
and
;
(i.6)
,
and
;
(i.7)
,
and
;
(ii) it is a source if one of the following conditions holds:
(ii.1)
,
and
;
(ii.2)
,
,
and
;
(ii.3)
,
and
;
(ii.4)
,
and
;
(ii.5)
,
and
;
(ii.6)
,
and
;
(ii.7)
,
and
;
(iii) it is a saddle if one of the following conditions holds:
(iii.1)
,
and
;
(iii.2)
and
;
(iii.3)
,
and
;
(iii.4)
,
and
;
(iv) it is a non-hyperbolic if one of the following conditions holds:
(iv.1)
,
and
;
(iv.2)
and
;
(iv.3)
and
.
Proof. The Jacobian matrix of the fixed point
is
(2.7)
By computation, we have the characteristic equation of
are follows:
where
,
,
and
.
Let
Obviously,
,
(2.8)
(2.9)
(2.10)
From Lemma 1.1, we can see that both
and
if and only if
and
. Then
If
, then
,
and
;
If
, then
and
;
If
, then we can divided d into there cases:
,
and
.
Case I: if
, then
and
;
Case II: if
, then
;
Case III: if
, then
and
.
If
, then
and
;
If
, then
and
;
If
, then
and
.
From
and
, we get conclusion (i). The propositions (ii) and (iii) are proved by using the above method.
Note that the fixed point
is non-hyperbolic if either
or
. There also are two cases:
and
. When
,
is non-hyperbolic if and only if
and
. Then
If
, then
; If
, then
. From
and
, we get conclusion (iv.1, 2).
When
,
is non-hyperbolic if and only if
. Then
From
and
, we get conclusion (iv.3). Thus, we finish the proof.
3. Fold Bifurcation
From Proposition 2.2 (iii), we know that the eigenvalues at
are
and
if both
and
.
Let
,
and
. We transform the fixed pointed
into the origin. After Taylor expansion, the system (1.2) becomes
(3.1)
where
and
When
, we construct an invertible matrix
Let
Then system (3.1) can be transformed into
(3.2)
where
By the center manifold theorem and the method of the references [13], we know that there exists a center manifold
of system (3.2) at the origin in a small neighborhood of
, which can be approximately represented as follows:
For
and
sufficiently small, we assume that a center manifold of the form
where
is a function with order at least three of its variables.
(3.3)
Also,
(3.4)
Compare the coefficients for (3.3) and (3.4), we have
System (3.2) restricted to the center manifold is given by
(3.5)
Let
then
Theorem 3.1 System (1.2) undergoes a fold bifurcation at
when parameter f varies in a small neighbourhood of
.
4. Flip Bifurcation
From Proposition 2.2 (iii), we know that the eigenvalues at
are
and
if both
and
(similarly for the case of
).
Let
,
and
. We transform the fixed pointed
into the origin. After Taylor expansion, the system (1.2) becomes
(4.1)
where
and
When
, we construct an invertible matrix
Let
Then system (4.1) can be transformed into
(4.2)
where
By the center manifold theorem and the method of the references [13], we know that there exists a center manifold
of system (4.2) at the origin in a small neighborhood of
, which can be approximately represented as follows:
For
and
sufficiently small, we assume that a center manifold of the form
where
is a function with order at least three of its variables.
(4.3)
Also,
(4.4)
Compare the coefficients for (4.3) and (4.4), we have
System (4.2) restricted to the center manifold is given by
(4.5)
Let
then
Based on the above analysis, we have the following results.
Theorem 4.1 System (1.2) undergoes a flip bifurcation at
when parameter a varies in a small neighbourhood of
.
5. Numerical Simulations
In this section, we use the bifurcation diagrams, phase portraits diagrams and the Maximum Lyapunov exponents diagrams of the system (1.2) to validate the above analytical result. Then, we consider bifurcation parameters in the following two cases.
Case I: Suppose that parameters
,
,
,
,
in system (1.2). We see that a fold bifurcation occurs at the fixed point
when
. This confirms Theorem 3.1.
Figure 1 shows the bifurcation diagram of system (1.2). We can see that the fixed point
is stable on the center manifold for
. When
the fixed point
on the center manifold becomes unstable. At the same time, there occurs a new fixed point
when
. Hence, the transcritical type of fold bifurcation takes place. The phase diagrams of system (1.2) which are associated with Figure 1 are displayed in Figures 2(a)-(h) with different f.
Case II: Suppose that parameters
,
,
,
,
in system (1.2). We see that a flip bifurcation occurs at the fixed point
when
. This confirms Theorem 4.1.
Figure 3(a) shows the bifurcation diagram of system (1.2). We can see that the fixed point
is stable for
. When
, system (1.2) occurs a subcritical flip bifurcation at the fixed point
. When a increases to 3.6, system (1.2) undergoes a cascade of inverse period-doubling bifurcations. Subsequently, the system (1.2) near the fixed point
is chaotic with narrow periodic windows. The Maximum Lyapunov exponents corresponding to Figure 3(a) are calculated in Figure 3(b), which is consistent with the bifurcation diagram. The phase diagrams of system (1.2) which are associated with Figure 3 are displayed in Figure 4(a)-(h) with different f.
(a) (b)
Figure 1. (a) Bifurcation diagram of the system (1.2) in (f, x)-plane when a = 2, b = 1, c = 1, d = 3, m = 1 and f varies in [3.5, 4.5], the initial value is (0.9, 0.1). (b) Bifurcation diagram of the system (1.2) in (f, y)-plane with same parameters in (a).
(a) (b)
Figure 3. (a) Supercritical flip bifurcation diagram of the system (1.2) on the (a, x)-plane with f = 2.9, b = 1, c = 1, d = 2, m = 1, x0 = 1.9, y0 = 0.1; (b) Maximum Lyapunov exponents corresponding to (a).
6. Conclusion
In this paper, a discrete prey-predator model with ratio-dependent functional response in the closed first quadrant is established. We obtain the equilibrium points of system (1.2) and the conditions for its existence. Moreover, we show that the boundary equilibrium point of system (1.2) can undergo flip bifurcation and Hopf bifurcation by using center manifold theorem and bifurcation theory. The system (1.2) displays interesting dynamical behaviors, including cascade of period-doubling, chaotic with narrow periodic windows and period-2 orbits, which implie that the boundary equilibrium points of system (1.2) lose stability via bifurcation. Nevertheless, the bifurcation of positive equilibrium point has not been discussed. It is a further research topic on how to set reasonable parameters to analyze the bifurcation at the positive equilibrium point.