Periodic solution of a stage-structured predator-prey model with Crowley-Martin type functional response

: In this paper, the existence of positive periodic solution of stage-structured predator-prey model with Crowley-Martin type functional response is investigated. The prey population fall into two categories: mature and immature prey. The predator population is dependent only on mature prey and is inﬂuenced by Crowley-Martin type functional response. Based on the Mawhin’s coincidence degree theory and nontrivial estimation techniques for a priori bounds of unknown solutions to the operator equation Fz = µ Nz , we prove the existence of positive periodic solution. Finally, the e ﬀ ectiveness of our result is veriﬁed by an example and numerical simulation.

On the other hand, in the study of predator-prey system, functional response also plays an important role, which can represent the quantity of prey killed by a predator per unit time and describe the amount of biological transfer between different nutritional levels. Holling type-II functional response [16,25,26] takes into account the average feeding rate of a predator to its prey. It can be expressed as: f (x, y) = a 1 x 1+a 2 x , where a 1 and a 2 are positive constants, which denote the capture rate and the influence of the processing time, respectively. In this functional response, competition between predators for food only occurs when prey is depleted. Another functional response function is the Beddinton-DeAnglis type [22,27,28], which is similar to the Holling type-II. However, it takes into account the interference between predators. Therefore, Beddinton-DeAnglis type functional response function describes that individuals from two or more predator groups not only take prey, but also meet and compete with other predators. It has the type: f (x, y) = a 1 x 1+a 2 x+a 3 y , where a 3 is a positive constant, describing the degree of disturbance between predators. Moreover, it is assumed that its influence on the predation rate can be ignored in the case of high prey density. In this paper, we consider the Crowley-Martin type functional response [29][30][31][32]: .
It also takes the interference between predators into account, but the biggest difference between it and Beddinton-DeAnglis type functional response is: the influence of predator disturbance on the predation rate is always an important factor, which cannot be ignored. Hence, Crowley-Martin type functional response is more consistent with the phenomenon in ecology and has more research value.
Maiti et al. [32] studied the global dynamics of an autonomous stage-structured predator-prey model with Crowley-Martin type functional response. However, they did not cosider the periodic behavior of this model with periodic parameters. The assumption of periodicity of the parameters is a way of incorporating the periodicity of the environment. The periodic oscillation of the parameters seems reasonable in view of seasonal factors, e.g. mating habits, availability of food, weather conditions, harvesting and hunting, etc. Cai et al. [30] presented the existence of positive periodic solutions of an Eco-Epidemic model with Crowley-Martin type functional response. Inspired by the above works, we study a stage-structured predator-prey system with Crowley-Martin type functional response: where x 1 (t) and x 2 (t) are the population density of immature and mature prey at time t, y(t) is the population density of predator at time t, and all the following parameters involved are continuous positive periodic functions: (I) for immature prey: (1) the ratio function s(t) represents the ratio of birth rate to available mature prey; (2) the ratio of the conversion of immature prey to mature prey to existing immature prey is denoted by r(t), and the ratio of the death rate of immature prey to existing immature prey is represented by d(t). (II) for mature prey: (1) α(t) is an internally specific interference function; (2) the ratio of the death rate to existing mature prey is denoted by d 1 (t); (3) the interaction between predator and mature prey is a Crowley-Martin type functional response with rate β(t). (III) for predator: (1) β 1 (t) denotes the intake of predators, and 0 < β 1 (t) < β(t); (2) d 2 (t) denotes the death rate of the predator; (3) γ(t) denotes the internal specific disturbance function for the predator.
In terms of the number of creatures, the initial conditions for model (1.1) are given by (x 1 (t), x 2 (t), y(t)) ∈ C + = C(0, R 3 The purpose in the present paper is to find some suitable conditions of the existence of positive periodic solution for system (1.1). The method is based on Mawhin's coincidence degree theory.

Existence of positive periodic solution
In this section, we establish the existence of positive periodic solution for the system (1.1). For this purpose, we first assume that the parameters in the system (1.1) are all ω-period. To obtain a positive periodic solution for system (1.1), we summarize the following lemmas. For the sake of simplicity, we use the notations: (H 3 )sp + > (r +d)e 2(d+d 1 )ω . From a biological viewpoint, the assumptions (H 1 ) and (H 2 ) imply that the intake of a mature predator is greater than its death rate, while the assumption (H 3 ) implies that the birth rate of the prey is influenced by the stage structure and the death rate of mature and immature prey. Now, we present a theorem on the existence. Proof: Firstly, replacing the variables by Then system (1.1) changes into (2.1) Z, W are both Banach space with the norm || · || as follows For any z = (z 1 , z 2 , z 3 ) ∈ Z, the periodicity of system (2.1) implies: and are ω-period functions. In fact, Obviously, Γ 2 (z, t), Γ 3 (z, t) are also both periodic functions by a similar way. Set From the above definition, we have Hence, F is a Fredholm map of index zero. Moreover, it is clear that P and Q are continuous projection operators with ker F = ImP, and ImF = ker Q = Im(I − Q).
Therefore, the inverse K p : F → DomF ∩ ker P exists and is given by Clearly, QN and K p (I − Q)N are continuous. Due to Z is a Banach space, using the Arzala-Ascoli theorem, we have that N is F-compact onŪ for any open bounded set U ⊂ Z.
Next, in order to apply the continuation theorem, we need construct an appropriate open bounded subset U. Therefore, the operator equation is defined by Fz = µNz, µ ∈ (0, 1), that is,
In view of (H 1 ) and (H 2 ), the product of two solutions of the above equation must be less than zero, thus there is only one positive solution z 3 (t). Since Q = ker F, we deduce that J = I. Thus where deg(·, ·, ·) is Brouwer's degree. Hence, the requirement (2) of Lemma 2.1 also holds. From Lemma 2.1 and the periodicity of the system, we know that system (1.1) has at least one positive periodic solution on DomF ∩Ū. This completes the proof of Theorem 2.1.
(3.1) where x 1 (t) and x 2 (t) are the population density of immature and mature prey at time t, y(t) is the population density of predator at time t. s(t) = 10 + sin t denotes the ratio of birth rate to available mature prey, r(t) = 7 + sin t represents the ratio of the conversion of immature prey to mature prey to existing immature prey, d(t) = 0.3 and d 1 (t) = 0.2 are the ratio of the death rate of immature prey and mature prey, respectively. α(t) = 0.1 is an internally specific interference coefficient. The term (2+sin t)x 2 (t)y(t) (1+0.5x 2 (t))(1+0.5y(t)) stands for the Crowley-Martin type functional response with rate β(t) = 2 + sin t. Moreover, β 1 (t) = 1.5 + sin t denotes the intake of predator; d 2 (t) = 0.1 is the death rate of the predator and γ(t) = 0.05 is the internal specific disturbance coefficient for the predator. Simple computation shows (H 1 )-(H 3 ) in Theorem 2.1 are satisfied, we conclude that system (3.1) has at least one positive periodic solution. Now we take the initial values x 1 (0) = 0.5, x 2 (0) = 0.5, y(0) = 5, then the following figure is obtained through Maple software. As can be seen from the above figure, system (3.1) has a periodic solution with a period 2π. Namely, immature prey, mature prey and predator all have a periodic solution with period 2π in the same periodic environment. On the other hand, it also shows that our results are feasible based on the method of Mawhins coincidence degree theory and some nontrivial estimation techniques.

Conclusion
In the present paper, a stage-structured predator-prey model with Crowley-Martin type functional response is considered. It is assumed that the prey population can divide into two parts: mature and immature prey. The predator population is only dependent on mature prey and is influenced by Crowley-Martin type functional response. Based on the method of Mawhin's coincidence degree theory and novel estimation techniques for a priori bounds of unknown solutions to the operator equation Fz = µNz, we obtain some interesting and novel sufficient conditions for the existence of positive periodic solution of the ecological model. Another interesting topic is the existence of analytic periodic solution. Kosov and Semenov [33,34] studied the existence of (analytic) exact periodic solutions of some nonlinear differential equations, and showed that these periodic solutions are analytic functions under some sufficiently conditions. However, it is impossible to get the analytic periodic solutions in this model due to its great complexity. The analytic periodic solutions can be obtained for some particular systems, it is difficult for the very complicated system. Moreover, our method is based on Mawhin's coincidence degree theory. Thus, we can prove the existence of the periodic solutions by this method.