predator-prey

In this paper, we study an impulsively controlled predator-prey model with MonodHaldane functional response. By using the Floquet theory, we prove that there exists a stable prey-free solution when the impulsive period is less than some critical value, and give the condition for the permanence of the system. In addition, we show the existence and stability of a positive periodic solution by u sing bifurcation theory.


Introduction
In population dynamics, one of central goals is to understand the dynamical relationship between predator and prey.One important component of the predator-prey relationship is the predator's rate of feeding on prey, i.e., the so-called predator's functional response.Functional response refers to the change in the density of prey attached per unit time per predator as the prey density changes.Holling [7] gave three different kinds of functional response for different kinds of species to model the phenomena of predation, which made the standard Lotka-Volterra system more realistic.These functional responses are monotonic in the first quadrant.But, some experiments and observations indicate that a non-monotonic response occurs at a level: when the nutrient concentration reaches a high level an inhibitory effect on the specific growth rate may occur.To model such an inhibitory effect, the authors in [1,14] where x(t) and y(t) represent population densities of prey and predator at time t.All parameters are positive constants.Usually, a is the intrinsic growth rate of the prey, K is the carrying capacity of the prey, the constant D is the death rate of the predator, m is the rate of conversion of a consumed prey to a predator and b measures the level of prey interference with predation.
As Cushing [5] pointed out that it is necessary and important to consider models with periodic ecological parameters or perturbations which might be quite naturally exposed (for example, those due to seasonal effects of weather, food supply, mating habits or harvesting seasons and so on).Such perturbations were often treated continually.But, there are still some other perturbations such as fire, flood, etc, that are not suitable to be considered continually.These impulsive perturbations bring sudden changes to the system.
In this paper, with the idea of impulsive perturbations, we consider the following predator-prey model with periodic constant impulsive immigration of the predator and periodic harvesting on the prey.
where T is the period of the impulsive immigration or stock of the predator, p i (0 ≤ p i < 1, = 1, 2) are the harvesting control parameters, q is the size of immigration or stock of the predator.This model is an example of impulsive differential equations whose theories and applications were greatly developed by the efforts of Bainov and Lakshmikantham et al. [4,8].
Recently, many researchers have intensively investigated systems with impulsive perturbations (cf, [2,3,10,11,12,13,15,16,17,18,19,20,21]).Most of such systems have dealt with impulsive harvesting and immigration of predators at different fixed times.On the contrary, here we consider the impulsive harvesting and immigration at the same time in our model which has not been studied well until now.EJQTDE, 2010 No. 19, p. 2 The main purpose of this paper is to study the dynamics of the system (2).The organization of this paper is as follows.In the next section, we introduce some notations and lemmas related to impulsive differential equations which are used in this paper.In Section 3, we show the stability of prey-free periodic solutions and give a sufficient condition for the permanence of system (2) by applying the Floquet theory and the comparison theorem.In Section 4, we show the existence of nontrivial periodic solutions via the bifurcation theorem.Finally, in conclusion, we give a bifurcation diagram that shows the system has various dynamic aspects including chaos.

Preliminaries
x, y ≥ 0}.Denote N the set of all of nonnegative integers and f = (f 1 , f 2 ) T the right hand of (2).Let V : R (2)V is locally Lipschitz in x.
The upper right derivatives of V (t, x) with respect to the impulsive differential system (2) is defined as Remarks 2.1.(1) The solution of the system (2) is a piecewise continuous function x : R + → R 2 + , x(t) is continuous on (nT, (n + 1)T ], n ∈ N and x(nT + ) = lim t→nT + x(t) exists.(2) The smoothness properties of f guarantee the global existence and uniqueness of solution of the system (2).(See [8] for the details).
We will use the following important comparison theorem on an impulsive differential equation [8].
Note that dx dt = dy dt = 0 whenever x(t) = y(t) = 0, t = nT .So, we can easily show the following lemma.

Lemma 2.4.
There is an M > 0 such that x(t), y(t) ≤ M for all t large enough, where (x(t), y(t)) is a solution of (2).
Proof.Let x(t) = (x(t), y(t)) be a solution of ( 2) and let Clearly, the right hand of (5), is bounded when From Lemma 2.2, we can obtain that is bounded by a constant for sufficiently large t.Hence there is an M > 0 such that x(t) ≤ M, y(t) ≤ M for a solution (x(t), y(t)) with all t large enough.Now, we consider the following impulsive differential equation.
We can easily obtain the following results.
It is from Lemma 2.5 that the general solution y(t) of ( 7) can be identified with the positive periodic solution y * (t) of ( 7) for sufficiently large t and we can obtain the complete expression for the prey-free periodic solution of (2) for t ∈ (nT, (n + 1)T ].

Extinction and permanence
Now, we present a condition which guarantees local stability of the prey-free periodic solution (0, y * (t)).
Proof.The local stability of the periodic solution (0, y * (t)) of ( 2) may be determined by considering the behavior of small amplitude perturbations of the solution.Let (x(t), y(t)) be any solution of (2).Define x(t) = u(t), y(t) = y * (t) + v(t).Then they may be written as where and Φ(0) = I, the identity matrix.The linearization of the third and fourth equation of (2) becomes u(nT EJQTDE, 2010 No. 19, p. 5 Note that all eigenvalues of S = we have . By Floquet Theory in [4], (0, y * (t)) is locally asymptotically stable if Definition 3.1.The system ( 2) is permanent if there exist M ≥ m > 0 such that, for any solution (x(t), y(t)) of ( 2) with x 0 > 0, .
Proof.Let (x(t), y(t)) be any solution of (2)with x 0 > 0. From Lemma 2.4, we may as- From Lemma 2.5, clearly we have y(t) ≥ m 2 for all t large enough.Now we shall find an m 1 > 0 such that x(t) ≥ m 1 for all t large enough.We will do this in the following two steps. and , t ∈ (nT, (n+1)T ].Then there exists T 1 > 0 such that y(t) ≤ u(t) ≤ u * (t) + ǫ 1 and Integrating ( 12) on (nT, (n + 1)T ](n ≥ N 1 ), we obtain x((n + 1)T ) ≥ x(nT + ) exp which is a contradiction.Hence there exists a t 1 > 0 such that x(t 1 ) ≥ m 3 .
(Step 2) If x(t) ≥ m 3 for all t ≥ t 1 , then we are done.If not, we may let t * = inf t>t 1 {x(t) < m 3 }.Then x(t) ≥ m 3 for t ∈ [t 1 , t * ] and, by the continuity of x(t), we have x(t * ) = m 3 .In this step, we have only to consider two possible cases.
. Since so we have that G(T ) = 0 has a unique positive solution T * .From Theorem 3.1 and Theorem 3.2, we know that the prey-free periodic solution is locally asymptotically stable if T < T * and otherwise, the prey and predator can coexist.Thus T * plays the role of a critical value that discriminates between stability and permanence.

Existence and stability of a positive periodic solution
Now, we deal with a problem of the bifurcation of the nontrivial periodic solution of the system (2), near (0, y * (t)).The following theorem establishes the existence of a positive periodic solution of the system (2) near (0, y * (t)).
Theorem 4.1.The system (2) has a positive periodic solution which is supercritical if T > T * .
Proof.We will use the results of [6,9,11] to prove this theorem.To use theorems of [6,9,11], it is convenient for the computation to exchange the variables of x and y and change the period T to τ .Thus the system (2) becomes as follows Let Φ be the flow associated with (15).We have X(t) = Φ(t, x 0 ), 0 < t ≤ τ , where x 0 = (x(0), y(0)).The flow Φ applies up to time τ .So, X(τ ) = Φ(τ, x 0 ).We will use EJQTDE, 2010 No. 19, p. 9 all notations in [9].Note that Then Now, we can compute where τ 0 is the root of d ′ 0 = 0. Actually, we know that To determine the sign of B, let φ(t) = a − c b y * (t).Then we obtain that and so φ(t) is strictly increasing.Since we have φ(τ 0 ) > 0. This implies that B < 0. Hence BC < 0. Thus, from Theorem 2 of [9], the statement follows.

Conclusions
In this paper, we have studied the influences of impulsive perturbations on a predatorprey system with the Monod-Haldane functional response.We have found out that there exists a threshold value that plays a key role on discriminating between the stability of the prey-free periodic solution and the permanence of the system via Floquet theory and the comparison theorem.Furthermore, we have shown that the system has a positive periodic solution which is supercritical under some conditions.In order to illustrate the dynamics of the system by a numerical example, we give bifurcation diagrams in Figures 1 and 2 when the parameters are fixed except the period T as follows: a = 4.0, K = 1.9, b = 0.6, c = 0.75, D = 0.25, m = 0.6, p 1 = 0.3, p 2 = 0.001 and q = 1.2.EJQTDE, 2010 No. 19, p. 12 These figures point out that system (2) has various dynamical behaviors such as quasiperiodic, periodic windows, strange attractors and periodic doubling and halving phenomena etc. (see Figure 2).In this case, we can obtain the critical value T * ≈ 1.58 suggested in Theorems 3.1 and 3.2.As mentioned in Theorem 4.1, the value T * plays an important part in the classification for the existence of a positive periodic solution as shown in Figure 1.
suggested a function p(x) = αx b + x 2 called the Monod-Haldane function, or Holling type-IV function.We can write a predatorprey model with Monod-Haldane type functional response as follows.