Next Article in Journal
Quadruple Integral Involving the Logarithm and Product of Bessel Functions Expressed in Terms of the Lerch Function
Next Article in Special Issue
Backward Stochastic Differential Equations (BSDEs) Using Infinite-Dimensional Martingales with Subdifferential Operator
Previous Article in Journal
Quadratic Lyapunov Functions for Stability of the Generalized Proportional Fractional Differential Equations with Applications to Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of an Impulsive Stochastic Predator–Prey System with the Beddington–DeAngelis Functional Response

College of Science, Guilin University of Technology, Guilin 541004, China
Axioms 2021, 10(4), 323; https://doi.org/10.3390/axioms10040323
Submission received: 28 October 2021 / Revised: 17 November 2021 / Accepted: 21 November 2021 / Published: 28 November 2021
(This article belongs to the Special Issue Recent Advances in Stochastic Differential Equations)

Abstract

:
Taking impulsive effects into account, an impulsive stochastic predator–prey system with the Beddington–DeAngelis functional response is proposed in this paper. First, the impulsive system is transformed into an equivalent system without pulses. Then, by constructing suitable functionals and applying the extreme-value theory of quadratic functions, sufficient conditions on the existence of periodic Markovian processes are provided. The uniform continuity and global attractivity of solutions are also investigated. Additionally, we investigate the extinction and permanence in the mean of all species with the help of comparison methods and inequality techniques. Sufficient conditions on the existence and ergodicity of the stationary distribution of solutions for the autonomous and non-impulsive case are given. Finally, numerical simulations are performed to illustrate the main results.

1. Introduction

Predator–prey systems with the Beddington–DeAngelis functional response have always attracted the attention of many researchers (to name a few, see [1,2,3] and the references therein). As we know, white noises are very common in the real world. They affect the growth rates of species to some extent. In order to understand the effects of white noises on the dynamics of ecological systems, random disturbance is introduced [4,5,6,7,8,9]. For example, Yagi and Ton [9] studied the following model:
d x 1 ( t ) = x 1 ( t ) a 11 ( t ) a 12 ( t ) x 1 ( t ) c 1 ( t ) x 2 ( t ) m 0 ( t ) + m 1 ( t ) x 1 ( t ) + m 2 ( t ) x 2 ( t ) d t + σ 1 ( t ) x 1 ( t ) d ω 1 ( t ) , d x 2 ( t ) = x 2 ( t ) a 21 ( t ) a 22 ( t ) x 2 ( t ) + c 2 ( t ) x 1 ( t ) m 0 ( t ) + m 1 ( t ) x 1 ( t ) + m 2 ( t ) x 2 ( t ) d t + σ 2 ( t ) x 2 ( t ) d ω 2 ( t ) ,
where x 1 ( t ) and x 2 ( t ) represent the quantities of prey and predator populations, respectively; ω 1 ( t ) and ω 2 ( t ) are independent standard Brownian motions defined on the complete probability space ( Ω , F , { F t } t 0 , P ) with a filtration { F t } t 0 satisfying the usual conditions—that is, it is increasing and right continuous and F 0 contains all P-null set. See [1,3] for the biological meanings of the parameter functions.
Moreover, the growths of species usually suffer rapid changes in relatively short time periods due to activities such as stocking and harvesting. To describe such phenomena, impulsive dynamical systems have been proposed and studied (e.g., [10,11,12,13,14,15] and references therein).
Furthermore, the life cycles of species and the variation of the environment play important roles in many biological and ecological dynamical systems. In particular, the effects of a periodically varying environment are important for evolutionary theory, as the selective forces on systems in a fluctuating environment differ from those in a stable environment. Therefore, the assumption of periodicity on parameters in the system (in a way) incorporates the periodicity of the environment [16,17,18].
Motivated by the above discussion and based on (1), we propose the following periodic impulsive stochastic predator–prey model with the Beddington–DeAngelis functional response:
d x 1 ( t ) = x 1 ( t ) a 11 ( t ) a 12 ( t ) x 1 ( t ) c 1 ( t ) x 2 ( t ) m 0 ( t ) + m 1 ( t ) x 1 ( t ) + m 2 ( t ) x 2 ( t ) d t + σ 1 ( t ) x 1 ( t ) d ω 1 ( t ) , d x 2 ( t ) = x 2 ( t ) a 21 ( t ) a 22 ( t ) x 2 ( t ) + c 2 ( t ) x 1 ( t ) m 0 ( t ) + m 1 ( t ) x 1 ( t ) + m 2 ( t ) x 2 ( t ) d t + σ 2 ( t ) x 2 ( t ) d ω 2 ( t ) , t t k , x 1 ( t k + ) = ( 1 + ξ 1 k ) x 1 ( t k ) , x 2 ( t k + ) = ( 1 + ξ 2 k ) x 2 ( t k ) , t = t k .
In the sequel, we always assume that the following conditions are satisfied by the constants and parameter functions of (2).
  • All functions a 11 ( · ) , a 12 ( · ) , a 21 ( · ) , a 22 ( · ) , c 1 ( · ) , c 2 ( · ) , m 0 ( · ) , m 1 ( · ) , m 2 ( · ) , σ 1 ( · ) , and σ 2 ( · ) are positive, bounded, continuous, and periodic with the same period T.
  • The impulsive points satisfy 0 < t 1 < t 2 < with lim k t k = and there exists a positive integer q such that t k + q = t k + T and ξ i , k + q = ξ i , k for i = 1 , 2 and k N = { 1 , 2 , } .
  • By the biological meanings, we assume ξ i k > 1 for i = 1 , 2 and k N = { 1 , 2 , } .
The motivations for investigating the dynamics of (2) are as follows. Firstly, for deterministic systems, the existence and attractivity of equilibria/periodic solutions are important topics for understanding the dynamics of ecological systems [19,20,21,22,23]. Similarly, for stochastic systems, it is crucial to study the existence of periodic Markovian processes. Secondly, extinction, permanence in mean for all species, and stochastic persistence are all important properties to understand the dynamics of biological systems [24,25,26,27]. Lastly, for the special case of (2) with parameter functions being constants, it is interesting to study the existence of a stationary distribution as the functional response is non-monotonic.
The rest of this paper is organized as follows. Section 2 begins with definitions, important lemmas, and some notations. Section 3 provides sufficient conditions on the existence of periodic Markovian processes for (2). Section 4 focuses on the extinction and permanence in the mean of all species of (2). Section 6 is devoted to the existence and uniqueness of distribution. Some numerical simulations are given in Section 7. The paper concludes with a brief discussion and conclusions.

2. Preliminaries

