Next Article in Journal
Numerical Simulation for Brinkman System with Varied Permeability Tensor
Next Article in Special Issue
Stochastic Modeling of Within-Host Dynamics of Plasmodium Falciparum
Previous Article in Journal
Meaningful Secret Image Sharing with Uniform Image Quality
Previous Article in Special Issue
Stochastic Epidemic Model for COVID-19 Transmission under Intervention Strategies in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of a Class of Predation-Predation Model Dynamics with Random Perturbations

Department of Mathematics, Yunnan Minzu University, No. 2929, Yuehua Street, Chenggong District, Kunming 650500, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(18), 3238; https://doi.org/10.3390/math10183238
Submission received: 2 July 2022 / Revised: 24 August 2022 / Accepted: 31 August 2022 / Published: 6 September 2022
(This article belongs to the Special Issue Stochastic Models in Mathematical Biology)

Abstract

:
In this paper, we study a class of predation–prey biological models with random perturbation. Firstly, the existence and uniqueness of systematic solutions can be proven according to Lipschitz conditions, and then we prove that the systematic solution exists globally. Moreover, the article discusses the long-term dynamical behavior of the model, which studies the stationary distribution and gradual properties of the system. Next, we use two different methods to give the conditions of population extinction. From what has been discussed above, we can safely draw the conclusion that our results are reasonable by using numerical simulation.

1. Introduction

Along with the remarkable development of science and technology, mathematics is applied everywhere around us. In this way, it is possible for us to use mathematics to solve various problems in our daily lives. The application of mathematics in biology is more significant, in 2015, J. Banasiak [1] solved systems biology problems by establishing a mathematical framework. He studied multi-scale problems in the complex domains in biology, which can help us to use mathematical methods to solve biological problems. By studying population models, we can also effectively predict and control biological populations. Grunert K [2], Luoyi Wu [3] and Àngel Calsina [4] have conducted research on this.
Research on biological populations has acquired many good results, which can be referred to [5,6,7,8,9,10], the study of biological populations is of great significance to us. By analyzing the rules of the biological population change and studying the factors affecting the population change, we can take appropriate measures to control the population number. In addition, the study of population models can make reasonable predictions for population numbers, which can provide correct suggestions for industrial and animal protection. Due to the fact that the number of populations is of great importance to ecosystems, the study of population models is reasonable.
Predators initially grow when they enter the ecosystem, but they are unlikely to grow indefinitely due to limiting factors, such as food and habitat space. When the population reaches a certain degree, the population reaches saturation, so the introduction of the saturation factor y ( t ) 2 + D y ( t ) is crucial. The non-monotonic functional predator-bait model is mainly explored in this paper, the specific mathematical expression is given below:
d x d t = r x ( t ) ( 1 x ( t ) k ) x ( t ) y ( t ) a + x ( t ) + y ( t ) d y d t = u x ( t ) y ( t ) a + x ( t ) + y ( t ) y ( t ) 2 + D y ( t )
where x ( t ) and y ( t ) are the populations of prey and predators at the time t, the parameter r represents the intrinsic growth rate of the prey, k represents the environmental accommodation, x ( t ) y ( t ) a + x ( t ) + y ( t ) is the functional response function for both predators and prey, u represents the maximum growth rate generated by the predation, the parameter D represents the natural growth rate of the predators, and all the parameters are positive numbers.
The model has been discussed a lot by plenty of scholars before. However, in spite of that, the population is influenced not only by the above parameters but also by environmental factors such as climate, pollution and disease, as well as various human factors, the model is determinstic, so the model ( 1 ) may have certain limitations, it is very necessary to add a certain random perturbation to the original model, this paper mainly analyzes the influence of white noise [11,12,13] on the model, so the model ( 1 ) becomes:
d x ( t ) = ( r x ( t ) ( 1 x ( t ) k ) x ( t ) y ( t ) a + x ( t ) + y ( t ) ) d t + σ 1 x ( t ) d B 1 ( t ) d y ( t ) = ( u x ( t ) y ( t ) a + x ( t ) + y ( t ) y ( t ) 2 + D y ( t ) ) d t + σ 2 y ( t ) d B 2 ( t )
where B 1 ( t ) and B 2 ( t ) are independent wiener processes, σ 1 and σ 2 are the strength of white noise, and they are all non-negative and bounded.
This model can be regarded as a stochastic process, so this system can be studied through the relevant laws of geometric Brownian motion. According to the modeling method of reference [14], this model can also be seen as the application of two stochastic resetting scenarios [15] in the heterogeneous diffusion process, one is a stochastic resetting with a diffusion coefficient σ 1 2 x 2 2 , and the other is a stochastic resetting with a diffusion coefficient σ 2 2 y 2 2 , so this model can also be considered with the resetting idea, by establishing the relationship between geometric Brownian motion and stochastic resetting, it facilitates us to study stochastic processes in different minds.
With further research, there are many models that focus on the role of resetting, especially in the geometric Brownian movement, related articles can be referred to [16]. Resetting dynamics plays an important role in quantum physics and biophysics, and specific information can refer to the literature [17,18,19,20], resetting also has important applications in computer fields [21] and economics [22].
White noise plays a crucial role in biological population research, and the introduction of white noise makes the model more reasonable in the actual situation. Moreover, the study of white noise has profound significance in population control and pest elimination. Through the control of noise intensity, it can achieve the purpose of population control, which is of great significance in ecological balance. Moreover, the purpose of eliminating pests can be achieved by reasonable control of noise intensity, correlational research is discussed in this paper.
In Section 2, firstly, the article explains that the solution of the system 2 exists locally and is unique according to the Lipshiz conditions [23,24,25], secondly, we show using proof by contradiction to prove that the solution exists globally.
In Section 3, the thesis mainly explores the properties of the solution. Firstly, the paper studies the long-term dynamical behavior of the system and gives a sufficient condition for the existence of a stationary distribution [26,27]. Secondly, the gradual stabilization of the solution is proved in the article, it proves that the solution of the system oscillates at the internal equilibrium state [28] of the corresponding deterministic system. Finally, the paper gives numerical simulations of the stationary distribution.
In Section 4, the paper explores the extinction of the populations with two methods. The first way is to achieve pest control purposes by controlling the intensity and size of the random perturbations, and another is to control the pests by natural predation.
In Section 5, we have summarized the full text. First, the advantages and disadvantages of the article are analyzed, and I put forward the relevant requirements and expectations for myself, hoping that I could make progress in future scientific research and study.

