Abstract and Applied Analysis

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

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

Research Article

Dynamical Behavior of a Stochastic Food-Chain System with Beddington-DeAngelis Functional Response

Yanming Ge and Zigen Song

College of Information Technology, Shanghai Ocean University, Shanghai 201306, China

Received 11 March 2014; Revised 13 May 2014; Accepted 18 May 2014; Published 11 June 2014

Academic Editor: Yongli Song

Copyright © 2014 Yanming Ge and Zigen Song. 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.

We investigate a stochastic Food-Chain System , , , where , = is a standard Brownian motion. Firstly, the existence, the uniqueness, and the positivity of the solution are proved. Secondly, the stochastically ultimate boundedness of the system is investigated. Thirdly, the boundedness of moments and upper-growth rate of the solution are obtained. Then the global attractivity of the system is discussed. Finally, the main results are illustrated by several examples.

1. Introduction

In recent decades, a number of authors have studied the predator-prey type systems [1–8], especially the simple Food-Chain Systems [9–14]. These deterministic systems assume that all parameters in the models are deterministic, and some beautiful and valuable conclusions are obtained. Then it is well known that biology systems are often subject to environmental noise [15]. It is therefore useful to reveal how the noise affects these systems. Though there are various types of environmental noise, May [16] have claimed that the birth rates in the system should be stochastic. At the same time, the natural growth of many species vary with time, for example, due to the seasonality.

On the other hand, basing on some natural phenomenon and statistics in laboratory, the predator-dependent functional response should be considered in the predator-prey type systems. As we know, there are three classical predator-dependent functional responses: Hassell and Varley [17], Beddington et al. [18, 19], and Crowley and Martin [20]. Recently, Skalski and Gilliam [21] pointed out that these functional responses can provide better descriptions, and the Beddington-DeAngelis (B-D) type functional response performed even better in some cases.

Based on this background, in this paper we will discuss the following stochastic system:

where , , is a Brownian motion defined on a complete probability space . All parametric functions are continuous and bounded functions on . , and denote the population density of resource, consumer, and predator at time , respectively. For convenience, the model, consisting of (1), together with , , , will be referred to as model SDE. On the other hand, for easy expression, if is a continuous bounded function on , we define and . Throughout this paper, we assume that , , and the others are similar. The notation used in this paper expresses the positive constant, for convenience, it denotes different constant in different places, and they only differ in positive constant times.

The rest of the paper is organized as follows. In Section 2, we show that the system has a global positive solution. In Section 3, the stochastically ultimate boundedness of the system is acquired. In Section 4, we obtain the asymptotic moment estimation and the upper-growth rate of the solution. In Section 5, the global attractivity of SDE is investigated. In Section 6, we illustrate our main results through several numerical examples. Last but not least, conclusions are drawn in Section 7.

The key method used in this paper is the analysis of Lyapunov functions. This Lyapunov function analysis for stochastic differential equations was developed by Khasminskii [22] and has been used by many authors [23–28].

2. Positive and Global Solutions

As the state of the SDE is the population density of species in the system at time , it should be nonnegative. Moreover, in order for a stochastic differential equation to have a unique global (i.e., no explosion in a finite time) solution for any given initial data, the coefficients of the equation are generally required to satisfy the linear growth condition and local Lipschitz condition [25]. However, the coefficients of each equation in SDE neither obey the linear growth condition nor local Lipschitz continuous. In this section, we show existence and uniqueness of the positive solution.

Lemma 1. For any initial value , , , system (SDE) has a unique positive local solution for almost surely (a.s.), where is the explosion time.

Proof. To begin with, consider the following equations:

on with initial value , , , where , , . Notice that the last equations’ coefficients satisfy the local Lipschitz condition; thus, there is a unique solution on . Therefore, it follows from Itô’s formula that , , is the unique positive local solution of SDE with initial value , , .

Theorem 2. Consider system SDE, for any given initial value , and there is a unique solution on and the solution will remain in with probability 1, where .

