Analysis of a Delayed Stochastic One-Predator Two-Prey Population Model in a Polluted Environment ()
1. Introduction
With the rapid development of economy, environmental pollution has gradually become the major social problem today. With a growing number of toxicant and contaminants entering into the ecosystem, the quality of our living environment has declined. Then many species have been extinct, and some of them are on the edge of extinction. Therefore, controlling environment pollution has become a major topic in many countries, which draws researchers to investigate the influence of environment pollution.
In the 1980s, Hallam et al. [1] [2] [3] firstly proposed the deterministic models to study the impact of environment toxicant on the survival of biological population. Their studies have provided useful bases about protecting species for us. However, population system is often affected by environmental noise, and there are many scholars who have studied the dynamics of stochastic models with toxicant [4] [5] [6] [7] .
On the other hand, more realistic models of population interactions should take the effects of time delay into account [8] [9] . Further, in the natural world, it is a common phenomenon that a predator feeds on some competing preys [10] [11] [12] . However, there is little research on the delayed stochastic one-predator two-prey model in a polluted environment. Thus we consider a stochastic delayed one-predator two-prey model with toxicant input in this article.
The rest is organized as follows. In Section 2, we show some notations and introduce a stochastic delayed one-predator two-prey model in polluted environment. In Section 3, we show that the system (SM) has a unique global positive solution. In Section 4, we give the main theorems and their proof. In Section 5, the attractively global system is investigated. In Section 6, we present numerical simulations to illustrate our mathematical findings.
2. The Model and Notations
In this section, we will give some notations on stochastic one-predator-two-prey system. The stochastic predator-prey system in a polluted environment takes the following form:
(2.1)
The above model does not incorporate the effect of time delay, but for a long time, it has been recognized that delays can have a complex effect on the dynamics of a system [9] [13] . In the same time, the natural growth of many populations is inevitably affected by many random disturbances. Considering the effects of random disturbances, we assume the growth rate of prey and the death rate of predator are perturbed with
where
is mutually independent one-dimensional standard Brownian motions with
and
being the intensities of white noises. Stochastic version corresponding to deterministic system with time delays can be rewritten as:
(2.2)
In past few decades, delay population systems with one predator and two competing preys have received great attention and have been investigated widely. However, as far as the authors concerned, no one has yet explored the predator-prey system with time delays and toxicant inputs in the same time. Therefore, on the basis of article [14] , we establish the following delayed stochastic one-predator two-prey model in a polluted environment:
(2.3)
with initial data
where
is the size of the the prey i,
, and
is the size of the predator;
is the growth rate of the i the species,
,
is the death rate of the predator;
is the intra-specific competition rate,
.
and
stand for the inter-specific competition rates between species 1 and 2,
and
stand for the capture rates,
and
are the efficiency of food conversion.
and
denote the concentrations of the toxicant in the organism of species and the environment at time t, respectively.
stand for dose-response of the prey and predator to the organismal toxicant, and
and
denote the excretion and depuration rates of the toxicant,
, respectively.
and
represent the absorption of toxicant per unit of mass by the environment and by food, respectively.
is the concentration of toxicant in the environment;
is the uptake rate of food per unit mass. Parameter h reflects the ability of the environment to clean up toxicant.
denotes the exogenous rate of toxicant input into the environment and it is supposed to be bounded and
. All coefficients mentioned above are positive constants.
represents the time delay.
, where
represents the space of all the continued functions from
to
.
Although the model is a five-dimensional system, because the explicit solutions of the latter two equations are easy to get, it is actually only necessary to study the first three stochastic differential equations of the model, which is called model
in this paper.
For the sake of simplification, we define some notations:
For a function g, we denote the following notations:
3. Existence and Uniqueness of the Global Positive Solution
In order to make the model be sense, we need to show the solution is non-negative and global.
Lemma 3.1 ( [15] ) For model (2.3), if
,
, then
,
.
This paper assumes that condition
,
is always true in model (2.3), then the solution process of model (SM) should be non-negative.
Lemma 3.2 For any initial value
, there is a unique global positive solution
. Moreover, there is a positive constants k such that
(3.1)
Proof. Consider the following system:
(3.2)
with initial value
. Since the coefficients of (3.2) obey the local Lipstchiz condition, then (3.2) has a unique local positive solution
on
, where
stands for the explosion time. Hence it follows from Itô’s formula that (SM) has the following unique positive global solution
Now we show that
is global, i.e.,
. Consider the following system:
(3.3)
with initial value
. By the stochastic comparison theorem [16] , one can see that for
,
Thanks to Theorem 4.2 in [17] , system (3.3) can be explicitly solved as follows
Note that
and
are existent on
, hence
.
Before we state the main theorem of this paper, we need to introduce several hypotheses.
Hypothesis 1.
. which imply that all the populations coexist if model (SM) frees from stochastic noises.
Hypothesis 2.
.
4. Permanence in Time Average
In this section, we study the permanent in time average of systems (2.3) and (SM). We firstly do some preparation.
Definition 4.1. System (2.3) is said to be permanent in time average if there are positive constants
and
(
) such that
holds for any solution
of system (2.3) with initial condition
Lemma 4.1 ( [18] ). Suppose
be a continuous function from
to
.
1) If there exist constants
and
such that
for
, where
are independent standard Brownian motions and
are constants,
, then we have:
2) If there exist positive constants
and T such that
for all
, where
are independent standard Brownian motions and
are constants,
, then if
, we have
.
Lemma 4.2 For arbitrary
we have
(4.1)
Proof. Consider the first two equations in model (3.3). By [19] , we have
Consequently, if
, then
If
, then
We now consider the following stochastic equation with delays
(4.2)
Lemma 4.3 Assume that
and
, then the solution of (4.2) has the following properties:
Proof. The solution of (4.2) has the property that ( [4] )
(4.3)
By Itô’s formula, we get
(4.4)
(4.5)
(4.6)
Dividing both sides of (4.4), (4.5) and (4.6) by t, one can see that
(4.7)
(4.8)
(4.9)
And it is well known that
(4.10)
Utilizing Lemma 4.1,
and
, it is easy to derive that
(4.11)
And
(4.12)
Plugging (4.11), (4.12) into (4.7), (4.8), respectively, then combining (4.10) leads to
(4.13)
and
(4.14)
Besides, computing
, one can derive that
(4.15)
In view of (4.1), (4.13), (4.14), (4.15) and Lemma 4.1, we can gain that
(4.16)
By (3.11)-(3.16), we have
This completes the proof.
Lemma 4.4 For arbitrarily
,
(4.17)
Proof. From Lemma 4.3, we can see that either
or
(a constant).
If
, then
If
, then it is easy to see that
This completes the proof.
Lemma 4.5 For any initial data
, the solution
of system (SM) has the property that
for
.
where
is the solution of system (4.2).
Proof. Let
. Then, by Itô’s formula, we have
This is
Then
By Lemma 3.2, we obtain that
is the solution of the following equation
Hence, we have
In the same way, we can get
On the other hand, let
. Then, by Itô’s formula, we derive that
We then have
Therefore
By Lemma 3.2, we know that
is a solution of the following equation
Hence the proof is completed.
Lemma 4.6 The solution
of (SM) satisfies
(4.18)
Proof. We know that either
or
constant,
. If
, then
If
constant, similar to the proof of [20] , we have
.
Hence
. The proof is completed.
Based on the above discussion, now we show the main results in this sectionm
Theorem 4.1 Assume that
and
, then for any initial data
, the solution
of (SM) has the following properties:
Proof. Let
be the solution of system (4.2). By Lemma 4.5, we obtain
Since
Then
(4.19)
Similarly,
(4.20)
It can be straightly shown by Lemma 4.3 that
(4.21)
Applying Itô’s formula to (SM) results in
(4.22)
(4.23)
(4.24)
Dividing both sides of (4.22), (4.23) and (4.24) by time t, one can obtain that
(4.25)
(4.26)
(4.27)
Let m and n satisfy the following system:
Then,
where
represents the complement minor of
in the determinant
.
After a simple calculation of
, we get
(4.28)
Substituting(4.18), (4.17), (4.1) into (4.28) gives
(4.29)
For sufficiently large t. By
, Lemma 4.1 and the arbitrariness of
, we have
(4.30)
Analogously, let
and
satisfy the following system
Then,
Multiplying both sides of (4.25), (4.26) and (4.27) by
, (−1) and
, respectively, then adding these three equations, we get
(4.31)
Plugging (4.1), (4.17), (4.30), (4.31) into (4.27) results in
(4.32)
Noting that
. Then an application of Lemma 4.1 to (4.32) yields that
(4.33)
Substituting (4.1), (4.17), (4.31) and (4.33) into (4.25) elicits that
for sufficiently large t. Hence, we can further get from Lemma 4.1 and the arbitrariness of
that
(4.34)
In the same way, we can show that
(4.35)
Substituting (4.34) and (4.35) into (4.27), and then using (4.1) and (4.17), Lemma 4.1 and the arbitrariness of
, we get
This together with (4.30), (4.31), (4.33), (4.34) and (4.35), yields
This completes the proof.
As to model (2.3), by utilizing similar techniques to those employed in the proof of Theorem 4.1, one can obtain the similar result.
5. Global Attractivity
Theorem 5.1 Let Hypothesis 2 hold, then system (SM) is globally attractive.
Proof. Let
be the cofactor of the ith diagonal element of
, where
Making use of Kirchhoff’s Matrix Tree Theorem (See, e.g., [21] ), one has
.
Define
Reckoning the right differential
deduces that
According to Theorem 2.3 in [22] , we have
As a consequence,
Therefore
(5.1)
According to (SM),
Clearly,
is differentiable. By virtue of (3.1),
where
is a constant. Consequently,
is uniformly continuous. In the same way,
and
are also uniformly continuous functions. By virtue of (5.1) and Barbalat’s conclusion [23] , we derive the required assertion.
6. Numerical Simulations
Now we introduce some numerical figures to support Theorem 4.1 by using the Euler scheme [24] .
Considering the parameters as following:
(6.1)
It is easy to get that:
Set
and
. Then by calculation, we have
,
,
and
. Furthermore,
. Therefore all conditions in Theorem 4.1 have been checked. Then we have
see Figure 1, which is obtained by applying the Milstein method [25] .
Figure 1 shows the simulations of the solutions of systems (SM), besides,
and
are shown in Figure 1. We see that
tend to constants, which is consistent with the results of Theorem 4.1. the solutions of systems (SM) fluctuate around a small zone. Thus, we think the system (SM) is permanent.
Since the parameters given above meet the hypothesis 2:
,
and
. According to Theorem 5.1, we can get the system (SM) is global attractively, see Figure 2.
Figure 2 shows the simulations of the solutions of systems (SM), From Figure 2(a) and Figure 2(b), we can see that the solution of the system is globally attractively, whether with or without random perturbations.
7. Conclusions and Discussions
The dynamic relationship between predator and their preys has been and will continue to be one of the major themes in ecology due to its importance and
Figure 1. System (6.1) with
and
. The red line and purplish red represent two prey population
, respectively, the pink line represents the predator population
. The yellow line represents
, the green line represent
and the blue line represent
. Figure 1(a) and Figure 1(b) represent system with and without random perturbations, respectively. (For interpretation of the references to color in this figure caption, the reader is referred to the web version of this article.)
Figure 2. System (6.1) with
,
,
. Time series diagram of
. The green line represents one prey population
The purplish red line represents another prey population
and the blue line represents one predator population
. Figure 2(a) and Figure 2(b) represent system with and without random perturbations, respectively.
universal existence. This paper is concerned with a delayed stochastic one-predator-two-prey population model in a polluted environment. Firstly, we show that there exists a unique positive solution in our system. Secondly, the permanence of this system is investigated. Conditions for the system to be permanent in time average are given. Our main result in part 4 reveals the impacts of stochastic perturbations on the persistence and extinction of every species. Finally, our results are confirmed by numerical simulation.
Some questions deserve further explorations. In the first place, it is significant to study the delay population system with other disturbance, such as Lévy jumps, or Markovian switching. Another problem is to consider population models with different functional responses, such as Holling II-IV type and Beddington-DeAngelis functional response. We leave these investigations for future work.
Acknowledgements
This work was supported by the National Natural Science Foundation of China (11861027). This work was supported by the Natural Science Foundation of Guangxi (2016GXNSFAA380194).