Mathematical Problems in Engineering
Volume 2016 (2016), Article ID 3407463, 9 pages
http://dx.doi.org/10.1155/2016/3407463

Research Article

On a Stochastic Lotka-Volterra Competitive System with Distributed Delay and General Lévy Jumps

Lijie Zhang,¹ Chun Lu,1,2 and Hui Liu¹

¹Department of Mathematics, Qingdao University of Technology, Qingdao 266520, China
²Department of Mathematics, Harbin Institute of Technology, Weihai 264209, China

Received 6 August 2016; Accepted 9 November 2016

Academic Editor: Ana Carpio

Copyright © 2016 Lijie Zhang et al. 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 considers a stochastic competitive system with distributed delay and general Lévy jumps. Almost sufficient and necessary conditions for stability in time average and extinction of each population are established under some assumptions. And two facts are revealed: both stability in time average and extinction have closer relationships with the general Lévy jumps, firstly; and secondly, the distributed delay has no effect on the stability in time average and extinction of the stochastic system. Some simulation figures, which are obtained by the split-step -method to discretize the stochastic model, are introduced to support the analytical findings.

1. Introduction

In recent years, delay differential equations has been used in the study of population dynamics. A famous competitive system with distributed delay can be expressed by

where denotes the size of the th population, , , and are all positive constants, and is a probability measure on . There is an extensive literature concerned with the dynamics of (1) and we here only mention Kuang and Smith [1], Faria [2], Freedman and Wu [3], Bereketoglu and Győri [4], and Gopalsamy [5] among many others. In particular, Kuang (see [6, p. 231]) claimed that if and , then model (1) has a positive equilibrium which is globally asymptotically stable, where , , and . It is important to point out that if and , then .

In the real world, the intrinsic growth rates of many species are always disturbed by environmental noises (see, e.g., [7–10]), which was recognized by many scholars in recent years (see, e.g., [11–14]). In particular, May [7] has pointed out that, due to environmental noises, the birth rates, carrying capacity, and other parameters involved in the system should be stochastic. In this paper, we assume that the parameters and are stochastic; then by the central limit theorem, we can replace and by

where, for , represents a standard Brownian motion defined on a complete probability space and is the intensity of the noise.

On the other hand, the population systems may suffer sudden environmental perturbations, that is, some jump type stochastic perturbations, for example, earthquakes, hurricanes, and epidemics. Some scholars have concentrated on the population systems with compensator jumps, and some significant and interesting results have been obtained (see, e.g., [15–19]). Bao et al. [15, 16] did pioneering work in this field. In addition, Zou et al. [20–22] introduce a general Lévy jumps, which is more reasonable and complicated than the compensator jumps from the viewpoint of biomathematics (see [20]), into population models for the first time. However, there are no articles introducing the general Lévy jumps into population models with distributed delay, to the best of our knowledge. Motivated by these, we consider the famous stochastic competitive system with distributed delay and general Lévy jumps:

where , is a real-valued Poisson counting measure with characteristic measure on a measurable subset of with , , is bounded function, and , ; furthermore, we assume that is independent of . Let the initial data , where represents the family of continuous functions from to with the norm , , . Other parameters are defined and required as before.

For convenience, we introduce the following notations:

Moreover, we impose the following assumptions in this paper.

Assumption 1. There exists a positive constant such that for .

In this paper, we consider a stochastic competitive system with distributed delay and general Lévy jumps. Unlike the deterministic system, the stochastic system does not have an interior equilibrium. Therefore, we cannot investigate the stability of the stochastic system. In Section 2, we show that the solution to system (3) will tend to a point in time average. Furthermore, we establish almost sufficient and necessary conditions for stability in time average and extinction of each population. In Section 3, we present an example to illustrate our mathematical findings. Section 4 gives the conclusions and future directions of the research.

2. Main Content

Lemma 2 (see Liu et al. [23]). Suppose that

(i)If there exist two positive constants and such that for all , where , , are constants, then

(ii)If there exist three positive constants , , and such that for all , then a.s.

In order for the model to be significant, we shall show that the solution is global and nonnegative. However, theorem of existence and uniqueness ([24–28]) is not satisfied in system (3). By using method established by Mao et al. [8], we will show existence and uniqueness of the global positive solution of system (3).