Proof. From Lemma 1, to complete the proof we only need to show that . Hence, let be sufficiently large for , , and lying within the interval . For each integer , define the stopping times:

Clearly, is increasing as . Set , from which . Now, we only need to prove . If this statement is false, there is a pair of constants and such that . Thus there exists an integer such that

Define . If , we have

Integrating both sides of the above inequality from 0 to and then taking the expectations lead to

Set ; then it follows from inequality (4) that . Clearly, for every , there is , , or which equals either or ; therefore, is no less than

Then an application of (6) results in

where is the indicator function of . Letting leads to the contradiction

3. Stochastically Ultimate Boundedness

In the previous section, we have showed that the solutions of SDE will remain in the positive cone forever. However, this nonexplosion property in a population dynamical system is often not good enough while the property of ultimate boundedness is more desired. Let us discuss it in this section. We begin with giving the definition of stochastically ultimate boundedness.

Definition 3. The SDE is said to be stochastically ultimate bounded if for any , there is a positive constant such that for any initial value, the solution of SDE has the property that

Theorem 4. For any there is a positive constant such that for any initial value, the solution of SDE has the property that

In particular, the SDE is stochastically ultimate bounded.

Proof. Define

Applying Itôs formula gives

Using Itôs formula again we have

Integrating the above inequality and then taking expectations on both sides, we can see that

This implies immediately that

On the other hand,

and then we have

which yields the required assertion by setting . Especially, let ; then

Now, for any , let , by Chebyshev’s inequality,

Hence,

In other word, SDE is stochastically ultimate bounded.

4. Boundedness of Moments and Upper-Growth Rate Estimation

Theorems 2 and 4 show that the solutions of SDE will remain in a subset of . Now let us construct different types of Lyapunov functions to investigate how the solutions behave in .

Theorem 5. Consider system SDE for any positive numbers , there exists a positive number such that, for any initial value , the solution of SDE obeys

Proof. Define

Applying Itôs formula gives

Integrating the above inequality and then taking expectations on both sides, we can see that

Together with , this inequality shows that

From this inequality, we obtain the desired assertion.

This theorem, together with Theorem 2, denotes , in other word, the average value of on time and “space” is bounded. It is similar to and .

Theorem 6. For any positive numbers , , and , solutions of SDE satisfy

Proof. Define

Applying Itôs formula gives

Let

where , and by the Itôs formula, we have

Therefore,

Define

Notice that the quadratic variation of is

Let and . In view of the exponential martingale inequality, for every and , we have

An application of Borel-Cantelli lemma yields that, for almost any , there is a random integer such that, for ,

So

for , a.s. Substituting the above inequality into (32), we have

In other words, we have shown that

Thus, for all , , we have

Then

Letting , , and , we can obtain the required assertion.

Based on this theorem, obviously, the upper-growth rate is linear.

5. Global Attractivity

In this section, we turn to establishing sufficient criteria for the global attractivity of system SDE. Firstly, we prepare some useful lemmas.

Lemma 7. Let be a solution to SDE with initial value . Then for all , there exist three constants related to : for , such that

Proof. Define for . Making use of Itôs formula we have

Applying Itôs formula again to and then taking expectations on both sides, we can see that

which is the desired assertion of , and the other two inequalities can be proved similarly.

Lemma 7 tells us that there is a , such that , , and for all . At the same time, it follows from the continuity of , , and that there exist , , and such that , , and for all . Denote , , and . Then for all ,

Lemma 8 (see [29]). Suppose that an -dimensional stochastic process on satisfies the condition

for some positive constants , , and . Then there exists a continuous modification of which has the property that for every there is a positive random variable such that

In other words, almost every sample path of is locally but uniformly Hölder continuous with exponent .

Lemma 9. Let be a solution of SDE on with initial value ; then almost every sample path of is uniformly continuous on .

Proof. The second equation of SDE is equivalent to the following stochastic integral equation:

We estimate

On the other hand, as an application of the moment inequality for stochastic integrals [30], we see that for and ,

Then for , , , we obtain that