We start with the definition of a periodic Markovian process.
Definition 1
([9,18]). A stochastic process ξ ( t ) = ξ ( t , ω ) ( < t < ) is said to be periodic with period T if for every finite sequence of numbers t 1 , t 2 , …, t n the joint distribution of the random variables ξ ( t 1 + h ) , …, ξ ( t n + h ) is independent of h, where h = k T with k = ± 1 , ± 2 , …
Definition 2
([13,15]). For the following impulsive stochastic system of differential equations:
d x ( t ) = f ( t , x ( t ) ) d t + g ( t , x ( t ) ) d ω ( t ) , t t k , x ( t k + ) = ( 1 + λ k ) x ( t k ) , k N ,
with the initial value x ( 0 ) = x 0 R n , a stochastic process x ( t ) = ( x 1 ( t ) , x 2 ( t ) , , x n ( t ) ) T , t R + = [ 0 , ) is said to be a solution if
(i) 
x ( t ) is f t -adapted and is continuous on ( 0 , t 1 ) and each interval ( t k , t k + 1 ) R + ,   k N , f ( t , x ( t ) ) L ( R + , R n ) ,   g ( t , x ( t ) ) L 2 ( R + , R n ) , where L p ( R + , R n ) is the set of all R n -valued measurable f t -adapted processes satisfying 0 T | f ( t ) | p d t < almost surely for all T > 0 , p = 1 , 2;
(ii) 
For every t k , k N , x ( t k + ) = lim t t k + x ( t ) and x ( t k ) = lim t t k x ( t ) exist, and x ( t k ) = x ( t k ) with the probability one;
(iii) 
For all t [ 0 , t 1 ] ,   x ( t ) obeys the integral equation
x ( t ) = x ( 0 ) + 0 t f ( s , x ( s ) ) d s + 0 t g ( s , x ( s ) ) d ω ( s )
and for all t [ t k , t k + 1 ] ,   k N ,   x ( t ) obeys the following integral equation:
x ( t ) = x ( t k + ) + t k t f ( s , x ( s ) ) d s + t k t g ( s , x ( s ) ) d ω ( s ) .
Definition 3
([27]). System (2) is globally attractive if
lim t | x 1 ( t ) x ˜ 1 ( t ) | = lim t | x 2 ( t ) x ˜ 2 ( t ) | = 0
for any two positive solutions ( x 1 ( t ) , x 2 ( t ) ) and ( x ˜ 1 ( t ) , x ˜ 2 ( t ) ) .
Lemma 1
([15]). For the following Itô’s differential equation:
d x ( t ) = b ( t , x ( t ) ) d t + σ ( t , x ( t ) ) d ω ( t ) ,
if all the coefficients are T-periodic in t and satisfy the linear growing condition and the Lipschitz condition in every cylinder U l × R + for l > 0 , where U l = { x : | | x | | l } and there exists a function v = v ( t , x ) which is twice continuously differentiable with respect to x and continuously differentiable with respect to t in R n × R + , T-periodic in t and satisfies the following conditions:
inf | | x | | > l v ( t , x ) a s l ,
L v ( t , x ) 1 o u t s i d e   s o m e   c o m p a c t   s e t ,
then there exists a solution of (4) which is a T-periodic Markovian process, where L is the generator of (4) defined by
L v ( t , x ) = v t + x v · b + 0.5 t r ( σ σ T v x x ) ( t , x ) .
Lemma 2
([13]). Suppose that Z C ( R + , ( 0 , ) ) and lim t F ( t ) / t = 0 a.s.
(a) 
If there exist two positive constants t 0 > 0 and λ 0 > 0 such that, for all t > t 0 ,
ln Z ( t ) λ t λ 0 0 t z ( s ) d s + F ( t ) a . s . ,
then
Z * λ / λ 0 a . s . i f λ 0 , lim t Z ( t ) = 0 a . s . i f λ < 0 .
(b) 
If there exist some constants t 0 > 0 , λ 0 > 0 , and λ such that, for all t > t 0 ,
ln Z ( t ) λ t λ 0 0 t z ( s ) d s + F ( t ) a . s . ,
then
Z * λ / λ 0 a . s .
Lemma 3
([28]). Suppose that an n-dimensional stochastic process x ( t ) on t 0 satisfies
E | x ( t ) x ( s ) | β 1 c | t s | 1 + β 2 , 0 s , t < ,
for some positive constants β 1 , β 2 , and c. Then there exists a continuous modification x ˜ ( t ) of x ( t ) , which has the property that, for every ϑ ( 0 , β 2 β 1 ) , there is a positive random variable h ( ω ) such that
P ω : sup 0 < | t s | < h ( ω ) , 0 s , t < | x ˜ ( t , ω ) x ( t , ω ) | | t s | ϑ 2 1 2 ϑ = 1 .
In other words, almost every sample path of x ˜ ( t ) is locally but uniformly Hölder continuous with exponent ϑ.
Lemma 4
([29]). Let f be a nonnegative function defined on R + such that f is integrable on R + and is uniformly continuous on R + , then lim t f ( t ) = 0 .
To investigate the dynamics of (2), we first consider the following non-impulsive system:
d y 1 ( t ) = y 1 ( t ) ( a 11 ( t ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 ( t ) A 1 ( t ) y 1 ( t ) c 1 ( t ) A 2 ( t ) y 2 ( t ) m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 ( t ) + m 2 ( t ) A 2 ( t ) y 2 ( t ) ) d t + σ 1 ( t ) y 1 ( t ) d ω 1 ( t ) , d y 2 ( t ) = y 2 ( t ) ( a 21 ( t ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 ( t ) A 2 ( t ) y 2 ( t ) + c 2 ( t ) A 1 ( t ) y 1 ( t ) m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 ( t ) + m 2 ( t ) A 2 ( t ) y 2 ( t ) ) d t + σ 2 ( t ) y 2 ( t ) d ω 2 ( t ) ,
where
A i ( t ) = j = 1 q ( 1 + ξ i j ) t T 0 t k < t ( 1 + ξ i k ) , i = 1 , 2 ,
with the convention that the product of an empty index set is 1. It is easy to show that both A 1 ( · ) and A 2 ( · ) are T-periodic. We refer the readers to [22] for details.
Lemma 5.
Let x 1 ( t ) = A 1 ( t ) y 1 ( t ) and x 2 ( t ) = A 2 ( t ) y 2 ( t ) . Then the following two statements hold.
(i) 
If y ( t ) = ( y 1 ( t ) , y 2 ( t ) ) is a solution of (5), then x ( t ) = ( x 1 ( t ) , x 2 ( t ) ) is a solution of (2);
(ii) 
If x ( t ) = ( x 1 ( t ) , x 2 ( t ) ) is a solution of (2), then y ( t ) = ( y 1 ( t ) , y 2 ( t ) ) is a solution of (5).
The assertions in Lemma 5, similar to Theorem 3.1 in [13], follow easily from Definition 2 and hence the proof is omitted here. Lemma 5 tells us that when the dynamics of (2) is considered, it suffices to focus on that of (5).
With similar arguments as those in the proof of Theorem 2.1 in [26], we can obtain the existence and non-negativeness of solutions to (5).
Theorem 1.
For any given initial value ( y 1 ( 0 ) , y 2 ( 0 ) ) R + 2 , (5) has a unique solution ( y 1 ( t ) , y 2 ( t ) ) on R + with ( y 1 ( t ) , y 2 ( t ) ) R + 2 for t R + with probability one.
Finally, for any positive, bounded, and continuous function f ( · ) on R + , we denote
f u = sup { f ( t ) , t 0 } , f l = inf { f ( t ) , t 0 } , f * = lim sup t f ( t ) , f * = lim inf t f ( t ) .
If f ( · ) is further integrable, then
f = lim t 1 t 0 t f ( s ) d s a n d f T = 1 T 0 T f ( s ) d s .

3. Existence of Periodic Markovian Processes

In this section, we first study the existence of periodic Markovian processes of (5).
Theorem 2.
Suppose the following condition ( H 1 ) holds,
( H 1 ) λ T : = m 2 l a 11 ( t ) c 1 ( t ) a 21 ( t ) m 2 l σ 1 2 ( t ) + σ 2 2 ( t ) 2 T + m 2 l T j = 1 q ln ( 1 + ξ 1 k ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) > max ( A 2 u a 22 u c 1 l A 2 l a 21 l ) 2 4 c 1 l ( A 2 l ) 2 a 22 l , ( A 1 u ( c 2 u a 11 u + m 2 l a 12 u ) ) 2 4 c 2 u a 12 l A 1 u A 1 l .
Then there exists at least one positive T-Markovian process for (5).
Proof. 
It suffices to find a C 2 -function V ( t , y 1 , y 2 ) and a closed set U R + 2 such that all conditions of Lemma 1 hold. Define
V 1 ( t , y 1 ) = c 2 u A 1 u y 1 m 2 l ln y 1 , V 2 ( y 2 ) = c 1 l A 2 l y 2 ln y 2 ,
and V ( t , y 1 , y 2 ) = V 1 ( t , y 1 ) + V 2 ( t , y 2 ) + W ( t ) , where W ( t ) is a positive differentiable function satisfying
W ( t ) = m 2 l a 11 ( t ) c 1 ( t ) a 21 ( t ) m 2 l σ 1 2 ( t ) + σ 2 2 ( t ) 2 m 2 l a 11 ( t ) c 1 ( t ) a 21 ( t ) m 2 l σ 1 2 ( t ) + σ 2 2 ( t ) 2 T .
It is not difficult to verify that W ( t ) is a T-periodic function on R + and
lim inf ( y 1 , y 2 ) R + 2 / U k V ( t , y 1 , y 2 ) a s k ,
where U k = { ( y 1 , y 2 ) : ( y 1 , y 2 ) ( 1 k , k ) × ( 1 k , k ) } . Hence V ( t , y 1 , y 2 ) is T-periodic and satisfies the first condition of Lemma 1. Next we verify the second condition of Lemma 1. Applying Itô’s formula to V 1 ( t , y 1 ) and V 2 ( t , y 2 ) gives
L V 1 ( t , y 1 ) = c 2 u A 1 u m 2 l y 1 y 1 ( a 11 ( t ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 ( t ) A 1 ( t ) y 1 c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 ) + m 2 l σ 1 2 ( t ) 2 m 2 l a 11 ( t ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) σ 1 2 ( t ) 2 + m 2 l c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 c 2 u A 1 u c 1 l A 2 l y 1 y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 + ( m 2 l a 12 ( t ) A 1 ( t ) + c 2 u A 1 u a 11 ( t ) ) y 1 c 2 u A 1 u a 12 ( t ) A 1 ( t ) y 1 2 ,
and
L V 2 ( t , y 2 ) = c 1 l A 2 l 1 y 2 y 2 ( a 21 ( t ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 ( t ) A 2 ( t ) y 2 + c 2 ( t ) A 1 ( t ) y 1 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 ) + σ 2 2 ( t ) 2 a 21 ( t ) 1 T j = 1 q ln ( 1 + ξ 2 k ) + σ 2 2 ( t ) 2 + c 1 l A 2 l c 2 u A 1 u y 1 y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 c 2 ( t ) A 1 ( t ) y 1 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 + ( a 22 ( t ) A 2 ( t ) c 1 l A 2 l a 21 ( t ) ) y 2 c 1 l A 2 l a 22 ( t ) A 2 ( t ) y 2 2 ,
respectively. Then
L V m 2 l a 11 ( t ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) σ 1 2 ( t ) 2 + a 21 ( t ) + σ 2 2 ( t ) 2 1 T j = 1 q ln ( 1 + ξ 2 k ) + m 2 l c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 + W ( t ) + ( c 2 u A 1 u a 11 ( t ) + m 2 l a 12 ( t ) A 1 ( t ) ) y 1 c 2 u A 1 u a 12 ( t ) A 1 ( t ) y 1 2 + ( a 22 ( t ) A 2 ( t ) c 1 l A 2 l a 21 ( t ) ) y 2 c 1 l A 2 l a 22 ( t ) A 2 ( t ) y 2 2 λ T c 2 u A 1 u a 12 l A 1 l y 1 2 + A 1 u ( c 2 u a 11 u + m 2 l a 12 u ) y 1 c 1 l A 2 l a 22 l A 2 l y 2 2 + ( A 2 u a 22 u c 1 l A 2 l a 21 l ) y 2 .
Choose ε ( 0 , 1 ) such that
1 2 λ T ( A 2 u a 22 u c 1 l A 2 l a 21 l ) 2 4 c 1 l ( A 2 l ) 2 a 22 l > A 1 u ( c 2 u a 11 u + m 2 l a 12 u ) ε , 1 2 λ T ( A 1 u ( c 2 u a 11 u + m 2 l a 12 u ) ) 2 4 c 2 u a 12 l A 1 u A 1 l > ( A 2 u a 22 u c 1 l A 2 l a 21 l ) ε ,
c 2 u A 1 u a 12 l A 1 l y 1 2 + A 1 u ( c 2 u a 11 u + m 2 l a 12 u ) y 1 < 1 2 for   y 1 > 1 ε , c 1 l A 2 l a 22 l A 2 l y 2 2 + ( A 2 u a 22 u c 1 l A 2 l a 21 l ) y 2 < 1 2 for   y 2 > 1 ε .
Denote
D ε = ( y 1 , y 2 ) : ( y 1 , y 2 ) [ ε , 1 ε ] × [ ε , 1 ε ] .
Then D ε is compact and its complement D ε C = R + 2 / D ε = i = 1 4 D ε i , where
D ε 1 = { ( y 1 , y 2 ) R + 2 | y 1 < ε } , D ε 2 = { ( y 1 , y 2 ) R + 2 | y 2 < ε } , D ε 3 = { ( y 1 , y 2 ) R + 2 | y 1 > 1 ε } , D ε 4 = { ( y 1 , y 2 ) R + 2 | y 2 > 1 ε } .
Accordingly, we distinguish four cases to finish the proof.
Case (i): ( y 1 , y 2 ) D ε 1 . Note that
c 1 l ( A 2 l ) 2 a 22 l y 2 2 + ( A 2 u a 22 u c 1 l A 2 l a 21 l ) y 2 ( A 2 u a 22 u c 1 l A 2 l a 21 l ) 2 4 c 1 l ( A 2 l ) 2 a 22 l .
By (6), we get
L V 1 2 λ T ( A 2 u a 22 u c 1 l A 2 l a 21 l ) 2 4 c 1 l ( A 2 l ) 2 a 22 l < 0 .
Case (ii): ( y 1 , y 2 ) D ε 2 . Similar to Case (i), we can get
L V 1 2 λ T ( A 1 u ( c 2 u a 11 u + m 2 l a 12 u ) ) 2 4 c 2 u a 12 l A 1 u A 1 l < 0 .
Case (iii): ( y 1 , y 2 ) D ε 3 . In this case, by (8) and the definition of λ T and the choice of ε , we have
L V λ T 1 2 + ( A 2 u a 22 u c 1 l A 2 l a 21 l ) 2 4 c 1 l ( A 2 l ) 2 a 22 l 1 2 .
Case (iv): ( y 1 , y 2 ) D ε 4 . Similar to case (iii), we can show that
L V 1 2 .
In summary, we have verified the second condition of Lemma 1. This completes the proof. □
Before proving the global attractivity of (5), we first study its boundedness and uniform continuity.
Lemma 6.
For any initial value y ( 0 ) = ( y 1 ( 0 ) , y 2 ( 0 ) ) T R + 2 , the solution ( y 1 ( t ) , y 2 ( t ) ) of (5) satisfies
lim sup t E ( y i p ) < H i ( p )
for any p > 0 , i = 1 , 2. Further
lim sup t P { | y i ( t ) | > η i } < ε , i = 1 , 2 ,
that is, the solution of (5) is stochastically ultimately bounded, where H i ( p ) and η i are positive constants and ε is a corresponding small positive constant.
Proof. 
Using Itô’s formula to y 1 p and y 2 p , we have
d y 1 p = [ p y 1 p a 11 ( t ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 ( t ) A 1 ( t ) y 1 c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) ( y 1 + m 2 ( t ) A 2 ( t ) y 2 + σ 1 2 ( t ) 2 p ( p 1 ) y 1 p ] d t + p y 1 p 1 σ 1 ( t ) y 1 d ω 1 ( t ) p y 1 p a 11 u + 1 T j = 1 q ln ( 1 + ξ 1 k ) + p ( σ 1 u ) 2 2 a 12 l A 1 l y 1 d t + p σ 1 ( t ) y 1 p d ω 1 ( t ) ,
and
d y 2 p = [ p y 2 p a 21 ( t ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 ( t ) A 2 ( t ) y 2 + c 2 ( t ) A 1 ( t ) y 1 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 + p ( p 1 ) σ 2 2 ( t ) 2 y 2 p ] d t + p y 2 p σ 2 ( t ) d ω 2 ( t ) p y 2 p a 21 l + 1 T j = 1 q ln ( 1 + ξ 2 k ) + p ( σ 2 u ) 2 2 + c 2 u m 1 l a 22 l A 2 l y 2 d t + p σ 2 ( t ) y 2 p d ω 1 ( t ) ,
respectively. By the properties of expectation, it follows that
d E ( y 1 p ) d t p a 11 u + 1 T j = 1 q ln ( 1 + ξ 1 k ) + p ( σ 1 u ) 2 2 E ( y 1 p ) p a 12 l A 1 l [ E ( y 1 p ) ] 1 + 1 p ,
and
d E ( y 2 p ) d t p a 21 l + 1 T j = 1 q ln ( 1 + ξ 2 k ) + p ( σ 2 u ) 2 2 + c 2 u m 1 l E ( y 2 p ) p a 22 l A 2 l [ E ( y 2 p ) ] 1 + 1 p .
Using the comparison theorem of differential equations, we obtain
lim sup t E ( y 1 p ( t ) ) a 11 u + 1 T j = 1 q ln ( 1 + ξ 1 k ) + p ( σ 1 u ) 2 2 a 12 l A 1 l p : = H ˘ 1 ( p ) ,
and
lim sup t E ( y 2 p ( t ) ) a 21 l + 1 T j = 1 q ln ( 1 + ξ 2 k ) + p ( σ 2 u ) 2 2 + c 2 u m 1 l a 22 l A 2 l p : = H ˘ 2 ( p ) .
Therefore, there exists a t ˜ > 0 such that E ( y 1 p ( t ) ) 3 2 H ˘ 1 ( p ) and E ( y 2 p ( t ) ) 3 2 H ˘ 2 ( p ) for any t t ˜ . Moreover, by the continuity of E ( y 1 p ( t ) ) and E ( y 2 p ( t ) ) , there exist H ˜ 1 ( p ) > 0 and H ˜ 2 ( p ) > 0 such that E ( y 1 p ( t ) ) H ˜ 1 ( p ) and E ( y 2 p ( t ) ) H ˜ 2 ( p ) for any t [ 0 , t ˜ ] . Then H 1 ( p ) = max { 3 2 H ˘ 1 ( p ) , H ˜ 1 ( p ) } and H 2 ( p ) = max { 3 2 H ˘ 2 ( p ) , H ˜ 2 ( p ) } are what we are looking for.
Next, applying Chebyshev inequality, we can directly derive that
lim sup t P { | y i ( t ) | > η i } < ε , i = 1 , 2 ,
which means that (5) is stochastically ultimately bounded. This completes the proof. □
Lemma 7.
Let ( y 1 ( t ) , y 2 ( t ) ) be a solution of (5) on R + with the initial value ( y 1 ( 0 ) , y 2 ( 0 ) ) R + 2 . Then almost every sample path of ( y 1 ( t ) , y 2 ( t ) ) is uniformly continuous on R + .
Proof. 
By Lemma 6 and the properties of expectation, we get
E y 1 a 11 ( t ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 ( t ) A 1 ( t ) y 1 c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 p 1 2 E a 11 ( t ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 ( t ) A 1 ( t ) y 1 c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 2 p + 1 2 E | y 1 | 2 p 1 2 H 1 ( 2 p ) + 3 2 p 1 ( a 11 u ) 2 p + a 12 u A 1 u H 1 ( 2 p ) + c 1 u m 2 l 2 p : = G 1 ( p ) .
By using stochastic integral inequality, for 0 < t 1 < t 2 and p > 2 , we have
E t 1 t 2 σ 1 ( s ) y 1 ( s ) d ω 1 ( s ) p ( σ 1 u ) p p ( p 1 ) 2 p 2 ( t 2 t 1 ) p 2 2 t 1 t 2 E ( | y 1 | p ) d s ( σ 1 u ) p p ( p 1 ) 2 p 2 ( t 2 t 1 ) p 2 H 1 ( p ) .
Thus when 0 < t 1 < t 2 < , t 2 t 1 1 , and 1 p + 1 q = 1 ,
E | y 1 ( t 2 ) y 1 ( t 1 ) | p = E t 1 t 2 y 1 a 11 ( s ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 ( s ) A 1 ( t ) y 1 c 1 ( s ) A 2 ( t ) y 2 m 0 ( s ) + m 1 ( s ) A 1 ( t ) y 1 + m 2 ( s ) A 2 ( t ) y 2 d s + t 1 t 2 σ 1 ( s ) y 1 ( s ) d ω 1 ( s ) p 2 p 1 E t 1 t 2 y 1 a 11 ( s ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 12 ( s ) A 1 ( t ) y 1 c 1 ( s ) A 2 ( t ) y 2 m 0 ( s ) + m 1 ( s ) A 1 ( t ) y 1 + m 2 ( s ) A 2 ( t ) y 2 d s p + 2 p 1 E t 1 t 2 σ 1 ( s ) y 1 ( s ) d ω 1 ( s ) p 2 p 1 ( t 2 t 1 ) p q t 1 t 2 E | y 1 ( a 11 ( s ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 12 ( s ) A 1 ( t ) y 1 c 1 ( s ) A 2 ( t ) y 2 m 0 ( s ) + m 1 ( s ) A 1 ( t ) y 1 + m 2 ( s ) A 2 ( t ) y 2 ) | p d s + 2 p 1 ( σ 1 u ) p p ( p 1 ) 2 p 2 ( t 2 t 1 ) p 2 H 1 ( p ) 2 p 1 ( t 2 t 1 ) p q + 1 G 1 ( P ) + 2 p 1 ( σ 1 u ) p p ( p 1 ) 2 p 2 ( t 2 t 1 ) p 2 H 1 ( p ) 2 p 1 ( t 2 t 1 ) p 2 1 + p ( p 1 ) 2 p 2 G 2 ( p ) ,
where G 2 ( p ) = max { G 1 ( p ) , σ 1 u H 1 ( p ) } . Then an application of Lemma 3 yields that almost every sample path of y 1 is locally but uniformly Hölder-continuous with exponent ϑ ( 0 , p 2 2 p ) and hence almost every sample path of y 1 is uniformly continuous on R + . The conclusion on y 2 is proved with similarity. This completes the proof. □
Now, we are ready to state and prove the global attractivity of (5).
Theorem 3.
Suppose the following condition holds,
( H 2 ) ι 1 = a 12 l c 1 u m 1 u m 0 l m 2 l c 2 u m 0 l > 0 and ι 2 = a 22 l c 1 u m 0 l c 2 u m 2 u m 0 l m 1 l > 0 .
Then all solutions of system (5) are globally attractive.
Proof. 
Let ( y 1 , y 2 ) and ( y ˜ 1 , y ˜ 2 ) be any two solutions of (5). Define V ( t ) = | ln y 1 ln y ˜ 1 | + | ln y 2 ln y ˜ 2 | . Then V ( · ) is continuous and positive on R + . Using Itô’s formula, we get
d + V = sgn ( y 1 y ˜ 1 ) { a 12 ( t ) A 1 ( t ) ( y 1 y ˜ 1 ) [ c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 c 1 ( t ) A 2 ( t ) y ˜ 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y ˜ 1 + m 2 ( t ) A 2 ( t ) y ˜ 2 ] } d t + sgn ( y 2 y ˜ 2 ) { a 22 ( t ) A 2 ( t ) ( y 2 y ˜ 2 ) + [ c 2 ( t ) A 1 ( t ) y 1 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 c 2 ( t ) A 1 ( t ) y ˜ 1 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y ˜ 1 + m 2 ( t ) A 2 ( t ) y ˜ 2 ] } d t a 12 ( t ) A 1 ( t ) | y 1 y ˜ 1 | d t a 22 ( t ) A 2 ( t ) | y 2 y ˜ 2 | d t + | c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 c 1 ( t ) A 2 ( t ) y ˜ 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y ˜ 1 + m 2 ( t ) A 2 ( t ) y ˜ 2 | d t + | c 2 ( t ) A 1 ( t ) y 1 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 c 2 ( t ) A 1 ( t ) y ˜ 1 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y ˜ 1 + m 2 ( t ) A 2 ( t ) y ˜ 2 | d t = a 12 ( t ) A 1 ( t ) | y 1 y ˜ 1 | d t a 22 ( t ) A 2 ( t ) | y 2 y ˜ 2 | d t + c 1 ( t ) A 2 ( t ) | ( m 0 ( t ) + m 1 ( t ) A 1 ( t ) y ˜ 1 ) ( y 2 y ˜ 2 ) + m 1 ( t ) A 1 ( t ) y ˜ 2 ( y 1 y ˜ 1 ) [ m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 ] [ m 0 ( t ) + m 1 ( t ) A 1 ( t ) y ˜ 1 + m 2 ( t ) A 2 ( t ) y ˜ 2 ] | d t + c 2 ( t ) A 1 ( t ) | m 2 ( t ) A 2 ( t ) y 1 ( y 2 y ˜ 2 ) + ( m 0 ( t ) + m 2 ( t ) A 2 ( t ) y 2 ) ( y 1 y ˜ 1 ) [ m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 ] [ m 0 ( t ) + m 1 ( t ) A 1 ( t ) y ˜ 1 + m 2 ( t ) A 2 ( t ) y ˜ 2 ] | d t a 12 ( t ) A 1 ( t ) | y 1 y ˜ 1 | d t + c 1 u m 1 u A 1 ( t ) m 0 l m 2 l | y 1 y ˜ 1 | d t + c 1 u A 2 ( t ) m 0 l | y 2 y ˜ 2 | d t a 22 ( t ) A 2 ( t ) | y 2 y ˜ 2 | d t + c 2 u A 1 ( t ) m 0 l | y 1 y ˜ 1 | d t + c 2 u m 2 u A 2 ( t ) m 0 l m 1 l | y 2 y ˜ 2 | d t A 1 l ι 1 | y 1 y ˜ 1 | d t A 2 l ι 2 | y 2 y ˜ 2 | d t .
Integrating both sides of the last inequality from 0 to t and using ( H 2 ), we get
V ( t ) V ( 0 ) 0 t α 1 ι 1 | y 1 y ˜ 1 | + α 2 ι 2 | y 2 y ˜ 2 | d s V ( 0 ) < .
Thus | y 1 y ˜ 1 | L 1 [ 0 , ) and | y 2 y ˜ 2 | L 1 [ 0 , ) . Now the result follows from Definition 3 and Lemma 4. This completes the proof. □
Remark 1.
Theorem 3 implies that if the pulses are bounded then they have no influence on the global attractivity of system (5), which agrees with ([27], Theorem 4) (stochastic but non-impulsive case).

4. Extinction and Permanence in Mean

For convenience, we denote
λ 1 : = a 11 ( t ) σ 1 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 1 k ) , λ 2 : = a 21 ( t ) σ 2 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 2 k ) , λ ˜ 1 : = a 11 ( t ) σ 1 2 ( t ) 2 c 1 ( t ) m 2 l T + 1 T j = 1 q ln ( 1 + ξ 1 k ) , λ ˜ 2 : = a 21 ( t ) σ 2 2 ( t ) 2 + c 2 ( t ) m 1 l T + 1 T j = 1 q ln ( 1 + ξ 2 k ) .
Obviously, they satisfy λ ˜ 1 λ 1 and λ 2 λ ˜ 2 by a direct verification.
Theorem 4.
For system (5), the following results hold.
(i) 
If λ 1 < 0 and λ ˜ 2 < 0 , then all species of (5) are extinct, i.e., lim t y i ( t ) = 0 , i = 1 , 2 .
(ii) 
If λ ˜ 1 > 0 and λ ˜ 2 < 0 , then y 1 is permanent in mean and y 2 is extinct, i.e., ς ˜ 1 y 1 ζ 1 for some positive numbers ς ˜ 1 and ζ 1 , lim t y 2 ( t ) = 0 .
(iii) 
If λ 1 < 0 and λ 2 > 0 , then y 1 is extinct and y 2 is permanent in mean, i.e., lim t y 1 ( t ) = 0 and ς 2 y 2 ζ ˜ 2 for some positive numbers ς 2 and ζ 2 ˜ .
(iv) 
If λ ˜ 1 > 0 and λ 2 > 0 , then y i is permanent in mean, i.e., ς i y i ζ i for some positive numbers ς i and ζ i , i = 1 , 2 .
Proof. 
Applying Itô’s formula to ln y 1 and ln y 2 gives
d ln y 1 = ( a 11 ( t ) + 1 T j = 1 q ln ( 1 + ξ 1 k ) σ 1 2 ( t ) 2 a 12 ( t ) A 1 ( t ) y 1 c 1 ( t ) A 2 ( t ) y 2 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 ) + σ 1 ( t ) d ω 1 ( t ) ,
and
d ln y 2 = ( a 21 ( t ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) σ 2 2 ( t ) 2 a 22 ( t ) A 2 ( t ) y 2 + c 2 ( t ) A 1 ( t ) y 1 m 0 ( t ) + m 1 ( t ) A 1 ( t ) y 1 + m 2 ( t ) A 2 ( t ) y 2 ) + σ 2 ( t ) d ω 2 ( t ) ,
respectively. Integrating (9) from 0 to t produces
ln y 1 ln y 1 ( 0 ) t = 1 t 0 t a 11 ( s ) σ 1 2 ( s ) 2 d t + 1 T j = 1 q ln ( 1 + ξ 1 k ) 1 t 0 t a 12 ( s ) A 1 ( s ) y 1 ( s ) d s 1 t 0 t c 1 ( s ) A 2 ( s ) y 2 m 0 ( s ) + m 1 ( s ) A 1 ( s ) y 1 + m 2 ( s ) A 2 ( s ) y 2 d s + 1 t 0 t σ 1 ( s ) d ω 1 ( s ) .
By virtue of the strong law of large numbers for martingales, lim sup t 1 t 0 t σ 1 ( s ) d ω 1 ( s ) = 0 . Then let t in the previous equality to obtain
ln y 1 t a 11 ( s ) σ 1 2 ( s ) 2 T + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 l A 1 l y 1 .
(i)
If λ 1 < 0 , applying Lemma 2, we obtain from (11) that lim t y 1 ( t ) = 0 . Then (10) reads
d ln y 2 = a 21 ( t ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) σ 2 2 ( t ) 2 a 22 ( t ) A 2 ( t ) y 2 d t + σ 2 ( t ) d ω 2 ( t ) .
Integrating both sides from 0 to t yields
ln y 2 ln y 2 ( 0 ) t 1 t 0 t a 21 ( s ) + σ 2 2 ( s ) 2 d s + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 l A 2 l y 2 + 1 t 0 t σ 2 ( t ) d ω 2 ( t ) .
With a similar argument as above, we can obtain lim t y 2 ( t ) = 0 . Therefore, all species are extinct.
(ii)
If λ ˜ 1 > 0 , then λ 1 > 0 . Lemma 2 and (11) imply
y 1 * a 11 ( t ) σ 1 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 l A 1 l : = ζ 1 .
By monotonicity, we can derive from (10) that
d ln y 2 a 21 ( t ) + 1 T j = 1 q ln ( 1 + ξ 2 k ) σ 2 2 ( t ) 2 a 22 ( t ) A 2 ( t ) y 2 + c 2 ( t ) A 1 ( t ) y 1 m 1 ( t ) A 1 ( t ) y 1 d t + σ 2 ( t ) d ω 2 ( t ) .
Thus
ln y 2 ln y 2 ( 0 ) t 1 t 0 t c 2 ( s ) m 1 l a 21 ( s ) σ 2 2 ( s ) 2 d s + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 l A 2 l y 2 + 1 t 0 t σ 2 ( t ) d ω 2 ( t ) .
Since λ ˜ 2 < 0 , it follows from Lemma 2 that lim t y 2 ( t ) = 0 . Using (9) again, we have
ln y 1 ln y 1 ( 0 ) t 1 t 0 t a 11 ( s ) σ 1 2 ( s ) 2 d s + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 u A 1 u y 1 + 1 t 0 t σ 1 ( s ) d s .
Letting t and using Lemma 2, we arrive at
y 1 * a 11 ( t ) σ 1 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 u A 1 u : = ς ˜ 1 .
In summary,
ς ˜ 1 = a 11 ( t ) σ 1 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 u A 1 u y 1 a 11 ( t ) σ 1 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 l A 1 l = ζ 1 ,
which means y 1 is permanent in mean.
(iii)
If λ 1 < 0 , case (i) implies that lim t y 1 ( t ) = 0 . It follows from (10) that
d ln y 2 a 21 ( t ) σ 2 2 ( t ) 2 d t + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 l A 2 l y 2 d t + σ 2 ( t ) d ω 2 ( t )
and
d ln y 2 a 21 ( t ) σ 2 2 ( t ) 2 d t + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 u A 2 u y 2 d t + σ 2 ( t ) d ω 2 ( t ) .
Since λ 2 > 0 , using Lemma 2 gives
ς 2 = a 21 ( t ) σ 2 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 u A 2 u y 2 a 21 ( t ) σ 2 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 l A 2 l = ζ ˜ 2 ,
which means y 2 is permanent in mean.
(iv)
Obviously, λ ˜ 1 > 0 and λ 2 > 0 imply that λ 1 > 0 and λ ˜ 2 > 0 , respectively. Therefore, species y 2 can not be extinct. It follows from (9) that
ln y 1 ln y 1 ( 0 ) t 1 t 0 t a 11 ( s ) c 1 ( s ) m 2 l σ 1 2 ( s ) 2 d s + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 u A 1 u y 1 + 1 t 0 t σ 1 ( s ) d s .
Apply Lemma 2 to get
y 1 * a 11 ( t ) c 1 ( t ) m 2 l σ 1 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 1 k ) a 12 u A 1 u : = ς 1 .
This, combined with (ii), produces ς 1 y 1 ζ 1 , i.e., y 1 is permanent in mean.
Since λ ˜ 2 > 0 , applying Lemma 2 to (12) yields
y 2 * c 2 ( t ) m 1 l a 21 ( t ) σ 2 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 l A 2 l : = ζ 2 .
Moreover, we get from (10) that
ln y 2 ln y 2 ( 0 ) t = 1 t 0 t a 21 ( s ) σ 2 2 ( s ) 2 d s a 22 ( t ) A 2 ( t ) y 2 + 1 T j = 1 q ln ( 1 + ξ 2 k ) + 1 t 0 t c 2 ( s ) A 1 ( s ) y 1 m 0 ( s ) + m 1 ( s ) A 1 ( s ) y 1 + m 2 ( s ) A 2 ( s ) y 2 d s + 1 t 0 t σ 2 ( s ) d ω 2 ( s ) c 2 ( t ) A 1 l ε 0 m 0 u + m 1 u A 1 u η 1 + m 2 u A 2 u η 2 a 21 ( t ) σ 2 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 u A 2 u y 2 + 1 t 0 t σ 2 ( t ) d ω 2 ( t ) .
In view of λ 2 > 0 and lim sup t 1 t 0 t σ 2 ( s ) d ω 2 ( s ) = 0 , applying Lemma 2 again yields
y 2 * c 2 ( t ) A 1 l ε 0 m 0 u + m 1 u A 1 u η 1 + m 2 u A 2 u η 2 a 21 ( t ) σ 2 2 ( t ) 2 T + 1 T j = 1 q ln ( 1 + ξ 2 k ) a 22 u A 2 u ς 2 .
Thus ς 2 y 2 ζ 2 , which means that y 2 is permanent in mean. The proof is completed. □
Remark 2.
Theorem 4 provides sufficient conditions on the extinction of all species. Obviously, for the autonomous case, they are in accordance with ([26], Theorem 4.1). Moreover, we have also established sufficient conditions ensuring the permanence in mean of every species for the non-autonomous case. Therefore, Theorem 4 not only improves but also generalizes the corresponding results in [26].

5. Stationary Distribution

In this section, we consider the distribution of the autonomous and non-impulsive case of (2). The model is as follows:
d x 1 ( t ) = x 1 ( t ) a 11 a 12 x 1 ( t ) c 1 x 2 ( t ) m 0 + m 1 x 1 ( t ) + m 2 x 2 ( t ) d t + σ 1 x 1 ( t ) d ω 1 ( t ) , d x 2 ( t ) = x 2 ( t ) a 21 a 22 x 2 ( t ) + c 2 x 1 ( t ) m 0 + m 1 x 1 ( t ) + m 2 x 2 ( t ) d t + σ 2 x 2 ( t ) d ω 2 ( t ) .
First, we introduce an important lemma for the existence of a stationary distribution. For the autonomous stochastic differential equation
d x ( t ) = b ( x ) d t + k = 1 r σ k ( x ) d ω k ( t ) ,
where x ( t ) is a homogeneous Markovian process in n-dimensional Euclidean space R n . The diffusion matrix A ( x ) = ( a i j ( x ) ) n × n ,   a i j ( x ) = k = 1 r σ k i ( x ) σ k j ( x ) .
Lemma 8
([16]). Assume there exists a bounded domain U R n with regular boundary Γ satisfying the following two conditions:
(i) 
F is uniformly elliptical in the domain U and some neighborhood thereof, where F u = b ( x ) u x + t r ( A ( x ) ) u x x 2 .
(ii) 
There is a non-negative C 2 -function V ( x ) and a positive constant C such that L V ( x ) C for any x R n / U .
Then the Markovian process x ( t ) of system (14) has a unique stationary distribution μ ( · ) , and for any integrable function f ( · ) with respect the measure μ, we have
P lim t 1 t 0 t f ( x ( t ) ) d t = R n f ( x ) μ ( d x ) = 1 .
Now we consider the stationary distribution of (13).
Theorem 5.
Suppose the following condition holds:
( H 3 ) λ 0 : = m 2 a 11 σ 1 2 2 a 21 + c 1 + σ 2 2 2 > max ( a 22 c 1 a 21 ) 2 4 c 1 a 22 , ( m 2 a 12 + c 2 a 11 ) 2 4 c 2 a 12 ,
then there exists a stationary distribution for system (13), which is ergodic.
Proof. 
Define
V 1 ( t ) = c 2 y 1 ( t ) m 2 ln y 1 ( t ) and V 2 ( t ) = c 1 y 2 ( t ) ln y 2 ( t ) .
Applying Itô’s formula to V 1 ( · ) and V 2 ( · ) , then
L V 1 = c 2 m 2 y 1 ( t ) y 1 ( t ) a 11 a 12 y 1 ( t ) c 1 y 2 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) + m 2 σ 1 2 2 = m 2 a 11 σ 1 2 2 + c 2 a 11 y 1 ( t ) c 2 a 12 y 1 2 ( t ) c 1 c 2 y 1 ( t ) y 2 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) + a 12 m 2 y 1 ( t ) + m 2 c 1 y 2 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) ,
and
L V 2 = c 1 1 y 2 ( t ) y 2 ( t ) a 21 a 22 y 2 ( t ) + c 2 y 1 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) + σ 2 2 2 = a 21 + σ 2 2 2 c 1 a 21 y 2 ( t ) c 1 a 22 y 2 2 ( t ) + c 1 c 2 y 1 ( t ) y 2 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) + a 22 y 2 ( t ) c 2 y 1 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) .
Let V ( t ) = V 1 ( t ) + V 2 ( t ) , then
L V = m 2 a 11 σ 1 2 2 c 2 a 12 y 1 2 ( t ) c 1 c 2 y 1 ( t ) y 2 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) + c 2 a 11 y 1 ( t ) + a 12 m 2 y 1 ( t ) + m 2 c 1 y 2 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) + a 21 + σ 2 2 2 + a 22 y 2 ( t ) c 1 a 21 y 2 ( t ) c 1 a 22 y 2 2 ( t ) + c 1 c 2 y 1 ( t ) y 2 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) c 2 y 1 ( t ) m 0 + m 1 y 1 ( t ) + m 2 y 2 ( t ) m 2 a 11 σ 1 2 2 ( a 21 + c 1 + σ 2 2 2 ) c 2 a 12 y 1 2 ( t ) + ( a 12 m 2 + c 2 a 11 ) y 1 ( t ) c 1 a 22 y 2 2 ( t ) + ( a 22 c 1 a 21 ) y 2 ( t ) = λ 0 c 2 a 12 y 1 2 ( t ) + ( a 12 m 2 + c 2 a 11 ) y 1 ( t ) c 1 a 22 y 2 2 ( t ) + ( a 22 c 1 a 21 ) y 2 ( t ) .
It is not difficult to verify that
lim inf ( y 1 , y 2 ) R + 2 / U k V ( t , y 1 , y 2 ) a s k ,
where U k = { ( y 1 , y 2 ) : ( y 1 , y 2 ) ( 1 k , k ) × ( 1 k , k ) } . Hence the first condition of Lemma 8 is satisfied. Next we verify the second condition of Lemma 8. Define a closed and compact set D ε as before, where ε > 0 is small enough such that
1 2 λ 0 ( a 22 c 1 a 21 ) 2 4 c 1 a 22 > ( a 12 m 2 + c 2 a 11 ) ε ,
and
1 2 λ 0 ( a 12 m 2 + c 2 a 11 ) 2 4 c 2 a 12 > a 22 ε .
Then the component D ε C = R + 2 / D ε = i = 1 4 D ε i is the same as in Theorem 2. We consider the following cases.
Case (i): ( y 1 , y 2 ) D ε 1 . It is clear that
c 1 a 22 y 2 2 ( t ) + ( a 22 c 1 a 21 ) y 2 ( t ) ( a 22 c 1 a 21 ) 2 4 c 1 a 22 .
By (15), then
L V 1 2 λ 0 ( a 22 c 1 a 21 ) 2 4 c 1 a 22 < 0 .
Case (ii): ( y 1 , y 2 ) D ε 2 . Similarly, we have
c 2 a 12 y 1 2 ( t ) + ( a 12 m 2 + c 2 a 11 ) y 1 ( t ) ( a 12 m 2 + c 2 a 11 ) 2 4 c 2 a 12 ,
and hence,
L V 1 2 λ 0 ( a 12 m 2 + c 2 a 11 ) 2 4 c 2 a 12 < 0 .
Obviously, if ( y 1 , y 2 ) D ε i , i = 3 , 4 , then L V . Lemma 8 shows system (13) has a unique stationary distribution, which is ergodic. This completes the proof. □
Remark 3.
For a stochastic but non-impulsive case, the stationary distribution is also studied in [26], but their result (Theorem 3.1 in [26]) is based on the existence of equilibrium, which is not necessary in Theorem 5, and the proof is also different.