Lemma 3. Let Assumption 1 hold. For any given initial value ; then system (3) has a unique positive solution on a.s. and the solution satisfies

Proof. The proof is similar to Han et al. [29] by defining

where

In addition, applying the inequality, for ,

So we omit it here. Now let us prove inequality (6).
Case  1 (). For any , applying the generalized Itô’s formula [30] to (3) results in

Thus

The rest of proof is analogous with Lemma  4.4 in [15]; we omitted it here.
Case  2 (). The proof is similar to Case  1; we left out it here. The proof is complete.

Theorem 4. For system (3), we suppose that Assumption 1, and , holds.

(I)If and , then both and are extinctive almost surely (a.s.); that is, a.s., .(II)If and , then is extinctive a.s. and is stable in time average a.s.; that is,

(III)If and , then is extinctive a.s. and is stable in time average a.s.; that is,

(IV)If , ,

(A)If and , then is extinctive a.s. and is stable in time average a.s.:

(B)If and , then is extinctive a.s. and is stable in time average a.s.:

(C)If and , then both and are stable in time average a.s.:

Proof. Applying Itô’s formula [30] to the first equality in (3), we get

Making use of the Fubini theorem and a substitution technique, we have

Therefore, we derive

Similarly,

(I) Assume that and From (9),

By the strong law of large numbers for martingales, we therefore have a.s., ; meanwhile, we define, for , , . Under Assumption 1,, making use of the strong law of large numbers for local martingales (see, e.g., [31]), we then derive

and that

Consequently, , a.s. Similarly, by (19), we can show that if , then , a.s.
(II) Suppose that and Since , then, by (I), , a.s. Hence, making use of (9), for arbitrary , there is such that, for ,

Substituting the above inequalities into (19), we can see that, for ,

Since , we can choose sufficiently small such that Applying (i) and (ii) in Lemma 2 to (25) and (26), respectively, we derive

Let . Then, we have , a.s.
The proof of (III) is homogeneous with (II) by symmetry; hence it is omitted.
Now let us prove (IV). For , consider the following equation:

In virtue of the classic stochastic comparison theorem [32], we can find that

Since , , similar to the proof of (II), we can show that

Thus

which, together with (29), implies that

On the other hand, calculating deduces

By virtue of (6) and (32), for arbitrary , there is such that, for ,

Substituting the above inequalities into (33) results in

for . Meanwhile, calculating yields

Using the same way, by (36) we can have that, for ,

(A) Suppose and Note that , and then let be sufficiently small such that Applying (i) in Lemma 2 to (35) gives , a.s. The proof of , a.s., is similar to (II) and hence is omitted.
The proof of (B) is similar to (A) by symmetry and hence is left out.
(C) Suppose that and Since , it then follows from (33) and Lemma 2 that

Making use of the arbitrariness of , we can see that

It follows from (37), Lemma 2, and the arbitrariness of that

likewise. Let be sufficiently small such that When (32) and (39) are used in (19), we get

for sufficiently large . In virtue of (ii) in Lemma 2 and the arbitrariness of , we get

Similarly, substituting (32) and (40) into (20) brings about , a.s. This, together with (39), (40), and (42), means and , a.s.

Remark 5. It is important to designate that if , , and , then and cannot hold simultaneously.

Remark 6. Theorem 4 implies an important fact that when , , the jump process can result in extinction of the population , for example, earthquakes and hurricanes, and when , , the jump process is always advantage for the population , for example, ocean red tide.

Remark 7. From the perspective of the condition in Theorem 4, the distributed delay does not influence some the properties including extinction and stability in time average.

3. Numerical Simulations

In this section, we employ the split-step -method, whose approximate solution is mean-square convergent with order (see [32, 33]), to discretize (3). Here, we choose the initial data , , , , , , , , , and . Then . The main difference between the conditions of the following Case  1–Case  4 is that the values of and are different. In Case  1, we let and . Then by virtue of Kuangs work [6], we have that the positive equilibrium is globally asymptotically stable. Figure 1(a) verifies this. In Case  2, we set , , and . Then and . In view of (A) in Theorem 4, goes to extinction and