Then it follows from Lemma 8 that almost every sample path of is locally but uniformly Hölder-continuous with exponent for every and therefore almost every sample path of is uniformly continuous on . Similarly, we can show that almost every sample path of and is uniformly continuous on .

Lemma 10 (see [31]). Let be a nonnegative function defined on such that is integrable on and is uniformly continuous on . Then .

Now, we give our main result of this section.

Definition 11 (see [32]). Let and be two arbitrary solutions of system SDE with initial values , . If

then we say that the system is globally attractive.

Theorem 12. If , , , then system SDE is globally attractive.

Proof. Define

then it is a continuous and positive function on . A direct calculation of the right differential of , and then applying Itö’s formula, we have

Integrating both sides gives

Therefore,

Making use of and the known conditions result in that

Then the required assertion follows from Lemmas 9 and 10 immediately.

6. Numerical Results

We present numerical experiments in this section to show how the proposed model works in the constructive examples. The results enhance the readers to understand the theoretical conclusions from the practical applications.

Here, we use the Milstein method [33] to construct the discretization equation of SDE. That is,

where , , and , , are the Gaussian random variables. We obtain the following results.

As we pointed out in Theorem 6, a.s., and Figure 1 confirms this. Here we choose , , , , , , , , , , , and .

Figure 1: The curve of for initial values and .

Figure 2 confirms that the system is globally attractive (Theorem 12). Here , , , , , , , , , . The initial values are and . Then , , , and the assumptions of Theorem 12 are satisfied.

Figure 2: The curves of , and for initial values and , .

7. Conclusions

This paper studies a stochastic nonautonomous food-chain system with Beddington-DeAngelis functional response. This system retains some conventional properties of stochastic systems. That is, it has a unique global positive solution for any initial conditions, the moments of the solutions are bounded, and the solutions are global attractive and stochastically ultimate bounded.

Nevertheless, there are rooms to continue working on this issue. It is interesting to study other parameters; for example, are stochastic. The permanence and extinction of the system and the stability in distribution deserve to be investigated too.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research is supported by the National Natural Science Foundation of China under Grant no. 11302126.

References

C. Azar, J. Holmberg, and K. Lindgren, “Stability analysis of harvesting in a predator-prey model,” Journal of Theoretical Biology, vol. 174, no. 1, pp. 13–19, 1995.

M. R. Myerscough, M. J. Darwen, and W. L. Hogarth, “Stability, persistence and structural stability in a classical predator-prey model,” Ecological Modelling, vol. 89, no. 1–3, pp. 31–42, 1996.

J. A. Axelsen, N. Holst, T. Hamers, and P. H. Krogh, “Simulations of the predator prey interactions in a two species ecotoxicological test system,” Ecological Modelling, vol. 101, no. 1, pp. 15–25, 1997.

H. Malchow, “Motional instabilities in prey-predator systems,” Journal of Theoretical Biology, vol. 204, no. 4, pp. 639–647, 2000.

D. V. Vayenas, G. Aggelis, V. Tsagou, and S. Pavlou, “Dynamics of a two-prey-one-predator system with predator switching regulated by a catabolic repression control-like mode,” Ecological Modelling, vol. 186, no. 3, pp. 345–357, 2005.

A. Verdy, “Modulation of predator-prey interactions by the Allee effect,” Ecological Modelling, vol. 221, no. 8, pp. 1098–1107, 2010.

M. M. Molina, M. A. Moreno-Armendáriz, and J. C. Seck Tuoh Mora, “On the spatial dynamics and oscillatory behavior of a predator-prey model based on cellular automata and local particle swarm optimization,” Journal of Theoretical Biology, vol. 336, pp. 173–184, 2013.

V. Křivan, “Behavioral refuges and predator-prey coexistence,” Journal of Theoretical Biology, vol. 339, pp. 112–121, 2013.

P. Turchin, Complex Population Dynamics: A Theoretical/Empirical Synthesis, Princeton University Press, Princeton, NJ, USA, 2003.

R. Xu, L. Chen, and F. Hao, “Periodic solutions of a discrete time Lotka-Volterra type food-chain model with delays,” Applied Mathematics and Computation, vol. 171, no. 1, pp. 91–103, 2005.

R. Xu and Z. Ma, “Stability and Hopf bifurcation in a ratio-dependent predator-prey system with stage structure,” Chaos, Solitons & Fractals, vol. 38, no. 3, pp. 669–684, 2008.

Z. Hu, G. Gao, and W. Ma, “Dynamics of a three-species ratio-dependent diffusive model,” Nonlinear Analysis: Real World Applications, vol. 11, no. 3, pp. 2106–2114, 2010.

W. Li and L. Wang, “Stability and bifurcation of a delayed three-level food chain model with Beddington-DeAngelis functional response,” Nonlinear Analysis: Real World Applications, vol. 10, no. 4, pp. 2471–2477, 2009.

Q. Gan, R. Xu, and P. Yang, “Bifurcation and chaos in a ratio-dependent predatory-prey system with time delay,” Chaos, Solitons and Fractals, vol. 39, no. 4, pp. 1883–1895, 2009.

T. C. Gard, Introduction to stochastic differential equations, Marcel Dekker, New York, NY, USA, 1988.

R. M. May, Stability and Complexity in Model Ecosystems, Princeton University Press, Princeton, NJ, USA, 2001.

M. P. Hassell and G. C. Varley, “New inductive population model for insect parasites and its bearing on biological control,” Nature, vol. 223, no. 5211, pp. 1133–1137, 1969.

J. R. Beddington, “Mutual interference between parasites or predators and its effect on searching efficiency,” Journal of Animal Ecology, vol. 44, pp. 331–341, 1975.

D. L. DeAngelis, R. A. Goldsten, and R. Neill, “A model for trophic interaction,” Ecology, vol. 56, pp. 881–892, 1975.

P. H. Crowley and E. K. Martin, “Functional response and interference within and between year classes of a dragonfly population,” Journal of the North American Benthological Society, vol. 8, pp. 211–221, 1989.

G. T. Skalski and J. F. Gilliam, “Functional responses with predator interference: viable alternatives to the Holling type II model,” Ecology, vol. 82, no. 11, pp. 3083–3092, 2001.

R. Z. Khasminskii, Stochastic Stability of Differential Equations, Sijthoff & Noordhoff, Germantown, Md, USA, 1980.

L. Arnold, Stochastic Differential Equations: Theory and Applications, John Wiley & Sons, New York, NY, USA, 1972.

A. Friedman, Stochastic Differential Equations and Their Applications, Academic Press, New York, NY, USA, 1976.

X. Mao and C. Yuan, Stochastic Differential Equations with Markovian Switching, Imperial College Press, London, UK, 2006.

X. Mao, Stochastic Differential Equations and Applications, Horwood Publishing, Chichester, UK, 2nd edition, 2007.

M. B. Nevelson and R. Z. Khasminskii, Stochastic Approximation and Recurrent Estimation, Nauka, Moscow, Russia, 1972.

H. Qiu, M. Liu, K. Wang, and Y. Wang, “Dynamics of a stochastic predator-prey system with Beddington-DeAngelis functional response,” Applied Mathematics and Computation, vol. 219, no. 4, pp. 2303–2312, 2012.

I. Karatzas and S. E. Shreve, Brownian Motion and Stochastic Calculus, Springer, Berlin, Germany, 1991.

X. Mao, Stochastic Differential Equations and Applications, Horwood Publishing, Chichester, UK, 2008.

I. Barbălat, “Systèmes d'équations différentielles d'oscillations non linéaires,” Revue Roumaine de Mathematiques Pures et Appliquees, vol. 4, pp. 267–270, 1959.

X. Li and X. Mao, “Population dynamical behavior of non-autonomous Lotka-Volterra competitive system with random perturbation,” Discrete and Continuous Dynamical Systems A, vol. 24, no. 2, pp. 523–545, 2009.

D. J. Higham, “An algorithmic introduction to numerical simulation of stochastic differential equations,” SIAM Review, vol. 43, no. 3, pp. 525–546, 2001.