6. Examples and Simulations

In this section, by applying the Milstein simulation method [30], we give some examples and simulations to illustrate our theoretical results and reveal the effects of random disturbance and impulsive factors.
Example 1.
In order to numerically examine the existence of a T-periodic solution of system (2), we choose the following parameters:
a 11 ( t ) = 5 + 0.2 sin 2 t , a 12 ( t ) = 2 + 0.1 cos 2 t , c 1 ( t ) = 1 + 0.3 cos 2 t , m 0 ( t ) = 1 + 0.3 sin 2 t , m 1 ( t ) = 0.2 + 0.1 sin 2 t , m 2 ( t ) = 1 + 0.5 sin 2 t , a 21 ( t ) = 1 + 0.2 sin 2 t , a 22 ( t ) = 2 + 0.1 sin 2 t , c 2 ( t ) = 0.8 + 0.3 cos 2 t , σ 1 ( t ) = 0.02 + 0.02 sin 2 t , σ 2 ( t ) = 0.01 + 0.01 sin 2 t , ξ 1 = 0.03 + 0.01 sin 2 t , ξ 2 = 0.02 + 0.01 sin 2 t , t k = k π , k N ,
such that (H 1 ) holds. That is to say, the conditions of Theorem 2 are satisfied. By Theorem 2, we can see that system (2) has a unique π-periodic solution, illustrated in Figure 1a–c.
Furthermore, if there are no impulsive or stochastic effects, i.e., σ i ( t ) = 0 or ξ i k = 0 , other parameters remain unchanged, then by verification, (H 1 ) holds and system (2) still has a π -periodic solution. The deterministic case ( σ i ( t ) = ξ i k = 0 ), the impulsive case ( σ i ( t ) = 0 ) and the stochastic case ( ξ i k = 0 ) are shown in Figure 1d–f, respectively.
By comparing (c) with (d–f), we can observe that small impulse or stochastic noise has little influence on the existence of the π -periodic solution.
Similarly, we can verify that all parameters meet (H 2 ) and Theorem 3 shows that the solution of (2) is globally attractive. Figure 2 confirms this.
Example 2.
We choose a set of parameters as follows in order to investigate the permanence and extinction of system (2). Let
a 11 ( t ) = 2 + 0.2 sin 2 t , a 12 ( t ) = 0.8 + 0.1 cos 2 t , c 1 ( t ) = 0.5 + 0.3 cos 2 t , m 0 ( t ) = 0.5 + 0.3 sin 2 t , m 1 ( t ) = 0.2 + 0.1 sin 2 t , m 2 ( t ) = 0.8 + 0.5 sin 2 t , a 21 ( t ) = 0.6 + 0.2 sin 2 t , a 22 ( t ) = 0.4 + 0.1 sin 2 t , c 2 ( t ) = 1 + 0.3 cos 2 t , σ 1 ( t ) = 0.02 + 0.02 sin 2 t , σ 2 ( t ) = 0.01 + 0.01 sin 2 t , ξ 1 = 0.03 + 0.01 sin 2 t , ξ 2 = 0.02 + 0.01 sin 2 t , t k = 0.1 k π , k N .
It is not difficult to derive that λ ˜ 1 > 0 and λ 2 > 0 . Then by Theorem 4, all species of (2) are permanent in the mean. Figure 3a confirms this.
Next we illustrate the cases of extinction in Theorem 4. In order to reveal the impacts of stochastic or impulsive factors, we only change the values of stochastic or impulsive parameters and keep all other parameters the same.
(i)
Let σ 2 ( t ) = 2 + 0.01 sin 2 t , then λ ˜ 1 > 0 and λ ˜ 2 < 0 . Theorem 4 implies x 1 is permanent in the mean and x 2 is extinct, see Figure 3b. It shows that too much white noise results in the extinction of the predator.
(ii)
If ξ 1 = 2 + 0.01 sin 2 t , ξ 2 = 1.8 + 0.01 sin 2 t , then λ 1 < 0 and λ 2 > 0 . Theorem 4 implies x 1 is extinct and x 2 is permanent in the mean, as illustrated in Figure 3c, which shows that too large a pulse leads to the extinction of the prey.
(iii)
If σ 1 ( t ) = 2 + 0.02 sin 2 t , σ 2 ( t ) = 1 + 0.01 sin 2 t , then λ 1 < 0 and λ ˜ 2 < 0 . Theorem 4 shows that both prey and predator are extinct (Figure 3d). This indicates that the white noise has a huge influence on the system permanence, and too much noise will make all species extinct.
Example 3.
For the autonomous and non-impulsive case, let a 11 = 2 ,   a 12 = 0.8 ,   c 1 = 1 ,   m 0 = 0.3 , m 1 = 0.2 ,   m 2 = 1 ,   a 21 = 0.8 ,   a 22 = 0.7 ,   c 2 = 0.8 ,   σ 1 = 0.2 ,   σ 2 = 0.1 . An easy computation shows that (H 3 ) holds, so the distribution of (2) is stable (see Figure 4). Figure 4a,b shows the time series graphs of x 1 ( t ) and x 2 ( t ) , respectively. Figure 4c is the density of distribution of x 1 ( t ) and x 2 ( t ) .