Figure 1(b) confirms this. In Case  3, we choose , , , and . Then and . It follows from (B) in Theorem 4 that goes to extinction and

See Figure 1(c). In Case  4, we let , , , and , that is, and . According to (C) in Theorem 4, we obtain

Figure 1(d) validates this.

(a)

(b)

(c)

(d)

Figure 1: The horizontal axis and the vertical axis in this and following figures represent the time and the populations size (step size ).

4. Conclusions and Remarks

This paper investigates a stochastic competitive system with distributed delay and general Lévy jumps. Under the assumption , the almost complete parameter analysis is fulfilled in detail. Our results imply that the general Lévy jumps can significantly change the properties of population models.

Some interesting and significant topics deserve our further engagement. One may put forward a more realistic and sophisticated model to integrate the colored noise into the model [10, 11, 34]. Another significant problem is devoted to stochastic model with infinite delays and general Lévy jumps. We will leave these for future investigation.

It should also be mentioned that “stability in time average” is not a good definition of persistence for stochastic population models. Some papers have introduced more appropriate definitions of permanence for stochastic population models, that is, stochastically persistent in probability or stochastic permanence (see, e.g., [35–37]). We will research these kinds of permanence of model (3) in detail in our following study.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (no. 11501150), a Project of Shandong Province Higher Educational Science and Technology Program of China (J16LI09), the National Social Science Foundation of China (16BJL087), and the instructional reform item of Higher Education of Shandong Province (2015M091).

References

Y. Kuang and H. L. Smith, “Global stability for infinite delay Lotka-Volterra type systems,” Journal of Differential Equations, vol. 103, no. 2, pp. 221–246, 1993.

T. Faria, “Asymptotic stability for delayed logistic type equations,” Mathematical and Computer Modelling, vol. 43, no. 3-4, pp. 433–445, 2006.

H. I. Freedman and J. H. Wu, “Periodic solutions of single-species models with periodic delay,” SIAM Journal on Mathematical Analysis, vol. 23, no. 3, pp. 689–701, 1992.

H. Bereketoglu and I. Győri, “Global asymptotic stability in a nonautonomous Lotka-Volterra type system with infinite delay,” Journal of Mathematical Analysis and Applications, vol. 210, no. 1, pp. 279–291, 1997.

K. Gopalsamy, Stability and Oscillations in Delay Differential Equations of Population Dynamics, vol. 74 of Mathematics and Its Applications, Kluwer Academic, Dordrecht, Netherlands, 1992.

Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, Academic Press, Boston, Mass, USA, 1993.

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

X. Mao, G. Marion, and E. Renshaw, “Environmental Brownian noise suppresses explosions in population dynamics,” Stochastic Processes and their Applications, vol. 97, no. 1, pp. 95–110, 2002.

L. Zu, D. Jiang, D. O'Regan, and B. Ge, “Periodic solution for a non-autonomous Lotka-Volterra predator-prey model with random perturbation,” Journal of Mathematical Analysis and Applications, vol. 430, no. 1, pp. 428–437, 2015.

X. Li and X. Mao, “A note on almost sure asymptotic stability of neutral stochastic delay differential equations with Markovian switching,” Automatica, vol. 48, no. 9, pp. 2329–2334, 2012.

M. Liu and C. Bai, “Dynamics of a stochastic one-prey two-predator model with Lévy jumps,” Applied Mathematics and Computation, vol. 284, pp. 308–321, 2016.

W. Li, H. Song, Y. Qu, and K. Wang, “Global exponential stability for stochastic coupled systems on networks with Markovian switching,” Systems & Control Letters, vol. 62, no. 6, pp. 468–474, 2013.

M. Liu and K. Wang, “A note on a delay Lotka-Volterra competitive system with random perturbations,” Applied Mathematics Letters, vol. 26, no. 6, pp. 589–594, 2013.

C. Lu and X. Ding, “Persistence and extinction in general non-autonomous logistic model with delays and stochastic perturbation,” Applied Mathematics and Computation, vol. 229, pp. 1–15, 2014.

J. Bao, X. Mao, G. Yin, and C. Yuan, “Competitive Lotka-Volterra population dynamics with jumps,” Nonlinear Analysis: Theory, Methods & Applications, vol. 74, no. 17, pp. 6601–6616, 2011.

