Existence of Periodic Solutions in a Discrete Predator-Prey System with Beddington-DeAngelis Functional Responses

A discrete predator-prey model with Holling II and Beddington-DeAngelis functional responses is investigated. With the aid of differential equations with piecewise constant arguments, a discrete version of continuous nonautonomous delayed predator-prey model with BeddingtonDeAngelis functional responses is proposed. By using Gaines andMawhin’s continuation theorem of coincidence degree theory, sufficient conditions for the existence of positive solutions of the model are established.


Introduction
In population dynamics, the functional response refers to the number of prey eaten per predator per unit time as a function of prey density.Based on a lot of experiments on mammals, Holling 1 proposed three kinds of functional responses as follows: 1 p 1 x ax, where x represents the density of prey.Functions p 1 x , p 2 x , and p 3 x refer to Holling types I, II, and III, respectively, and a > 0 is the predation rate of the predator.After that, the dynamical properties of predator-prey systems with functional response have received great attention from both theoretical and mathematical biologists; for example, Liu et   Recently, Song and Li 9 investigated the following two-prey one-predator model, where two prey are competitive and the predator has Holling functional II functional response: International Journal of Mathematics and Mathematical Sciences 3 where x 1 t and x 2 t stand for the population size of prey pest species and x 3 t is the population size of the predator natural species species; b i > 0 i 1, 2, 3 are intrinsic rates of increase or decrease, α > 0 and β > 0 are parameters representing competitive effects between two prey, η > 0, μ > 0 are positive constants, ηx 1 t / 1 ω 1 x 1 t and μx 2 t / 1 ω 2 x 2 t are the Holling II functional responses.d > 0 is the rate of conversing prey into predator.
To model mutual interference among predators, Beddington 10 andDeAngelis et al. 11 argued that the well-known Holling type II functional response will be replaced by the Beddington-DeAngelis functional responses which is similar to the Holling type II functional response but has an extra term ν i x 3 i 1, 2 in the denominator, then system 1.3 may be modified as the following predator-prey system with Beddington-DeAngelis functional responses: where x 1 t and x 2 t stand for the population size of prey pest species and x 3 t is the population size of the predator natural species species; b i > 0 i 1, 2, 3 are intrinsic rates of increase or decrease, α > 0 and β > 0 are parameters representing competitive effects between two prey, η > 0, μ > 0 are positive constants, ηx 3 t / 1 ω 1 x 1 t ν 1 x 3 t , μx 3 t / 1 ω 2 x 2 t ν 2 x 3 t , ηx 1 t / 1 ω 1 x 1 t ν 1 x 3 t , and μx 2 t / 1 ω 2 x 2 t ν 2 x 3 t are the Beddington-DeAngelis functional responses.d > 0 is the rate of conversing prey into predator.
In the natural word, any biological and environmental parameters are naturally subject to fluctuation in time.The effect of a periodically varying environment is important for evolutionary theory as the selective forces on systems in a fluctuating environment differ from those in a stable environment.Thus, assumptions of periodicity of parameters are a way of incorporating periodicity of the environment, such as seasonal effect of weather, food supplies, and mating habits 12 .Based on the point of view, the modification of 1.4 according to the environmental variation is the nonautonomous differential equations As is known to us, discrete time models governed by difference equations are more appropriate to describe the dynamics relationship among populations than continuous ones when the populations have nonoverlapping generations.Moreover, discrete time models can International Journal of Mathematics and Mathematical Sciences also provide efficient models of continuous ones for numerical simulations.Therefore, it is reasonable and interesting to study discrete time systems governed by difference equations.
Recently, there are some papers which deal with these topics see 13-22 .The principle object of this article is to propose a discrete analogue system 1.5 and explore its dynamics.The remainder of the paper is organized as follows: in Section 2, with the help of differential equations with piecewise constant arguments, we first propose a discrete analogue of system 1.5 , modelling the dynamics of time nonautonomous predator-prey system with Beddington-DeAngelis functional responses, where populations have nonoverlapping generations.In Section 3, based on the coincidence degree and the related continuation theorem, sufficient conditions for the existence of positive solutions of the model are obtained.

