This paper concerns the dynamics of a stochastic, hybrid delay, one-predator-two-prey model with harvesting and Lévy jumps in a polluted environment. Under some basic assumptions, sufficient conditions of stochastic persistence in the mean and extinction of each species are obtained, as well as the existence of optimal harvesting strategy (OHS). Our results show that both time delays and environmental noises affect the survival state of the species. Moreover, the accurate expressions for the optimal harvesting effort (OHE) and the maximum of expectation of sustainable yield (MESY) are given. Finally, some numerical simulations are provided to support our results.
Citation: Sheng Wang, Baoli Lei. Dynamics of a stochastic hybrid delay one-predator-two-prey model with harvesting and jumps in a polluted environment[J]. Mathematical Modelling and Control, 2025, 5(1): 85-102. doi: 10.3934/mmc.2025007
[1] | Mohammed Ahmed Alomair, Ali Muhib . On the oscillation of fourth-order canonical differential equation with several delays. AIMS Mathematics, 2024, 9(8): 19997-20013. doi: 10.3934/math.2024975 |
[2] | H. Salah, M. Anis, C. Cesarano, S. S. Askar, A. M. Alshamrani, E. M. Elabbasy . Fourth-order differential equations with neutral delay: Investigation of monotonic and oscillatory features. AIMS Mathematics, 2024, 9(12): 34224-34247. doi: 10.3934/math.20241630 |
[3] | Osama Moaaz, Wedad Albalawi . Differential equations of the neutral delay type: More efficient conditions for oscillation. AIMS Mathematics, 2023, 8(6): 12729-12750. doi: 10.3934/math.2023641 |
[4] | Abdelkader Moumen, Amin Benaissa Cherif, Fatima Zohra Ladrani, Keltoum Bouhali, Mohamed Bouye . Fourth-order neutral dynamic equations oscillate on timescales with different arguments. AIMS Mathematics, 2024, 9(9): 24576-24589. doi: 10.3934/math.20241197 |
[5] | Fahd Masood, Osama Moaaz, Shyam Sundar Santra, U. Fernandez-Gamiz, Hamdy A. El-Metwally . Oscillation theorems for fourth-order quasi-linear delay differential equations. AIMS Mathematics, 2023, 8(7): 16291-16307. doi: 10.3934/math.2023834 |
[6] | Clemente Cesarano, Osama Moaaz, Belgees Qaraad, Ali Muhib . Oscillatory and asymptotic properties of higher-order quasilinear neutral differential equations. AIMS Mathematics, 2021, 6(10): 11124-11138. doi: 10.3934/math.2021646 |
[7] | Maryam AlKandari . Nonlinear differential equations with neutral term: Asymptotic behavior of solutions. AIMS Mathematics, 2024, 9(12): 33649-33661. doi: 10.3934/math.20241606 |
[8] | Maged Alkilayh . Nonlinear neutral differential equations of second-order: Oscillatory properties. AIMS Mathematics, 2025, 10(1): 1589-1601. doi: 10.3934/math.2025073 |
[9] | Bouharket Bendouma, Fatima Zohra Ladrani, Keltoum Bouhali, Ahmed Hammoudi, Loay Alkhalifa . Solution-tube and existence results for fourth-order differential equations system. AIMS Mathematics, 2024, 9(11): 32831-32848. doi: 10.3934/math.20241571 |
[10] | Ali Muhib, Hammad Alotaibi, Omar Bazighifan, Kamsing Nonlaopon . Oscillation theorems of solution of second-order neutral differential equations. AIMS Mathematics, 2021, 6(11): 12771-12779. doi: 10.3934/math.2021737 |
This paper concerns the dynamics of a stochastic, hybrid delay, one-predator-two-prey model with harvesting and Lévy jumps in a polluted environment. Under some basic assumptions, sufficient conditions of stochastic persistence in the mean and extinction of each species are obtained, as well as the existence of optimal harvesting strategy (OHS). Our results show that both time delays and environmental noises affect the survival state of the species. Moreover, the accurate expressions for the optimal harvesting effort (OHE) and the maximum of expectation of sustainable yield (MESY) are given. Finally, some numerical simulations are provided to support our results.
In this article, we study the oscillatory behavior of the fourth-order neutral nonlinear differential equation of the form
{(r(t)Φp1[w′′′(t)])′+q(t)Φp2(u(ϑ(t)))=0,r(t)>0, r′(t)≥0, t≥t0>0, | (1.1) |
where w(t):=u(t)+a(t)u(τ(t)) and the first term means the p-Laplace type operator (1<p<∞). The main results are obtained under the following conditions:
L1: Φpi[s]=|s|pi−2s, i=1,2,
L2: r∈C[t0,∞) and under the condition
∫∞t01r1/(p1−1)(s)ds=∞. | (1.2) |
L3: a,q∈C[t0,∞), q(t)>0, 0≤a(t)<a0<∞, τ,ϑ∈C[t0,∞), τ(t)≤t, limt→∞τ(t)=limt→∞ϑ(t)=∞
By a solution of (1.1) we mean a function u ∈C3[tu,∞), tu≥t0, which has the property r(t)(w′′′(t))p1−1∈C1[tu,∞), and satisfies (1.1) on [tu,∞). We assume that (1.1) possesses such a solution. A solution of (1.1) is called oscillatory if it has arbitrarily large zeros on [tu,∞), and otherwise it is called to be nonoscillatory. (1.1) is said to be oscillatory if all its solutions are oscillatory.
We point out that delay differential equations have applications in dynamical systems, optimization, and in the mathematical modeling of engineering problems, such as electrical power systems, control systems, networks, materials, see [1]. The p-Laplace equations have some significant applications in elasticity theory and continuum mechanics.
During the past few years, there has been constant interest to study the asymptotic properties for oscillation of differential equations with p-Laplacian like operator in the canonical case and the noncanonical case, see [2,3,4,11] and the numerical solution of the neutral delay differential equations, see [5,6,7]. The oscillatory properties of differential equations are fairly well studied by authors in [16,17,18,19,20,21,22,23,24,25,26,27]. We collect some relevant facts and auxiliary results from the existing literature.
Liu et al. [4] studied the oscillation of even-order half-linear functional differential equations with damping of the form
{(r(t)Φ(y(n−1)(t)))′+a(t)Φ(y(n−1)(t))+q(t)Φ(y(g(t)))=0,Φ=|s|p−2s, t≥t0>0, |
where n is even. This time, the authors used comparison method with second order equations.
The authors in [9,10] have established sufficient conditions for the oscillation of the solutions of
{(r(t)|y(n−1)(t)|p−2y(n−1)(t))′+∑ji=1qi(t)g(y(ϑi(t)))=0,j≥1, t≥t0>0, |
where n is even and p>1 is a real number, in the case where ϑi(t)≥υ (with r∈C1((0,∞),R), qi∈C([0,∞),R), i=1,2,..,j).
We point out that Li et al. [3] using the Riccati transformation together with integral averaging technique, focuses on the oscillation of equation
{(r(t)|w′′′(t)|p−2w′′′(t))′+∑ji=1qi(t)|y(δi(t))|p−2y(δi(t))=0,1<p<∞, , t≥t0>0. |
Park et al. [8] have obtained sufficient conditions for oscillation of solutions of
{(r(t)|y(n−1)(t)|p−2y(n−1)(t))′+q(t)g(y(δ(t)))=0,1<p<∞, , t≥t0>0. |
As we already mentioned in the Introduction, our aim here is complement results in [8,9,10]. For this purpose we discussed briefly these results.
In this paper, we obtain some new oscillation criteria for (1.1). The paper is organized as follows. In the next sections, we will mention some auxiliary lemmas, also, we will use the generalized Riccati transformation technique to give some sufficient conditions for the oscillation of (1.1), and we will give some examples to illustrate the main results.
For convenience, we denote
A(t)=q(t)(1−a0)p2−1Mp1−p2(ϑ(t)), B(t)=(p1−1)εϑ2(t)ζϑ′(t)r1/(p1−1)(t), ϕ1(t)=∫∞tA(s)ds,R1(t):=(p1−1)μt22r1/(p1−1)(t),ξ(t):=q(t)(1−a0)p2−1Mp2−p11ε1(ϑ(t)t)3(p2−1),η(t):=(1−a0)p2/p1Mp2/(p1−2)2∫∞t(1r(δ)∫∞δq(s)ϑp2−1(s)sp2−1ds)1/(p1−1)dδ,ξ∗(t)=∫∞tξ(s)ds, η∗(t)=∫∞tη(s)ds, |
for some μ∈(0,1) and every M1,M2 are positive constants.
Definition 1. A sequence of functions {δn(t)}∞n=0 and {σn(t)}∞n=0 as
δ0(t)=ξ∗(t), and σ0(t)=η∗(t),δn(t)=δ0(t)+∫∞tR1(t)δp1/(p1−1)n−1(s)ds, n>1σn(t)=σ0(t)+∫∞tσp1/(p1−1)n−1(s)ds, n>1. | (2.1) |
We see by induction that δn(t)≤δn+1(t) and σn(t)≤σn+1(t) for t≥t0, n>1.
In order to discuss our main results, we need the following lemmas:
Lemma 2.1. [12] If the function w satisfies w(i)(ν)>0, i=0,1,...,n, and w(n+1)(ν)<0 eventually. Then, for every ε1∈(0,1), w(ν)/w′(ν)≥ε1ν/n eventually.
Lemma 2.2. [13] Let u(t) be a positive and n-times differentiable function on an interval [T,∞) with its nth derivative u(n)(t) non-positive on [T,∞) and not identically zero on any interval of the form [T′,∞), T′≥T and u(n−1)(t)u(n)(t)≤0, t≥tu then there exist constants θ, 0<θ<1 and ε>0 such that
u′(θt)≥εtn−2u(n−1)(t), |
for all sufficient large t.
Lemma 2.3 [14] Let u∈Cn([t0,∞),(0,∞)). Assume that u(n)(t) is of fixed sign and not identically zero on [t0,∞) and that there exists a t1≥t0 such that u(n−1)(t)u(n)(t)≤0 for all t≥t1. If limt→∞u(t)≠0, then for every μ∈(0,1) there exists tμ≥t1 such that
u(t)≥μ(n−1)!tn−1|u(n−1)(t)| for t≥tμ. |
Lemma 2.4. [15] Assume that (1.2) holds and u is an eventually positive solution of (1.1). Then, (r(t)(w′′′(t))p1−1)′<0 and there are the following two possible cases eventually:
(G1) w(k)(t)>0, k=1,2,3,(G2) w(k)(t)>0, k=1,3, and w′′(t)<0. |
Theorem 2.1. Assume that
liminft→∞1ϕ1(t)∫∞tB(s)ϕp1(p1−1)1(s)ds>p1−1pp1(p1−1)1. | (2.2) |
Then (1.1) is oscillatory.
proof. Assume that u be an eventually positive solution of (1.1). Then, there exists a t1≥t0 such that u(t)>0, u(τ(t))>0 and u(ϑ(t))>0 for t≥t1. Since r′(t)>0, we have
w(t)>0, w′(t)>0, w′′′(t)>0, w(4)(t)<0 and (r(t)(w′′′(t))p1−1)′≤0, | (2.3) |
for t≥t1. From definition of w, we get
u(t)≥w(t)−a0u(τ(t))≥w(t)−a0w(τ(t))≥(1−a0)w(t), |
which with (1.1) gives
(r(t)(w′′′(t))p1−1)′≤−q(t)(1−a0)p2−1wp2−1(ϑ(t)). | (2.4) |
Define
ϖ(t):=r(t)(w′′′(t))p1−1wp1−1(ζϑ(t)). | (2.5) |
for some a constant ζ∈(0,1). By differentiating and using (2.4), we obtain
ϖ′(t)≤−q(t)(1−a0)p2−1wp2−1(ϑ(t)).wp1−1(ζϑ(t))−(p1−1)r(t)(w′′′(t))p1−1w′(ζϑ(t))ζϑ′(t)wp1(ζϑ(t)). |
From Lemma 2.2, there exist constant ε>0, we have
ϖ′(t)≤−q(t)(1−a0)p2−1wp2−p1(ϑ(t))−(p1−1)r(t)(w′′′(t))p1−1εϑ2(t)w′′′(ϑ(t))ζϑ′(t)wp1(ζϑ(t)). |
Which is
ϖ′(t)≤−q(t)(1−a0)p2−1wp2−p1(ϑ(t))−(p1−1)εr(t)ϑ2(t)ζϑ′(t)(w′′′(t))p1wp1(ζϑ(t)), |
by using (2.5) we have
ϖ′(t)≤−q(t)(1−a0)p2−1wp2−p1(ϑ(t))−(p1−1)εϑ2(t)ζϑ′(t)r1/(p1−1)(t)ϖp1/(p1−1)(t). | (2.6) |
Since w′(t)>0, there exist a t2≥t1 and a constant M>0 such that
w(t)>M. |
Then, (2.6), turns to
ϖ′(t)≤−q(t)(1−a0)p2−1Mp2−p1(ϑ(t))−(p1−1)εϑ2(t)ζϑ′(t)r1/(p1−1)(t)ϖp1/(p1−1)(t), |
that is
ϖ′(t)+A(t)+B(t)ϖp1/(p1−1)(t)≤0. |
Integrating the above inequality from t to l, we get
ϖ(l)−ϖ(t)+∫ltA(s)ds+∫ltB(s)ϖp1/(p1−1)(s)ds≤0. |
Letting l→∞ and using ϖ>0 and ϖ′<0, we have
ϖ(t)≥ϕ1(t)+∫∞tB(s)ϖp1/(p1−1)(s)ds. |
This implies
ϖ(t)ϕ1(t)≥1+1ϕ1(t)∫∞tB(s)ϕp1/(p1−1)1(s)(ϖ(s)ϕ1(s))p1/(p1−1)ds. | (2.7) |
Let λ=inft≥Tϖ(t)/ϕ1(t) then obviously λ≥1. Thus, from (2.2) and (2.7) we see that
λ≥1+(p1−1)(λp1)p1/(p1−1) |
or
λp1≥1p1+(p1−1)p1(λp1)p1/(p1−1), |
which contradicts the admissible value of λ≥1 and (p1−1)>0.
Therefore, the proof is complete.
Theorem 2.2. Assume that
liminft→∞1ξ∗(t)∫∞tR1(s)ξp1/(p1−1)∗(s)ds>(p1−1)pp1/(p1−1)1 | (2.8) |
and
liminft→∞1η∗(t)∫∞t0η2∗(s)ds>14. | (2.9) |
Then (1.1) is oscillatory.
proof. Assume to the contrary that (1.1) has a nonoscillatory solution in [t0,∞). Without loss of generality, we let u be an eventually positive solution of (1.1). Then, there exists a t1≥t0 such that u(t)>0, u(τ(t))>0 and u(ϑ(t))>0 for t≥t1. From Lemma 2.4 there is two cases (G1) and (G2).
For case (G1). Define
ω(t):=r(t)(w′′′(t))p1−1wp1−1(t). |
By differentiating ω and using (2.4), we obtain
ω′(t)≤−q(t)(1−a0)p2−1wp2−1(ϑ(t))wp1−1(t)−(p1−1)r(t)(w′′′(t))p1−1wp1(t)w′(t). | (2.10) |
From Lemma 2.1, we get
w′(t)w(t)≤3ε1t. |
Integrating again from t to ϑ(t), we find
w(ϑ(t))w(t)≥ε1ϑ3(t)t3. | (2.11) |
It follows from Lemma 2.3 that
w′(t)≥μ12t2w′′′(t), | (2.12) |
for all μ1∈(0,1) and every sufficiently large t. Since w′(t)>0, there exist a t2≥t1 and a constant M>0 such that
w(t)>M, | (2.13) |
for t≥t2. Thus, by (2.10), (2.11), (2.12) and (2.13), we get
ω′(t)+q(t)(1−a0)p2−1Mp2−p11ε1(ϑ(t)t)3(p2−1)+(p1−1)μt22r1/(p1−1)(t)ωp1/(p1−1)(t)≤0, |
that is
ω′(t)+ξ(t)+R1(t)ωp1/(p1−1)(t)≤0. | (2.14) |
Integrating (2.14) from t to l, we get
ω(l)−ω(t)+∫ltξ(s)ds+∫ltR1(s)ωp1/(p1−1)(s)ds≤0. |
Letting l→∞ and using ω>0 and ω′<0, we have
ω(t)≥ξ∗(t)+∫∞tR1(s)ωp1/(p1−1)(s)ds. | (2.15) |
This implies
ω(t)ξ∗(t)≥1+1ξ∗(t)∫∞tR1(s)ξp1/(p1−1)∗(s)(ω(s)ξ∗(s))p1/(p1−1)ds. | (2.16) |
Let λ=inft≥Tω(t)/ξ∗(t) then obviously λ≥1. Thus, from (2.8) and (2.16) we see that
λ≥1+(p1−1)(λp1)p1/(p1−1) |
or
λp1≥1p1+(p1−1)p1(λp1)p1/(p1−1), |
which contradicts the admissible value of λ≥1 and (p1−1)>0.
For case (G2). Integrating (2.4) from t to m, we obtain
r(m)(w′′′(m))p1−1−r(t)(w′′′(t))p1−1≤−∫mtq(s)(1−a0)p2−1wp2−1(ϑ(s))ds. | (2.17) |
From Lemma 2.1, we get that
w(t)≥ε1tw′(t) and hence w(ϑ(t))≥ε1ϑ(t)tw(t). | (2.18) |
For (2.17), letting m→∞and using (2.18), we see that
r(t)(w′′′(t))p1−1≥ε1(1−a0)p2−1wp2−1(t)∫∞tq(s)ϑp2−1(s)sp2−1ds. |
Integrating this inequality again from t to ∞, we get
w′′(t)≤−ε1(1−a0)p2/p1wp2/p1(t)∫∞t(1r(δ)∫∞δq(s)ϑp2−1(s)sp2−1ds)1/(p1−1)dδ, | (2.19) |
for all ε1∈(0,1). Define
y(t)=w′(t)w(t). |
By differentiating y and using (2.13) and (2.19), we find
y′(t)=w′′(t)w(t)−(w′(t)w(t))2≤−y2(t)−(1−a0)p2/p1M(p2/p1)−1∫∞t(1r(δ)∫∞δq(s)ϑp2−1(s)sp2−1ds)1/(p1−1)dδ, | (2.20) |
hence
y′(t)+η(t)+y2(t)≤0. | (2.21) |
The proof of the case where (G2) holds is the same as that of case (G1). Therefore, the proof is complete.
Theorem 2.3. Let δn(t) and σn(t) be defined as in (2.1). If
limsupt→∞(μ1t36r1/(p1−1)(t))p1−1δn(t)>1 | (2.22) |
and
limsupt→∞λtσn(t)>1, | (2.23) |
for some n, then (1.1)is oscillatory.
proof. Assume to the contrary that (1.1) has a nonoscillatory solution in [t0,∞). Without loss of generality, we let u be an eventually positive solution of (1.1). Then, there exists a t1≥t0 such that u(t)>0, u(τ(t))>0 and u(ϑ(t))>0 for t≥t1. From Lemma 2.4 there is two cases.
In the case (G1), proceeding as in the proof of Theorem 2.2, we get that (2.12) holds. It follows from Lemma 2.3 that
w(t)≥μ16t3w′′′(t). | (2.24) |
From definition of ω(t) and (2.24), we have
1ω(t)=1r(t)(w(t)w′′′(t))p1−1≥1r(t)(μ16t3)p1−1. |
Thus,
ω(t)(μ1t36r1/(p1−1)(t))p1−1≤1. |
Therefore,
limsupt→∞ω(t)(μ1t36r1/(p1−1)(t))p1−1≤1, |
which contradicts (2.22).
The proof of the case where (G2) holds is the same as that of case (G1). Therefore, the proof is complete.
Corollary 2.1. Let δn(t) and σn(t) be defined as in (2.1). If
∫∞t0ξ(t)exp(∫tt0R1(s)δ1/(p1−1)n(s)ds)dt=∞ | (2.25) |
and
∫∞t0η(t)exp(∫tt0σ1/(p1−1)n(s)ds)dt=∞, | (2.26) |
for some n, then (1.1) is oscillatory.
proof. Assume to the contrary that (1.1) has a nonoscillatory solution in [t0,∞). Without loss of generality, we let u be an eventually positive solution of (1.1). Then, there exists a t1≥t0 such that u(t)>0, u(τ(t))>0 and u(ϑ(t))>0 for t≥t1. From Lemma 2.4 there is two cases (G1) and (G2).
In the case (G1), proceeding as in the proof of Theorem 2, we get that (2.15) holds. It follows from (2.15) that ω(t)≥δ0(t). Moreover, by induction we can also see that ω(t)≥δn(t) for t≥t0, n>1. Since the sequence {δn(t)}∞n=0 monotone increasing and bounded above, it converges to δ(t). Thus, by using Lebesgue's monotone convergence theorem, we see that
δ(t)=limn→∞δn(t)=∫∞tR1(t)δp1/(p1−1)(s)ds+δ0(t) |
and
δ′(t)=−R1(t)δp1/(p1−1)(t)−ξ(t). | (2.27) |
Since δn(t)≤δ(t), it follows from (2.27) that
δ′(t)≤−R1(t)δ1/(p1−1)n(t)δ(t)−ξ(t). |
Hence, we get
δ(t)≤exp(−∫tTR1(s)δ1/(p1−1)n(s)ds)(δ(T)−∫tTξ(s)exp(∫sTR1(δ)δ1/(p1−1)n(δ)dδ)ds). |
This implies
∫tTξ(s)exp(∫sTR1(δ)δ1/(p1−1)n(δ)dδ)ds≤δ(T)<∞, |
which contradicts (2.25). The proof of the case where (G2) holds is the same as that of case (G1). Therefore, the proof is complete.
Example 2.1. Consider the differential equation
(u(t)+12u(t2))(4)+q0t4u(t3)=0, | (2.28) |
where q0>0 is a constant. Let p1=p2=2, r(t)=1, a(t)=1/2, τ(t)=t/2, ϑ(t)=t/3 and q(t)=q0/t4. Hence, it is easy to see that
A(t)=q(t)(1−a0)(p2−1)Mp2−p1(ϑ(t))=q02t4, B(t)=(p1−1)εϑ2(t)ζϑ′(t)r1/(p1−1)(t)=εt227 |
and
ϕ1(t)=q06t3, |
also, for some ε>0, we find
liminft→∞1ϕ1(t)∫∞tB(s)ϕp1/(p1−1)1(s)ds>(p1−1)pp1/(p1−1)1.liminft→∞6εq0t3972∫∞tdss4>14q0>121.5ε. |
Hence, by Theorem 2.1, every solution of Eq (2.28) is oscillatory if q0>121.5ε.
Example 2.2. Consider a differential equation
(u(t)+a0u(τ0t))(n)+q0tnu(ϑ0t)=0, | (2.29) |
where q0>0 is a constant. Note that p=2, t0=1, r(t)=1, a(t)=a0, τ(t)=τ0t, ϑ(t)=ϑ0t and q(t)=q0/tn.
Easily, we see that condition (2.8) holds and condition (2.9) satisfied.
Hence, by Theorem 2.2, every solution of Eq (2.29) is oscillatory.
Remark 2.1. Finally, we point out that continuing this line of work, we can have oscillatory results for a fourth order equation of the type:
{(r(t)|y′′′(t)|p1−2y′′′(t))′+a(t)f(y′′′(t))+∑ji=1qi(t)|y(σi(t))|p2−2y(σi(t))=0,t≥t0, σi(t)≤t, j≥1,, 1<p2≤p1<∞. |
The paper is devoted to the study of oscillation of fourth-order differential equations with p-Laplacian like operators. New oscillation criteria are established by using a Riccati transformations, and they essentially improves the related contributions to the subject.
Further, in the future work we get some Hille and Nehari type and Philos type oscillation criteria of (1.1) under the condition ∫∞υ01r1/(p1−1)(s)ds<∞.
The authors express their debt of gratitude to the editors and the anonymous referee for accurate reading of the manuscript and beneficial comments.
The author declares that there is no competing interest.
[1] |
M. Liu, X. He, J. Yu, Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays, Nonlinear Anal., 28 (2018), 87–104. https://doi.org/10.1016/j.nahs.2017.10.004 doi: 10.1016/j.nahs.2017.10.004
![]() |
[2] |
G. Liu, X. Meng, Optimal harvesting strategy for a stochastic mutualism system in a polluted environment with regime switching, Phys. A, 536 (2019), 120893. https://doi.org/10.1016/j.physa.2019.04.129 doi: 10.1016/j.physa.2019.04.129
![]() |
[3] |
N. Cao, X. L. Fu, Stationary distribution and extinction of a Lotka-Volterra model with distribute delay and nonlinear stochastic perturbations, Chaos Solitons Fract., 169 (2023), 113246. https://doi.org/10.1016/j.chaos.2023.113246 doi: 10.1016/j.chaos.2023.113246
![]() |
[4] | Y. Kuang, Delay differential equations: with applications in population dynamics, Academic Press, 1993. |
[5] |
W. Zuo, D. Jiang, X. Sun, T. Hayat, A. Alsaedi, Long-time behaviors of a stochastic cooperative Lotka-Volterra system with distributed delay, Phys. A, 506 (2018), 542–559. https://doi.org/10.1016/j.physa.2018.03.071 doi: 10.1016/j.physa.2018.03.071
![]() |
[6] |
H. J. Alsakaji, S. Kundu, F. A. Rihan, Delay differential model of one-predator two-prey system with Monod-Haldane and holling type II functional responses, Appl. Math. Comput., 397 (2021), 125919. https://doi.org/10.1016/j.amc.2020.125919 doi: 10.1016/j.amc.2020.125919
![]() |
[7] |
S. Wang, G. Hu, T. Wei, On a three-species stochastic hybrid Lotka-Volterra system with distributed delay and Lévy noise, Filomat, 36 (2022), 4737–4750. https://doi.org/10.2298/FIL2214737W doi: 10.2298/FIL2214737W
![]() |
[8] |
Q. Liu, D. Jiang, Influence of the fear factor on the dynamics of a stochastic predator-prey model, Appl. Math. Lett., 112 (2021), 106756. https://doi.org/10.1016/j.aml.2020.106756 doi: 10.1016/j.aml.2020.106756
![]() |
[9] |
Q. Zhang, D. Jiang, Z. Liu, D. O'Regan, Asymptotic behavior of a three species eco-epidemiological model perturbed by white noise, J. Math. Anal. Appl., 433 (2016), 121–148. https://doi.org/10.1016/j.jmaa.2015.07.025 doi: 10.1016/j.jmaa.2015.07.025
![]() |
[10] |
Y. Zhao, L. You, D. Burkow, S. Yuan, Optimal harvesting strategy of a stochastic inshore-offshore hairtail fishery model driven by Lévy jumps in a polluted environment, Nonlinear Dyn., 95 (2019), 1529–1548. https://doi.org/10.1007/s11071-018-4642-y doi: 10.1007/s11071-018-4642-y
![]() |
[11] |
N. Tuerxun, Z. Teng, Global dynamics in stochastic n-species food chain systems with white noise and Lévy jumps, Math. Methods Appl. Sci., 45 (2022), 5184–5214. https://doi.org/10.1002/mma.8101 doi: 10.1002/mma.8101
![]() |
[12] |
Q. Yang, X. Zhang, D. Jiang, Dynamical behaviors of a stochastic food chain system with Ornstein-Uhlenbeck process, J. Nonlinear Sci., 32 (2022), 34. https://doi.org/10.1007/s00332-022-09796-8 doi: 10.1007/s00332-022-09796-8
![]() |
[13] |
X. Yu, S. Yuan, T. Zhang, Persistence and ergodicity of a stochastic single species model with Allee effect under regime switching, Commun. Nonlinear Sci., 59 (2018), 359–374. https://doi.org/10.1016/j.cnsns.2017.11.028 doi: 10.1016/j.cnsns.2017.11.028
![]() |
[14] |
M. Liu, X. He, J. Yu, Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays, Nonlinear Anal., 28 (2018), 87–104. https://doi.org/10.1016/j.nahs.2017.10.004 doi: 10.1016/j.nahs.2017.10.004
![]() |
[15] |
M. Liu, Y. Zhu, Stationary distribution and ergodicity of a stochastic hybrid competition model with Lévy jumps, Nonlinear Anal., 30 (2018), 225–239. https://doi.org/10.1016/j.nahs.2018.05.002 doi: 10.1016/j.nahs.2018.05.002
![]() |
[16] |
J. Bao, C. Yuan, Stochastic population dynamics driven by Lévy noise, J. Math. Anal. Appl., 391 (2012), 363–375. https://doi.org/10.1016/j.jmaa.2012.02.043 doi: 10.1016/j.jmaa.2012.02.043
![]() |
[17] |
M. Liu, K. Wang, Stochastic Lotka-Volterra systems with Lévy noise, J. Math. Anal. Appl., 410 (2014), 750–763. https://doi.org/10.1016/j.jmaa.2013.07.078 doi: 10.1016/j.jmaa.2013.07.078
![]() |
[18] |
Y. Zhao, S. Yuan, Optimal harvesting policy of a stochastic two-species competitive model with Lévy noise in a polluted environment, Phys. A, 477 (2017), 20–33. https://doi.org/10.1016/j.physa.2017.02.019 doi: 10.1016/j.physa.2017.02.019
![]() |
[19] |
Q. Liu, D. Jiang, N. Shi, T. Hayat, A. Alsaedi, Stochastic mutualism model with Lévy jumps, Commun. Nonlinear Sci., 43 (2017), 78–90. https://doi.org/10.1016/j.cnsns.2016.05.003 doi: 10.1016/j.cnsns.2016.05.003
![]() |
[20] |
M. Liu, Dynamics of a stochastic regime-switching predator-prey model with modified Leslie-Gower Holling-type II schemes and prey harvesting, Nonlinear Dyn., 96 (2019), 417–442. https://doi.org/10.1007/s11071-019-04797-x doi: 10.1007/s11071-019-04797-x
![]() |
[21] |
Q. Han, D. Jiang, C. Ji, Analysis of a delayed stochastic predator-prey model in a polluted environment, Appl. Math. Model., 38 (2014), 3067–3080. https://doi.org/10.1016/j.apm.2013.11.014 doi: 10.1016/j.apm.2013.11.014
![]() |
[22] |
Q. Liu, Q. Chen, Analysis of a stochastic delay predator-prey system with jumps in a polluted environment, Appl. Math. Comput., 242 (2014), 90–100. https://doi.org/10.1016/j.amc.2014.05.033 doi: 10.1016/j.amc.2014.05.033
![]() |
[23] |
H. Qiu, W. Deng, Optimal harvesting of a stochastic delay logistic model with Lévy jumps, J. Phys. A, 49 (2016), 405601. https://doi.org/10.1088/1751-8113/49/40/405601 doi: 10.1088/1751-8113/49/40/405601
![]() |
[24] |
M. Liu, C. Bai, Optimal harvesting of a stochastic mutualism model with Lévy jumps, Appl. Math. Comput., 276 (2016), 301–309. https://doi.org/10.1016/j.amc.2015.11.089 doi: 10.1016/j.amc.2015.11.089
![]() |
[25] |
M. Liu, X. He, J. Yu, Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays, Nonlinear Anal., 28 (2018), 87–104. https://doi.org/10.1016/j.nahs.2017.10.004 doi: 10.1016/j.nahs.2017.10.004
![]() |
[26] |
S. K. Mandal, S. Poria, A study of Michaelis-Menten type harvesting effects on a population in stochastic environment, Phys. A, 611 (2023), 128469. https://doi.org/10.1016/j.physa.2023.128469 doi: 10.1016/j.physa.2023.128469
![]() |
[27] |
S. Wang, L. Dong, Stochastic dynamics of a hybrid delay food chain model with harvesting and jumps in a polluted environment, Methodol. Comput. Appl., 25 (2023), 94. https://doi.org/10.1007/s11009-023-10064-9 doi: 10.1007/s11009-023-10064-9
![]() |
[28] |
J. Yu, M. Liu, Stationary distribution and ergodicity of a stochastic food-chain model with Lévy jumps, Phys. A, 482 (2017), 14–28. https://doi.org/10.1016/j.physa.2017.04.067 doi: 10.1016/j.physa.2017.04.067
![]() |
[29] |
M. Liu, J. Yu, P. Mandal, Dynamics of a stochastic delay competitive model with harvesting and Markovian switching, Appl. Math. Comput., 337 (2018), 335–349. https://doi.org/10.1016/j.amc.2018.03.044 doi: 10.1016/j.amc.2018.03.044
![]() |
[30] |
Y. Deng, M. Liu, Analysis of a stochastic tumor-immune model with regime switching and impulsive perturbations, Appl. Math. Model., 78 (2020), 482–504. https://doi.org/10.1016/j.apm.2019.10.010 doi: 10.1016/j.apm.2019.10.010
![]() |
[31] |
M. Liu, C. Bai, On a stochastic delayed predator-prey model with Lévy jumps, Appl. Math. Comput., 228 (2014), 563–570. https://doi.org/10.1016/j.amc.2013.12.026 doi: 10.1016/j.amc.2013.12.026
![]() |
[32] |
N. Tuerxun, Z. Teng, Global dynamics in stochastic n-species food chain systems with white noise and Lévy jumps, Math. Methods Appl. Sci., 45 (2022), 5184–5214. https://doi.org/10.1002/mma.8101 doi: 10.1002/mma.8101
![]() |
[33] |
M. Liu, K. Wang, Survival analysis of stochastic single-species population models in polluted environments, Ecol. Model., 220 (2009), 1347–1357. https://doi.org/10.1016/j.ecolmodel.2009.03.001 doi: 10.1016/j.ecolmodel.2009.03.001
![]() |
[34] |
S. Wang, L. Wang, T. Wei, Optimal harvesting for a stochastic logistic model with S-type distributed time delay, J. Differ. Equations Appl., 23 (2017), 618–632. https://doi.org/10.1080/10236198.2016.1269761 doi: 10.1080/10236198.2016.1269761
![]() |
[35] |
M. Liu, K. Wang, Q. Wu, Survival analysis of stochastic competitive models in a polluted environment and stochastic competitive exclusion principle, Bull. Math. Biol., 73 (2011), 1969–2012. https://doi.org/10.1007/s11538-010-9569-5 doi: 10.1007/s11538-010-9569-5
![]() |
[36] |
M. Kinnally, R. Williams, On existence and uniqueness of stationary distributions for stochastic delay differential equations with positivity constraints, Electron. J. Probab., 15 (2010), 409–451. https://doi.org/10.1214/ejp.v15-756 doi: 10.1214/ejp.v15-756
![]() |
[37] |
M. Hairer, J. C. Mattingly, M. Scheutzow, Asymptotic coupling and a general form of Harris' theorem with applications to stochastic delay equations, Probab. Theory Related Fields, 149 (2011), 223–259. https://doi.org/10.1007/s00440-009-0250-6 doi: 10.1007/s00440-009-0250-6
![]() |
[38] | G. Prato, J. Zabczyk, Ergodicity for infinite dimensional systems, Cambridge University Press, 1996. https://doi.org/10.1017/CBO9780511662829 |
1. | Omar Saber Qasim, Ahmed Entesar, Waleed Al-Hayani, Solving nonlinear differential equations using hybrid method between Lyapunov's artificial small parameter and continuous particle swarm optimization, 2021, 0, 2155-3297, 0, 10.3934/naco.2021001 | |
2. | Fahd Masood, Osama Moaaz, Shyam Sundar Santra, U. Fernandez-Gamiz, Hamdy A. El-Metwally, Oscillation theorems for fourth-order quasi-linear delay differential equations, 2023, 8, 2473-6988, 16291, 10.3934/math.2023834 | |
3. | Peng E, Tingting Xu, Linhua Deng, Yulin Shan, Miao Wan, Weihong Zhou, Solutions of a class of higher order variable coefficient homogeneous differential equations, 2025, 20, 1556-1801, 213, 10.3934/nhm.2025011 |