7. Discussion and Conclusions

In this paper, a stochastic predator–prey system with impulsive effects and Beddington–DeAngelis functional responses is studied. The complicated dynamical behaviors are revealed in our main results. Sufficient conditions of the existence of periodic Markovian process are established in Theorem 2. The attractivity of solutions is shown in Theorem 3. The extinction and permanence in the mean of predator and prey species are represented in Theorem 4. The existence and ergodicity of the stationary distribution for the autonomous and non-impulsive cases of system (2) are given in Theorem 5. Finally, using Matlab, some simulations (Figure 1, Figure 2, Figure 3 and Figure 4) are given to verify the main results.
By giving some remarks (Remarks 1, 2 and 3), the main difference from some existing results are shown in detail. Particularly, we establish the sufficient conditions assuring the existence of a periodic Markovian process of (2), which has not been studied previously. In the process of our analysis, the Beddington–DeAngelis functional responses brought some difficulties and inequality techniques were applied to overcome them. When there are too many kinds of functional responses, what is the best way to deal with other functional responses? Further, time delays often appear in biological models; how should the effect of time delays be discussed? These are necessary and interesting questions for us to study in the future.

Funding

This research was supported by the National Natural Science Foundation of China (grant number 11861027).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author would like to thank the editor and the anonymous reviewers for their valuable comments, which helped to improve the presentation of the paper.

