A singular prey-predator model with time delays is formulated and analyzed. Allee effect is considered on the growth of the prey population. The singular prey-predator model is transformed into its normal form by using differential-algebraic system theory. We study its dynamics in terms of local analysis and Hopf bifurcation. The existence of periodic solutions via Hopf bifurcation with respect to two delays is established. In particular, we study the direction of Hopf bifurcation and the stability of bifurcated periodic solutions by applying the normal form theory and the center manifold argument. Finally, numerical simulations are included supporting the theoretical analysis and displaying the complex dynamical behavior of the model outside the domain of stability.

1. Introduction

The economic theory proposed by Gordon [1] in 1954 was described as follows:net economic revenue (NER) = total revenue (TR) − total cost (TC).

This provides theoretical evidence for the establishment of singular bioeconomic model, which is described by differential-algebraic equations. Recently, there has been extensive literature dealing with such systems, describing the interactions between the different species and harvesting effort regarding activity, stabilities of equilibrium, bifurcations, and other dynamics (see, e.g., [25] and the references therein).

In most of ecosystems, since one species does not respond instantaneously to interactions with other species, some delays due to several reasons, such as gestation, hunting, and maturation, are required. To incorporate this idea in a modelling approach, time delays have been introduced into ecosystems. For a long time, it has been recognized that delays can have very complicated impact on the dynamics of a system (see, e.g., monographs by Hale [6], Hale and Verduyn Lunel [7], Kuang [8], Yuan and Song [9], and Wu [10]). Generally speaking, delays can cause the loss of stability and lead to periodic solutions. Some authors like in [11] described the effects of time delay in a prey-predator model incorporating parasite infection for the prey population. The literature [12] described the dynamics of a stage structured population model with time delay in fluctuating environment. The literature [13] investigated dynamics of delayed prey-predator model with harvesting. Chakraborty et al. [14] introduced a single discrete gestation delay in a differential-algebraic biological economic system and showed Hopf bifurcation in the neighborhood of coexisting equilibrium point through considering the delay as a bifurcation parameter. Zhang et al. [15] studied a ratio-dependent prey-predator singular model and analyzed the direction and stability of periodic solutions. Some other authors have discussed dynamical models with multiple delays. While analyzing, the multiply delayed models are mostly simplified by taking equal magnitude for two delays [16, 17] or choosing the sum of two delays as a bifurcation parameter [18, 19]. However, the delays appearing in different terms of an ecological system are not always equal. Therefore, it is necessary to discuss dynamical system with different delays.

In addition, numerous examples demonstrate that the growth of natural populations can exhibit Allee effect [20], which describes a positive relation between population density and the per capita growth rate. When populations get larger, the effect will saturate or disappear. Allee effect may result from some causes, such as mate finding, social dysfunction, inbreeding depression, fool exploitation, and predator avoidance of defense. Therefore, in this paper, we consider a differential-algebraic prey-predator model with two time delays and the Allee effect as follows: where is the Allee effect constant of the prey species, respectively. and are the intrinsic mortality rate of the prey and predator population, respectively. is the handling time. denotes food utilization efficiency. Prey dynamics is delayed by due to crowing and the predator takes time to convert the food into its growth; is harvesting effort for predator, and , , and are harvesting reward per unit harvesting effort for unit weight of predator, harvesting cost per unit harvesting effort for predator, and the net economic revenue per unit harvesting effort, respectively. All the parameters are positive constants.

In this paper, we investigate the effects of two delays on the dynamics of singular bioeconomic model with Allee effect. The existence of periodic solution has been explored through Hopf bifurcation. The domain of stability is defined. Then, the formula for determining the properties of Hopf bifurcation is derived by using the normal form method and center manifold theorem. Finally, numerical simulations show the bifurcation plot with respect to time delay and give the effectiveness of the result mentioned above.

2. Stability and Existence of Hopf Bifurcation

There exists a positive equilibrium for model (1), where coordinates , , and satisfy the following equations:

For the sake of the simplicity, denote the two differential equations by and and the algebraic equation by .

Let , , and . Therefore, the system (1) can be rewritten in the following form: Consider the following local parametric of the third equation of the system (3): where , , , , and is a smooth mapping and has the following expression: Then, the linearization part of the system (3) is For the sake of the simplicity, let

According to the Jacobian matrix , the characteristic equation of (1) at can be written as follows:

When the time delays are not considered, the characteristic polynomial for model (1) can be expressed as the following form: It is clear that if the following conditions are satisfied: then is locally asymptotically stable for model (1).

Next, we consider the effects of the two time delays on the stability of model (1) in the following cases.

Case 1 ( and ). Due to , (8) becomes

Assume that a purely imaginary solution of the form exists in (11). Substituting it into (11) and separating the real and imaginary parts, we have Taking square on both sides of (12) and summing them up, we obtain From condition (10), (13) has unique positive real root if The critical value of the delay corresponding to is given by If , (13) has no real root and model (1) is asymptotically stable for any time delay .

Now, differentiating (8) with respect to and substituting the eigenvalue , it follows that

Theorem 1. Assume that condition (14) holds; then, there exists a given by such that the positive equilibrium of model (1) is locally asymptotically stable for and undergoes a Hopf bifurcation when . That is, model (1) has a branch of periodic solutions bifurcating from the positive equilibrium near .

Case 2 ( and ). According to a similar discussion as in Case 1, we can obtain the following result.

Theorem 2. There exists a such that the positive equilibrium of model (1) is locally asymptotically stable for and undergoes a Hopf bifurcation when given by where is the root of corresponding characteristic equation.

Case 3 ( and ). This case states that is regarded as a parameter and is in its stable interval. Assume that (8) has purely imaginary solution of the form . Substituting it into (8) and separating the real and imaginary parts, we have Eliminating leads to It can be seen that there exists at least one real positive root for (20) if the condition holds. Equation (19) can also be written as Equation (21) can be solved as where and .

Now, differentiating (8) with respect to and substituting the eigenvalue and time delay , it follows that where

Solving (23), we obtain

Theorem 3. Assume and . The positive equilibrium of model (1) is asymptotically stable for and undergoes Hopf bifurcation at , where where and .

Case 4 ( and ). According to a similar discussion as in Case 3, we can obtain the following result.

Theorem 4. The positive equilibrium of model (1) is asymptotically stable for and undergoes Hopf bifurcation at , where where .

3. Direction and the Stability of Hopf Bifurcation

In the previous section, we obtain the conditions under which a family of periodic solutions bifurcate from the positive equilibrium at the critical values of time delays and . As pointed out by Hassard et al., by employing the normal form and center manifold theory, the formulae for determining the directions, stability, and period of Hopf bifurcation can be presented. Following the ideas of Hassard et al., this section discusses the directions, stability, and period of Hopf bifurcation. These properties are studied with respect to for fixed . Throughout this section, it is considered that the system (1) undergoes Hopf bifurcation at , at the equilibrium . Without loss of generality, this section assumes that .

Let , , , , and and dropping the bars for simplification of notations, the system (1) can be written as a functional differential equation in where , and , are given, respectively, by where . By the Riesz representation theorem, there exists a matrix whose components are bounded variation for , such that where we choose For , define Then, system (29) is equivalent to where for .

For , define and bilinear form where , , and are adjoint operators. Since are eigenvalues of , they are also eigenvalues of . Next, we need to compute the eigenvector of and corresponding to and , respectively.

In terms of the discussion mentioned above, we see that are eigenvalues of and . Next, we calculate the eigenvector of corresponding to and the eigenvector of belonging to , respectively. By the definition of and , it is not difficult to obtain that Similarly, let be the eigenvector of corresponding to , and we can obtain It follows from (36) that Then, we obtain such that , .

Based on the ideas in Hassard et al. [21], we can obtain and the following expressions: where those important parameters can be computed as follows: According to the computation process similar to Hassard’s, parameters of can be written as follows:

Therefore, we can compute by the parameters and delays. Then, the following results can be expressed as These parameters give a description of the bifurcating periodic solutions in the center manifold of system (1) at critical values , which can be expressed as follows.(i) determines the direction of the Hopf bifurcation: if , the Hopf bifurcation is supercritical (subcritical).(ii) gives the stability of periodic solution: if , the periodic solution is stable (unstable).(iii) expresses the period of the bifurcating periodic solutions: if , the period of the bifurcating periodic solution increases (decreases).

4. Numerical Simulation

With the help of MATLAB, we present some numerical simulations to verify and extend our theoretical analysis proved in Sections 2 and 3. The parameters of model (1) are chosen as , , , , , , , and .

According to the parameters given above, we can obtain the positive equilibrium point of model (1). It is easy to obtain that is asymptotically stable in the absence of delay (see Figure 1). Moreover, since , there is a positive root for (13). Hence, when , the critical value of the delay corresponding to is , which is shown in Figure 2. When , by a similar computation, we can obtain that the critical delay for is .

Based on the results of Cases 3 and 4 in Section 2, we draw a critical curve with respect to two parameters and (see Figure 3). The main algorithm for drawing the critical curve is searching the values () for fixed positive values of   . From Figure 3, it is clear that the domain surrounded by the critical curve and the two axes, that is, the gray region, is stable for model (1).

In addition, according to the algorithm derived in Section 3, we can obtain the following values: , , , , , , and . By the discussion in Section 3, we know that model (1) can undergo a supercritical Hopf bifurcation at the positive equilibrium and the bifurcating periodic solution occurs when crosses to the right with fixed and the bifurcating periodic solution is stable. Figure 4 plots the bifurcating periodic solution of the model (1) at and .

To explore the possibility of occurrence of chaos, bifurcation diagrams are plotted for the key parameter . Bifurcation diagram in Figure 5 of model (1) in () plane is given for the fixed . We can see that model (1) experiences the processes of periodic, periodic doubling cascade and chaos in the region . When , model (1) exhibits a stable limit cycle. When is increasing, the limit cycle is unstable and there is a cascade of periodic doubling bifurcation (see Figure 6) leading to chaos.

5. Conclusions

In this paper, dynamical complexity of a multiple-delayed predator-prey model with Allee effect is analyzed by using differential-algebraic system theory. The prey-predator model with single delay has been investigated by many researchers (see, e.g., [2225]). However, to our best knowledge, there are few papers on the bifurcations of differential-algebraic population dynamics with two or multiple delays. On the other hand, it should be noted that almost the existing bioeconomic models (see [25, 14, 15]) only investigate the simplest case of a logistic prey growth function. Compared with these works, the introduction of Allee effect makes the work studied in this paper novel.

This paper provides a new and efficient method for the qualitative analysis of the Hopf bifurcation of differential-algebraic biological economic system with multiple delays. From the analysis of the proposed model, we have obtained some interesting and useful results. The main contribution of this paper is that the stability domain with respect to two delays is defined and plotted. By choosing the time delay as a bifurcation parameter, the direction of Hopf bifurcation and stability of the bifurcating periodic orbit are discussed by using normal form and center manifold. Moreover, numerical simulations are carried out which substantiate the theoretical results. At the same time, bifurcation diagrams and attracting sets are plotted, which validate the existence of chaos due to periodic doubling bifurcation. However, those theoretical analyses, such as the existence of chaos and conditions of periodic doubling bifurcation, will be the topics of the future research.

Conflict of Interests

The authors do not have conflict of interests with any commercial identities.


This work was supported by the National Science Foundation of China (61273008), Doctor Startup Fund of Liaoning Province (20131026), Fundamental Research Funds for the Central University (N120405009), and Soft Science Research Project of Hubei Province (2012GDA01309). The authors also gratefully thank those anonymous authors whose work largely constitutes this sample paper.