A Discrete Version of Model 1.5
There are several different ways of deriving discrete time version of dynamical systems corresponding to continuous time formulations.One of the ways of deriving difference equations modelling the dynamics of populations with nonoverlapping generations that we will use in the following is based on appropriate modifications of models with overlapping generations.For more details about the approach, we refer to 15, 18 .
In the following, we will discrete the system 1.5 .Assume that the average growth rates in system 1.5 change at regular intervals of time, then we can obtain the following modified system: where t denotes the integer part of t, t ∈ 0, ∞ and t / 0, 1, 2, . . . .Equations of type 2.1 are known as differential equations with piecewise constant arguments, and these equations occupy a position midway between differential equations and difference equations.By a solution of 2.1 , we mean a function x x 1 , x 2 , x 3 T , which is defined for t ∈ 0, ∞ and has the following properties: 2 the derivative dx/dt dx 1 t /dt, dx 2 t /dt, dx 3 t /dt T exists at each point t ∈ 0, ∞ with the possible exception of the points t ∈ {0, 1, 2, . ..},where left-sided derivative exists; 3 The equations in 2.1 are satisfied on each interval k, k 1 with k 0, 1, 2, . . . .We integrate 2.1 on any interval of the form k, k 1 , k 0, 1, 2, . .., and obtain for k ≤ t < k 1, k 0, 1, 2, . . .,

Existence of Positive Periodic Solutions
For convenience and simplicity in the following discussion, we always use the notations below throughout the paper: where f k is an ω-periodic sequence of real numbers defined for k ∈ Z.For system 2.3 , we always assume that In order to explore the existence of positive periodic solutions of 2.3 and for the reader's convenience, we will first summarize below a few concepts and results without proof, borrowing from 23 .
Let X, Y be normed vector spaces, L : Dom L ⊂ X → Y is a linear mapping, N : X → Y is a continuous mapping.The mapping L will be called a Fredholm mapping of index zero if dim Ker L codim Im L < ∞ and Im L is closed in Y .If L is a Fredholm mapping of index zero and there exist continuous projectors P : X → X and Q : Y → Y such that International Journal of Mathematics and Mathematical Sciences Im P Ker L, Im L Ker Q Im I − Q , it follows that L| Dom L ∩ Ker P : I − P X → Im L is invertible.We denote the inverse of that map by K P .If Ω is an open bounded subset of X, the mapping N will be called L compact on Ω if QN Ω is bounded and Let l ω ⊂ l 3 denote the subspace of all ω-periodic sequences equipped with the usual supremum norm • , that is, u max k∈I ω |u k |, for any u {u k : k ∈ Z} ∈ l ω .It is easy to show that l ω is a finite-dimensional Banach space. Let

3.4
Then it follows that l ω 0 and l ω c are both closed linear subspaces of l ω and l ω l ω 0 l ω c , dim l ω c 3.

3.5
In the following, we will be ready to establish our result.
Theorem 3.3.Let S 1 , S 1 , S 2 , S 2 , k 0 , and h 0 be defined by 3.30 , 3.50 , 3.34 , 3.46 , 3.39 , and 3.55 , respectively.In addition to the condition (H1), assume further that the following conditions: , and u 3 k ln x 3 k .Then 2.3 takes the form where where u ∈ X, k ∈ Z. Then it is trivial to see that L is a bounded linear operator and Ker L l ω c , Im L l ω 0 , dim Ker L 3 codim Im L.

International Journal of Mathematics and Mathematical Sciences
Then it follows that L is a Fredholm mapping of index zero.Define

3.11
It is not difficult to show that P and Q are continuous projectors such that Furthermore, the generalized inverse to L K P : Im L → Ker P Dom L exists and is given by Obviously, QN and K P I −Q N are continuous.Since X is a finite-dimensional Banach space, using the Ascoli-Arzela theorem, it is not difficult to show that Now we are at the point to search for an appropriate open, bounded subset Ω for the application of the continuation theorem.Corresponding to the operator equation Lu λNu, λ ∈ 0, 1 , we have