2. Global Existence and Uniqueness of the Solution

Before discussing the various properties of the system, first judges whether the system solution exists, if so, whether the system solution exists is global. In this section, the article focuses on the global existence and uniqueness of the solution by constructing an auxiliary function.
Theorem 1.
For any initial condition ( x ( 0 ) , y ( 0 ) ) , system ( 2 ) must exists a unique solution ( x ( t ) , y ( t ) ) , and the solution exists globally.
Proof of Theorem 1.
We observe that parameters of the system do not satisfy the Lipschitz condition, first, we discuss the following system of equations:
d p ( t ) = [ r σ 1 2 2 r e p ( t ) k e q ( t ) a + e p ( t ) + e q ( t ) ] d t + σ 1 d B 1 t d q ( t ) = [ D σ 2 2 2 e q ( t ) + u e p ( t ) a + e p ( t ) + e q ( t ) ] d t + σ 2 d B 2 t
For any initial value p ( 0 ) = ln x ( 0 ) ) , q ( 0 ) = ln y ( 0 ) , system ( 3 ) has a unique solution ( p ( t ) , q ( t ) ) at t [ 0 , τ e ) , τ e is blasting time, according to the Ito formula, system 2 has a unique solution ( x ( t ) , y ( t ) ) = ( e p ( t ) , e q ( t ) ) . The next step is to prove that the solution exists globally, that is τ e a.s. Let m 0 be big enough, the initial condition ( x ( 0 ) , y ( 0 ) ) belongs to the interval [ 1 m 0 , m 0 ] , for any positive number m m 0 , define the stopping time [29]:
τ m = i n f { t [ 0 , τ e ) : x ( t ) [ 1 m 0 , m 0 ] o r y ( t ) [ 1 m 0 , m 0 ] }
It is easy to see that τ m is rising monotonically, stipulate τ = lim t + τ m , in the following proof, we prove that τ = + a.s.
Proof by contradiction: Suppose τ = + does not hold, there must be a finite positive value T and 0 < ϵ < 1 , the following results hold true
P ( τ m T ) ϵ
Construct a Lyapunov function [30,31,32] V 1 ( x , y ) = ( x + 1 ln x ) + c ( y + 1 ln y ) , and c is a positive parameter that is determined in the proof, apparently V 1 is a map from R 2 + to R + , according to the Ito formula:
d V 1 ( x , y ) = L V 1 ( x , y ) d t + ( 1 1 x ) σ 1 x d B 1 ( t ) + ( c c y ) σ 1 y d B 2 ( t )
where
L V 1 ( x , y ) = ( 1 1 x ) ( r x r x 2 k x y a + x + y ) + σ 1 2 x 2 2 x 2 + ( c c y ) ( u x y a + x + y + D y y 2 ) + c σ 2 2 y 2 2 y 2 = r x r k x 2 x y a + x + y r + r k x + y a + x + y + σ 1 2 2 + c u x y a + x + y + c D y c y 2 c u x a + x + y c D + c y + c σ 2 2 2 ( r + r k ) x + ( 1 a + c D + c ) y + c u 1 a + x + y x y + λ
where λ = σ 1 2 2 + c σ 2 2 2 , take c = 1 u , obviously, z < 2 ( z + 1 ln z ) is established by the relation of derivatives and extremvalues in the interval ( 0 , + ) , so above formula becomes
L V 1 ( x , y ) ( r + r k ) x + ( 1 a + c D + c ) y + λ β ( x + c y ) + λ 2 β [ x + 1 ln x + c ( y + 1 ln y ) ] + λ = 2 β V 1 ( x , y ) + λ
where β = m a x { r + r k , 1 a c + D + 1 } , substitute ( 7 ) to Formula ( 6 ) , next both sides of the upper formula simultaneously integrate in the interval [ 0 , τ m ] , and take the expectation.
E 0 τ m d V 1 ( x , y ) E 0 τ m [ λ + 2 β V ( x , y ) ] d t
that is
E V 1 ( x ( τ m ) , y ( τ m ) ) V 1 ( x 0 , y 0 ) + λ T + 2 β 0 τ m E V 1 ( x , y ) d t
according to Gronwall’s inequation
E V 1 ( x ( τ m ) , y ( τ m ) ) [ V 1 ( x 0 , y 0 ) + λ T ] e 0 τ m 2 β d t [ V 1 ( x 0 , y 0 ) + λ T ] e 2 β T = γ
When m > m 0 , take event Ω m , according to the hypothesis (5), then P ( Ω m ) ϵ > 0 , while ω Ω m , according to the expectation definition has the following results:
γ E [ I Ω m V 1 ( x ( τ m ) , y ( τ m ) ) ] ϵ [ m + 1 ln m + 1 u ( m + 1 ln m ) ]
then take the limit on both sides
γ
there is a contradiction, so γ = + a.s., so T = + , so solutions exist globally. □

