Journal of Applied Mathematics

Volume 2014 (2014), Article ID 107565, 11 pages

http://dx.doi.org/10.1155/2014/107565

Research Article

Bifurcation Analysis and Control of a Differential-Algebraic Predator-Prey Model with Allee Effect and Time Delay

Xue Zhang and Qing-ling Zhang

Department of Science, Northeastern University, No. 3-11, Wenhua Road, Heping District, Shenyang, Liaoning 110004, China

Received 3 August 2013; Accepted 22 January 2014; Published 3 March 2014

Academic Editor: Junjie Wei

Copyright © 2014 Xue Zhang and Qing-ling Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper studies systematically a differential-algebraic prey-predator model with time delay and Allee effect. It shows that transcritical bifurcation appears when a variation of predator handling time is taken into account. This model also exhibits singular induced bifurcation as the economic revenue increases through zero, which causes impulsive phenomenon. It can be noted that the impulsive phenomenon can be much weaker by strengthening Allee effect in numerical simulation. On the other hand, at a critical value of time delay, the model undergoes a Hopf bifurcation; that is, the increase of time delay destabilizes the model and bifurcates into small amplitude periodic solution. Moreover, a state delayed feedback control method, which can be implemented by adjusting the harvesting effort for biological populations, is proposed to drive the differential-algebraic system to a steady state. Finally, by using Matlab software, numerical simulations illustrate the effectiveness of the results.

1. Introduction

In recent years, the growing human needs for more food and more energy have led to increased exploitation of these resources. The problems related to many fields like fishery, forestry, and wildlife. Therefore, mankind is facing the dual problems of resource shortages and environmental degradation. Concerning the conservation for the long-term benefits of humanity, there is a wide-range of interest in analysis and modelling of bioeconomic systems. In many earlier studies, it has been shown that harvesting has a strong impact on population dynamics, ranging from rapid depletion to complete preservation of biological populations. Two main kinds of harvesting were focused on nonzero constant harvesting [1–3] and constant harvesting effort [4–7]. With constant harvesting, the generalized Gause prey-predator model is found to exhibit saddle-node bifurcations, Hopf bifurcation, heteroclinic bifurcation, and nilpotent saddle bifurcation [1]. Xiao et al. [2] have investigated the dynamical properties of a ratio-dependent predator-prey model with nonzero constant rate predator harvesting. These results reveal far richer dynamics compared to the model with no harvesting. The literature [3] shows that harvesting effort as control parameter is not only possible to control the cyclic behavior of populations leading to the persistence of all species, but other desired stable equilibrium including disease-free can also be obtained. Das et al. [4] discussed the bioeconomic harvesting of a prey-predator fishery in which both species are infected by some toxicant. Ji and Wu [5] studied a predator-prey model with a constant-rate prey harvesting incorporating a constant prey refuge, where the influence of harvesting effort on the density of two species was discussed. Chakraborty et al. [6] describes a prey-predator model with stage structure for predator and selective harvesting effort on predator population. Geometric approach is used to drive the sufficient conditions for global stability of the system, and fishing effort is used to investigate the optimal utilization of the resource. Xiang et al. [7] consider a Lotka-Volterra model with impulsive harvest for the prey and investigate globally attractive periodic solution. However, most of these discussions are only based on differential equations or difference equations.

In power systems, neural networks, and genetic networks [8–12], differential-algebraic equations have been studied widely and a lot of results have been obtained, such as local stability, optimal control, singularity induced bifurcation, and feasibility regions. From 2009, several differential-algebraic biological models were reported [13–16]. Surprisingly enough, all the existing differential-algebraic biological modelling literature considers only the simplest case of a logistic prey growth function. However, numerous examples demonstrate that the growth of natural populations can exhibit Allee effect, which is a phenomenon in biology named after Allee [17]. Allee effect describes a positive relation between population density and the per capita growth rate. In other words, for smaller populations, the reproduction and survival of individuals decrease. This effect usually saturates or disappears as populations get larger. The effect may be due to any number of causes, for example, mate finding, social dysfunction, inbreeding depression, food exploitation, and predator avoidance or defense.