Conflicts of Interest

The author declares no conflict of interest. I have no financial and personal relationships with other people or organizations that could inappropriately influence this work. The funders had no role in the design of the study.

References

  1. Cantrell, R.S.; Cosner, C. On the dynamics of predator-prey models with the Bedding-DeAngelis functional response. J. Math. Anal. Appl. 2001, 257, 206–222. [Google Scholar] [CrossRef]
  2. Fan, M.; Kuang, Y. Dynamics of a non-autonomous predator-prey system with the Beddington-DeAngelis functional response. J. Math. Anal. Appl. 2004, 295, 15–39. [Google Scholar] [CrossRef] [Green Version]
  3. Wang, H. Global analysis of the predator-prey system with Beddington-DeAngelis functional response. J. Math. Anal. Appl. 2003, 281, 395–401. [Google Scholar] [CrossRef] [Green Version]
  4. Friedman, A. Stochastic Differential Equations and Their Applications; Academic Press: New York, NY, USA, 1976. [Google Scholar]
  5. Du, N.H.; Sam, V.H. Dynamics of a stochastic Lotka-Volterra model perturbed by white noise. J. Math. Anal. Appl. 2006, 324, 82–97. [Google Scholar] [CrossRef] [Green Version]
  6. Mao, X.; Marion, G.; Renshaw, E. Environmental Brownian noise suppresses explosions in populations dynamics. Stoch. Process. Appl. 2002, 97, 95–110. [Google Scholar] [CrossRef]
  7. Li, X.; Mao, X. Population dynamical behavior of non-autonomous Lotka-Volterra competitive system with random perturbation. Discret. Contin. Dyn. Syst. 2009, 24, 523–545. [Google Scholar] [CrossRef] [Green Version]
  8. Zuo, W.; Jiang, D.; Sun, X.; Hayat, T.; Alsaedib, A. Long-time behaviors of a stochastic cooperative Lotka-Volterra system with distributed delay. Physica A 2018, 506, 542–559. [Google Scholar] [CrossRef]
  9. Yagi, A.; Ton, T.V. Dynamic of a stochastic predator-prey population. Appl. Math. Comput. 2011, 218, 3100–3109. [Google Scholar] [CrossRef]
  10. Shao, Y.; Li, Y. Dynamical analysis of a stage structured predator-prey system with impulsive diffusion and generic functional response. Appl. Math. Comput. 2013, 220, 472–481. [Google Scholar] [CrossRef]
  11. Wu, R.H.; Zou, X.L.; Wang, K. Asymptotic behavior of a stochastic non-autonomous predator-prey model with impulsive perturbations. Commun. Nonlinear Sci. Numer. Simul. 2015, 20, 965–974. [Google Scholar] [CrossRef]
  12. Liu, M.; Wang, K. On a stochastic logistic equation with impulsive perturbations. Comput. Math. Appl. 2012, 63, 871–886. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, S.W.; Tan, D.J. Dynamics of a stochastic predator-prey system in a polluted environment with pulse toxicant input and impulsive perturbations. Appl. Math. Model. 2015, 39, 6319–6331. [Google Scholar] [CrossRef]
  14. Chen, L.J.; Chen, F.D. Dynamic behaviors of the periodic predator-prey system with distributed time delays and impulsive effect. Nonlinear Anal. Real World Appl. 2011, 12, 2467–2473. [Google Scholar] [CrossRef]
  15. Zhang, S.Q.; Meng, X.Z.; Feng, T.; Zhang, T.H. Dynamics analysis and numerical simulations of a stochastic non-autonomous predator-prey system with impulsive effects. Nonlinear Anal. Hybrid Syst. 2017, 26, 19–37. [Google Scholar] [CrossRef]
  16. Jiang, D.; Zuo, W.J.; Hayat, T.; Alsaedi, A. Stationary distribution and periodic solutions for stochastic Holling-Leslie predator-prey systems. Physica A 2016, 460, 16–28. [Google Scholar] [CrossRef]
  17. Zhao, Y.; Yuan, S.; Zhang, T.H. Stochastic periodic solution of a non-autonomous toxic-producing phytoplankton allelopathy model with environmental fluctuation. Commun. Nonlinear Sci. Numer. Simul. 2017, 44, 266–276. [Google Scholar] [CrossRef]
  18. Jiang, D.; Zhang, Q.M.; Hayat, T.; Alsaedi, A. Periodic solution for a stochastic non-autonomous competitive Lotka-Volterra model in a polluted environment. Physica A 2017, 471, 276–287. [Google Scholar] [CrossRef]
  19. Shao, Y. Globally asymptotical stability and periodicity for a nonautonomous two-species system with diffusion and impulses. Appl. Math. Model. 2012, 36, 288–300. [Google Scholar] [CrossRef]
  20. Li, D.; Xu, D. Periodic solutions of stochastic delay differential equations and applications to Logistic equation and neural networks. J. Korean Math. Soc. 2013, 50, 1165–1181. [Google Scholar] [CrossRef] [Green Version]
  21. Zhang, B.; Gopalsamy, K. On the periodic solution of N-dimensional stochastic population models. Stoch. Anal. Appl. 2000, 18, 323–331. [Google Scholar] [CrossRef]
  22. Zhang, Y.; Chen, S.H.; Gao, S.J.; Wei, X. Stochastic periodic solution for a perturbed non-autonomous predator-prey model with generalized nonlinear harvesting and impulses. Physica A 2017, 486, 347–366. [Google Scholar] [CrossRef]
  23. Zuo, W.J.; Jiang, D. Periodic solutions for a stochastic non-autonomous Holling-Tanner predator-prey system with impulses. Nonlin. Anal. Hybrid Syst. 2016, 22, 191–201. [Google Scholar] [CrossRef]
  24. Liu, M.; Wang, K. Asymptotic behavior of a stochastic nonautonomous Lotka-Volterra competitive system with impulsive perturbations. Math. Comput. Model. 2013, 57, 909–925. [Google Scholar] [CrossRef]
  25. Khas’minskii, R.Z. Stochastic Stability of Differential Equations; Sijthoff Noordhoff: Alphen aan den Rijn, The Netherlands; Springer Science & Business Media: Dordrecht, The Netherlands, 1980. [Google Scholar]
  26. Ji, C.Y.; Jiang, D. Dynamics of a stochastic density dependent predator-prey system with Beddington-DeAngelis functional response. J. Math. Anal. Appl. 2011, 381, 441–453. [Google Scholar] [CrossRef] [Green Version]
  27. Qiu, H.; Liu, M.; Wang, K.; Wang, Y. Dynamics of a stochastic predator-prey system with Beddington-DeAngelis functional response. Appl. Math. Comput. 2012, 219, 2303–2312. [Google Scholar] [CrossRef]
  28. Karatzas, I.; Shreve, S.E. Brownian Motion and Stochastic Calculus; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  29. Barbalat, I. Systems dequations differentielles d’osci d’oscillations nonlineaires. Rev. Roum. Math. Pures Appl. 1959, 4, 267–270. [Google Scholar]
  30. Higham, D.J. An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Rev. 2001, 43, 525–546. [Google Scholar] [CrossRef]
