We analyze a three species predator-prey chain model with stochastic perturbation. First, we show that this system has a unique positive solution and its pth moment is bounded. Then, we deduce conditions that the system is persistent in time average. After that, conditions for the system going to be extinction in probability are established. At last, numerical simulations are carried out to support our results.
1. Introduction
Recently, the dynamical relationship between predator-prey has been one of the dominant themes in both ecology and mathematical ecology due to its universal importance. Especially, the predator-prey chain model is the typical representative. Thereby it significantly changed the biology, the understanding of the existence, and development of the basic law and has made the model become a research hot spot. One of the most famous models for population dynamics is the Lotka-Volterra predator-prey system which has received plenty of attention and has been studied extensively; see [1–4]. Specially persistence and extinction of this model are interesting topics.
The three species predator-prey chain model is described as follows:
(1)x1˙(t)=x1(t)(a1-b11x1(t)-b12x2(t)),x2˙(t)=x2(t)(-a2+b21x1(t)-b22x2(t)-b23x3(t)),x3˙(t)=x3(t)(-a3+b32x2(t)-b33x3(t)),
where xi(t)(i=1,2,3) denotes the population densities of the species at time t. The parameters a1,a2,a3,bii(i=1,2,3) are positive constants that stand for intrinsic growth rate, predator death rate of the second species, predator death rate of the third species, coefficient of internal competition, respectively. b21, b32 represent saturated rate of the second and the third predator and b12, b23 represent the decrement rate of predator to prey.
System (1) describes a three species predator-prey chain model in which the latter preys on the former. From a biological viewpoint, we not only require the positive solution of the system but also require its unexploded property in any finite time and stability.
We know that the global asymptotic stability of a positive equilibrium x*=(x1*,x2*,x3*) holds and is global stability if the following condition holds:
(2)a1-b11b21a2-b11b22+b12b21b21b32a3>0,
which could refer to [5]. However, population dynamics in the real world is inevitably affected by environmental noise (see, e.g., [6, 7]). Parameters involved in the system are not absolute constants, they always fluctuate around some average values. The deterministic models assume that parameters in the systems are deterministic irrespective of environmental fluctuations which impose some limitations in mathematical modeling of ecological systems. So we cannot omit the influence of the noise on the system. Recently many authors have discussed population systems subject to white noise (see, e.g., [8–15]). May (see, e.g., [16]) pointed out that due to continuous fluctuation in the environment, the birth rates, death rates, saturated rate, competition coefficients, and all other parameters involved in the model exhibit random fluctuation to some extent, and as a result the equilibrium population distribution never attains a steady value but fluctuates randomly around some average value. Sometimes, large amplitude fluctuation in population will lead to the extinction of certain species, which does not happen in deterministic models.
Therefore, Lotka-Volterra predator-prey chain models in random environments are becoming more and more popular. Ji et al. [14, 15] investigated the asymptotic behavior of the stochastic predator-prey system with perturbation. Liu and Chen introduced periodic constant impulsive immigration of predator into predator-prey system and gave conditions for the system to be extinct and permanence. Polansky [17] and Barra et al. [18] have given some special systems of their invariant distribution. After that, Gard [5] analysed that under some conditions the stochastic food chain model exists an invariant distribution. However, seldom people study the persistent and nonpersistent of the food chain model with stochastic perturbation.
In this paper, we introduce the white noise into the intrinsic growth rate of system (1), and suppose ai→ai+σiB˙i(t)(i=1,2,3); then we obtain the following stochastic system:
(3)x˙1(t)=x1(t)(a1-b11x1(t)-b12x2(t))+σ1x1(t)B˙1(t),x˙2(t)=x2(t)(-a2+b21x1(t)-b22x2(t)-b23x3(t))-σ2x2(t)B˙2(t),x˙3(t)=x3(t)(-a3+b32x2(t)-b33x3(t))-σ3x3(t)B˙3(t),
where Bi(t)(i=1,2,3) are independent white noises with Bi(0)=0, σi2>0(i=1,2,3) representing the intensities of the noise.
The aim of this paper is to discuss the long time behavior of system (3). We have mentioned that x*=(x1*,x2*,x3*) is the positive equilibrium of system (1). But, when it suffers stochastic perturbations, there is no positive equilibrium. Hence, it is impossible that the solution of system (3) will tend to a fixed point. In this paper, we show that system (3) is persistent in time average. Furthermore, under certain conditions, we prove that the population of system (3) will die out in probability which will not happen in deterministic system and could reveal that large white noise may lead to extinction.
The rest of this paper is organized as follows. In Section 2, we show that there is a unique nonnegative solution of system (3), and its pth moment is bounded. In Section 3, we show that system (3) is persistent in time average. While in Section 4, we consider three situations when the population of the system will be extinction. In Section 5, numerical simulations are carried out to support our results.
Throughout this paper, unless otherwise specified, let (Ω,{ℱt}t≥0,P) be a complete probability space with a filtration {ℱt}t≥0 satisfying the usual conditions (i.e., it is right continuous and ℱ0 contains all P-null sets). Let R+3 denote the positive cone of R3; namely, R+3={x∈R3:xi>0,1≤i≤3}, R-+3={x∈R3:xi≥0,1≤i≤3}.
2. Existence and Uniqueness of the Nonnegative Solution
To investigate the dynamical behavior, the first concern thing is whether the solution is global existence. Moreover, for a population model, whether the solution is nonnegative is also considered. Hence, in this section, we show that the solution of system (3) is global and nonnegative. As we have known, in order for a stochastic differential equation to have a unique global (i.e., no explosion at a finite time) solution with any given initial value, the coefficients of the equation are generally required to satisfy the linear growth condition and local Lipschitz condition (see, e.g., [19]). It is easy to see that the coefficients of system (3) are locally Lipschitz continuous, so system (3) has a local solution. By Lyapunov analysis method, we show the global existence of this solution.
Theorem 1.
For any given initial value x(0)=x0∈R+3, system (3) has a unique global positive solution x(t)=(x1(t),x2(t),x3(t)) for all t≥0 with probability one.
Proof.
It is clear that the coefficients of system (3) are locally Lipschitz continuous for the given initial value x(0)=x0∈R+3. So there is a unique local solution x(t) on t∈[0,τe), where τe is the explosion time (see, e.g., [19]). To show that this solution is global, we need to show that τe=∞ a.s. Let m0≥1 be sufficiently large so that each component of x0 all lies within the interval [1/m0,m0]. For each integer m≥m0, define the stopping time:
(4)τm=inf{1mt∈[0,τe):min{x1(t),x2(t),x3(t)}≤1mormax{x1(t),x2(t),x3(t)}≥m1m}.
Throughout this paper, we set inf∅=∞ (as usual ∅ denotes the empty set). Clearly, τm is increasing as m→∞. Set τ∞=limm→∞τm; then τ∞≤τe a.s. If we can show that τ∞=∞ a.s., then τe=∞ and (x1(t),x2(t),x3(t))∈R+3 a.s. for all t≥0. In other words, to complete the proof all we need to show is that τ∞=∞ a.s. For if this statement is false, then there is a pair of constants T>0 and ϵ∈(0,1) such that
(5)P{τ∞≤T}>ϵ.
Hence there is an integer m1≥m0 such that
(6)P{τm≤T}≥ϵ∀m≥m1.
Define a C2-function V:R+3→R-+ by
(7)V(x1,x2,x3)=b32[b21(x1-1-logx1)+b12(x2-1-logx2)]+b23b12(x3-1-logx3);
the nonnegativity of this function can be seen from u-1-(1/2)logu≥0, for all u>0. Using Itô's formula, we get
(8)dV∶=LVdt+b32b21σ1(x1-1)dB1(t)+b32b12σ2(x2-1)dB2(t)+b23b12σ3(x3-1)dB3(t),
where
(9)LV=b32b21(x1-1)(a1-b11x1-b12x2)+b32b12(x2-1)(-a2+b21x1-b22x2-b23x3)+b23b12(x3-1)(-a3+b32x2-b33x3)+b32b21σ122+b32b12σ222+b23b12σ122=b32b21(-a1+σ122)+b32b12(a2+σ222)+b23b12(a3+σ322)+(b32b21a1+b32b21b11-b32b12b21)x1+(b32b21b12+b32b12b22-b23b12b32-b32b12a2)x2+(b32b12b23+b23b12b33-b23b12a3)x3-b32b21b11x12-b32b12b22x22-b23b12b33x32≤M^,
where M^ is a constant. Therefore
(10)∫0τm∧TdV(x(t))≤∫0τm∧TM^dt+∫0τm∧Tb32b21σ1(x1-1)dB1(t)+∫0τm∧Tb32b12σ2(x2-1)dB2(t)+∫0τm∧Tb23b12σ3(x3-1)dB3(t),
which implies that
(11)E[V(x(τm∧T))]≤V(x0)+M^T.
Set Ωm={τm≤T} for m≥m1. By (6), we know P(Ωm)≥ϵ. Notice that for every ω∈Ωm, there is at least one of xi(τm,ω) equals either m or 1/m; then
(12)V(x(τm))≥(m-1-logm)∧(m-1-1+logm).
It then follows from (11) that
(13)V(x0)+M^T≥E[1Ωm(ω)V(x(τm))]≥ϵ(m-1-logm)∧(m-1-1+logm),
where 1Ωm(ω) is the indicator function of Ωm. Letting m→∞ leads to the contradiction that ∞>V(x0)+M^T=∞. So we must have τ∞=∞ a.s.
Theorem 2.
Let x(t)=(x1(t),x2(t), and let x3(t)) be the solution of system (3) with any given initial value x(0)=x0∈R+3, then there exists a positive constant K(p) such that
(14)E[(b32b21x1(t)+b32b12x2(t)+b23b12x3(t))p]≤K(p)iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii∀t∈(0,∞),p>1.
Proof.
Let y(t)=b32b21x1(t)+b32b12x2(t)+b23b12x3(t); then
(15)dy(t)=(x32a1b21b32x1-a2b12b32x2-a3b12b23x3-b11b21b32x12-b12b22b32x22-b12b23b33x32)dt+σ1b21b32x1dB1(t)+σ2b12b32x2dB2(t)+σ3b12b23x3dB3(t),
and so
(16)dyp=pyp-1(x32a1b21b32x1-a2b12b32x2-a3b12b23x3-b11b21b32x12-b12b22b32x22-b12b23b33x32)dt+pyp-1(σ1b21b32x1dB1(t)+σ2b12b32x2dB2(t)+σ3b12b23x3dB3(t))+12p(p-1)yp-2×(σ12b212b322x12+σ22b122b322x22+σ32b122b232x32)dt.
Note that
(17)b11b21b32x12+b12b22b32x22+b12b23b33x32≥min{b11,b22,b33}b32b21+b32b12+b23b12×(b32b21x1+b32b12x2+b23b12x3)2=min{b11,b22,b33}b32b21+b32b12+b23b12y2.
Then
(18)dyp=pyp-1(a1y-min{b11,b22,b33}b32b21+b32b12+b23b12y2)dt+12p(p-1)yp-2max{σ12,σ22,σ32}y2dt+pyp-1(σ1b21b32x1dB1(t)+σ2b12b32x2dB2(t)+σ3b12b23x3dB3(t))≤[pmin{b11,b22,b33}b32b21+b32b12+b23b12p(a1+p2max{σ12,σ22,σ32})yp-pmin{b11,b22,b33}b32b21+b32b12+b23b12yp+1]dt+pyp-1(σ1b21b32x1dB1(t)+σ2b12b32x2dB2(t)+σ3b12b23x3dB3(t)).
Hence
(19)dE[yp(t)]dt≤p(a1+pmax{σ12,σ22,σ32}2)E[yp(t)]-pmin{b11,b22,b33}b32b21+b32b12+b23b12E[yp+1(t)]≤p(a1+pmax{σ12,σ22,σ32}2)E[yp(t)]-pmin{b11,b22,b33}b32b21+b32b12+b23b12E[yp(t)](p+1)/p.
Therefore, by comparison theorem, we get
(20)limsupt→∞E[yp(t)]≤[(a1+pmax{σ12,σ22,σ32})(b32b21+b32b12+b23b12)min{b11,b22,b33}]p.
Besides, note that E[yp(t)] is continuous; then there is a positive constant K(p) such that
(21)E[yp(t)]≤K(p),∀t∈[0,∞).
3. Persistent in Time Average
There is no equilibrium of system (3). Hence we cannot show the permanence of the system by proving the stability of the positive equilibrium as the deterministic system. In this section we first show that this system is persistent in mean. Before we give the result, we should do some prepared work.
L. S. Chen and J. Chen in [20] proposed the definition of persistence in mean for the deterministic system. Here, we also use this definition for the stochastic system.
Definition 3.
System (3) is said to be persistent in mean, if
(22)liminft→∞1t∫0tx3(s)ds>0,a.s.
Lemma 4 (see [21, Lemma 17]).
Let f∈C([0,+∞)×Ω,(0,+∞)) and F∈C([0,+∞)×Ω,R). If there exist positive constants λ0, λ, such that
(23)logf(t)≥λt-λ0∫0tf(s)ds+F(t),t≥0a.s.,
and limt→∞(F(t)/t)=0a.s., then
(24)liminft→∞1t∫0tf(s)ds≥λλ0,a.s.
From Lemma 4, it is easy to see that we could get Lemmas 5 and 6 with the same method.
Lemma 5.
Let f∈C([0,+∞)×Ω,(0,+∞)) and F∈C([0,+∞)×Ω,R). If there exist positive constants λ0, λ, such that
(25)logf(t)≤λt-λ0∫0tf(s)ds+F(t),t≥0a.s.,
and limt→∞(F(t)/t)=0a.s., then
(26)limsupt→∞1t∫0tf(s)ds≤λλ0,a.s.
Lemma 6.
Let f∈C([0,+∞)×Ω,(0,+∞)) and F∈C([0,+∞)×Ω,R). If there exist positive constants λ0, λ, such that
(27)logf(t)=λt-λ0∫0tf(s)ds+F(t),t≥0a.s.,
and limt→∞(F(t)/t)=0a.s., then
(28)limt→∞1t∫0tf(s)ds=λλ0,a.s.
From the stochastic comparison theorem [11], it is easy to get the following result.
Lemma 7.
Let x(t)∈R+3 be a solution of system (3) with x(0)=(x1(0),x2(0),x3(0)). Then one has
(29)x(t)≤Φ(t);
that is,
(30)xi(t)≤Φi(t),i=1,2,3,
where
(31)Φ(t)=(Φ1(t),Φ2(t),Φ3(t))⊤,Φi(t) is solutions of the following stochastic differential equations:
(32)dΦ1(t)=Φ1(t)(a1-b11Φ1(t))dt+σ1Φ1(t)dB1(t),Φ1(0)=x1(0),dΦ2(t)=Φ2(t)(-a2+b21Φ1(t)-b22Φ2(t))dt-σ2Φ2(t)dB2(t),Φ2(0)=x2(0),dΦ3(t)=Φ3(t)(-a3+b32Φ2(t)-b33Φ3(t))dt-σ3Φ3(t)dB3(t),Φ3(0)=x3(0).
If Assumption 8 is satisfied, the solution Φ(t) of system (32) with any initial value Φ(0)∈R+3 has the following property:
(34)limt→∞logΦi(t)t=0,limt→∞1t∫0tΦi(s)ds=Mi,a.s.,
where
(35)M1=r1b11,M2=r1b21-r2b11b11,M3=r1b21b32-r2b11b32-r3b11b22b11b22b33.
Proof.
From the result in [14] and Assumption 8 being satisfied, we know
(36)limt→∞logΦ1(t)t=0,liminft→∞1t∫0tΦ1(s)ds=a1-σ12/2b11=r1b11,
Besides, according to Itô's formula, the second population of system (32) is changed into
(37)dlogΦ2(t)=(-r2+b21Φ1(t)-b22Φ2(t))dt-σ2dB2(t).
It then follows
(38)logΦ2(t)=logΦ2(0)-r2t+b21∫0tΦ1(s)ds-b22∫0tΦ2(s)ds-σ2B2(t),
With Lemma 6 and Assumption 8, we could get
(39)limt→∞1t∫0tΦ2(s)ds=-r2+b21(r1/b11)b22=r1b21-r2b11b11b22>0.
Let (38) divide t, and t→∞, together with (36) and (39), consequently
(40)limt→∞logΦ2(t)t=0.
Similarly, according to Itô's formula, the third population of system (25) is changed into
(41)dlogΦ3(t)=(-r3+b32Φ2(t)-b33Φ3(t))dt-σ3dB3(t);
it then follows
(42)logΦ3(t)=logΦ3(0)-r3t+b32∫0tΦ2(s)ds-b33∫0tΦ3(s)ds-σ3B3(t),limt→∞1t∫0tΦ3(s)ds=-r3+b32((r1b21-r2b11)/b11b22)b33>0,limt→∞logΦ3(t)t=0.
From this, together with Lemmas 7 and 9, the following result is obviously true.
Theorem 10.
If Assumption 8 is satisfied, the solution x(t) of system (3) with any initial value x(0)∈R+3 has the following property:
(43)limsupt→∞logxi(t)t≤0,i=1,2,3.
Above all, we could get.
Theorem 11.
If Assumption 8 is satisfied, the the solution x(t) of system (3) with any initial value x(0)∈R+3 has the following property:
(44)liminft→∞1t∫0tx3(s)ds≥x~3*,a.s.,
where x~*=(x~1*,x~2*,x~3*) is the only nonnegative solution of the following equation:
(45)r1-b11x1-b12x2=0,-r2+b21x1-b22x2-b23x3=0,-r3+b32x2-b33x3=0.
Proof.
From system (3), such that
(46)d(c1logx1(t)+c2logx2(t)+c3logx3(t))=[(r1c1-r2c2-r3c3)+(-b11c1+b21c2)x1+(-b12c1-b22c2+b32c3)x2-(b23c2+b33c3)x3]dt+c1σ1dB1(t)-c2σ2dB2(t)-c3σ3dB3(t).
Let c1=b21, c2=b11, and c3=(b11b22+b12b21)/b32, together with Assumption 8, we know
(47)r1c1-r2c2-r3c3>0;
hence
(48)(c1(logx1(t)-logx1(0))+c2(logx2(t)-logx2(0))+c3(logx3(t)-logx3(0)))×(t)-1=(r1c1-r2c2-r3c3)-(c2b23+c3b33)1t∫0tx3(s)ds+c1σ1B1(t)-c2σ2B2(t)-c3σ3B3(t)t.
According to Theorem 10, where
(49)limsupt→∞logxi(t)t≤0,i=1,2,3,
and limt→∞(Bi(t)/t)=0, i=1,2,3,
(50)liminft→∞1t∫0tx3(s)ds≥r1c1-r2c2-r3c3c2b23+c3b33=x~3*,
where x~*=(x~1*,x~2*,x~3*) is the only nonnegative solution of the following equation when Assumption 8 is satisfied:
(51)r1-b11x1-b12x2=0,-r2+b21x1-b22x2-b23x3=0,-r3+b32x2-b33x3=0.
4. Nonpersistence
In this section, we show the situation when the population of system (3) will be extinction in three cases.
Case 1 (r1<0).
According to Itô's formula, the first population of system (25) is changed into
(52)dlogΦ1(t)≤(r1-b11Φ1(t))dt-σ1dB1(t).
If r1<0, we could get
(53)limsupt→∞logΦ1(t)t≤r1b11<0a.s.
From the stochastic comparison theorem, we have
(54)limsupt→∞logx1(t)t≤r1b11<0a.s.,
hence
(55)limt→∞x1(t)=0,a.s.
From the second population of system (25), we have
(56)limsupt→∞logΦ2(t)t≤-a2+b21limsupt→∞1t∫0tΦ1(s)ds≤-a2a.s.;
similarly
(57)limsupt→∞logΦ3(t)t≤-a3a.s.,limt→∞xi(t)=0a.s.i=2,3.
Case 2 (r1>0,r1-(b11/b21)r2<0).
It is clear that from the proof section of Case 1, we get
(58)logΦ2(t)-logΦ2(0)t≤-r2+b211t∫0tΦ1(s)ds-σ2dB2(t)ta.s.,
hence
(59)limsupt→∞logΦ2(t)t≤-r2+b21M1=-r2+b21r1b11<0a.s.
Similarly
(60)limsupt→∞logΦ3(t)t≤-r3+b32limsupt→∞1t∫0tΦ2(s)ds≤-a3<0a.s.;
thus,
(61)limt→∞xi(t)=0a.s.,i=2,3.
Above all, and from the conclusion in [22], we could easily know that the distribution of x1(t) converges weekly to the probability measure with density:
(62)f*(ζ)=C0ζ2r1/σ12-1e-2b11ζ/σ12,
where C0=(2b11/σ12)2r1/σ12/Γ(2r1/σ12) and
(63)limt→∞1t∫0tx1(s)ds=r1b11,a.s.
Case 3 (r1-(b11/b21)r2-((b11b22+b12b21)/b21b32)r3<0).
It is clear that
(64)d(c1logx1(t)+c2logx2(t)+c3logx3(t))=[(r1c1-r2c2-r3c3)+(-b11c1+b21c2)x1+(-b12c1-b22c2+b32c3)x2-(b23c2+b33c3)x3]dt+c1σ1dB1(t)-c2σ2dB2(t)-c3σ3dB3(t).
Since c1=b21, c2=b11, c3=(b11b22+b12b21)/b32, we get
(65)c1logx1(t)+c2logx2(t)+c3logx3(t)≤(r1c1-r2c2-r3c3)t+c1logx1(0)+c2logx2(0)+c3logx3(0)+c1σ1B1(t)-c2σ2B2(t)-c3σ3B3(t),
Moreover,
(66)logx1c1(t)x2c2(t)x3c3(t)t≤(r1c1-r2c2-r3c3)+c1logx1(0)+c2logx2(0)+c3logx3(0)t+c1σ1B1(t)t-c2σ2B2(t)t-c3σ3B3(t)t.
And limt→∞(Bi(t)/t)=0, i=1,2,3, implies
(67)limsupt→∞logx1c1(t)x2c2(t)x3c3(t)t≤r1c1-r2c2-r3c3<0;
then
(68)limt→∞x1c1(t)x2c2(t)x3c3(t)=0a.s.
Therefore, by the above arguments, we get the following conclusion.
Theorem 12.
Let x(t) be the solution of system (3) with any initial value x(0)∈R+3. Then
(1) if r1<0, then
(69)limt→∞xi(t)=0a.s.,i=1,2,3,
(2) if r1>0, r1-(b11/b21)r2<0, then
(70)limt→∞xi(t)=0a.s.,i=2,3,
and the distribution of x1(t) converges weekly to the probability measure with density:
(71)f*(ζ)=C0ζ2r1/σ12-1e-2b11ζ/σ12,
where C0=(2b11/σ12)2r1/σ12/Γ(2r1/σ12) and
(72)limt→∞1t∫0tx1(s)ds=r1b11,a.s;
(3) if r1-(b11/b21)r2-((b11b22+b12b21)/b21b32)r3<0, then
(73)limt→∞x1c1(t)x2c2(t)x3c3(t)=0,a.s.,
where c1=b21, c2=b11, and c3=(b11b22+b12b21)/b32.
5. Numerical Simulation
In this section, we give out the numerical experiment to support our results. Consider
(74)x˙1(t)=x1(t)(a1-b11x1(t)-b12x2(t))+σ1x1(t)B˙1(t),x˙2(t)=x2(t)(-a2+b21x1(t)-b22x2(t)-b23x3(t))-σ2x2(t)B˙2(t),x˙3(t)=x3(t)(-a3+b32x2(t)-b33x3(t))-σ3x3(t)B˙3(t).
By the Milstein method in [23], we have the difference equation:
(75)x1,k+1=x1,k+x1,k×[σ122(a1-b11x1,k-b12x2,k)Δt+σ1ϵ1,kΔt+σ122(ϵ1,k2Δt-Δt)],x2,k+1=x2,k+x2,k×[σ122(-a2+b21x1,k-b22x2,k-b23x3,k)Δt-σ2ϵ2,kΔt+σ222(ϵ2,k2Δt-Δt)],x3,k+1=x3,k+x3,k×[σ122(-a3+b32x2,k-b33x3,k)Δt-σ3ϵ3,kΔt+σ322(ϵ3,k2Δt-Δt)],
where ϵ1,k, ϵ2,k, and ϵ3,k, i=1,2,3, are the Gaussian random variables N(0,1), r1=a1-σ12/2>0, and ri=ai+σi2/2, i=2,3. Choosing (x1(0),x2(0),x3(0))∈R+3, and suitable parameters, by Matlab, we get Figures 1, 2, and 3.
The solution of system (1) and system (3) with (x1(0),x2(0),x3(0))=(1,0.8,0.5), a1=0.3, a2=0.4, a3=0.1, b11=0.1, b12=0.1, b21=0.6, b22=0.6, b23=0.6, b32=0.8, and b33=0.6. The red lines represent the solution of system (1), while the blue lines represent the solution of system (3) with σ1=0.02, σ2=0.01, and σ3=0.01.
Two of the species will die out in probability. The red lines represent the solution of system (1), while the blue lines represent the solution of system (3) with (x1(0),x2(0),x3(0))=(1,0.8,0.5), a1=0.4, a2=0.4, a3=0.1, b11=0.1, b12=0.1, b21=0.6, b22=0.6, b23=0.6, b32=0.8, and b33=0.6. The red lines represent the solution of system (1), while the blue lines represent the solution of system (3) with σ1=0.02, σ2=3, and σ3=0.01.
One of the species or both species will die out in probability. The red lines represent the solution of system (1), while the blue lines represent the solution of system (3) with (x1(0),x2(0),x3(0))=(1,0.8,0.5), a1=-0.1, a2=0.4, a3=0.1, b11=0.1, b12=0.1, b21=0.6, b22=0.6, b23=0.6, b32=0.8, and b33=0.6. The red lines represent the solution of system (1), while the blue lines represent the solution of system (3) with σ1=0.02, σ2=3, and σ3=0.01.
In Figure 1, when the noise is small, choosing parameters satisfying the condition of Theorem 10, the solution of system (3) will persist in time average.
In Figure 2, we observe case (3) in Theorem 12 and choose parameters r1>0, r1-(b11/b21)r2<0. As Theorem 12 indicated that two predators will die out in probability. The prey solution of system (3) will persist in time average.
In Figure 3, we observe case (1) in Theorem 12 and choose parameters r1<0. As Theorem 12 indicated that not only predators but also prey will die out in probability when the noise of the prey is large, and it does not happen in the deterministic system.
Acknowledgments
The work was supported by the Ministry of Education of China (no. 109051), the Ph.D. Programs Foundation of Ministry of China (no. 200918), and NSFC of China (no. 11371085), Natural Science Foundation of Jilin Province of China (no. 201115133).
GohB. S.Global stability in many species system1997111135143FreedmanH. I.198057New York, NY, USAMarcel Dekkerx+254MR586941KuangY.SmithH. L.Global stability for infinite delay Lotka-Volterra type systems1993103222124610.1006/jdeq.1993.1048MR1221904ZBL0786.34077GopalsamyK.Global asymptotic stability in a periodic Lotka-Volterra system198524160170GardT. C.1988114New York, NY, USAMarcel Dekkerxiv+234MR917064GardT. C.Persistence in stochastic food web models198446335737010.1016/S0092-8240(84)80044-0MR748544ZBL0533.92028GardT. C.Stability for multispecies population models in random environments198610121411141910.1016/0362-546X(86)90111-2MR869549ZBL0598.92017MaoX.MarionG.RenshawE.Environmental Brownian noise suppresses explosions in population dynamics20029719511010.1016/S0304-4149(01)00126-0MR1870962ZBL1058.60046MaoX.Delay population dynamics and environmental noise20055214916210.1142/S021949370500133XMR2147279ZBL1093.60033HuY.WuF.HuangC.Stochastic Lotka-Volterra models with multiple delays20113751425710.1016/j.jmaa.2010.08.017MR2735693ZBL1245.92063WantanabeN. I.1981Amsterdam, The NetherlandsNorth-HollandJiC.JiangD.ShiN.O'ReganD.Existence, uniqueness, stochastic persistence and global stability of positive solutions of the logistic equation with random perturbation2007301778910.1002/mma.778MR2277928ZBL1148.34040ZhuC.YinG.On competitive Lotka-Volterra model in random environments2009357115417010.1016/j.jmaa.2009.03.066MR2526815ZBL1182.34078JiC.JiangD.ShiN.Analysis of a predator-prey model with modified Leslie-Gower and Holling-type II schemes with stochastic perturbation2009359248249810.1016/j.jmaa.2009.05.039MR2546763ZBL1190.34064JiC.JiangD.LiX.Qualitative analysis of a stochastic ratio-dependent predator-prey system201123551326134110.1016/j.cam.2010.08.021MR2728069ZBL1229.92076MayR. M.2001Princeton, NJ, USAPrinceton University PressPolanskyP.Invariant distributions for multipopulation models in random environments1979161253410.1016/0040-5809(79)90004-2MR545478BarraM.GrossoD. G.GerardiA.KochG.MarchettiF.1979Berlin, GermanySpringerMaoX. R.1997Chichester, UKHorwoodChenL. S.ChenJ.1993Beijing, USAScience PressXiaP. Y.ZhengX. K.JiangD. Q.Persistence and nonpersistence of a nonautonomous stochastic mutualism system2013201313256249JiC. Y.JiangD. Q.Dynamics of a stochastic density dependent predator-prey system with Beddington-DeAngelis functional response2011381144145310.1016/j.jmaa.2011.02.037MR2796222ZBL1232.34072HighamD. J.An algorithmic introduction to numerical simulation of stochastic differential equations200143352554610.1137/S0036144500378302MR1872387ZBL0979.65007