3.14
Suppose that u k u 1 k , u 2 k , u 3 k T ∈ X is an arbitrary solution of system 3.14 for a certain λ ∈ 0, 1 .Summing both sides of 3.14 from 0 to ω − 1 with respect to k, respectively, we obtain

3.15
International Journal of Mathematics and Mathematical Sciences 9 It follows from 3.14 and 3.15 that In view of the hypothesis that u {u k } ∈ X, there exist ξ i , η i ∈ I ω such that Then it is obvious that where ∇ denotes the forward difference operator ∇u k u k 1 − u k .In view of 3.8 , we get Then we have which leads to International Journal of Mathematics and Mathematical Sciences From 3.16 , we have In the sequel, we consider two cases.

3.29
It follows from 3.29 that max

3.33
It follows from 3.33 that max

3.36
From 3.16 , we get

3.37
International Journal of Mathematics and Mathematical Sciences Then, where

3.45
Combining both equation of

1. 2
By using dynamical system technique and the geometrical singular perturbation theory, Cantrell et al. 3 made a discussion on the coexistence of predators and prey of system 1.2 which occurred along a stable positive equilibrium.Ko and Ryu 4 discussed the existence, stability, and uniqueness of coexistence states and the extinction and permanence of a diffusive two-competing-prey and one-predator system with Beddington-DeAngelis functional response.Chen et al. 5 analyzed the extinction of the predator and the global asymptotic stability of the boundary solution of a nonautonomous predator-prey system with the Beddington-DeAngelis functional response.Huo et al. 6  studied the global existence of positive periodic solutions for a delayed predator-prey model with the Beddington-DeAngelis functional response.For more knowledge about this topic, one can see 7, 8 .

Lemma 3 .
1 see Continuation Theorem 23 .Let L be a Fredholm mapping of index zero, and let N be L compact on Ω. Suppose that, a for each λ ∈ 0, 1 , every solution x of Lx λNx is such that x / ∈ ∂Ω; b QNx / 0 for each x ∈ Ker L ∩ ∂Ω, and deg{JQN, Ω ∩ Ker L, 0} / 0, then the equation Lx Nx has at least one solution lying in Dom L ∩ Ω. Lemma 3.2 see 15 .Let g : Z → R be ω-periodic, that is, g k ω g k .Then for any fixed k 1 , k 2 ∈ I ω and any k ∈ Z, one has

− exp u 1 b 2 − exp u 2 −
k exp u 1 1 ω 1 k exp u 1 ν 1 k exp u 3 J be the identity mapping.According to the definition of topology, direct calculation yieldsdeg JQN u 1 , u 2 , u 3 T ; Ω ∩ ker L; 0 deg QN u 1 , u 2 , u 3 T ; Ω ∩ Ker L; 0 deg QN u 1 , u 2 , u 3 , 1 T ; Ω ∩ Ker L; 0 deg QN u 1 , u 2 , u 3 , 0 T ; Ω ∩ Ker L; al. 2International Journal of Mathematics and Mathematical Sciences investigated the coexistence of predators and preys of the following predator-prey system with Holling II functional response: stands for the density of the prey, y t and z t are the densities of the predators, respectively,and r, a, A, d, e, D, E, K, b, and B are positive constants.Assuming that one predator consumes prey according to the Holling II functional response and the other predator consumes prey according to the Beddington-DeAngelis functional response, Cantrell et al. 3 proposed the revised version of system 1.1 as follows: By now, we have proved that Ω verifies all requirements of Lemma 3.1.Then it follows that Lu Nu has at least one solution in Dom L ∩ Ω, that is to say, 3.6 has at least one ω-periodic solution in Dom L ∩ Ω, say u * * 3 k T is an ω-periodic solution of system 2.3 with strictly positive components.We complete the proof.