Figure 1. The periodic Markovian process for (2) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 0.02 . (a) The time series graph of x 1 ( t ) for Example 1; (b) the time series graph of x 2 ( t ) for Example 1; (c) the phase graph of Example 1 with σ i ( t ) 0 , ξ i k 0 ; (d) the phase graph of deterministic system ( ξ i k = σ i ( t ) = 0 ) ; (e) the phase graph of impulsive system ( σ i ( t ) = 0 , ξ i k 0 ) ; (f) the phase graph of stochastic system ( σ i ( t ) 0 , ξ i k = 0 ) i = 1 , 2 , k N .
Figure 1. The periodic Markovian process for (2) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 0.02 . (a) The time series graph of x 1 ( t ) for Example 1; (b) the time series graph of x 2 ( t ) for Example 1; (c) the phase graph of Example 1 with σ i ( t ) 0 , ξ i k 0 ; (d) the phase graph of deterministic system ( ξ i k = σ i ( t ) = 0 ) ; (e) the phase graph of impulsive system ( σ i ( t ) = 0 , ξ i k 0 ) ; (f) the phase graph of stochastic system ( σ i ( t ) 0 , ξ i k = 0 ) i = 1 , 2 , k N .
Axioms 10 00323 g001
Figure 2. The attractivity of (2). The red lines are the time series of x 1 ( t ) and x 2 ( t ) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 0.02 , and the blue lines are the time series of x 1 ( t ) and x 2 ( t ) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 2 . (a) The attractivity of x 1 ( t ) ; (b) The attractivity of x 2 ( t ) .
Figure 2. The attractivity of (2). The red lines are the time series of x 1 ( t ) and x 2 ( t ) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 0.02 , and the blue lines are the time series of x 1 ( t ) and x 2 ( t ) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 2 . (a) The attractivity of x 1 ( t ) ; (b) The attractivity of x 2 ( t ) .
Axioms 10 00323 g002
Figure 3. The permanence in the mean and extinction of (2) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 0.02 . (a) the time series graph of x 1 ( t ) and x 2 ( t ) ; (b) the permanence in the mean of x 1 ( t ) and extinction of x 2 ( t ) with σ 1 ( t ) = 0.02 + 0.02 sin 2 t , σ 2 ( t ) = 2 + 0.01 sin 2 t ; (c) the extinction of x 1 ( t ) and permanence in the mean of x 2 ( t ) with ξ 1 = 2 + 0.01 sin 2 t , ξ 2 = 1.8 + 0.01 sin 2 t ; (d) the extinction of all species with σ 1 ( t ) = 2 + 0.02 sin 2 t and σ 2 ( t ) = 1 + 0.01 sin 2 t .
Figure 3. The permanence in the mean and extinction of (2) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 0.02 . (a) the time series graph of x 1 ( t ) and x 2 ( t ) ; (b) the permanence in the mean of x 1 ( t ) and extinction of x 2 ( t ) with σ 1 ( t ) = 0.02 + 0.02 sin 2 t , σ 2 ( t ) = 2 + 0.01 sin 2 t ; (c) the extinction of x 1 ( t ) and permanence in the mean of x 2 ( t ) with ξ 1 = 2 + 0.01 sin 2 t , ξ 2 = 1.8 + 0.01 sin 2 t ; (d) the extinction of all species with σ 1 ( t ) = 2 + 0.02 sin 2 t and σ 2 ( t ) = 1 + 0.01 sin 2 t .
Axioms 10 00323 g003
Figure 4. The distribution of (13) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 0.02 . (a) the time series graph of x 1 ( t ) ; (b) the time series graph of x 2 ( t ) ; (c) the density of distribution of x 1 ( t ) and x 2 ( t ) .
Figure 4. The distribution of (13) with initial data x 1 ( 0 ) = x 2 ( 0 ) = 0.02 . (a) the time series graph of x 1 ( t ) ; (b) the time series graph of x 2 ( t ) ; (c) the density of distribution of x 1 ( t ) and x 2 ( t ) .
Axioms 10 00323 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shao, Y. Dynamics of an Impulsive Stochastic Predator–Prey System with the Beddington–DeAngelis Functional Response. Axioms 2021, 10, 323. https://doi.org/10.3390/axioms10040323

AMA Style

Shao Y. Dynamics of an Impulsive Stochastic Predator–Prey System with the Beddington–DeAngelis Functional Response. Axioms. 2021; 10(4):323. https://doi.org/10.3390/axioms10040323

Chicago/Turabian Style

Shao, Yuanfu. 2021. "Dynamics of an Impulsive Stochastic Predator–Prey System with the Beddington–DeAngelis Functional Response" Axioms 10, no. 4: 323. https://doi.org/10.3390/axioms10040323

APA Style

Shao, Y. (2021). Dynamics of an Impulsive Stochastic Predator–Prey System with the Beddington–DeAngelis Functional Response. Axioms, 10(4), 323. https://doi.org/10.3390/axioms10040323

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop