Next Article in Journal
Time-Fractional Diffusion-Wave Equation with Mass Absorption in a Sphere under Harmonic Impact
Next Article in Special Issue
On the Oscillation of Non-Linear Fractional Difference Equations with Damping
Previous Article in Journal
Truncated Fubini Polynomials
Previous Article in Special Issue
Global Dynamics of Leslie-Gower Competitive Systems in the Plane
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability and Bifurcation Analysis on a Predator–Prey System with the Weak Allee Effect

1
Department of Mathematics, School of Science, Zhejiang Sci-Tech University, Hangzhou 310018, China
2
College of Mathematics and Systems Science, Shandong University of Science and Technology Qingdao 266590, China
3
School of Mathematical Sciences, Qufu Normal University, Qufu 273165, China
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(5), 432; https://doi.org/10.3390/math7050432
Submission received: 9 April 2019 / Revised: 26 April 2019 / Accepted: 6 May 2019 / Published: 16 May 2019

Abstract

:
In this paper, the dynamics of a predator-prey system with the weak Allee effect is considered. The sufficient conditions for the existence of Hopf bifurcation and stability switches induced by delay are investigated. By using the theory of normal form and center manifold, an explicit expression, which can be applied to determine the direction of the Hopf bifurcations and the stability of the bifurcating periodic solutions, are obtained. Numerical simulations are performed to illustrate the theoretical analysis results.

1. Introduction

The Allee effect refers to the phenomenon that individuals of many species benefit from the presence of conspecifics, which implies when the per capita growth rate becomes bigger with the increase of the density and gets smaller as long as density passes through a critical value [1]. When the growth rate decreases but still keeps positive at a low population density, the Allee effect is called weak. It is easy to see that this situation is totally different from the logistic growth. During the last decade this phenomenon attracted a lot of attention with the rise of conservation biology.
The dynamics with the Allee effect, such as the stability of positive equilibrium [2,3,4,5,6,7], the existence of periodic solution [8,9,10,11,12], or bifurcation analysis [13,14,15], on nonlinear systems have been extensively investigated. Especially, Boukal and L. Berec [16] critically reviewed and classified models of single-species population dynamics with the demographic Allee effect. Meanwhile in consideration of the carrying capacity of environment with respect to the prey in the per capita growth rate of the population, they also proposed the following differential equation with the weak Allee effect
N ˙ ( t ) = ε K N 2 ( t ) 1 N ( t ) K ,
where N ( t ) denotes the prey population at moment t, 0 < ε < 1 is the per capita growth rate of the population and K > 0 is the carrying capacity of the environment. In biology, the group defense enriches the ability of the prey to defend or escape from the predators, and thus it helps decrease the predation. Tabares et al. [17] proposed the coupled system of equations
N ˙ ( t ) = 4 ε K N 2 ( t ) 1 N ( t ) K α N ( t ) P ( t ) , P ˙ ( t ) = γ P ( t ) + a β P ( t ) t N ( τ ) exp ( a ( t τ ) ) d τ .
as a model to describe the dynamical behaviour of the predator-prey system incorporating a weak Allee effect in the prey population N ( t ) and a distributed delay in the predator populations P ( t ) at time t. They analyzed the influence of weak Allee effect in the system. Here ε , α , β , γ , a and K are all positive parameters. α is the rate of predation, β is the conversion rate, γ is the predator mortality in the absence of prey and a is the delay parameter. K denotes the carrying capacity of the environment with respect to the prey. Equation (2) is regarded as a predator-prey system with group defense [18]. J. D. Ferreira et al. [1] proposed the system
N ˙ ( t ) = 4 ε K N 2 ( t ) 1 N ( t ) K α N ( t ) P ( t ) , P ˙ ( t ) = γ P ( t ) + β P ( t ) t N ( τ ) G ( t τ ) d τ ,
which describes the dynamics of interactions between a predator and a prey species with weak Allee effect. They studied the influence of the weak Allee effect with the density function G : R + R + defined by
G ( t ) = a 2 t e a t , a > 0
and got analogous stability results as (2). For all biological processes in nature, the time delay is an inherent phenomenon. Hence, it is practical to assume that the growth rate of the prey population is determined by its recent history. So taking into account the maturation time of the prey, we define the density function G ( t ) by
G ( t ) = a e a t , a > 0 ,
and it makes system (2) the following delayed predator-prey model with weak Allee effect.
N ˙ ( t ) = 4 ε K N 2 ( t ) 1 N ( t τ ) K α N ( t ) P ( t ) , P ˙ ( t ) = γ P ( t ) + a β P ( t ) t N ( τ ) exp ( a ( t τ ) ) d τ .
In this paper, we study the effect of the delay τ on the positive equilibrium of system (4), and investigate the normal form for the Hopf bifurcation. In Section 2, stability of the equilibrium and existence of Hopf bifurcations are given. In Section 3, the direction and stability of the bifurcating periodic solutions of the system will be analyzed. Finally, in Section 4, numerical simulations are performed to illustrate our theoretical results.

2. Stability of the Equilibrium and Existence of Hopf Bifurcations

System (4) can be rewritten as
N ˙ ( t ) = 4 ε K N 2 ( t ) 1 N ( t τ ) K α N ( t ) P ( t ) , P ˙ ( t ) = γ P ( t ) + β P ( t ) Q ( t ) , Q ˙ ( t ) = a ( N ( t ) Q ( t ) ) ,
where
Q ( t ) = a t N ( τ ) exp ( a ( t τ ) ) d τ .
Introducing N = K x 1 , P = K x 2 , Q = K x 3 , t = t ¯ ε and then dropping the bar, system (5) becomes
x ˙ 1 ( t ) = 4 x 1 2 ( t ) ( 1 x 1 ( t τ ) ) K α ε x 1 ( t ) x 2 ( t ) , x ˙ 2 ( t ) = γ ε x 2 ( t ) + K β ε x 2 ( t ) x 3 ( t ) , x ˙ 3 ( t ) = a ε ( x 1 ( t ) x 3 ( t ) ) .
When τ = 0 , system (6) is system (3) in [17]. System (6) has only one positive equilibrium E γ β K , 4 ε γ ( β K γ ) α β 2 K 3 , γ β K . In [17], sufficient conditions, which guarantee the positive equilibrium is locally asymptotically stable for system (6) with τ = 0 , have been given. We now recall some of their conclusions.
Lemma 1.
Equilibrium E of system (6) is locally asymptotically stable provided the following inequalities are satisfied
β K γ > 0 ,
2 γ β K > 0 ,
a β 2 K 2 + 4 γ ε ( 2 γ β K ) β 2 ε K 2 × 4 a γ ( 2 γ β K ) β 2 ε K 2 4 a γ 2 ( β K γ ) β 2 ε 2 K 2 > 0 .
Let b = 1 β K , then Equations (7)–(9) are equivalent to
1 2 < b γ < 1 ,
a > γ [ 1 ( γ + 4 ε ) b + 16 ε γ b 2 16 ε γ 2 b 3 ] 2 b γ 1 .
Denote the three curves on γ b plane b γ = 1 2 , b γ = 1 and a = γ [ 1 ( γ + 4 ε ) b + 16 ε γ b 2 16 ε γ 2 b 3 ] 2 b γ 1 by l 1 , l 2 and l 3 , respectively. For the case when a = 0.2 , ε = 0.15 , the three curves divide the first quadrant of the γ b plane into four regions D i , i = 1 , 2 , 3 , 4 (see Figure 1a). We know from Lemma 1 that the positive equilibrium of system (6) is stable when τ = 0 .
The characteristic equation of system (6) at E is
λ 4 b γ ( 1 b γ b γ e λ τ ) α K b γ ε 0 0 λ 4 γ ( 1 b γ ) α K a ε 0 λ + a ε = 0 ,
i.e.,
λ 3 + a ε 4 b γ ( 1 b γ ) λ 2 4 a b γ ε ( 1 b γ ) λ + 4 a b γ 2 ( 1 b γ ) ε 2 + 4 b 2 γ 2 λ 2 + a ε λ e λ τ = 0 .
Clearly, λ = 0 is not a root of (13). Let λ = i ω ( ω > 0 ) be a root of (13), then
i ω 3 r 1 ω 2 + i r 2 ω + r 3 + r 4 ( ω 2 + i a ω ε ) e i ω τ = 0 ,
where r 1 = a ε 4 b γ ( 1 b γ ) , r 2 = 4 a b γ ε ( 1 b γ ) , r 3 = 4 a b γ 2 ( 1 b γ ) ε 2 and r 4 = 4 b 2 γ 2 . Separating the real and imaginary parts leads to
r 1 ω 2 + r 3 + r 4 ω ( ω cos ( ω τ ) + a ε sin ( ω τ ) ) = 0 , ω 3 + r 2 ω + r 4 ω ( a ε cos ( ω τ ) + ω sin ( ω τ ) ) = 0 .
It follows directly from (15) that
ω 6 + ( r 1 2 2 r 2 r 4 2 ) ω 4 + ( r 2 2 2 r 1 r 3 a 2 r 4 2 ε 2 ) ω 2 + r 3 2 = 0 .
In fact, (15) can be rewritten as
r 1 ω 2 + r 3 = r 4 ω ( ω cos ( ω τ ) + a ε sin ( ω τ ) ) , ω 3 + r 2 ω = r 4 ω ( a ε cos ( ω τ ) + ω sin ( ω τ ) ) .
Now, taking the square of two sides of the equations and adding we obtain (16) directly.
Let ω 2 = z , then (16) can be rewritten as
f ( z ) z 3 + ( r 1 2 2 r 2 r 4 2 ) z 2 + ( r 2 2 2 r 1 r 3 a 2 r 4 2 ε 2 ) z + r 3 2 = 0 .
Since the constant term in (17) is positive, f ( z ) = 0 has at least one positive root if and only if there exists z * > 0 such that
f ( z * ) = 0 a n d f ( z * ) 0 .
In fact,
f ( z ) = 3 z 2 + 2 A z + B ,
where A = r 1 2 2 r 2 r 4 2 and B = r 2 2 2 r 1 r 3 a 2 r 4 2 ε 2 . When B < 0 , it is clearly that f ( z ) = 0 has a positive root z * = A + A 2 3 B 3 . Through direct calculation, we get
B = 8 a b γ 2 ε 2 2 a b ( 1 b γ ) 2 2 a b 3 γ 2 a ε ( 1 b γ ) + 4 b γ ( 1 b γ ) 2 .
Thus, B < 0 if
2 a b ( 1 b γ ) 2 2 a b 3 γ 2 a ε ( 1 b γ ) + 4 b γ ( 1 b γ ) 2 < 0 .
We rewrite (19) as
a > 4 b ε γ ( 1 b γ ) 2 1 b γ + 2 ε b 3 γ 2 2 b ε ( 1 b γ ) 2 .
Denote l 4 : a = 4 b ε γ ( 1 b γ ) 2 1 b γ + 2 ε b 3 γ 2 2 b ε ( 1 b γ ) 2 and l 5 : f ( z * ) = 0 , then curve l 4 divides D 4 into two regions D 5 and D 6 , curve l 5 divides D 3 into two regions D 7 and D 8 (see Figure 1b). This means that Equation (17) has two positive roots in region D 7 . Thus, we can draw the following conclusions.
Lemma 2.
Assume (10), (11), (18) and (20) hold,
(i) when f ( z * ) = 0 , f ( z ) = 0 has one positive root;
(ii) when f ( z * ) < 0 , f ( z ) = 0 has two positive roots.
A positive root of f ( z ) = 0 corresponds to a pair of purely imaginary roots of Equation (3). Thus we get the following results.
Lemma 3.
Assume (10), (11), (18) and (20) hold.
(i) if f ( z ) = 0 has one positive root z * , then (16) has one pair of purely imaginary roots ± i ω 0 with double, where ω 0 2 = z * ;
(ii) if f ( z ) = 0 has two positive roots z + and z , then (16) has two pair of purely imaginary roots ± i ω + and ± i ω , where ω + 2 = z + and ω 2 = z . Furthermore f ( z + ) > 0 and f ( z ) < 0 .
Denote ω * as a positive root of (16). By (15), we have
sin ω * τ = ω * 4 + a ε r 1 r 2 ω * 2 a ε r 3 r 4 ω * ω * 2 + a 2 ε 2 .
cos ω * τ = a ε r 1 ω * 3 + r 3 a r 2 ε ω * r 4 ω * ω * 2 + a 2 ε 2 G ( ω * ) .
Since a ε r 1 = 4 b γ ( 1 b γ ) > 0 and r 3 a r 2 ε = 4 a b γ ε 2 ( 1 b γ ) ( γ + a ) > 0 , we known cos ω * τ > 0 .
Let
h ( z ) = z 2 + a ε r 1 r 2 z a ε r 3
Since ( a ε r 1 r 2 ) 2 + 4 a ε r 3 > 0 and a ε r 3 < 0 , then h ( z ) = 0 has one positive root and one negative root. Denote the positive root by z 0 . According to the discussions above and Ruan. S et al. [19], we get the following conclusion.
Lemma 4.
(i) If (10), (11), (18) and (20) are not all satisfied, then Equation (13) has at least one root with positive real part for τ > 0 ;
(ii) If (10), (11), (20) and f ( z * ) = 0 hold, then there exists τ 0 + < τ 1 + < < τ j + < , such that Equation (13) has a pair of purely imaginary roots ± i ω + with double when τ = τ j + ;
(iii) If (10), (11), (20) and f ( z * ) < 0 hold, then there exists τ 0 + < τ 1 + < < τ j + < and τ 0 < τ 1 < < τ j < such that Equation (13) has a pair of purely imaginary roots ± i ω ± when τ = τ j ± , respectively, where
τ j + = 1 ω + [ arccos G ( ω + ) + 2 j π ] , z 0 < z + , 1 ω + [ arccos G ( ω + ) + 2 ( j + 1 ) π ] , z 0 > z + ,
and
τ j = 1 ω [ arccos G ( ω ) + 2 ( j + 1 ) π ] , z 0 > z , 1 ω [ arccos G ( ω ) + 2 j π ] , z 0 < z .
Let λ ( τ ) = r ( τ ) + i ω ( τ ) be the root of Equation (13) satisfying r ( τ j ± ) = 0 , ω ( τ j ± ) = ω ± , respectively. Substituting λ ( τ ) = r ( τ ) + i ω ( τ ) into Equation (13) and taking derivative with respect to τ , we get
d λ d τ = λ r 4 λ 2 + a ε λ e λ τ 3 λ 2 + 2 r 1 λ + r 2 + r 4 2 λ + a ε e λ τ r 4 τ λ 2 + a ε λ e λ τ = r 4 ω 2 a ε cos ω τ + ω sin ω τ + i ω cos ω τ a ε sin ω τ H ( ω , τ ) ,
where
H ( ω , τ ) = 3 ω 2 + r 2 + r 4 a ε ω 2 τ cos ω τ + 2 ω + a ω τ ε sin ω τ + i [ 2 r 1 ω + r 4 2 ω + a ω τ ε cos ω τ + ω 2 τ a ε sin ω τ ,
It follows that
R e ( d λ d τ ) = ω 2 r 4 Δ × { { 3 ω 2 + r 2 + r 4 a ε ω 2 τ cos ω τ + 2 ω + a ω τ ε sin ω τ } ( a ε cos ω τ + ω sin ω τ ) + [ 2 r 1 ω + r 4 2 ω + a ω τ ε cos ω τ + ω 2 τ a ε sin ω τ ω cos ω τ a ε sin ω τ } = ω 2 Δ 3 ω 4 2 ( r 1 2 2 r 2 r 4 ) ω 2 ( r 2 2 2 r 1 r 3 r 4 2 a 2 ε 2 ) = ω 2 Δ 3 z ± 2 2 ( r 1 2 2 r 2 r 4 ) z ± ( r 2 2 2 r 1 r 3 r 4 2 a 2 ε 2 ) = ω 2 Δ f ( z ± ) ,
namely, r ( τ j + ) < 0 and r ( τ j ) > 0 ( j = 0 , 1 , 2 , ), where = | H ( ω , τ ) | . Therefore, the transversality condition holds, and thus Hopf-bifurcation occurs at τ = τ j ± . Furthermore, we have the following conclusion.
Lemma 5.
Assume (10), (11), (20) and (18) hold.
(i) For f ( z * ) < 0 , there exists a integer m such that
τ 0 + < τ 0 < τ 1 + < τ 1 < < τ m 1 < τ m + < τ m + 1 + < τ m + 1 .
Equation (13) has at least a pair of positive real parts when τ j = 0 m ( τ j + , τ j ) ;
(ii) All roots of Equation (13) have negative real parts when τ j = 1 m ( τ j 1 , τ j + ) ; Equation (13) has at least a pair of roots with positive real parts when τ > τ m + .
(iii) when τ = τ j + and τ = τ j ( j = 0 , 1 , 2 , ) , all roots of Equation (13) have negative real parts except the purely imaginary roots ± i ω + and ± i ω .
From (24) and (25), it is easy to verify that
τ 0 + < τ 0 , and τ j + τ j 1 + < τ j τ j 1 .
Therefore, together with the define of τ j ± and ω ± , means that the lemma is true.
Now we have found the stability switches. In another words, the stability of the equilibrium switching from stability to instability when τ passes through τ j + ( j = 0 , 1 , 2 , , m ) and back to stability when τ passes through τ j ( j = 0 , 1 , 2 , , m 1 ) . When τ > τ m + , the equilibrium is instability for ever.

3. Direction and Stability of Hopf Bifurcations

In Section 2, some conditions which guarantee system (5) and its modification (6) undergo Hopf bifurcation at some critical values of τ are obtained. We now apply the normal form and center manifold theory introduced by Hassard et al. [20] to study the direction of these Hopf bifurcations and stability of the bifurcated periodic solutions.
We rescale the time by t t τ to normalize the delay. Let τ = τ * + μ , μ R and τ * { τ j + } { τ j } , then μ = 0 is the Hopf bifurcation value and system (6) is equivalent to the following functional differential equation in C = C ( [ 1 , 0 ] , R 3 )
x ˙ ( t ) = L μ ( x t ) + F ( μ , x t ) .
where x t ( θ ) = x ( t + θ ) C , and L μ : C R 3 , F : R × C R 3 are given respectively, by
L μ ( ϕ ) = ( τ * + μ ) a 11 a 12 0 0 0 a 23 a 31 0 a 33 ϕ 1 ( 0 ) ϕ 2 ( 0 ) ϕ 3 ( 0 ) + ( τ * + μ ) M 0 0 0 0 0 0 0 0 ϕ 1 ( 1 ) ϕ 2 ( 1 ) ϕ 3 ( 1 )
and
F ( μ , ϕ ) = ( τ * + μ ) 4 ( 1 b γ ) ϕ 1 2 ( 0 ) 8 ϕ 1 ( 0 ) ϕ 1 ( 1 ) K α ε ϕ 1 ( 0 ) ϕ 2 ( 0 ) + h . o . t K β ε ϕ 2 ( 0 ) ϕ 3 ( 0 ) 0 ,
where a 11 = 4 b γ ( 1 b γ ) , a 12 = K α b γ ε , a 13 = 4 γ ( 1 b γ ) K α , a 31 = a ε , a 33 = a ε , and M = 4 b 2 γ 2 .
By the Reisz representation theorem, there exist a function η ( θ , μ ) of bounded variation for θ [ 1 , 0 ] such that
L μ ϕ = 1 0 d η ( θ , μ ) ϕ ( θ ) f o r ϕ C .
In fact, we can choose
η ( θ , μ ) = ( τ * + μ ) a 11 a 12 0 0 0 a 23 a 31 0 a 33 δ ( θ ) ( τ * + μ ) M 0 0 0 0 0 0 0 0 δ ( θ + 1 ) ,
where
δ ( θ ) = 0 , θ 0 , 1 , θ = 0 .
For ϕ C 1 ( [ 1 , 0 ] , R 3 ) , define
A ( μ ) ϕ = d ϕ d θ , θ [ 1 , 0 ) , 1 0 d η ( s , μ ) ϕ ( s ) , θ = 0 ,
and
R ( μ ) ϕ = 0 , θ [ 1 , 0 ) , F ( μ , ϕ ) , θ = 0 .
Then (28) can be rewritten as
x ˙ ( t ) = A ( μ ) x t + R ( μ ) x t .
For ψ C 1 ( [ 0 , 1 ] , ( R 3 ) * ) , define
A * ψ ( s ) = d ψ d s , s ( 0 , 1 ] , 1 0 d η T ( t , 0 ) ψ ( t ) , s = 0 .
and a bilinear inner product
< ψ ( s ) , ϕ ( θ ) > = ψ ¯ ( 0 ) ϕ ( 0 ) 1 0 ξ θ ψ ¯ ( ξ θ ) d η ( θ ) ϕ ( ξ ) d ξ ,
where η ( θ ) = η ( θ , 0 ) . Then A * and A ( 0 ) are adjoint operators, and ± i τ * ω * are eigenvalues of A ( 0 ) . Of course, they are also eigenvalues of A * . Suppose that q ( θ ) = q ( 0 ) e i ω * τ * θ is an eigenvector of A ( 0 ) corresponding to i ω * τ * , where q ( 0 ) = ( q 1 , q 2 , q 3 ) T . From (31)–(33), we get
τ * i ω * a 11 M e i ω * τ * a 12 0 0 i ω * a 23 a 31 0 i ω * a 33 q 1 q 2 q 3 = 0 0 0 .
Direct computation leads to
q 1 = 1 , q 2 = i ω * a 11 M e i ω * τ * a 12 , q 3 = i ω * a 33 a 31 .
It is not difficult to verify that q * = D ( q 1 * , q 2 * , q 3 * ) e i ω * τ * θ is an eigenvector of A * corresponding to i ω * τ * such that < q * , q > = 1 , < q * , q ¯ > = 0 , where
q 1 * = 1 , q 2 * = a 12 i ω * , q 3 * = i ω * a 11 M e i ω * τ * a 31 , and D = [ 1 + q ¯ 2 * q 2 + q ¯ 3 * q 3 + τ M e i ω * τ * ] 1 .
Using the same notations as in literatures [20,21,22,23], we shall compute the coordinates to describe the center manifold C 0 at μ = 0 .
Let x t be the solution of (28) when μ = 0 and define
z ( t ) = < q * , x t > , W ( t , θ ) = x t z ( t ) q ( θ ) z ¯ ( t ) q ¯ ( θ ) = x t ( θ ) 2 R e { z ( t ) q ( θ ) } .
On the center manifold C 0 , we have
W ( t , θ ) = W ( z ( t ) , z ¯ ( t ) , θ ) ,
where
W ( z ( t ) , z ¯ ( t ) , θ ) = W 20 ( θ ) z 2 2 + W 11 ( θ ) z z ¯ + W 02 ( θ ) z ¯ 2 2 + ,
z ( t ) and z ¯ ( t ) are local coordinates of center manifold C 0 in the direction of the q and q * , respectively. According to the center manifold theory, we have
W ( 0 , 0 , θ ) = 0 , W ( 0 , 0 , θ ) = 0 .
For solution x t in C 0 of Equation (28),
< ψ , A ϕ > = < A * ψ , ϕ > , for ( ϕ , ψ ) D ( A ) × D ( A * ) since μ = 0 ,
so we have
z ˙ ( t ) = < q * , x ˙ ( t ) > = < q * , L μ ( x t ) + F ( μ , x t ) > = < q * , i ω * τ * x t > + < q * , F ( 0 , W + 2 R e { ( z ( t ) q ( θ ) } ) > = i ω * τ * z t + q ¯ * ( 0 ) F ( 0 , W ( z , z ¯ , 0 ) ) + 2 R e { ( z ( t ) q ( 0 ) } ) = def i ω * τ * z t + q ¯ * ( 0 ) F 0 ( z , z ¯ ) .
Consider the formal Taylor expansion of the F 0 ( z , z ¯ ) as follows
F 0 ( z , z ¯ ) = F z 2 z 2 2 + F z ¯ z ¯ 2 2 + F z z ¯ z z ¯ + F z 2 z ¯ z 2 z ¯ 2 +
We rewrite (40) as
z ˙ ( t ) = i ω * τ * z ( t ) + g ( z , z ¯ ) ,
where
g ( z , z ¯ ) = q ¯ * ( 0 ) F 0 ( z , z ¯ ) = g 20 z 2 2 + g 11 z z ¯ + g 02 z ¯ 2 2 + g 21 z 2 z ¯ 2 +
According to (30), we have
g ( z , z ¯ ) = D ¯ τ * ( 1 , q 2 * ¯ , q 3 * ¯ ) l 11 ϕ 1 2 ( 0 ) + l 12 ϕ 1 ( 0 ) ϕ 2 ( 0 ) + l 13 ϕ 1 ( 0 ) ϕ 1 ( 1 ) + l 14 ϕ 1 2 ( 0 ) ϕ 1 ( 1 ) l 21 ϕ 2 ( 0 ) ϕ 3 ( 0 ) 0 = D ¯ τ * { ( l 11 + l 12 q 2 + l 13 e i ω * τ * + l 21 q 2 q 3 q 2 * ¯ ) z 2 + [ 2 l 11 + l 12 ( q 2 + q ¯ 2 ) + l 21 q ¯ 2 * ( q 2 q ¯ 3 + q ¯ 2 q 3 ) ] z z ¯ + ( l 11 + l 12 q ¯ 2 + l 13 e i ω * τ * + l 21 q ¯ 2 * q ¯ 2 q ¯ 3 ) z ¯ 2 + [ l 11 ( 2 W 11 ( 1 ) ( 0 ) + W 20 ( 1 ) ( 0 ) ) + l 12 ( W 11 ( 2 ) ( 0 ) + W 20 ( 1 ) ( 0 ) q 12 2 + W 20 ( 2 ) ( 0 ) 2 + W 11 ( 1 ) ( 0 ) q 12 ) + l 13 ( W 11 ( 1 ) ( 1 ) + W 20 ( 1 ) ( 1 ) 2 + W 20 ( 1 ) ( 0 ) e i ω * τ * + W 11 ( 1 ) ( 0 ) e i ω * τ * ) + l 14 ( e i ω * τ * + 2 e i ω * τ * ) + q ¯ 2 * l 21 ( q 2 W 11 ( 3 ) ( 0 ) + q ¯ 2 W 20 ( 3 ) ( 0 ) + q ¯ 3 W 20 ( 2 ) ( 0 ) + q 3 W 11 ( 2 ) ( 0 ) ) ] z 2 z ¯ } .
Comparing the coefficients with (43), we have
g 20 = 2 D ¯ τ * ( l 11 + l 12 q 2 + l 13 e i ω * τ * + l 21 q 2 q 3 q 2 * ¯ ) , g 11 = D ¯ τ * ( 2 l 11 + l 12 ( q 2 + q ¯ 2 ) + l 13 ( e i ω * τ * + e i ω * τ * ) + l 21 q ¯ 2 * ( q 2 q ¯ 3 + q ¯ 2 q 3 ) ) , g 02 = 2 D ¯ τ * ( l 11 + l 12 q ¯ 2 + l 13 e i ω * τ * + l 21 q ¯ 2 * q ¯ 2 q ¯ 3 ) , g 21 = 2 D ¯ τ * [ l 11 ( 2 W 11 ( 1 ) ( 0 ) + W 20 ( 1 ) ( 0 ) ) + l 12 ( W 11 ( 2 ) ( 0 ) + W 20 ( 1 ) ( 0 ) q 12 2 + W 20 ( 2 ) ( 0 ) 2 + W 11 ( 1 ) ( 0 ) q 12 ) + l 13 ( W 11 ( 1 ) ( 1 ) + W 20 ( 1 ) ( 1 ) 2 + W 20 ( 1 ) ( 0 ) e i ω * τ * + W 11 ( 1 ) ( 0 ) e i ω * τ * ) + l 14 ( e i ω * τ * + 2 e i ω * τ * ) + q ¯ 2 * l 21 ( q 2 W 11 ( 3 ) ( 0 ) + q ¯ 2 W 20 ( 3 ) ( 0 ) + q ¯ 3 W 20 ( 2 ) ( 0 ) + q 3 W 11 ( 2 ) ( 0 ) ) ] .
In order to determine g 21 , we focus on the computation of W 20 ( θ ) and W 11 ( θ ) . By (35), we have
W ˙ = x ˙ t z ˙ ( t ) q ( θ ) z ¯ ˙ ( t ) q ¯ ( θ ) = A W g q ( θ ) g ¯ q ¯ ( θ ) , θ [ 1 , 0 ) , A W g q ( θ ) g ¯ q ¯ ( θ ) + F ( μ , ϕ ) , θ = 0 .
Differentiating (39) with respect to t, we obtain
W ˙ = W z z ˙ + W z ¯ z ¯ ˙ = [ W 20 ( θ ) z + W 11 ( θ ) z ¯ ] z ˙ + [ W 11 ( θ ) z + W 02 ( θ ) z ¯ ] z ¯ ˙ + = [ W 20 ( θ ) z + W 11 ( θ ) z ¯ ] ( i ω * τ * z + g ( z , z ¯ ) ) + [ W 11 ( θ ) z + W 02 ( θ ) z ¯ ] ( i ω * τ * z ¯ + g ¯ ( z , z ¯ ) ) +
Substituting (43) and (46) into (45), and comparing the coefficients, we have
( 2 i ω * τ * A ) W 20 ( θ ) = g 20 q ( θ ) g ¯ 02 q ¯ ( θ ) , θ [ 1 , 0 ) , g 20 q ( 0 ) g ¯ 02 q ¯ ( 0 ) + F z 2 , θ = 0 .
and
A W 11 ( θ ) = g 11 q ( θ ) g ¯ 11 q ¯ ( θ ) , θ [ 1 , 0 ) , g 11 q ( 0 ) g ¯ 11 q ¯ ( 0 ) + F z z ¯ , θ = 0 .
According to (33) and (47), for θ [ 1 , 0 ) ,
W 20 ( θ ) = 2 i ω * τ * W 20 ( θ ) + g 20 q ( θ ) + g ¯ 02 q ¯ ( θ ) .
It is clear that the solution of above equation is
W 20 ( θ ) = i g 20 ω * τ * q ( 0 ) e i ω * τ * θ + i g ¯ 02 3 ω * τ * q ¯ ( 0 ) e i ω * τ * θ + E 1 e 2 i ω * θ .
For θ = 0 , we see from (47) that
1 0 d θ η ( 0 , θ ) W 20 = 2 i ω * τ * W 20 ( 0 ) + g 20 q ( 0 ) + g ¯ 02 q ¯ ( 0 ) F z 2 .
Substituting (49) into (50), we have
( 2 i ω * τ * I 1 0 d θ η ( θ ) e 2 i ω * τ * θ E 1 = 2 τ * l 11 + l 12 q 12 + l 13 e i ω * τ * l 21 q ¯ 2 * q 2 q 3 0 ,
namely,
2 i ω * a 11 M e 2 i ω * τ * θ a 12 0 0 2 i ω * a 13 a 31 0 2 i ω * a 33 E 1 = 2 l 11 + l 12 q 12 + l 13 e i ω * τ * l 21 q ¯ 2 * q 2 q 3 0 .
So we have
E 1 ( 1 ) = 2 M 1 * l 11 + l 12 q 2 + l 13 e i ω * τ * a 12 0 l 21 q ¯ 2 * q 2 q 3 2 i ω * a 13 0 0 2 i ω * a 33 ,
E 2 ( 1 ) = 2 M 1 * 2 i ω * a 11 M e 2 i ω * τ * θ l 11 + l 12 q 2 + l 13 e i ω * τ * 0 0 l 21 q ¯ 2 * q 2 q 3 a 13 a 31 0 2 i ω * a 33 ,
E 3 ( 1 ) = 2 M 1 * 2 i ω * a 11 M e 2 i ω * τ * θ a 12 l 11 + l 12 q 2 + l 13 e i ω * τ * 0 2 i ω * l 21 q ¯ 2 * q 2 q 3 a 31 0 0 ,
where
M 1 * = 2 i ω * a 11 M e 2 i ω * τ * θ a 12 0 0 2 i ω * a 13 a 31 0 2 i ω * a 33 .
Similarly, we have
E 1 ( 2 ) = 2 M 2 * 2 l 11 + l 12 ( q 2 + q ¯ 2 ) + l 13 ( e i ω * τ * + e i ω * τ * ) a 12 0 l 21 q ¯ 2 * ( q 2 q ¯ 3 + q ¯ 2 q 3 ) 0 a 13 0 0 a 33 ,
E 2 ( 2 ) = 2 M 2 * a 11 M 2 l 11 + l 12 ( q 2 + q ¯ 2 ) + l 13 ( e i ω * τ * + e i ω * τ * ) 0 0 l 21 q ¯ 2 * ( q 2 q ¯ 3 + q ¯ 2 q 3 ) a 13 0 0 a 33 ,
E 3 ( 2 ) = 2 M 2 * a 11 M a 12 2 l 11 + l 12 ( q 2 + q ¯ 2 ) + l 13 ( e i ω * τ * + e i ω * τ * ) 0 0 l 21 q ¯ 2 * ( q 2 q ¯ 3 + q ¯ 2 q 3 ) a 31 0 0 ,
where
M 2 * = a 11 M a 12 0 a 31 0 a 13 0 0 a 33 .
Thus, g 12 can be determined by the parameters and delay. So the following quantities can be derived:
C 1 ( 0 ) = i 2 ω * τ * ( g 20 g 11 2 | g 11 | 2 1 3 | g 02 | 2 ) + g 12 2 .
μ 2 = R e C 1 ( 0 ) R e λ ( τ * ) , β 2 = 2 R e C 1 ( 0 ) .
T 2 = I m ( C 1 ( 0 ) ) + μ 2 I m λ ( τ * ) τ * ω * .
It is well-known from [20] that the sign of μ 2 decides the direction of Hopf bifurcation, the sign of β 2 determines the stability of bifurcating periodic solutions and T 2 defines the period of bifurcating periodic solutions. Thus we have the following conclusion.
Theorem 1.
(1) The Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exist for τ > τ * ( τ < τ * ) if μ 2 > 0 ( μ 2 < 0 ) ;
(2) The bifurcating periodic solutions are stable (unstable) if β 2 < 0 ( β 2 > 0 ) ;
(3) The periodic of the bifurcating periodic solutions increase (decrease) if T 2 > 0 ( T 2 < 0 ) .

4. Numerical Simulations

In this section, we perform some numerical simulations to verify our theoretical analysis proved in previous sections. To investigate the effect of delay, we focus on the case when ( γ , b ) D 7 . Set a = 0.2 , ε = 0.15 , γ = 0.55 , b = 1.6 , K = 12.5 , β = 0.05 and α = 0.1 . By direct computation, we get τ 0 + 0.4414 and the equilibrium E is ( 0.88 , 0.0507 , 0.88 ) .
(i) Take τ = 0.4 . Starting from the initial value ( x 1 ( 0 ) , x 2 ( 0 ) , x 3 ( 0 ) ) = ( 0.2 , 0.05 , 0.2 ) , we get Figure 2. According the previous analysis, the equilibrium should be stable. The simulation results coincide exactly with the our conclusion.
(ii) Take τ = 0.45 . Starting from the initial value ( x 1 ( 0 ) , x 2 ( 0 ) , x 3 ( 0 ) ) = ( 0.8 , 0.05 , 0.8 ) , we get Figure 3. This implies that the Hopf bifurcation associated with the critical value τ 0 + 0.4414 is supercritical, the bifurcating periodic solution is stable, and the equilibrium E is unstable. All these results are consistent with our analysis.
(iii) Take τ = 0.55 . Starting from the initial value ( x 1 ( 0 ) , x 2 ( 0 ) , x 3 ( 0 ) ) = ( 0.8 , 0.05 , 0.8 ) , we get Figure 4. In this case, the bifurcating periodic solution vanishes, but the equilibrium is still unstable.

5. Conclusions

In this paper, we have studied the dynamics of a predator-prey system with the weak Allee effect. Firstly, we investigated the effect of the time delay τ on the stability of the positive equilibrium of system (5). Then, by using the normal form and center manifold theorems, we explored the existence of the stability switches, the Hopf bifurcation, the bifurcating direction and the stability of the bifurcating periodic solutions when the characteristic equation of system (6) has two pairs of purely imaginary roots. Meanwhile, we found the characteristic equation of system (6) has one pair of purely imaginary roots which are double roots when f ( z ) = 0 has one positive root. We believe system (6) perhaps has more interesting dynamics in this case which is more complicated than our present work and will be considered in the future. Furthermore, we can incorporate other time delays, such as predator maturation time, etc., into the mathematical model to investigate their dynamics by other methods, from which we may be able to obtain more interesting results.

Author Contributions

Conceptualization, J.Z. and L.Z.; methodology, J.Z.; validation, J.Z., L.Z. and Y.B.; investigation, J.Z.; resources, L.Z.; data curation, Y.B.; writing—original draft preparation, L.Z.; writing—review and editing, Y.B.

Funding

This research was funded by the National Nature Science Foundation of China (No. 11672270, No. 11501511), Zhejiang Provincial Natural Science Foundation of China under Grant No. LY15A010021 and Science Foundation of Zhejiang Sci-Tech University (ZSTU) under Grant No. 13062176-Y.

Acknowledgments

The authors are grateful to the reviewers for the careful reading of the paper with several productive suggestions and their valuable constructive comments, which certainly improved the quality of the original paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ferreira, J.D.; Salazarb, C.A.T.; Tabares, P.C.C. Weak Allee effect in a predator-prey model involving memory with a hump. Nonlinear Anal.-Real World Appl. 2013, 14, 536–548. [Google Scholar] [CrossRef]
  2. Yang, K.; Miao, Z.; Chen, F.; Xie, X. Influence of single feedback control variable on an autonomous Holling-II type cooperative system. J. Math. Anal. Appl. 2016, 435, 874–888. [Google Scholar] [CrossRef]
  3. Lin, Q. Dynamic behaviors of a commensal symbiosis model with non-monotonic functional response and non-selective harvesting in a partial closure. Commun. Math. Biol. Neurosci. 2018, 2018, 4. [Google Scholar]
  4. Chen, F.; Wu, H.; Xie, X. Global attractivity of a discrete cooperative system incorporating harvesting. Adv. Differ. Equ. 2016, 2016, 268. [Google Scholar] [CrossRef]
  5. Manna, D.; Maiti, A.; Samanta, G.P. A Michaelis-Menten type food chain model with strong Allee effect on the prey. Appl. Math. Comput. 2017, 311, 390–409. [Google Scholar] [CrossRef]
  6. Li, Y.N.; Li, Y.; Liu, Y.; Cheng, H. Stability Analysis and Control Optimization of a Prey-Predator Model with Linear Feedback Control. Discrete Dyn. Nat. Soc. 2018, 2018, 4945728. [Google Scholar] [CrossRef]
  7. Bai, Y.; Mu, X. Global asymptotic stability of a generalized SIRS epidemic model with transfer from infectious to susceptible. J. Appl. Anal. Comput. 2018, 8, 402–412. [Google Scholar]
  8. Han, M.; Zhang, L.; Wang, Y.; Khalique, C.M. The effects of the singular lines on the traveling wave solutions of modified dispersive water wave equations. Nonlinear Anal.-Real World Appl. 2019, 47, 236–250. [Google Scholar] [CrossRef]
  9. Ding, H.; Lim, C.W.; Chen, L.Q. Nonlinear vibration of a traveling belt with non-homogeneous boundaries. J. Sound Vib. 2018, 424, 78–93. [Google Scholar] [CrossRef]
  10. Zhang, L.; Wang, Y.; Khlique, C.M.; Bai, Y. Peakon and cuspon solutions of a generalized Camassa-Holm-Novikov equation. J. Appl. Anal. Comput. 2018, 8, 1938–1958. [Google Scholar]
  11. Li, T.; Lin, Q.; Chen, J. Positive periodic solution of a discrete commensal symbiosis model with Holling II functional response. Commun. Math. Biol. Neurosci. 2016, 2016, 22. [Google Scholar]
  12. Li, Y.; Cheng, H.; Wang, J.; Wang, Y. Dynamic analysis of unilateral diffusion Gompertz model with impulsive control strategy. Adv. Differ. Equ. 2018, 1, 32. [Google Scholar] [CrossRef]
  13. Bai, Y.; Li, Y. Stability and Hopf bifurcation for a stage-structured predator-prey model incorporating refuge for prey and additional food for predator. Adv. Differ. Equ. 2019, 2019, 42. [Google Scholar] [CrossRef]
  14. Shi, Z.; Wang, J.; Li, Q.; Cheng, H. Control optimization and homoclinic bifurcation of a prey-predator model with ratio-dependent. Adv. Differ. Equ. 2019, 2019, 2. [Google Scholar] [CrossRef]
  15. Huang, J.; Gong, Y.; Ruan, S. Bifurcation Analysis in a Predator-prey Model with Constant-yield Predator Harvesting. Disc. Cont. Dyn. Syst.-Ser. B 2013, 18, 2101–2121. [Google Scholar] [CrossRef]
  16. Boukal, D.; Berec, L. Single-species models of the Allee effect: Extinction boundaries, sex ratios and mate enconunters. J. Theor. Biol. 2002, 218, 375–394. [Google Scholar] [CrossRef]
  17. Tabares, P.C.C.; Ferreira, J.D.; Rao, V.S.H. Weak Allee effect in a predator-prey system involving distributed delays. Comput. Appl. Math. 2011, 30, 675–699. [Google Scholar] [CrossRef] [Green Version]
  18. Freedman, H.I.; Wolkowicz, G.S.K. Predator-Prey Systems with Group Defence: The Paradox of Enrichment Revisited. Bull. Math. Biol. 1986, 48, 493–508. [Google Scholar] [CrossRef]
  19. Ruan, W.; Wei, J. On the zeros of transcendental functions with applications to stability of delay differential equations with two delays. Dyn. Cont. Discret. Impuls. Syst. Ser. A Math. Anal. 2003, 10, 863–874. [Google Scholar]
  20. Hassard, B.D.; Kazarinoff, N.D.; Wan, Y.H. Theory and Applications of Hopf Bifurcation; Cambridge University Press: Cambridge, UK, 1981. [Google Scholar]
  21. Wu, S.; Song, Y. Stability and spatiotemporal dynamics in a diffusive predator-prey model with nonlocal prey competition. Nonlinear Anal.-Real World Appl. 2019, 48, 12–39. [Google Scholar] [CrossRef]
  22. Song, Y.; Jiang, H.; Liu, Q.; Yuan, Y. Spatiotemporal dynamics of the diffusive Mussel-Algae model near Turing-Hopf bifurcation. SIAM J. Appl. Dyn. Syst. 2017, 16, 2030–2062. [Google Scholar] [CrossRef]
  23. Zhang, J.; Zhang, L.; Khalique, C.M. Stability and Hopf Bifurcation Analysis on a Bazykin Model with Delay. Abstr. Appl. Anal. 2014, 11, 1–7. [Google Scholar] [CrossRef]
Figure 1. Bifurcation diagram of positive equilibrium in γ –b plane for a = 0.2 , ε = 0.15 .
Figure 1. Bifurcation diagram of positive equilibrium in γ –b plane for a = 0.2 , ε = 0.15 .
Mathematics 07 00432 g001
Figure 2. Behavior of system (6) with τ = 0.4 . The positive equilibrium E is asymptotically stable when τ < τ 0 + 0.4414 .
Figure 2. Behavior of system (6) with τ = 0.4 . The positive equilibrium E is asymptotically stable when τ < τ 0 + 0.4414 .
Mathematics 07 00432 g002
Figure 3. Behavior of system (6) with τ = 0.45 . The positive equilibrium E is unstable. Furthermore, a periodic solution bifurcates from the positive equilibrium and is orbitally asymptotically stable.
Figure 3. Behavior of system (6) with τ = 0.45 . The positive equilibrium E is unstable. Furthermore, a periodic solution bifurcates from the positive equilibrium and is orbitally asymptotically stable.
Mathematics 07 00432 g003aMathematics 07 00432 g003b
Figure 4. Behavior of system (6) with τ = 0.55 . The positive equilibrium E is unstable when τ > τ 0 + 0.4414 .
Figure 4. Behavior of system (6) with τ = 0.55 . The positive equilibrium E is unstable when τ > τ 0 + 0.4414 .
Mathematics 07 00432 g004

Share and Cite

MDPI and ACS Style

Zhang, J.; Zhang, L.; Bai, Y. Stability and Bifurcation Analysis on a Predator–Prey System with the Weak Allee Effect. Mathematics 2019, 7, 432. https://doi.org/10.3390/math7050432

AMA Style

Zhang J, Zhang L, Bai Y. Stability and Bifurcation Analysis on a Predator–Prey System with the Weak Allee Effect. Mathematics. 2019; 7(5):432. https://doi.org/10.3390/math7050432

Chicago/Turabian Style

Zhang, Jianming, Lijun Zhang, and Yuzhen Bai. 2019. "Stability and Bifurcation Analysis on a Predator–Prey System with the Weak Allee Effect" Mathematics 7, no. 5: 432. https://doi.org/10.3390/math7050432

APA Style

Zhang, J., Zhang, L., & Bai, Y. (2019). Stability and Bifurcation Analysis on a Predator–Prey System with the Weak Allee Effect. Mathematics, 7(5), 432. https://doi.org/10.3390/math7050432

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