3. The Nature of the Solution

3.1. Existence and Ergodicity of the Stationary Distribution

In this section, we mainly discuss the long-term kinetic behavior of predators and prey, we study the weak stability of the system, when the parameters of the system meet certain conditions, the populations are long-term lasting, and the population number will be stable around a determined value, the existence of the stationary distribution of the model ( 2 ) is discussed below.
Theorem 2.
Any initial value condition ( x 0 , y 0 ) , if r 1 σ 1 2 > 0 and D σ 2 2 > 0 hold true, a stationary distribution μ ( · ) exists in system (2), and the stationary distribution μ ( · ) is ergodic.
Proof of Theorem 2.
Define a Lyapunov function V 2 ( x , y ) = x + x δ + y + y δ , the δ is a constant to be determined later, obviously, V 3 ( x , y ) is a positive definite function, according to the Ito formula
d V 2 ( x , y ) = L V 2 ( x , y ) d t + ( 1 + δ x δ 1 ) σ 1 x d B 1 ( t ) + ( 1 + δ y δ 1 ) σ 2 y d B 2 ( t )
where
L V 2 ( x , y ) = ( x + δ x δ ) ( r r x k y a + x + y ) + 1 2 δ ( δ 1 ) σ 1 2 x δ + ( y + δ y δ ) ( u x a + x + y y + D ) + 1 2 δ ( δ 1 ) σ 2 2 y δ = r x r k x 2 x y a + x + y + r δ x δ r δ k x δ + 1 δ x δ y a + x + y + 1 2 δ ( δ 1 ) σ 1 2 x δ + u x y a + x + y y 2 + D y + δ u x y δ a + x + y δ y δ + 1 + D δ y δ + 1 2 δ ( δ 1 ) σ 2 2 y δ
So just make δ < 0 , let δ = 1 , the above formula sub changes to
L V 2 ( x , y ) r k x 2 + r x r x + r k + 1 x + σ 1 2 x y 2 + ( u + D ) y D y + 1 + σ 2 2 y = r k x 2 + r x r 1 σ 1 2 x y 2 + ( u + D ) y D σ 2 2 y + ( r k + 1 )
For any small amount ϵ 1 and ϵ 2 , consider the following set
U = { ( x , y ) R + 2 | ϵ 1 < x < 1 ϵ 1 , ϵ 2 < y < 1 ϵ 2 }
then for the set E = R + 2 \ U , L V equals to negative infinite constant holds, so there is a positive constant M, make the L V < M established at set E, so system ( 2 ) possesses a stationary distribution μ ( · ) . In addition, there is a constant N
N = m i n { α 1 2 x 2 , α 2 2 y 2 , ( x , y ) E } > 0
for any η , there is the following result
i , j = 1 2 a i j η i η j = α 1 2 x 2 η 1 2 + α 2 2 y 2 η 2 2 N | | η | | 2
so the stationary distribution μ ( · ) is ergodic. □
Brief summary: According to the proof of Theorem 2, when the system parameters meet the above condition, the system ( 2 ) has a stationary distribution, so both of the populations have always been positive, it indicates that all populations of the system ( 2 ) will persist.

3.2. The Progressive Nature of the System (2) at the Internal Equilibrium Point ( x 1 , y 1 ) in the System (1)