J. Bao and C. Yuan, “Stochastic population dynamics driven by Lévy noise,” Journal of Mathematical Analysis and Applications, vol. 391, no. 2, pp. 363–375, 2012.

M. Liu and K. Wang, “Dynamics of a Leslie-Gower Holling-type II predatorprey system with Lévy jumps,” Nonlinear Analysis: Theory, Methods & Applications, vol. 85, pp. 204–213, 2013.

M. Liu and C. Bai, “Optimal harvesting of a stochastic mutualism model with Lévy jumps,” Applied Mathematics and Computation, vol. 276, pp. 301–309, 2016.

Q. Ma, D. Ding, and X. Ding, “Mean-square dissipativity of several numerical methods for stochastic differential equations with jumps,” Applied Numerical Mathematics, vol. 82, pp. 44–50, 2014.

X. Zou and K. Wang, “Numerical simulations and modeling for stochastic biological systems with jumps,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 5, pp. 1557–1568, 2014.

R. Wu, X. Zou, and K. Wang, “Asymptotic properties of stochastic hybrid Gilpin-Ayala system with jumps,” Applied Mathematics and Computation, vol. 249, pp. 53–66, 2014.

X. Zhang and K. Wang, “Stability analysis of a stochastic Gilpin-Ayala model driven by Lévy noise,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 5, pp. 1391–1399, 2014.

M. Liu, K. Wang, and Q. Wu, “Survival analysis of stochastic competitive models in a polluted environment and stochastic competitive exclusion principle,” Bulletin of Mathematical Biology, vol. 73, no. 9, pp. 1969–2012, 2011.

B. Øksendal and A. Sulem, Applied Stochastic Control of Jump Diffusions, Springer, Berlin, Germany, 2005.

N. Ikeda and S. Watanabe, Stochastic Differential Equations and Diffusion Processes, North-Holland, Amsterdam, The Netherlands, 1981.

N. Jacob, Y. Wang, and C. Yuan, “Stochastic differential delay equations with jumps, under nonlinear growth condition,” Stochastics, vol. 81, no. 6, pp. 571–588, 2009.

J. Bao and C. Yuan, “Comparison theorem for stochastic differential delay equations with jumps,” Acta Applicandae Mathematicae, vol. 116, no. 2, pp. 119–132, 2011.

C. Yuan and X. Mao, “Stability of stochastic delay hybrid systems with jumps,” European Journal of Control, vol. 16, no. 6, pp. 595–608, 2010.

Q. Han, D. Jiang, and C. Ji, “Analysis of a delayed stochastic predator-prey model in a polluted environment,” Applied Mathematical Modelling, vol. 38, no. 13, pp. 3067–3080, 2014.

H. Kunita, “Itô’s stochastic calculus: its surprising power for applications,” Stochastic Processes and Their Applications, vol. 120, no. 5, pp. 622–652, 2010.

R. S. Liptser, “A strong law of large numbers for local martingales,” Stochastics, vol. 3, no. 3, pp. 217–228, 1980.

Q. Ma, D. Ding, and X. Ding, “Mean-square dissipativity of several numerical methods for stochastic differential equations with jumps,” Applied Numerical Mathematics. An IMACS Journal, vol. 82, pp. 44–50, 2014.

X. Ding, Q. Ma, and L. Zhang, “Convergence and stability of the split-step θ-method for stochastic differential equations,” Computers & Mathematics with Applications, vol. 60, no. 5, pp. 1310–1321, 2010.

X. Li and G. Yin, “Logistic models with regime switching: permanence and ergodicity,” Journal of Mathematical Analysis and Applications, vol. 441, no. 2, pp. 593–611, 2016.

S. J. Schreiber, M. Benaïm, and K. A. S. Atchadé, “Persistence in fluctuating environments,” Journal of Mathematical Biology, vol. 62, no. 5, pp. 655–683, 2011.

M. Liu and C. Bai, “Analysis of a stochastic tri-trophic food-chain model with harvesting,” Journal of Mathematical Biology, vol. 73, no. 3, pp. 597–625, 2016.

M. Liu and M. Fan, “Permanence of stochastic Lotka-Volterra systems,” Journal of Nonlinear Science, 2016.