On the other hand, since reproduction of predator after consuming prey is not instantaneous in most cases, some time lag for gestation is required. Therefore, in this paper, we consider a differential-algebraic prey-predator model with time delay and the Allee effect on the growth of the prey population. We analyze the stability properties and bifurcation behavior of this model. A state delayed feedback control method is also proposed, which can eliminate Hopf bifurcation and drive the differential-algebraic prey-predator model to stay at a steady state.

2. Model Equations

The general predator-prey model in its classical form is represented by

where and represent the prey density and predator density at time , respectively; is the per capita growth rate of prey in absence of predator; is the intrinsic mortality rate of predator in the absence of food; is predator’s functional response, defined as the amount of prey catch per predator per unit of time;    is the rate of conversion of nutrients into the reproduction rate of the predator.

In 1954, Gordon [18] studied the effect of harvest effort on ecosystem from an economic perspective and proposed the following economic theory:

Based on this theory, this paper studies a class of delayed differential-algebraic predator-prey model with Allee effect on prey species and Holling-II functional response, which is written in the following form:

where is the fertility rate of prey species. and are the intrinsic mortality rate of the prey and predator species, respectively. is the strength of intracompetition of prey population. is the attack coefficient and is the handling time. denotes food utilization efficiency. The predator takes time to convert the food into its growth; is harvesting effort for predator, , , 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.

Let the fertility rate increase with population density and be described by

where and are the per capita maximum fertility rate and the Allee effect constant of the prey species, respectively. If , the fertility of the species is zero when is zero and approaches to when becomes very large. The increasing of depends on the parameter . The larger is, the stronger Allee effect will be. In particular, the fertility rate is density independent when ; that is, .

When the prey population is subject to Allee effect above, the predator-prey model (3) becomes

Assume that the per capita maximum fertility rate of prey must exceed its death rate, that is, ; otherwise, both prey and predator will become extinct. In addition, the maximum growth rate of predator population must exceed its death rate, that is, . If not, prey population will never be able to sustain predator population.

3. Qualitative Analysis

We nondimensionalize the model (5) with the following scaling:

and then obtain the following form:

where the nondimensional parameters are defined as

From the assumption mentioned above, the following reasonable condition must be satisfied:

For simplicity of computation, we consider the above model (7) instead of the model (5). Hence, we will perform a qualitative analysis of the model (7).

Let

where is the state and denotes bifurcation parameter.

3.1. The Model (7) with Zero Economic Profit

Considering zero economic revenue, the model (7) can be reduced as follows:

By the analysis of roots for the model (11), we obtain the following result.

Theorem 1. (1) The model (11) has a trivial equilibrium point for any positive parameters.

(2) There exist two boundary equilibrium points and if , where

(3) There exists another boundary equilibrium point , where , , , if the following conditions are satisfied:

(4) The model (11) has a positive equilibrium point , where , , and is the root of the following equation:

The Jacobian matrix of the model (11) takes the following form:

Hence, the Jacobian matrix of the model (11) at is

Theorem 2. is always a locally stable node.

From Theorem 2, it can be seen that both prey and predator populations will become extinct when their population densities lie in the attraction region of .

Note that , then the Jacobian matrix of the model (11) at is

Since

and the other eigenvalue of is , we obtain the following results on the stability of .

Theorem 3. Assume that and . Then, for any time delay , one has the following.

(1) is a saddle point if

(2) is an unstable node if

Next, we analyze the stability of the equilibrium point . The Jacobian matrix of the model (11) at is

For , one eigenvalue is , and the other is given by the equation

By simple analysis, the root of (22) is negative if the time delay satisfies , where

Then, we have the following results on the stability of .

Theorem 4. Assume that , , and . Then, one has the following.

(1) is a saddle point if

(2) is locally asymptotically stable if

Based on Theorem 4, the following bifurcation result can be obtained.

Theorem 5. Assume that , , and . Then, the model (11) undergoes transcritical bifurcation at the equilibrium point when bifurcation parameter increases through .

Proof. When and the bifurcation parameter , the characteristic polynomial at the equilibrium point has a simple zero eigenvalue with left null vector

and right null vector

According to the literature [19], the model (11) undergoes transcritical bifurcation at the equilibrium point . This completes the proof.

Remark 6. Transcritical bifurcation implies that the equilibrium point remains stable if the handling time is longer than the critical point , otherwise the stability of is lost. From the view of biological explanation, stable equilibrium point means that predator population is to be extinct, which results from longer handling time reducing the amount of prey catch per predator per unit of time.

The Jacobian matrix of the model (11) at is

The corresponding characteristic equation is

In the absence of delay, since the constant term of the above equation is always positive, the sign of the following expression:

determines the stability of equilibrium . Moveover, if the following condition:

is satisfied, it is clear that is always locally asymptotically stable for the model (11) in the absence of delay.

In the presence of delay, assume that a purely imaginary solution of the form exists for the above characteristic equation, where is the root of the following equation:

Obviously, this equation has a positive solution . Therefore, the system undergoes Hopf bifurcation at the equilibrium point when

Then the following result is obtained.

Theorem 7. Assume that the condition (31) holds. The model (11) is stable for the time delay and undergoes Hopf bifurcation at the equilibrium point when .

In succession, we discuss the bifurcation behavior regarding as the bifurcation parameter; that is, .

Theorem 8. If , the model (11) undergoes singular induced bifurcation at the equilibrium point when the bifurcation parameter increases through . Moreover, the stability of the equilibrium point changes, that is, from stable to unstable.

Proof. We define a new variable as . Then, has a simple zero eigenvalue at the equilibrium point as follows:

On the other hand, we get

According to [20], all conditions of singular induced bifurcation are satisfied. Hence, the model (11) undergoes singular induced bifurcation at the equilibrium point if the economic revenue is zero. When economic revenue increases through , one eigenvalue of the model (11) moves from (the open complex left half plane) to (the open complex right half plane) along the real axis by diverging through , which causes impulsive phenomenon of differential-algebraic system, that is, rapid expansion of the population from the view of biological explanation. Therefore, the stability of the equilibrium point changes, that is, from stable to unstable. This completes the proof.

For the model (11) without time delay, eigenvalues of Jacobian matrix at the equilibrium point are the roots of the following equation:

where

Remark 9. After simple computation, it can be seen that the Jacobian matrix at the equilibrium point has two eigenvalues. One is

and the other is

since as .

3.2. The Model (7) with Positive Economic Profit

When the economic profit is positive, the model (7) has positive equilibrium point , where coordinates , , and are the root of the following equation:

where

The Jacobian matrix of the model (7) at the equilibrium point is

For the sake of simplicity, let

The characteristic polynomial for the model (7) without time delay at the equilibrium point takes the following form:

By using the Routh-Hurwitz criteria, is locally asymptotically stable for the model (7) without time delay if

According to the Jacobian matrix , we can obtain the characteristic equation of the differential-algebraic model (7) at , which can be expressed as follows:

where

Assume that a purely imaginary solution of the form exists in (46). Substituting it into (46) and separating the real and imaginary parts, we have

Taking square on both sides of (48) and summing them up, we obtain

From the condition (45), we have and . Therefore, (49) has at least one real root if

The critical value of the delay corresponding to is given by

Now, differentiating (46) with respect to , we obtain

Substituting the eigenvalue and noticing that the derivative of (49) at is positive, it follows that

Theorem 10. Assume that conditions (45) and (50) hold. For the model (7), there exists a such that is locally asymptotically stable when and is unstable when . Furthermore, the model (7) undergoes Hopf bifurcation at the equilibrium point when , .

From Theorem 10, it can be seen that Hopf bifurcation results in oscillation of population density; that is, the differential-algebraic model (7) becomes unstable when the time delay exceeds the critical value. From the view of ecological managers, it may be desirable to have stable population densities in order to keep sustainable development of ecosystem. However, food conversion delay is generally regarded as an inherence of biological populations. It is difficult to adjust the food conversion delay. Therefore, we propose the following state delayed feedback control method:

where is the feedback gain.

Under the state delayed feedback control, the Jacobian matrix at the positive equilibrium point is

Then the character equation corresponding to the matrix is

For the simplicity of computation, the feedback gain term is denoted as . Then, the character equation (56) is written as

Assume that (57) has a purely imaginary root , and substitute it into the equation, then we obtain the following equation:

By simple analysis, we can obtain the following result.

Theorem 11. Under the conditions of Theorem 10, for any time delay , the equilibrium point of the model (54) is locally asymptotically stable if the feedback gain satisfies the following conditions:

In fact, the two differential equations of the differential-algebraic system (54) can be rewritten as

where and .

Remark 12. From the practical point of view, and can be regarded as harvesting efforts for the prey and predator species. From the expression above, it is clear that and are related to the equilibrium, the past, and present population density. Hence, population densities can be kept stable by changing the strength of harvesting efforts; that is, the state delayed feedback control can be completed by adjusting harvesting efforts for two species. Hence, the state feedback controller can be used to eliminate Hopf bifurcation and drive the model system to stabilize at the interior equilibrium point.

4. Numerical Simulation

In this section, we firstly assign some parameter values of the model (7) and provide some numerical simulations to illustrate the effectiveness of the results which have been established in the previous sections. For the purpose of simulation experiments, we mainly use the software MATLAB 7.0.

The differential-algebraic model (7) with given values takes the following form:

4.1. Existence of Singular Induced Bifurcation

The existence of singularity induced bifurcation of the model (61) is clearly shown in Table 1. From Table 1, it can be seen that as economic revenue increases through zero, one eigenvalue remains almost constant and the other moves from to along the real axis by diverging through for the fixed Allee effect constant. Moreover, for the fixed Allee effect , we draw a figure to illustrate the movement of maximum eigenvalue with respect to the economic revenue (see Figure 1). From Figure 1, it can be seen that singularity induced bifurcation for the model (61) occurs at , which is corresponding to the theoretical analysis. Hence, the stability of the model (61) at positive equilibrium point changes from stable to unstable. Furthermore, the stronger Allee effect is, the weaker impulsive phenomenon will be; that is, Allee effect has an impact on the dynamical behavior of the proposed model.

Table 1: Equilibrium points and eigenvalues of the model (7) for different Allee effect constants and economic revenue.

Allee effect
Economic profit
Equilibrium point
Eigenvalues

(0.957, 0.0284, 1.424)

(0.956, 0.0288, 1.264)
−0.930, −201.571

−0.929, 161.025

(0.874, 0.0285, 2.367)

(0.874, 0.0288, 1.142)
−0.775, −559.217

−0.774, 131.598

(0.442, 0.0265, 0.136)

(0.435, 0.0311, 0.115)
−0.203, −1.684

−0.162, 1.420

Figure 1: The maximum eigenvalue of the model (61) w.r.t. economic revenue .

4.2. Existence and Control of Hopf Bifurcation

For the fixed parameter , we draw a critical curve w.r.t. parameter (see Figure 2). From Figure 2, it is clear that the domain surrounded by the critical curve and the two axes is stable for the model (7). Based on the analysis in Section 3.2, the equilibrium point is stable for the model (7) without time delay. Furthermore, since , there is a positive root for (49). The critical value of the delay corresponding to is . The interior equilibrium point remains stable for . In Figure 3, a random time delay is selected in the interval , which is enough to merit the above mathematical study. When the time delay passes through the critical value , the equilibrium point loses its stability and Hopf bifurcation occurs. The bifurcating period solution from at is depicted in Figure 4. Next, a state feedback controller is applied to the model (7). We choose the feedback gain , where the gain terms satisfy and . Then, the model (7) is stable at and the Hopf bifurcation is also eliminated. Figure 5 shows the dynamical responses of the differential-algebraic model (7).

Figure 2: The critical curve of the delay w.r.t. economic revenue .

Figure 3: When , the equilibrium point is stable.

Figure 4: When , bifurcation period solutions from the equilibrium point occur.

Figure 5: Dynamical responses of the controlled system (54) with the feedback gain .

5. Conclusions

Nowadays, much attention has been paid to preserving biological resources with the aim of stemming the damage and ensuring the balance of ecosystems, which inspires the introduction of harvesting in the biological system. In this paper, we analyze the dynamical behavior of a delayed predator-prey model with Allee effect and harvesting by using differential algebraic systems theory. From the analysis of the proposed model, we have obtained some interesting and useful results, which extend the work done in [21]. This paper is mainly divided into two parts.

In the first part, we consider a delayed differential-algebraic predator-prey model with zero economic revenue. It is observed that transcritical bifurcation and singular induced bifurcation phenomena take place when handling time of predator and economic revenue are regarded as bifurcation parameters, respectively.

As the handling time decreases through the critical point , the differential-algebraic model undergoes transcritical bifurcation and the stability of is lost. Due to longer handling time reducing the amount of prey catch per predator per unit of time, prey population density can stay at a positive value and predator population is to be extinct; that is, is stable.

Singular induced bifurcation may cause impulsive phenomenon due to the variation of the economic revenues of harvesting. From a biological point of view, singular induced bifurcation implies rapid expansion of biological population, which may cause ecosystem unbalance and hamper the sustainable development. Hence, it is necessary to investigate the singular induced bifurcation in the presence of a reserve depending on the variation of economic revenue. The analysis of singular induced bifurcation can provide more information of forecasting so that ecological managers can lay down better management strategy. Furthermore, numerical simulation shows that the impulsive phenomenon can be much weaker by increasing Allee effect constant, which implies that Allee effect has an impact on the dynamical behavior of the proposed model. Therefore, ecological managers need to consider the inherent characters of the biological population and some external factors comprehensively.

In general, an individual prey killed does not contribute instantaneously to the growth of predator population. And differential equations with time delay always exhibit much more complicated dynamics than ordinary differential equations. Hence, in the second part of this paper, the effect of time delay on dynamical behavior of the differential-algebraic model is discussed. It shows that time delay plays an important role in the dynamical behavior of the differential-algebraic model. Hopf bifurcation occurs as time delay increases through a certain threshold, and time delay switches the stability of the proposed model. Furthermore, a state delayed feedback controller is designed to eliminate bifurcational phenomenon and keep population density at steady state. From numerical simulation, it is noted that the stronger Allee effect is not beneficial to the stability of biological species.

It should be noted that almost the existing bioeconomic models (see [13–16]) 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.

Conflict of Interests

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

Acknowledgments

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

References

R. M. Etoua and C. Rousseau, “Bifurcation analysis of a generalized Gause model with prey harvesting and a generalized Holling response function of type III,” Journal of Differential Equations, vol. 249, no. 9, pp. 2316–2356, 2010.

D. Xiao, W. Li, and M. Han, “Dynamics in a ratio-dependent predator-prey model with predator harvesting,” Journal of Mathematical Analysis and Applications, vol. 324, no. 1, pp. 14–29, 2006.

N. Bairagi, S. Chaudhuri, and J. Chattopadhyay, “Harvesting as a disease control measure in an eco-epidemiological system—a theoretical study,” Mathematical Biosciences, vol. 217, no. 2, pp. 134–144, 2009.

T. Das, R. N. Mukherjee, and K. S. Chaudhuri, “Harvesting of a prey-predator fishery in the presence of toxicity,” Applied Mathematical Modelling, vol. 33, no. 5, pp. 2282–2292, 2009.

L. Ji and C. Wu, “Qualitative analysis of a predator-prey model with constant-rate prey harvesting incorporating a constant prey refuge,” Nonlinear Analysis: Real World Applications, vol. 11, no. 4, pp. 2285–2295, 2010.

K. Chakraborty, S. Jana, and T. K. Kar, “Global dynamics and bifurcation in a stage structured prey-predator fishery model with harvesting,” Applied Mathematics and Computation, vol. 218, no. 18, pp. 9271–9290, 2012.

Z. Xiang, D. Long, and X. Song, “A delayed Lotka-Volterra model with birth pulse and impulsive effect at different moment on the prey,” Applied Mathematics and Computation, vol. 219, no. 20, pp. 10263–10270, 2013.

S. Ayasun, C. O. Nwankpa, and H. G. Kwatny, “Computation of singular and singularity induced bifurcation points of differential-algebraic power system model,” IEEE Transactions on Circuits and Systems. I, vol. 51, no. 8, pp. 1525–1538, 2004.

R. Riaza, “Singularity-induced bifurcations in lumped circuits,” IEEE Transactions on Circuits and Systems. I, vol. 52, no. 7, pp. 1442–1450, 2005.

V. M. Marchenko, O. N. Poddubnaya, and Z. Zaczkiewicz, “On the observability of linear differential-algebraic systems with delays,” IEEE Transactions on Automatic Control, vol. 51, no. 8, pp. 1387–1392, 2006.

L. Chen and K. Aihara, “Global searching ability of chaotic neural networks,” IEEE Transactions on Circuits and Systems. I, vol. 46, no. 8, pp. 974–993, 1999.

R. Somogyi and C. A. Sniegoski, “Modeling the complexity of genetic networks: understanding multigenic and pleiotropic regulation,” Complexity, vol. 1, no. 6, pp. 45–63, 1995/96.

X. Zhang, Q.-l. Zhang, and Y. Zhang, “Bifurcations of a class of singular biological economic models,” Chaos, Solitons & Fractals, vol. 40, no. 3, pp. 1309–1318, 2009.

W. Liu, C. Fu, and B. Chen, “Hopf bifurcation for a predator-prey biological economic system with Holling type II functional response,” Journal of the Franklin Institute, vol. 348, no. 6, pp. 1114–1127, 2011.

K. Chakraborty, M. Chakraborty, and T. K. Kar, “Bifurcation and control of a bioeconomic model of a prey-predator system with a time delay,” Nonlinear Analysis: Hybrid Systems, vol. 5, no. 4, pp. 613–625, 2011.

C. Liu, Q. Zhang, X. Zhang, and X. Duan, “Dynamical behavior in a stage-structured differential-algebraic prey-predator model with discrete time delay and harvesting,” Journal of Computational and Applied Mathematics, vol. 231, no. 2, pp. 612–625, 2009.

W. C. Allee, Animal Aggretions: A Study in General Sociology, University of Chicago Press, Chicago, Ill, USA, 1931.

H. S. Gordon, “Economic theory of a common property resource: the fishery,” Journal of Political Economy, vol. 63, pp. 116–124, 1954.

J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, vol. 42, Springer, New York, NY, USA, 1983.

V. Venkatasubramanian, H. Schättler, and J. Zaborszky, “Local bifurcations and feasibility regions in differential-algebraic systems,” IEEE Transactions on Automatic Control, vol. 40, no. 12, pp. 1992–2013, 1995.

J. Zu and M. Mimura, “The impact of Allee effect on a predator-prey system with Holling type II functional response,” Applied Mathematics and Computation, vol. 217, no. 7, pp. 3542–3556, 2010.