If ( x 1 , y 1 ) is the internal equilibrium state of the system ( 1 ) , we can clearly see that ( x 1 , y 1 ) is not the solution of a system ( 2 ) , so what is the relationship between ( x 1 , y 1 ) and the solution of a stochastic system ( 2 ) ? Next, the solution of the system ( 2 ) proceeds near the internal equilibrium state ( x 1 , y 1 ) of the deterministic system ( 1 ) is studied.
Theorem 3.
For any initial value condition ( x 0 , y 0 ) , if r k y 1 a ( a + x 1 + y 1 ) > 0 , then the solution of the stochastic system ( 2 ) has the following properties
lim t s u p 1 t · E 0 t [ ( x x 1 ) 2 + ( y y 1 ) 2 ] d s w h
both w and h are positive constants, they are given specifically later.
Proof of Theorem 3.
Define a function V 3 ( x , y ) = ( x x 1 x 1 ln x x 1 ) + λ 1 ( y y 1 y 1 ln y y 1 ) = V 31 ( x , y ) + λ 1 V 32 ( x , y ) , where λ 1 is a positive constant to be determined, obviously R 2 + R + .
Paramount consideration V 31 ( x , y ) , according to the Ito formula.
L V 31 ( x , y ) = ( x x 1 ) ( r r k x y a + x + y ) + x 1 σ 1 2 2 = ( x x 1 ) [ r k x 1 + y 1 a + x 1 + y 1 r k x y a + x + y ] + x 1 σ 1 2 2 = ( x x 1 ) [ r k ( x x 1 ) + a y 1 + y 1 x a y x 1 y ( a + x 1 + y 1 ) ( a + x + y ) ] + x 1 σ 1 2 2 = a ( y y 1 ) ( x x 1 ) + y 1 ( x x 1 ) 2 x 1 ( y y 1 ) ( x x 1 ) ( a + x 1 + y 1 ) ( a + x + y ) r k ( x x 1 ) 2 + x 1 σ 1 2 2
Then discuss V 32 ( x , y ) , according to the Ito formula
L V 32 ( x , y ) = ( y y 1 ) ( u x a + x + y y u x 1 a + x + y y 1 ) + y 1 σ 2 2 2 = ( y y 1 ) ( u x ( a + x 1 + y 1 ) u x 1 ( a + x + y ) a + x + y ( y y 1 ) ) + y 1 σ 2 2 2 = ( y y 1 ) ( a u ( x x 1 ) u ( x 1 y y 1 x ) ( a + x + y ) ( a + x 1 + y 1 ) ( y y 1 ) ) + y 1 σ 2 2 2 = ( y y 1 ) ( a u ( x x 1 ) u ( x 1 ( y y 1 ) y 1 ( x x 1 ) ) ( a + x + y ) ( a + x 1 + y 1 ) ( y y 1 ) ) + y 1 σ 2 2 2 = ( a u + u y 1 ) ( x x 1 ) ( y y 1 ) u x 1 ( y y 1 ) 2 ( a + x + y ) ( a + x 1 + y 1 ) ( y y 1 ) 2 + y 1 σ 2 2 2
Thus
L V 3 ( x , y ) ( r k y 1 a ( a + x 1 + y 1 ) ) ( x x 1 ) 2 λ 1 ( y y 1 ) 2 + a x 1 + λ 1 a u + λ 1 u y 1 ( a + x + y ) ( a + x 1 + y 1 ) ( x x 1 ) ( y y 1 ) + x 1 σ 1 2 2 + y 1 σ 2 2 2
so just take λ 1 = a + x 1 a u + u y 1 , h = m i n { r k y 1 a ( a + x 1 + y 1 ) , λ 1 } > 0 , bring them into the Formula ( 12 ) , both sides are integrated in the interval t ( 0 , t ) , and then take the expectation.
E [ V ( x , y ) ] V ( 0 , 0 ) h E 0 t [ ( x x 1 ) 2 + ( y y 1 ) 2 ] d s + w t
where w = x 1 σ 1 2 2 + y 1 σ 2 2 2 , so
lim t s u p 1 t · E 0 t [ ( x x 1 ) 2 + ( y y 1 ) 2 ] d s w h
so the Theorem 3 is proved. □
Brief summary: If the system ( 2 ) parameters satisfy the above condition, both predators and prey are persistent, and the solution of the system ( 2 ) will fluctuate up and down near the internal equilibrium state of its corresponding deterministic system. We can also get a conclusion, the system ( 2 ) is deterministic when the perturbation coefficient σ i are zero i = 1 , 2 , and the system is randomly stable.

3.3. Numerical Simulation

According to the milsteins discretization method [33,34,35], the following discretization equations are obtained
x i + 1 = x i + ( r x i ( 1 x i k ) x i y i a + x i + y 1 ) Δ t + σ 1 x i Δ t ξ 1 + 1 2 σ 1 2 x i ( ξ 1 2 1 ) Δ t y i + 1 = y i + ( u x i y i a + x i + y i y i 2 + D y ) Δ t + σ 2 x i Δ t ξ 2 + 1 2 σ 2 2 x i ( ξ 2 2 1 ) Δ t
Let initial value x 0 = 1 , y 0 = 12 , and take the parameter r = 10 , a = 10 , σ 1 = 1 , σ 2 = 1 , u = 0.3 , D = 5 , Δ = 200 2 16 , k = 1 , then 10 1 1 2 = 8 > 0 , 5 1 2 = 4 > 0 , thus the values satisfy the conditions of Theorem 2, the numerical simulation are as follows, see Figure 1.
By observing Figure 1, it shows that the system does have a stationary distribution, populations x ( t ) and y ( t ) are persistent, and they fluctuate up and down around the determined values, so the conclusion of the Theorem 2 is correct.

4. The Final Behavior of the Population

4.1. Environmental Forces to Eliminate Insect Pests

In nature, the number of pests will be affected by various factors such as climate, infectious diseases, so environmental factors can affect the number of pests, the next subsection discusses adjusting the stochastic disturbance coefficient to eliminate pests.
Theorem 4.
For any initial value condition x 0 > 0 , if r σ 1 2 2 < 0 , then lim t x ( t ) = 0 , the pests can be eliminated.
Proof of Theorem 4.
Observe the first formula of the system ( 2 ) , we have following results
d x ( t ) ( r x ( t ) ( 1 x ( t ) k ) ) d t + σ 1 x ( t ) d B 1 ( t )
consider the following formula
d ψ ( t ) = ( r ψ ( t ) ( 1 ψ ( t ) k ) d t + σ 1 ψ ( t ) d B 1 ( t ) ψ ( 0 ) = x ( 0 )
according to the principle of comparison
x ( t ) ψ ( t ) a . s .
the solution of Equation (15) is that
ψ ( t ) = e ( r σ 1 2 2 ) t + σ 1 B 1 ( t ) 1 ψ ( 0 ) + r k 0 t e ( r σ 1 2 2 ) s + σ 1 B 1 ( s ) d s
when r σ 1 2 2 < 0 , both sides are simultaneously divided t, and take the limit
lim t x ( t ) lim t ψ ( 0 ) e ( r σ 1 2 2 ) t + σ 1 B 1 ( t ) = 0
so lim t x ( t ) = 0 . □
Brief summary: According to the conclusion x ( t ) ψ ( t ) = ψ ( 0 ) e ( r σ 1 2 2 ) t + σ 1 B 1 ( t ) , it shows that the number of pests can be controlled by reasonably controlling the strength of the stochastic disturbance, and if σ 1 is large enough, the purpose of pest eradication will be achieved.

4.2. The Extinction of the Populations

Theorem 2 shows that both predators and bait can survive in the system (2) when the parameters meet certain conditions. Theorem 4 shows that it can control and even eliminate pests by controlling the size of the σ 1 . We should also focus on the conditions in which both populations become extinct, and specific proof is given in this section.
Theorem 5.
For any initial value condition ( x 0 , y 0 ) , if r σ 1 2 2 < 0 and ( r σ 1 2 2 ) ( u + D σ 2 2 2 ) ( r + u + D ) 2 4 > 0 , then both predators and prey become extinct, that is lim t x ( t ) = 0 , lim t y ( t ) = 0 .
Proof of Theorem 5.
Define a Lyapunov function V 4 ( x , y ) = ln ( x + y ) , according to the Ito formula
d V 4 ( x , y ) = L V 4 ( x , y ) d t + 1 x + y σ 1 x d B 1 ( t ) + 1 x + y σ 2 y d B 2 ( t )
where
L V 4 ( x , y ) = x x + y ( r r x k y a + x + y ) σ 1 2 x 2 2 ( x + y ) 2 + y x + y ( u x a + x + y + D y ) σ 2 2 y 2 2 ( x + y ) 2 = x ( x + y ) ( r r x k y a + x + y ) + y ( x + y ) ( u x a + x + y + D y ) σ 1 2 x 2 2 σ 2 2 y 2 2 ( x + y ) 2 = r x 2 r x 3 k x 2 y a + x + y + r x y r x 2 y k x y 2 a + x + y + u x 2 y a + x + y + D x y x y 2 + u x y 2 a + x + y + D y 2 y 3 ( x + y ) 2 + σ 1 2 x 2 2 σ 2 2 y 2 2 ( x + y ) 2 = g ( x , y ) ( x + y ) 2
where
g ( x , y ) = r x 2 r x 3 k x 2 y a + x + y + r x y r x 2 y k x y 2 a + x + y + u x 2 y a + x + y + D x y x y 2 + u x y 2 a + x + y + D y 2 y 3 σ 1 2 x 2 2 σ 2 2 y 2 2 ( r σ 1 2 2 ) x 2 + ( r + u + D ) x y + ( u + D σ 2 2 2 ) y 2 = f ( x , y )
where f ( x , y ) = x y B x y , it is a quadratic type, and B is a 2 × 2 matrix as follows
B 2 × 2 = r σ 1 2 2 r + u + D 2 r + u + D 2 u + D σ 2 2 2
because the question set conditions
r σ 1 2 2 < 0 ( r σ 1 2 2 ) ( u + D σ 2 2 2 ) ( r + u + D ) 2 4 > 0
thus the B is a negative definite matrix, and the eigenvalues of the matrix B are all negative, then there must be a reversible matrix P, makes the B = P 1 λ 1 0 0 λ 2 P established, so the f ( x , y ) can also be represented by the following way
f ( x , y ) = λ 1 x 1 2 + λ 2 y 1 2
where λ 1 , λ 2 are the two eigenvalues of the matrix B, and λ 1 < 0 , λ 2 < 0 , in that way
f ( x , y ) = λ 1 x 1 2 + λ 2 y 1 2 λ i x 2 + λ j y 2 λ m a x ( x 2 + y 2 )
where λ m a x = m a x { λ i , λ j } < 0 , i, j = 1 or 2, combined with conclusions (18)–(20), the following result holds
L V 4 ( x , y ) λ m a x ( x 2 + y 2 ) ( x + y ) 2 λ m a x ( x 2 + y 2 ) 2 ( x 2 + y 2 ) = λ m a x 2 < 0
bring (21) into the Formula (17), and then both sides simultaneously integrate in the interval ( 0 , t ) .
ln ( x + y ) ln ( x 0 , y 0 ) + λ m a x 2 t + 0 t 1 x + y σ 1 x d B 1 ( t ) + 0 t 1 x + y σ 2 d B 2 ( t )
then the both sides are simultaneously divided by the time t, and take the limit on t.
lim t 1 t ln ( x + y ) λ m a x 2 + lim t 1 t 0 t 1 x + y σ 1 x d B 1 ( t ) + lim t 1 t 0 t 1 x + y σ 2 y d B 2 ( t )
according to the strong number law [36] of the martingale
lim t 1 t 0 t 1 x + y σ 1 x d B 1 ( t ) = 0 a . s . lim t 1 t 0 t 1 x + y σ 1 y d B 1 ( t ) = 0 a . s .
bring (23) into Formula (22)
lim t s u p 1 t ln ( x + y ) λ m a x 2 < 0 a . s .
since both x ( t ) and y ( t ) are non-negative, so
lim t x ( t ) = 0 a . s . lim t y ( t ) = 0 a . s .
Brief summary: Through the discussion and study of this subsection, when the parameters meet the above conditions, both the predator and the predation become extinct with probability 1. We observe that when other parameters are determined if the two-parameter perturbation coefficient σ 1 and σ 2 are large enough, both groups become extinct, which is the same as the conclusion of Theorem 4.

4.3. Numerical Simulation

Theorem 5 shows that the system populations perish when the parameters meet the above conditions. In order to verify the correctness of the theorem, take the parameter set x 0 = 8 , y 0 = 5 , r = 0.1 , k = 10 , σ 1 = 0.8 , u = 0.1 , a = 10 , D = 0.2 , σ 2 = 1 , Δ = 200 2 16 , it has 0.1 0 . 8 2 2 = 0.22 < 0 and ( 0.4 0 . 8 2 2 ) × ( 0.1 + 0.2 0 . 8 2 2 ) ( 0.1 + 0.1 + 0.2 ) 2 4 = 0.004 > 0 , conditions satisfy the Theorem 5, numerical simulations are shown in Figure 2.

5. Conclusions

In the course of the research on biological mathematics, population size does play an important role in ecosystems. The paper studies a non-linear predation–prey model, in which the number of populations will be affected by mortality rate, growth rate and many other factors, and the population will be affected by the environmental factors, because the environment is always changing, in this paper, we consider the influence of white noise on model.
The paper studies a predation–prey model with Brownian motion, and the article is carried out in three parts. In the first part, the thesis proves that the positive solution is present and unique according to the Lipshiz condition, then we use the proof by contradiction to show that the solution of the system exists globally. In the second part, the article studies the long-term properties of the systematic solution. Firstly, we study the existing conditions of stationary distribution, then we prove that the stationary distribution is ergodic. Secondly, we study the progressive properties of the stochastic systems, we study the progressive behavior of the solution of a stochastic system near the internal equilibrium of its corresponding deterministic system. In the third part, using two methods explores the extinction behavior of the system populations. We also verify the correctness of the conclusions through numerical simulation.
This paper proves the properties of the system through the construction of the Lyapunov function and mathematical derivation according to Ito’s formula, and the conclusions are reached by numerical simulation, making the conclusions scientific and rigorous. Moreover, two different methods are used to prove the population extinction conditions. However, the paper discusses the effect of white noise on the populations except for the delay factors and the effect of Levy noise. I sincerely hope that I can make great progress in the future work.

Author Contributions

Formal analysis, P.L. and W.L.; Funding acquisition, X.T.; Software, H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China Tianyuan Mathematical Foundation (no. 12126363).

Data Availability Statement

Data openly available in a public repository.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Banasiak, J. Multi-scale problems in complex domains—A mathematical framework for systems biology: Comment on “On the interplay between mathematics and biology, hallmarks toward a new systems biology” by N. Bellomo, A. Elaiw, AM Althiabi and MA Alghamdi. Phys. Life Rev. 2015, 12, 65–67. [Google Scholar] [CrossRef] [PubMed]
  2. Grunert, K.; Holden, H.; Jakobsen, E.R.; Stenseth, N.C. Evolutionarily stable strategies in stable and periodically fluctuating populations: The Rosenzweig–MacArthur predator—Prey model. Proc. Natl. Acad. Sci. USA 2021, 118, e2017463118. [Google Scholar] [CrossRef] [PubMed]
  3. Wu, L.; Zheng, H.; Zhang, S. Dynamics of a non-autonomous predator-prey system with Hassell-Varley-Holling? function response and mutual interference. AIMS Math. 2021, 6, 6033–6049. [Google Scholar] [CrossRef]
  4. Calsina, À.; Sanchón, M. Stability and instability of equilibria of an equation of size structured population dynamics. J. Math. Anal. Appl. 2003, 286, 435–452. [Google Scholar] [CrossRef]
  5. Shen, S.H.; Yuan, S.L.; Luo, J.M.; Zhao, Y. Qualitative analysis of a ratio-dependent predator-prey system with Holling-(n+1) functional response. J. Univ. Shanghai Sci. Technol. 2009, 31, 11–13. [Google Scholar]
  6. Bezabih, A.F.; Edessa, G.K.; Rao, K.P. Eco-epidemiological modelling and analysis of prey-predator population. J. Appl. Math. Stat. 2021, 2021, 6679686. [Google Scholar]
  7. Karev, G.P.; Kareva, I.G. Replicator Equations and Models of Biological Populations and Communities. Math. Model. Nat. Phenom. 2014, 9, 68–95. [Google Scholar] [CrossRef]
  8. Lv, J.; Zou, X.; Li, Y. Dynamical properties of a stochastic predator-prey model with functional response. J. Appl. Anal. Comput. 2020, 10, 1242–1255. [Google Scholar] [CrossRef]
  9. Liu, G.; Wang, S.; Yan, J. Positive Periodic Solutions for Neutral Delay Ratio-Dependent Predator-Prey Model with Holling-Tanner Functional Response. Int. J. Math. Math. Sci. 2011, 2011, 376862. [Google Scholar] [CrossRef]
  10. Moaaz, O.; Awrejcewicz, J.; Mahjoub, H. Global behavior and the periodic character of some biological models. Adv. Differ. Equ. 2020, 2020, 410. [Google Scholar] [CrossRef]
  11. Ettaieb, A.; Ouerdiane, H.; Rguigui, H. Quantum White Noise Stochastic Analysis Based on Nuclear Algebras of Entire Functions. Bull. Malays. Math. Sci. Soc. 2021, 44, 599–623. [Google Scholar] [CrossRef]
  12. Wu, Z.; Huang, N.E. A study of the characteristics of white noise using the empirical mode decomposition method. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2004, 460, 1597–1611. [Google Scholar] [CrossRef]
  13. Tang, L.; Zhou, S.; Han, Z. Random exponential attractor for a non-autonomous Zakharov lattice system with multiplicative white noise. J. Differ. Equ. Appl. 2021, 27, 902–921. [Google Scholar] [CrossRef]
  14. dos Santos, M.A.F.; Menon, L., Jr.; Anteneodo, C. Efficiency of random search with space-dependent diffusivity. arXiv 2022, arXiv:2206.14229. [Google Scholar]
  15. Radice, M. Diffusion processes with Gamma-distributed resetting and non-instantaneous returns. J. Phys. A Math. Theor. 2022, 55, 224002. [Google Scholar] [CrossRef]
  16. Toledo-Marin, J.Q.; Boyer, D.; Sevilla, F.J. Predator-prey dynamics: Chasing by stochastic resetting. arXiv 2019, arXiv:1912.02141. [Google Scholar]
  17. Navascués, M. Resetting uncontrolled quantum systems. Phys. Rev. X 2018, 8, 031008. [Google Scholar] [CrossRef]
  18. Wald, S.; Böttcher, L. From classical to quantum walks with stochastic resetting on networks. Phys. Rev. E 2021, 103, 012122. [Google Scholar] [CrossRef]
  19. Reuveni, S.; Urbakh, M.; Klafter, J. Role of substrate unbinding in Michaelis–Menten enzymatic reactions. Proc. Natl. Acad. Sci. USA 2014, 111, 4391–4396. [Google Scholar] [CrossRef]
  20. Roldán, É.; Lisica, A.; Sánchez-Taltavull, D.; Grill, S.W. Stochastic resetting in backtrack recovery by RNA polymerases. Phys. Rev. E 2016, 93, 062411. [Google Scholar] [CrossRef]
  21. Luby, M.; Sinclair, A.; Zuckerman, D. Optimal speedup of Las Vegas algorithms. Inf. Process. Lett. 1993, 47, 173–180. [Google Scholar] [CrossRef]
  22. Stojkoski, V.; Sandev, T.; Kocarev, L.; Pal, A. Geometric Brownian motion under stochastic resetting: A stationary yet nonergodic process. Phys. Rev. E 2021, 104, 014121. [Google Scholar] [CrossRef] [PubMed]
  23. Jiang, G.; Wang, X. Regularity property of solution to two-parameter stochastic volterra equation with non-lipschitz coefficients. Acta Math. Sci. 2013, 33, 872–882. [Google Scholar] [CrossRef]
  24. Chen, P.; Xin, Z.; Zhang, X. Lipschitz stability of nonlinear ordinary differential equations with non-instantaneous impulses in ordered Banach spaces. Int. J. Nonlinear Sci. Numer. Simul. 2021, 22, 657–663. [Google Scholar] [CrossRef]
  25. Li, C.; Liu, S. Stochastic invariance for hybrid stochastic differential equation with non-Lipschitz coefficients. AIMS Math. 2020, 5, 3612–3633. [Google Scholar] [CrossRef]
  26. Liu, Q.; Jiang, D. Dynamics of a stochastic multigroup SEI epidemic model. Stoch. Anal. Appl. 2022, 40, 623–656. [Google Scholar] [CrossRef]
  27. Zhang, Y.; Ma, X.; Din, A. Stationary distribution and extinction of a stochastic SEIQ epidemic model with a general incidence function and temporary immunity. AIMS Math. 2021, 6, 12359–12378. [Google Scholar] [CrossRef]
  28. Duong, M.H. Analysis of the expected density of internal equilibria in random evolutionary multi-player multi-strategy games. J. Math. Biol. 2016, 73, 1727–1760. [Google Scholar] [CrossRef]
  29. Baxter, J.R.; Hernandez, M.N. Stopping time convergence for processes associated with Dirichlet forms. Potential Anal. 2019, 50, 245–277. [Google Scholar] [CrossRef]
  30. Agarwal, R.; Hristova, S.; O’Regan, D. Stability of generalized proportional Caputo fractional differential equations by Lyapunov functions. Fractal Fract. 2022, 6, 34. [Google Scholar] [CrossRef]
  31. Ignatyev, A.O. Lyapunov Function Method for Systems of Difference Equations: Stability with Respect to Part of the Variables. Differ. Equ. 2022, 58, 405–414. [Google Scholar] [CrossRef]
  32. Zhang, S.; Hu, L. The existence of solutions and generalized Lyapunov-type inequalities to boundary value problems of differential equations of variable order. AIMS Math. 2020, 5, 2923–2943. [Google Scholar] [CrossRef]
  33. Alnafisah, Y. A New Approach to Compare the Strong Convergence of the Milstein Scheme with the Approximate Coupling Method. Fractal Fract. 2022, 6, 339. [Google Scholar] [CrossRef]
  34. Ren, Q.; Tian, H. Generalized two-step Milstein methods for stochastic differential equations. Int. J. Comput. Math. 2020, 97, 1363–1379. [Google Scholar] [CrossRef]
  35. Zhang, W.; Yin, X.; Song, M.H.; Liu, M.Z. Convergence rate of the truncated Milstein method of stochastic differential delay equations. Appl. Math. Comput. 2019, 357, 263–281. [Google Scholar] [CrossRef]
  36. Hao, S.; Liu, Q. Convergence rates in the law of large numbers for arrays of martingale differences. J. Math. Anal. Appl. 2014, 417, 733–773. [Google Scholar] [CrossRef]
Figure 1. Stationary distribution of the model (2).
Figure 1. Stationary distribution of the model (2).
Mathematics 10 03238 g001
Figure 2. Model (2) population extinction behavior.
Figure 2. Model (2) population extinction behavior.
Mathematics 10 03238 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tan, X.; Liu, P.; Luo, W.; Chen, H. Analysis of a Class of Predation-Predation Model Dynamics with Random Perturbations. Mathematics 2022, 10, 3238. https://doi.org/10.3390/math10183238

AMA Style

Tan X, Liu P, Luo W, Chen H. Analysis of a Class of Predation-Predation Model Dynamics with Random Perturbations. Mathematics. 2022; 10(18):3238. https://doi.org/10.3390/math10183238

Chicago/Turabian Style

Tan, Xuewen, Pengpeng Liu, Wenhui Luo, and Hui Chen. 2022. "Analysis of a Class of Predation-Predation Model Dynamics with Random Perturbations" Mathematics 10, no. 18: 3238. https://doi.org/10.3390/math10183238

APA Style

Tan, X., Liu, P., Luo, W., & Chen, H. (2022). Analysis of a Class of Predation-Predation Model Dynamics with Random Perturbations. Mathematics, 10(18), 3238. https://doi.org/10.3390/math10183238

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