Optimal harvesting strategy for stochastic hybrid delay Lotka-Volterra systems with Lévy noise in a polluted environment

Sheng Wang; Lijuan Dong; Zeyan Yue; Sheng Wang; Lijuan Dong; Zeyan Yue

doi:10.3934/mbe.2023263

Mathematical Biosciences and Engineering

2023, Volume 20, Issue 4: 6084-6109. doi: 10.3934/mbe.2023263

Previous Article Next Article

Research article Special Issues

Optimal harvesting strategy for stochastic hybrid delay Lotka-Volterra systems with Lévy noise in a polluted environment

School of Mathematics and Information Science, Henan Polytechnic University, Jiaozuo 454003, China

Academic Editor: Xiaodi Li

Received: 21 September 2022 Revised: 13 December 2022 Accepted: 27 December 2022 Published: 31 January 2023

This paper concerns the dynamics of two stochastic hybrid delay Lotka-Volterra systems with harvesting and Lévy noise in a polluted environment (i.e., predator-prey system and competitive system). For every system, sufficient and necessary conditions for persistence in mean and extinction of each species are established. Then, sufficient conditions for global attractivity of the systems are obtained. Finally, sufficient and necessary conditions for the existence of optimal harvesting strategy are provided. The accurate expressions for the optimal harvesting effort (OHE) and the maximum of expectation of sustainable yield (MESY) are given. Our results show that the dynamic behaviors and optimal harvesting strategy are closely correlated with both time delays and three types of environmental noises (namely white Gaussian noises, telephone noises and Lévy noises).

Keywords:

Citation: Sheng Wang, Lijuan Dong, Zeyan Yue. Optimal harvesting strategy for stochastic hybrid delay Lotka-Volterra systems with Lévy noise in a polluted environment[J]. Mathematical Biosciences and Engineering, 2023, 20(4): 6084-6109. doi: 10.3934/mbe.2023263

Related Papers:

[1]	Mati ur Rahman, Mehmet Yavuz, Muhammad Arfan, Adnan Sami . Theoretical and numerical investigation of a modified ABC fractional operator for the spread of polio under the effect of vaccination. AIMS Biophysics, 2024, 11(1): 97-120. doi: 10.3934/biophy.2024007
[2]	Shaimaa A. M. Abdelmohsen, D. Sh. Mohamed, Haifa A. Alyousef, M. R. Gorji, Amr M. S. Mahdy . Mathematical modeling for solving fractional model cancer bosom malignant growth. AIMS Biophysics, 2023, 10(3): 263-280. doi: 10.3934/biophy.2023018
[3]	Mohammed Alabedalhadi, Mohammed Shqair, Ibrahim Saleh . Analysis and analytical simulation for a biophysical fractional diffusive cancer model with virotherapy using the Caputo operator. AIMS Biophysics, 2023, 10(4): 503-522. doi: 10.3934/biophy.2023028
[4]	Ken Takahashi, Takayuki Oda, Keiji Naruse . Coarse-grained molecular dynamics simulations of biomolecules. AIMS Biophysics, 2014, 1(1): 1-15. doi: 10.3934/biophy.2014.1.1
[5]	Nor Afiqah Mohd Aris, Siti Suhana Jamaian . Dynamical analysis of fractional-order chemostat model. AIMS Biophysics, 2021, 8(2): 182-197. doi: 10.3934/biophy.2021014
[6]	Alexander Galashev . Computer implementation of the method for electrolytic production of thin films for biomedical applications: short review. AIMS Biophysics, 2024, 11(1): 39-65. doi: 10.3934/biophy.2024004
[7]	Dorota Kantor, Gabriela Kanclerz, Grzegorz Tatoń . Is the occurrence of electromagnetic hypersensitivity related to sensitivity to air pollution and weather factors?. AIMS Biophysics, 2025, 12(2): 259-272. doi: 10.3934/biophy.2025014
[8]	Soumaya Eltifi-Ghanmi, Samiha Amara, Bessem Mkaouer . Kinetic and kinematic analysis of three kicks in Sanda Wushu. AIMS Biophysics, 2025, 12(2): 174-196. doi: 10.3934/biophy.2025011
[9]	Mehmet Yavuz, Fuat Usta . Importance of modelling and simulation in biophysical applications. AIMS Biophysics, 2023, 10(3): 258-262. doi: 10.3934/biophy.2023017
[10]	Carlo Bianca . Mathematical and computational modeling of biological systems: advances and perspectives. AIMS Biophysics, 2021, 8(4): 318-321. doi: 10.3934/biophy.2021025

Abstract

1. Introduction

In line with the Environmental Health Association of America's committee definition, the environment encompasses the habitats where humans reside, work, and engage in recreational activities. It includes the air we breathe, the food we consume, the water we drink, and the shelter we seek for protection against the elements. Additionally, it comprises harmful pollutants and other environmental factors that significantly impact our well-being and overall health. Our survival depends on access to nutritious food, clean air, and safe drinking water, which necessitates the promotion of a healthy atmosphere. The World Health Organization (WHO) emphasizes the importance of maintaining a clean environment that adheres to health standards. Regrettably, our water supply has become increasingly contaminated, with numerous pollutants affecting the quality of our drinking water sources. Thus, it is imperative that we develop a conscientious awareness of our surroundings. Environmental science encompasses scientific methodologies, economic considerations, and political interventions. Moreover, toxic substances introduced into lakes, streams, rivers, seas, and other bodies of water can dissolve, remain suspended, or settle at the bottom. Consequently, water contamination ensues, leading to a decline in groundwater quality and the deterioration of marine habitats [1].

Pollution can also seep into the groundwater and alter the sediment composition. Water pollution has various causes, with sewage and industrial waste being the most significant contributors that enter rivers. In many developed countries, the resources allocated for waste management are insufficient [2]. Currently, only around 10% of the generated wastewater is properly treated, while the rest is discharged and deposited into our water bodies. Consequently, impurities infiltrate groundwater, tributaries, and other water sources. The water that ultimately reaches our homes is often highly contaminated, carrying viruses that can harm microorganisms. Another major source of water pollution is agricultural runoff, which involves the drainage of water from fields containing pesticides and fertilizers. Domestic sewage, also known as sanitary sewage, refers to wastewater discharged from households. This type of water contains a wide array of dissolved and suspended contaminants.

Mathematical frameworks play a crucial role in predicting and addressing social issues, and they have been increasingly utilized in recent decades to mitigate the impact of these issues. By applying mathematical models, we can effectively control the spread of real-world problems. As many global challenges exhibit quasi-linear characteristics, relying solely on linear models can often lead to idealistic and unrealistic outcomes. Therefore, non-linear mathematical models provide a more accurate depiction of real-world issues. In the context of modeling lake contamination, a collocation approach has been proposed, which utilizes Bessel polynomials and collocation points to generate an updated matrix problem [3]. This approach allows for a more comprehensive understanding of the contamination dynamics within a lake system. Additionally, the homotopy perturbation approach has been employed to offer approximate and analytical solutions for nonlinear ordinary differential equation systems, such as those encountered in modeling lake contamination [4]. Furthermore, the differential transform method (DTM) has been utilized to analyze pollution models involving interconnected systems of three lakes connected by waterways [5]. These mathematical methods enhance our ability to study and address the complex dynamics of pollution in various environmental systems.

A fractional differential equation system has been investigated by Haq EU [6]. He designed the system to simulate the pollution of a group of lakes and assess the results using other well-known analytical approaches, such as iterative variation and Adomian decomposition. The Laplace Adomian decomposition method is widely applied to solve many real-life problems (see [7]–[9]). Several researchers have used the technique of LADM to investigate different types of problems. For example, in [10], the authors examined the free vibrations of a non-uniform Bernoulli beam using LADM. In [11], a fractional-order sterile insect technology model was studied with the help of LADM. Other studies include the investigation of MHD flow of incompressible fluid between two parallel plates [12], and the study of the effect of a magnetic field on the heat transfer of carbon nanotubes [13]. Different approaches have been used for disease models, studying them with various fractional operator techniques [14]–[23].

Aguirre and Tully [24] employed simple mixture problem methodologies to develop a differential equation that characterizes the pollution concentration in a lake. Prakasha and Veeresha [25] explored three alternative scenarios of the model and established an approximate analytical solution for the system of three fractional differential equations using the q-homotopy analysis transform technique (q-HATM). The results obtained demonstrated the high effectiveness and systematic nature of the proposed strategy. Khalid and Sultana [26] utilized the perturbation-iteration approach to estimate the solution of three input models: periodic, linear step, and exponentially decaying. As a reference, the analytical simulation of the lake system problem was conducted using the fourth-order Runge-Kutta technique (RK4). Biazar and Farrokhi [27] employed compartmental modeling to predict the contamination of a system of lakes through a set of mathematically-researched equations. Bazar and Shahbala [28] utilized the variational iteration method (VIM), and the results demonstrated that the variational iteration approach was easier to implement compared to the Adomian decomposition method for pollution monitoring in lakes. Three distinct types of input models were considered: sinusoidal, impulse, and step. Furthermore, Toufik and Atangana's novel numerical technique was applied [29] to examine the updated model of the contaminated lake system.

Fractional calculus, despite its name, deals with integrals and derivatives of any positive real order and can be considered a branch of mathematical modeling. It focuses on integrodifferential operators and equations with convolution-type integrals that involve weakly singular kernels of power-law form. It is closely related to the theory of pseudo-differential operators. In recent years, fractional calculus has garnered significant attention from researchers, and various aspects of this subject are being explored in research [30], [31]. This is because fractional derivatives serve as an important tool for describing the dynamic behavior of diverse physical systems [32]–[34]. The distinct characteristic of these differential operators lies in their non-local nature, which is absent in integer-order differential operators [35]. In fact, fractional order models are more accurate and practical than classical integer order models, and techniques such as the Legendre operational matrix can be extended to incorporate fractional calculus [36], [37]. The application of fractional differential equations in applied sciences is another area of interest [38], [39].

Fractional order derivatives offer increased flexibility in modeling various biomaterials and systems [40]. They provide powerful tools for characterizing the dynamic behavior of these systems. There are different types of fractional derivatives, such as Riemann and Liouville, and Caputo operators, which are commonly used in practice. Caputo [41] introduced a fractional derivative that allows for the incorporation of conventional initial and boundary conditions relevant to real-world problems. While these fractional derivatives offer improved accuracy in describing real phenomena compared to integer-order derivatives, their kernel functions may result in singularities that lead to computational challenges. To explore fractional mathematical models beyond the traditional Caputo derivative, various methodologies, including iterative and numerical methods, have been employed [42]–[44]. In this study, we focus on a system of three equations that describes lake contamination. The model represents the contamination of a three-lake system [27], as illustrated in Figure 1.

Figure 1. System of three lakes with interconnecting channels.

DownLoad: Full-Size Img PowerPoint

The arrows in the figure indicate the direction of flow in the channels or pipes. In this system, a contaminant is introduced into the first lake through a factory, and the rate of entrance of the pollutant into the lake per unit of time is represented by the function p(t). The function p(t) can be constant or vary over time. The goal is to determine the level of pollution in each lake at any given time. At time t ≥ 0, the contamination level c_i(t) in Lake i is proportional to the volume of water v_i and the concentration of pollutant y_i(t) in the following manner:

$c_{i} (t) = \frac{y_{i} (t)}{v_{i}} .$

In the lake system, it is assumed that each lake is initially free of any contaminants, so the concentration of pollutant in each lake at time t = 0 is given by y_i(0) = 0 for i = 1, 2, 3. To simulate the dynamic behavior of the lake system, a constant g_ji is used to represent the flow rate from lake i to lake j. The flux of pollutant flowing from lake i into lake j at any time t is denoted by r_ji(t) and is defined as follows:

$r_{j i} (t) = g_{j i} c_{i} (t) = \frac{g_{j i}}{v_{i}} y_{i} (t) .$

Thus, r_ji(t) represents the proportion at which the contamination level in lake i flows into lake j at time t. It is important to note this relationship.

Rate of change of pollutant = Input rate - Output rate.

Applying this principle to each lake, we obtain the following system of differential equations:

$(\begin{matrix} \frac{d y_{1}}{d t} = \frac{g_{13}}{v_{3}} y_{3} (t) + p (t) - \frac{g_{31}}{v_{1}} y_{1} (t) - \frac{g_{21}}{v_{1}} y_{1} (t), \\ \frac{d y_{2}}{d t} = \frac{g_{21}}{v_{1}} y_{1} (t) - \frac{g_{32}}{v_{2}} y_{2} (t), \\ \frac{d y_{3}}{d t} = \frac{g_{31}}{v_{1}} y_{1} (t) + \frac{g_{32}}{v_{2}} y_{2} (t) - \frac{g_{13}}{v_{3}} y_{3} (t) . \end{matrix}$ (1.1)

It is apparent from this that the incoming and outgoing flow rates in each lake are constant, and as a result, the water volume in each lake remains fixed. Therefore, we can establish the following conditions:

$(\begin{matrix} L a k e 1 : g_{13} = g_{21} + g_{31}, \\ L a k e 2 : g_{21} = g_{32}, \\ L a k e 3 : g_{31} + g_{32} = g_{13} . \end{matrix}$

To evaluate the performance of the proposed model, we consider specific predicted values as listed in [27], which can be expressed as follows:

$g_{21} = g_{32} = 18 m i^{3} / y e a r, g_{31} = 20 m i^{3} / y e a r, g_{13} = 38 m i^{3} / y e a r, v_{1} = 2900 m i^{3}, v_{2} = 850 m i^{3}, v_{3} = 1180 m i^{3} .$

In this study, to obtain a more accurate qualitative and numerical iterative analysis of the proposed model, we consider the system (1.1) in terms of the Caputo fractional order differential operator. The system can be expressed as follows:

$(\begin{matrix} ^{c} D_{t}^{α} y_{1} (t) = \frac{g_{13}^{α}}{v_{3}^{α}} y_{3} (t) + p (t) - \frac{g_{31}^{α}}{v_{1}^{α}} y_{1} (t) - \frac{g_{21}^{α}}{v_{1}^{α}} y_{1} (t), \\ ^{c} D_{t}^{α} y_{2} (t) = \frac{g_{21}^{α}}{v_{1}^{α}} y_{1} (t) - \frac{g_{32}^{α}}{v_{2}^{α}} y_{2} (t), \\ ^{c} D_{t}^{α} y_{3} (t) = \frac{g_{31}^{α}}{v_{1}^{α}} y_{1} (t) + \frac{g_{32}^{α}}{v_{2}^{α}} y_{2} (t) - \frac{g_{13}^{α}}{v_{3}^{α}} y_{3} (t), \end{matrix}$ (1.2)

where α denotes the fractional order.

The objective of this study is to analyze a system of fractional differential equations that describe lake pollution and provide meaningful insights through a mathematical model that can explain real-world scenarios using a simple and efficient approach. It is important to highlight that the investigated model heavily relies on time and its history, which can be systematically incorporated and represented using the newly developed fractional operator. To obtain the solution of the model, we utilize the Laplace transform in conjunction with the Adomian decomposition method and the Homotopy perturbation method. These techniques enable us to obtain an approximate solution for the considered model. Additionally, we compare the results obtained from both techniques, demonstrating their identical nature. Furthermore, we employ MATLAB to facilitate the numerical solutions and depict the pollution levels in each lake graphically, considering arbitrary fractional orders. This visual representation aids in understanding the behavior of the lake system under different scenarios.

The structure of this paper is organized as follows: Section 2 provides fundamental results and definitions related to the topic. In Section 3, the general solutions of the suggested model are obtained using the Laplace Adomian decomposition technique. This section focuses on the qualitative analysis of the model. Section 4 presents the general solutions of the considered model using the Homotopy perturbation technique. This section discusses the application of the Homotopy perturbation method for obtaining the solutions. In Section 5, a graphical comparison of the numerical solutions obtained from both the Laplace Adomian decomposition technique and the Homotopy perturbation technique is presented. The numerical simulation results are briefly discussed in this section. Finally, in Section 6, the concluding remarks of the study are provided.

2. Preliminaries and basic definitions

In this section, some essential preliminaries and fractional calculus results are presented in line with [45], [46].

Definition 2.1. For α > 0 with $n - 1 < α < n$ , $n \in ℕ$ , the Caputo fractional order derivative of the function f(x) is described as follows:

$^{c} D_{x}^{α} f (x) = \frac{d^{α}}{d t^{α}} f (t) = \frac{1}{Γ (n - α)} \int_{0}^{x} {(x - t)}^{n - α - 1} \frac{d^{n}}{d t^{n}} f (t) d t,$

such that the integral part on the right exists and n = [α] + 1, and the symbol Γ denotes the Gamma function which is described as

$Γ (n) = \int_{0}^{\infty} e^{- t} t^{n - 1} d t .$

Definition 2.2. The Laplace transform of g(t) in the general framework of Caputo definition is described as follows:

$ℒ [^{c} D_{t}^{α} g (t)] = s^{α} G (s) - \sum_{j = 0}^{m - 1} s^{α - j - 1} g^{k} (0), m - 1 < α, m \in ℕ .$

Definition 2.3. The Homotopy perturbation method (HPM) is a semi-analytical technique used to solve linear and nonlinear ordinary and partial differential equations, as well as systems of equations. It is also applicable to systems consisting of both linear and nonlinear differential equations. One notable advantage of using the HPM over decomposition methods is that it can handle nonlinear problems without requiring the use of Adomian polynomials. The HPM was first proposed by the Chinese mathematician He [47], who introduced the concept of forming a homotopy $ν (l, p) : Ω \times [0,1] \to ℝ$ for an equation involving both linear and nonlinear components.

$F (ν, p) = (1 - p) [L (ν) - L (u_{0})] + p [L (ν) + N (ν) - f (l)] = 0,$

where L is used for the linear part, N for the nonlinear part, l ∈ Ω, Ω is a topological space and p ∈ [0,1] is the embedding parameter. Furthermore, u₀ is an initial approximation that satisfies the boundary conditions.

3. General solution of the model (1.2) by (LADM)

In this section, we will discuss the general technique of the Laplace Adomian decomposition method (LADM) to construct the solution for the considered model (1.2), along with the initial conditions. We will also demonstrate how to integrate the LADM with the Caputo differential operator.

$(\begin{matrix} ^{c} D_{t}^{α} y_{1} (t) = \frac{g_{13}}{v_{3}} y_{3} (t) + p (t) - \frac{g_{31}}{v_{1}} y_{1} (t) - \frac{g_{21}}{v_{1}} y_{1} (t), \\ ^{c} D_{t}^{α} y_{2} (t) = \frac{g_{21}}{v_{1}} y_{1} (t) - \frac{g_{32}}{v_{2}} y_{2} (t), \\ ^{c} D_{t}^{α} y_{3} (t) = \frac{g_{31}}{v_{1}} y_{1} (t) + \frac{g_{32}}{v_{2}} y_{2} (t) - \frac{g_{13}}{v_{3}} y_{3} (t), \end{matrix}$ (3.1)

with subject to the initial conditions:

$y_{1} (0) = n_{1}, y_{2} (0) = n_{2}, y_{3} (0) = n_{3} .$

Now, taking the Laplace transform of (3.1) in the Caputo sense, one may get

$(\begin{matrix} ℒ [^{c} D_{t}^{α} y_{1} (t)] = ℒ [\frac{g_{13}}{v_{3}} y_{3} (t) + p (t) - \frac{g_{31}}{v_{1}} y_{1} (t) - \frac{g_{21}}{v_{1}} y_{1} (t)], \\ ℒ [^{c} D_{t}^{α} y_{2} (t)] = ℒ [\frac{g_{21}}{v_{1}} y_{1} (t) - \frac{g_{32}}{v_{2}} y_{2} (t)], \\ ℒ [^{c} D_{t}^{α} y_{3} (t)] = ℒ [\frac{g_{31}}{v_{1}} y_{1} (t) + \frac{g_{32}}{v_{2}} y_{2} (t) - \frac{g_{13}}{v_{3}} y_{3} (t)], \end{matrix}$ (3.2)

using the initial conditions yields

$(\begin{matrix} ℒ [y_{1} (t)] = \frac{n_{1}}{s} + \frac{1}{s^{α}} ℒ [\frac{g_{13}}{v_{3}} y_{3} (t) + p (t) - \frac{g_{31}}{v_{1}} y_{1} (t) - \frac{g_{21}}{v_{1}} y_{1} (t)], \\ ℒ [y_{2} (t)] = \frac{n_{2}}{s} + \frac{1}{s^{α}} ℒ [\frac{g_{21}}{v_{1}} y_{1} (t) - \frac{g_{32}}{v_{2}} y_{2} (t)], \\ ℒ [y_{3} (t)] = \frac{n_{3}}{s} + \frac{1}{s^{α}} ℒ [\frac{g_{31}}{v_{1}} y_{1} (t) + \frac{g_{32}}{v_{2}} y_{2} (t) - \frac{g_{13}}{v_{3}} y_{3} (t)] . \end{matrix}$ (3.3)

Assuming the solution for y₁(t), y₂(t) and y₃(t) in an infinite series is given below:

$\begin{matrix} (\begin{matrix} y_{1} (t) = \sum_{n = 0}^{\infty} y_{1 (n)} (t), \\ y_{2} (t) = \sum_{n = 0}^{\infty} y_{2 (n)} (t), \\ y_{3} (t) = \sum_{n = 0}^{\infty} y_{3 (n)} (t), \end{matrix} \end{matrix}$

using the above series in (3.3), and comparing like terms on both sides, one can get

$(\begin{matrix} ℒ [y_{1 (0)} (t)] = \frac{n_{1}}{s}, ℒ [y_{2 (0)} (t)] = \frac{n_{2}}{s}, ℒ [y_{3 (0)} (t)] = \frac{n_{3}}{s}, \\ ℒ [y_{1 (1)} (t)] = \frac{1}{s^{α}} ℒ [\frac{g_{13}}{v_{3}} y_{3 (0)} (t) + p (t) - \frac{g_{31}}{v_{1}} y_{1 (0)} (t) - \frac{g_{21}}{v_{1}} y_{1 (0)} (t)], \\ ℒ [y_{2 (1)} (t)] = \frac{1}{s^{α}} ℒ [\frac{g_{21}}{v_{1}} y_{1 (0)} (t) - \frac{g_{32}}{v_{2}} y_{2 (0)} (t)], \\ ℒ [y_{3 (1)} (t)] = \frac{1}{s^{α}} ℒ [\frac{g_{31}}{v_{1}} y_{1 (0)} (t) + \frac{g_{32}}{v_{2}} y_{2 (0)} (t) - \frac{g_{13}}{v_{3}} y_{3 (0)} (t)], \\ ⋮ \\ ℒ [y_{1 (n + 1)} (t)] = \frac{1}{s^{α}} ℒ [\frac{g_{13}}{v_{3}} y_{3 (n)} (t) + p (t) - \frac{g_{31}}{v_{1}} y_{1 (n)} (t) - \frac{g_{21}}{v_{1}} y_{1 (n)} (t)], \\ ℒ [y_{2 (n + 1)} (t)] = \frac{1}{s^{α}} ℒ [\frac{g_{21}}{v_{1}} y_{1 (n)} (t) - \frac{g_{32}}{v_{2}} y_{2 (n)} (t)], \\ ℒ [y_{3 (n + 1)} (t)] = \frac{1}{s^{α}} ℒ [\frac{g_{31}}{v_{1}} y_{1 (n)} (t) + \frac{g_{32}}{v_{2}} y_{2 (n)} (t) - \frac{g_{13}}{v_{3}} y_{3 (n)} (t)] . \end{matrix}$ (3.4)

Further, utilizing the inverse Laplace transform to equation (3.4), we have

$(\begin{matrix} y_{1 (0)} (t) = n_{1}, y_{2 (0)} (t) = n_{2}, y_{3 (0)} (t) = n_{3}, \\ y_{1 (1)} (t) = [\frac{g_{13}}{v_{3}} n_{3} + p (t) - \frac{g_{31}}{v_{1}} n_{1} - \frac{g_{21}}{v_{1}} n_{1}] \frac{t^{α}}{Γ (α + 1)}, \\ y_{2 (1)} (t) = [\frac{g_{21}}{v_{1}} n_{1} - \frac{g_{32}}{v_{2}} n_{2}] \frac{t^{α}}{Γ (α + 1)}, \\ y_{3 (1)} (t) = [\frac{g_{31}}{v_{1}} n_{1} + \frac{g_{32}}{v_{2}} n_{2} - \frac{g_{13}}{v_{3}} n_{3}] \frac{t^{α}}{Γ (α + 1)}, \\ y_{1 (2)} (t) = [\frac{g_{13}}{v_{3}} x_{11} + p (t) - \frac{g_{31}}{v_{1}} w_{11} - \frac{g_{21}}{v_{1}} w_{11}] \frac{t^{2 α}}{Γ (2 α + 1)}, \\ y_{2 (2)} (t) = [\frac{g_{21}}{v_{1}} w_{11} - \frac{g_{32}}{v_{2}} u_{11}] \frac{t^{2 α}}{Γ (2 α + 1)}, \\ y_{3 (2)} (t) = [\frac{g_{31}}{v_{1}} w_{11} + \frac{g_{32}}{v_{2}} u_{11} - \frac{g_{13}}{v_{3}} x_{11}] \frac{t^{2 α}}{Γ (2 α + 1)}, \\ y_{1 (3)} (t) = [\frac{g_{13}}{v_{3}} x_{111} + p (t) - \frac{g_{31}}{v_{1}} w_{111} - \frac{g_{21}}{v_{1}} w_{111}] \frac{t^{3 α}}{Γ (3 α + 1)}, \\ y_{2 (3)} (t) = [\frac{g_{21}}{v_{1}} w_{111} - \frac{g_{32}}{v_{2}} u_{111}] \frac{t^{3 α}}{Γ (3 α + 1)}, \\ y_{3 (3)} (t) = [\frac{g_{31}}{v_{1}} w_{111} + \frac{g_{32}}{v_{2}} u_{111} - \frac{g_{13}}{v_{3}} x_{111}] \frac{t^{3 α}}{Γ (3 α + 1)} . \end{matrix}$

Additionally, the remaining terms can be derived in a similar fashion. The unknown values in the aforementioned equations are listed below:

$(\begin{matrix} w_{11} = \frac{g_{13}}{v_{3}} n_{3} + p (t) - \frac{g_{31}}{v_{1}} n_{1} - \frac{g_{21}}{v_{1}} n_{1}, \\ u_{11} = \frac{g_{21}}{v_{1}} n_{1} - \frac{g_{32}}{v_{2}} n_{2}, \\ x_{11} = \frac{g_{31}}{v_{1}} n_{1} + \frac{g_{32}}{v_{2}} n_{2} - \frac{g_{13}}{v_{3}} n_{3}, \\ w_{111} = \frac{g_{13}}{v_{3}} x_{11} + p (t) - \frac{g_{31}}{v_{1}} w_{11} - \frac{g_{21}}{v_{1}} w_{11}, \\ u_{111} = \frac{g_{21}}{v_{1}} w_{11} - \frac{g_{32}}{v_{2}} u_{11}, \\ x_{111} = \frac{g_{31}}{v_{1}} w_{11} + \frac{g_{32}}{v_{2}} u_{11} - \frac{g_{13}}{v_{3}} x_{11} . \end{matrix}$

4. General solution for model (1.2) with (HPM)

We will now apply the Homotopy Perturbation Method (HPM) to derive the general solution of Model (1.2) as:

$(\begin{matrix} (1 - q {) [}^{c} D_{t}^{α} (y_{1} (t)) -^{c} D_{t}^{α} (y_{1 (0)} (t))] + q [^{c} D_{t}^{α} (y_{1} (t)) - \frac{g_{13}}{v_{3}} y_{3} (t) - p (t) + \frac{g_{31}}{v_{1}} y_{1} (t) + \frac{g_{21}}{v_{1}} y_{1} (t)] = 0, \\ (1 - q {) [}^{c} D_{t}^{α} (y_{2} (t)) -^{c} D_{t}^{α} (y_{2 (0)} (t))] + q [^{c} D_{t}^{α} (y_{2} (t)) - \frac{g_{21}}{v_{1}} y_{1} (t) + \frac{g_{32}}{v_{2}} y_{2} (t)] = 0, \\ (1 - q {) [}^{c} D_{t}^{α} (y_{3} (t)) -^{c} D_{t}^{α} (y_{3 (0)} (t))] + q [^{c} D_{t}^{α} (y_{3} (t)) - \frac{g_{31}}{v_{1}} y_{1} (t) - \frac{g_{32}}{v_{2}} y_{2} (t) + \frac{g_{13}}{v_{3}} y_{3} (t)] = 0. \end{matrix}$ (4.1)

By substituting q = 0 into equation (4.1), we obtain the following system of fractional differential equations:

$(\begin{matrix} ^{c} D_{t}^{α} (y_{1} (t)) -^{c} D_{t}^{α} (y_{1 (0)} (t)) = 0, \\ ^{c} D_{t}^{α} (y_{2} (t)) -^{c} D_{t}^{α} (y_{2 (0)} (t)) = 0, \\ ^{c} D_{t}^{α} (y_{3} (t)) -^{c} D_{t}^{α} (y_{3 (0)} (t)) = 0. \end{matrix}$ (4.2)

The solution to the above equation is straightforward. Next, setting q = 1 in equation (4.1) results in a similar model to equation (1.2). We assume the solution takes the form of an infinite series as:

$(\begin{matrix} y_{1} (t) = \sum_{n = 0}^{\infty} q^{n} y_{1 (n)} (t), \\ y_{2} (t) = \sum_{n = 0}^{\infty} q^{n} y_{2 (n)} (t), \\ y_{3} (t) = \sum_{n = 0}^{\infty} q^{n} y_{3 (n)} (t) . \end{matrix}$ (4.3)

Furthermore, the original system can be obtained by substituting q = 1 into equation (4.1). By substituting equation (4.3) into equation (4.1) and comparing the terms with respect to the powers of q, we obtain:

$q^{0} : (\begin{matrix} y_{1 (0)} (t) = y_{1} (0) = n_{1}, \\ y_{2 (0)} (t) = y_{2} (0) = n_{2}, \\ y_{3 (0)} (t) = y_{3} (0) = n_{3} . \end{matrix}$ (4.4)

Similarly,

$q^{1} : (\begin{matrix} y_{1 (1)} (t) = [\frac{g_{13}}{v_{3}} n_{3} + p (t) - \frac{g_{31}}{v_{1}} n_{1} - \frac{g_{21}}{v_{1}} n_{1}] \frac{t^{α}}{Γ (α + 1)}, \\ y_{2 (1)} (t) = [\frac{g_{21}}{v_{1}} n_{1} - \frac{g_{32}}{v_{2}} n_{2}] \frac{t^{α}}{Γ (α + 1)}, \\ y_{3 (1)} (t) = [\frac{g_{31}}{v_{1}} n_{1} + \frac{g_{32}}{v_{2}} n_{2} - \frac{g_{13}}{v_{3}} n_{3}] \frac{t^{α}}{Γ (α + 1)} . \end{matrix}$ (4.5)

$q^{2} : (\begin{matrix} y_{1 (2)} (t) = [\frac{g_{13}}{v_{3}} x_{11} + p (t) - \frac{g_{31}}{v_{1}} w_{11} - \frac{g_{21}}{v_{1}} w_{11}] \frac{t^{2 α}}{Γ (2 α + 1)}, \\ y_{2 (2)} (t) = [\frac{g_{21}}{v_{1}} w_{11} - \frac{g_{32}}{v_{2}} u_{11}] \frac{t^{2 α}}{Γ (2 α + 1)}, \\ y_{3 (2)} (t) = [\frac{g_{31}}{v_{1}} w_{11} + \frac{g_{32}}{v_{2}} u_{11} - \frac{g_{13}}{v_{3}} x_{11}] \frac{t^{2 α}}{Γ (2 α + 1)} . \end{matrix}$ (4.6)

The approximate series solution is thus acquired. In the next section, we will do simulation for the aforementioned methods.

5. Numerical results and simulations

In this section, we conduct numerical simulations to complement the analytical findings of our proposed model. The simulations involve qualitative point analysis and consider the parameters from a biological feasibility perspective. By using the parametric values, we determine the following terms of the proposed model:

$(\begin{matrix} y_{1 (0)} (t) = 0, y_{2 (0)} (t) = 0, y_{3 (0)} (t) = 0, \\ y_{1 (1)} (t) = p (t) \frac{t^{α}}{Γ (α + 1)}, y_{2 (1)} (t) = 0, y_{3 (1)} (t) = 0, \\ y_{1 (2)} (t) = (0.999) p (t) \frac{t^{2 α}}{Γ (2 α + 1)}, y_{2 (2)} (t) = (0.006206) p (t) \frac{t^{2 α}}{Γ (2 α + 1)}, \\ y_{3 (2)} (t) = (0.00689) p (t) \frac{t^{2 α}}{Γ (2 α + 1)}, \\ y_{1 (3)} (t) = (0.98713) p (t) \frac{t^{3 α}}{Γ (3 α + 1)}, \\ y_{2 (3)} (t) = (0.006069) p (t) \frac{t^{3 α}}{Γ (3 α + 1)}, \\ y_{3 (3)} (t) = (0.00679) p (t) \frac{t^{3 α}}{Γ (3 α + 1)} . \end{matrix}$ (5.1)

Furthermore, the solutions to the first few terms are given as:

$(\begin{matrix} y_{1} (t) = p (t) \frac{t^{α}}{Γ (α + 1)} + (0.999) p (t) \frac{t^{2 α}}{Γ (2 α + 1)} + (0.98713) p (t) \frac{t^{3 α}}{Γ (3 α + 1)}, \\ y_{2} (t) = (0.006206) p (t) \frac{t^{2 α}}{Γ (2 α + 1)} + (0.006069) p (t) \frac{t^{3 α}}{Γ (3 α + 1)}, \\ y_{3} (t) = (0.00689) p (t) \frac{t^{2 α}}{Γ (2 α + 1)} + (0.00679) p (t) \frac{t^{3 α}}{Γ (3 α + 1)} . \end{matrix}$ (5.2)

Now, evaluating (5.2) for α = 0.97, one get

$(\begin{matrix} y_{1} (t) = 1.012468567 p (t) t^{0.97} + 0.5275576582 p (t) t^{1.94} + 0.1839995904 p (t) t^{2.91}, \\ y_{2} (t) = 0.006283379927 p (t) t^{1.94} + 0.003204952380 p (t) t^{2.91}, \\ y_{3} (t) = 0.006975908427 p (t) t^{1.94} + 0.003585702201 p (t) t^{2.91} . \end{matrix}$ (5.3)

Similarly, for α = 0.98, the approximations are:

$(\begin{matrix} y_{1} (t) = 1.008360917 p (t) t^{0.98} + 0.5181171462 p (t) t^{1.96} + 0.1773092790 p (t) t^{2.94}, \\ y_{2} (t) = 0.006257887851 p (t) t^{1.96} + 0.003147600561 p (t) t^{2.94}, \\ y_{3} (t) = 0.006947606718 p (t) t^{1.96} + 0.003521536960 p (t) t^{2.94} . \end{matrix}$ (5.4)

And for α = 0.99, one can obtain the approximate solutions are:

$(\begin{matrix} y_{1} (t) = 1.004204343 p (t) t^{0.99} + 0.5087639171 p (t) t^{1.98} + 0.1708178394 p (t) t^{2.97}, \\ y_{2} (t) = 0.006232092153 p (t) t^{1.98} + 0.003090778992 p (t) t^{2.97}, \\ y_{3} (t) = 0.006918967923 p (t) t^{1.98} + 0.003457964962 p (t) t^{2.97} . \end{matrix}$ (5.5)

Furthermore, we have considered three different types of input functions to represent the duration of the contaminant in each lake. These input functions include periodic, exponentially decaying, and linear inputs.

Case 1: Periodic Input Model: In this scenario, the model is evaluated when pollution is periodically introduced into Lake 1. We have chosen the input function $p (t) = a s i n (ω t) + c$ , where a and ω represent the amplitude and frequency of the fluctuations, respectively, and c represents the average pollutant concentration input. For this example, let us consider a = ω = c = 1 as an illustration. This could represent a manufacturing facility that discharges waste during the day and generates more pollutants during the night due to increased production during that time. As a result, the pollutant input exhibits a periodic pattern. By applying this periodic input to the model, the concentration levels in the lakes eventually converge to the average input level of the contaminant.

Figure 2. Adaptive nature of the approximated LADM solutions for y₁, y₂, and y₃ at different arbitrary fractional orders for the Case 1.

DownLoad: Full-Size Img PowerPoint

Figure 3. Dynamical behavior of the considered model with LADM solutions for y₁, y₂, and y₃ at fractional orders 0.80, 0.50, 0.20 for the Case 1.

DownLoad: Full-Size Img PowerPoint

Based on Figure 2, it is evident that the pollution concentration in the lakes increases over time. Moreover, for smaller fractional orders, the increase is more rapid. It is worth noting that the timescale in all the figures is measured in days. In the next step, we will provide a comparison of the lakes using both the Laplace Adomian decomposition method (LADM) and the Homotopy perturbation method (HPM). Specifically, we will focus on the first few terms of the simulations, which will demonstrate the similarities between the two methods.

Figure 4. Adaptive nature of the approximated LADM and HPM solutions for y₁, y₂, and y₃ at different arbitrary fractional orders for the Case 1.

DownLoad: Full-Size Img PowerPoint

Figure 5. Dynamical behavior of the considered model with LADM and HPM solutions for y₁, y₂, and y₃ at fractional orders 0.80, 0.50, 0.20 for the Case 1.

DownLoad: Full-Size Img PowerPoint

Case 2: Exponentially Decaying Input Model: In this situation, the model is examined when pollutants with huge dumping. We chose $p (t) = a e^{- b t}$ , where a = 200 and b = 10.

Figure 6. Adaptive nature of the approximated LADM solutions for y₁, y₂, and y₃ at different arbitrary fractional orders for the Case 2.

DownLoad: Full-Size Img PowerPoint

Figure 7. Dynamical behavior of the considered model with LADM solutions for y₁, y₂, and y₃ at fractional orders 0.80, 0.50, 0.20 for the Case 2.

DownLoad: Full-Size Img PowerPoint

Figure 6 illustrates that the pollution concentration in the lakes initially increases, then decreases, and eventually stabilizes at a certain level over time. Similarly, we will proceed with the comparison of the lakes using both the Laplace Adomian decomposition method (LADM) and the Homotopy perturbation method (HPM). By examining the first few terms of the simulations, we will demonstrate the similarities between the two methods, thus validating their simulation results.

Figure 8. Adaptive nature of the approximated LADM and HPM solutions for y₁, y₂, and y₃ at different arbitrary fractional orders for the Case 2.

DownLoad: Full-Size Img PowerPoint

Figure 9. Adaptive nature of the approximated LADM and HPM solutions for y₁, y₂, and y₃ at different arbitrary fractional orders for the Case 3.

DownLoad: Full-Size Img PowerPoint

Figure 10. Dynamical behavior of the considered model with LADM and HPM solutions for y₁, y₂, and y₃ at fractional orders 0.80, 0.50, 0.20 for the Case 3.

DownLoad: Full-Size Img PowerPoint

Case 3: Linear Input Model: In this scenario, the model considers the case where Lake 1 is initially contaminated with a pollutant with a linear proportion. We have chosen the input function p(t) = a, t, where a = 100. For this linear input, the pollutant starts flowing into the lake at time zero, and the amount of pollutant before time zero is assumed to be zero. In the case of a step input, the key characteristic is that the input abruptly increases at time zero and remains relatively constant thereafter. As an example, consider a manufacturing plant that starts its operations at time zero and immediately begins discharging untreated sewage at a consistent rate and intensity. This linear input represents the continuous influx of pollutants into the lake.

Figure 9 shows the concentration of pollution in lakes rising rapidly with passing time. Next, the comparison of these lakes using both proposed methods of LADM and HPM for the first few terms that demonstrate simulation similarities are presented.

Figure 11. Adaptive nature of the approximated LADM and HPM solutions for y₁, y₂, and y₃ at different arbitrary fractional orders for the Case 3.

DownLoad: Full-Size Img PowerPoint

6. Conclusion

In this research, we have analyzed the dynamic behavior of a lake pollution model using the Caputo differential operator and the utilities of fractional calculus. The model has been numerically explored using the Laplace transform with the Adomian decomposition method (LADM) and the Homotopy perturbation method (HPM). The numerical results obtained from both methods are highly similar and provide strong confirmation for the considered model in arbitrary order derivatives. The outcomes are influenced by various parameters used in the model, and both methods converge effectively for solving fractional-order differential equations. The graphical results, generated using MATLAB, demonstrate the dependency of the model on the fractional operator and the parameters utilized in the proposed methods. The present study highlights the significance of the fractional concept in understanding and analyzing the suggested lakes pollution model, which is highly dependent on time and its history. It also provides a foundation for future research in this field. For instance, the mathematical model can be further enhanced by considering various dynamic structures and investigating different types of derivatives. The advantages of the LADM and HPM methods, including simplicity, accuracy, flexibility, and efficiency, make them valuable tools for analyzing nonlinear problems in various scientific and engineering disciplines. These methods offer promising avenues for further exploration and application in fractional calculus and related fields.

References

[1]	X. Zou, K. Wang, Optimal harvesting for a stochastic regime-switching logistic diffusion system with jumps, Nonlinear Anal. Hybrid Syst., 13 (2014), 32–44. https://doi.org/10.1016/j.nahs.2014.01.001 doi: 10.1016/j.nahs.2014.01.001
[2]	J. Roy, D. Barman, S. Alam, Role of fear in a predator-prey system with ratio-dependent functional response in deterministic and stochastic environment, Biosystems, 197 (2020), 104176. https://doi.org/10.1016/j.biosystems.2020.104176 doi: 10.1016/j.biosystems.2020.104176
[3]	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
[4]	Q. Yang, X. Zhang, D. Jiang, Dynamical behaviors of a stochastic food chain system with Ornstein-Uhlenbeck process, J. Nonlinear Sci., 32 (2022), 1–40. https://doi.org/10.1007/s00332-021-09760-y doi: 10.1007/s00332-021-09760-y
[5]	L. Wang, D. Jiang, Ergodicity and threshold behaviors of a predator-prey model in stochastic chemostat driven by regime switching, Math. Meth. Appl. Sci., 44 (2021), 325–344. https://doi.org/10.1002/mma.6738 doi: 10.1002/mma.6738
[6]	Q. Luo, X. Mao, Stochastic population dynamics under regime switching, J. Math. Anal. Appl., 334 (2007), 69–84. https://doi.org/10.1016/j.jmaa.2006.12.032 doi: 10.1016/j.jmaa.2006.12.032
[7]	Q. Luo, X. Mao, Stochastic population dynamics under regime switching Ⅱ, J. Math. Anal. Appl., 355 (2009), 577–593. https://doi.org/10.1016/j.jmaa.2009.02.010 doi: 10.1016/j.jmaa.2009.02.010
[8]	C. Zhu, G. Yin, On hybrid competitive Lotka-Volterra ecosystems, Nonlinear Anal., 71 (2009), 1370–1379. https://doi.org/10.1016/j.na.2009.01.166 doi: 10.1016/j.na.2009.01.166
[9]	X. Li, A. Gray, D. Jiang, X. Mao, Sufficient and necessary conditions of stochastic permanence and extinction for stochastic logistic populations under regime switching, J. Math. Anal. Appl., 376 (2011), 11–28. https://doi.org/10.1016/j.jmaa.2010.10.053 doi: 10.1016/j.jmaa.2010.10.053
[10]	M. Ouyang, X. Li, Permanence and asymptotical behavior of stochastic prey-predator system with Markovian switching, Appl. Math. Comput., 266 (2015), 539–559. https://doi.org/10.1016/j.amc.2015.05.083 doi: 10.1016/j.amc.2015.05.083
[11]	J. Bao, J. Shao, Permanence and extinction of regime-switching predator-prey models, SIAM J. Math. Anal., 48 (2016), 725–739. https://doi.org/10.1137/15M1024512 doi: 10.1137/15M1024512
[12]	M. Liu, X. He, J. Yu, Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays, Nonlinear Anal. Hybrid Syst., 28 (2018), 87–104. https://doi.org/10.1016/j.nahs.2017.10.004 doi: 10.1016/j.nahs.2017.10.004
[13]	Y. Cai, S. Cai, X. Mao, Stochastic delay foraging arena predator-prey system with Markov switching, Stoch. Anal. Appl., 38 (2020), 191–212. https://doi.org/10.1080/07362994.2019.1679645 doi: 10.1080/07362994.2019.1679645
[14]	J. Bao, X. Mao, G. Yin, C. Yuan, Competitive Lotka-Volterra population dynamics with jumps, Nonlinear Anal., 74 (2011), 6601–6616. https://doi.org/10.1016/j.na.2011.06.043 doi: 10.1016/j.na.2011.06.043
[15]	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
[16]	M. Liu, K. Wang, Dynamics of a Leslie-Gower Holling-type Ⅱ predator-prey system with Lévy jumps, Nonlinear Anal., 85 (2013), 204–213. https://doi.org/10.1016/j.na.2013.02.018 doi: 10.1016/j.na.2013.02.018
[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]	M. Liu, M. Deng, B. Du, Analysis of a stochastic logistic model with diffusion, Appl. Math. Comput., 266 (2015), 169–182. https://doi.org/10.1016/j.amc.2015.05.050 doi: 10.1016/j.amc.2015.05.050
[19]	X. Zhang, W. Li, M. Liu, K. Wang, Dynamics of a stochastic Holling Ⅱ one-predator two-prey system with jumps, Phys. A, 421 (2015), 571–582. https://doi.org/10.1016/j.physa.2014.11.060 doi: 10.1016/j.physa.2014.11.060
[20]	D. Valenti, G. Denaro, A. Cognata, B. La Spagnolo, A. Bonanno, G. Basilone, et al., Picophytoplankton dynamics in noisy marine environment, Acta Phys. Pol. B, 43 (2012), 1227–1240. https://doi.org/10.5506/APhysPolB.43.1227 doi: 10.5506/APhysPolB.43.1227
[21]	C. Guarcello, D. Valenti, G. Augello, B. Spagnolo, The role of non-Gaussian sources in the transient dynamics of long Josephson junctions, Acta Phys. Pol. B, 44 (2013), 997–1005. https://doi.org/10.5506/APhysPolB.44.997 doi: 10.5506/APhysPolB.44.997
[22]	C. Guarcello, D. Valenti, B. Spagnolo, V. Pierro, G. Filatrella, Josephson-based threshold detector for Lévy-distributed current fluctuations, Phys. Rev. Appl., 11 (2019), 044078. https://doi.org/10.1103/PhysRevApplied.11.044078 doi: 10.1103/PhysRevApplied.11.044078
[23]	A. A. Dubkov, A. La Cognata, B. Spagnolo, The problem of analytical calculation of barrier crossing characteristics for Lévy flights, J. Stat. Mech. Theory Exp., 2009 (2019), P01002. https://doi.org/10.1088/1742-5468/2009/01/P01002 doi: 10.1088/1742-5468/2009/01/P01002
[24]	B. Lisowski, D. Valenti, B. Spagnolo, M. Bier, E. Gudowska-Nowak, Stepping molecular motor amid Lévy white noise, Phys. Rev. E, 91 (2015), 042713. https://doi.org/10.1103/PhysRevE.91.042713 doi: 10.1103/PhysRevE.91.042713
[25]	I. A. Surazhevsky, V. A. Demin, A. I. Ilyasov, A. V. Emelyanov, K. E. Nikiruy, V. V. Rylkov, et al., Noise-assisted persistence and recovery of memory state in a memristive spiking neuromorphic network, Chaos Solitons Fractals, 146 (2021), 110890. https://doi.org/10.1016/j.chaos.2021.110890 doi: 10.1016/j.chaos.2021.110890
[26]	A. N. Mikhaylov, D. V. Guseinov, A. I. Belov, D. S. Korolev, V. A. Shishmakova, M. N. Koryazhkina, et al., Stochastic resonance in a metal-oxide memristive device, Chaos Solitons Fractal, 144 (2021), 110723. https://doi.org/10.1016/j.chaos.2021.110723 doi: 10.1016/j.chaos.2021.110723
[27]	Y. V. Ushakov, A. A. Dubkov, B. Spagnolo, Spike train statistics for consonant and dissonant musical accords in a simple auditory sensory model, Phys. Rev. E, 81 (2010), 041911. https://doi.org/10.1103/PhysRevE.81.041911 doi: 10.1103/PhysRevE.81.041911
[28]	N. V. Agudov, A. V. Safonov, A. V. Krichigin, A. A. Kharcheva, A. A. Dubkov, D. Valenti, et al., Nonstationary distributions and relaxation times in a stochastic model of memristor, J. Stat. Mech. Theory Exp., 2020 (2020), 024003. https://doi.org/10.1088/1742-5468/ab684a doi: 10.1088/1742-5468/ab684a
[29]	D. O. Filatov, D. V. Vrzheshch, O. V. Tabakov, A. S. Novikov, A. I. Belov, I. N. Antonov, et al., Noise-induced resistive switching in a memristor based on $\mathrm{ZrO_{2}(Y)/Ta_{2}O_{5}}$ stack, J. Stat. Mech. Theory Exp., 2019 (2019), 124026. https://doi.org/10.1088/1742-5468/ab5704 doi: 10.1088/1742-5468/ab5704
[30]	A. Carollo, B. Spagnolo, A. A. Dubkov, D. Valenti, On quantumness in multi-parameter quantum estimation, J. Stat. Mech. Theory Exp., 2019 (2019), 094010. https://doi.org/10.1088/1742-5468/ab3ccb doi: 10.1088/1742-5468/ab3ccb
[31]	R. Stassi, S. Savasta, L. Garziano, B. Spagnolo, F. Nori, Output field-quadrature measurements and squeezing in ultrastrong cavity-QED, New J. Phys., 18 (2016), 123005. https://doi.org/10.1088/1367-2630/18/12/123005 doi: 10.1088/1367-2630/18/12/123005
[32]	S. Ciuchi, F. De Pasquale, B. Spagnolo, Nonlinear relaxation in the presence of an absorbing barrier, Phys. Rev. E, 47 (1993), 3915. https://doi.org/10.1103/PhysRevE.47.3915 doi: 10.1103/PhysRevE.47.3915
[33]	Y. Kuang, Delay Differential Equations: With Applications in Population Dynamics, Academic Press, Boston, 1993.
[34]	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
[35]	F. A. Rihan, H. J. Alsakaji, Stochastic delay differential equations of three-species prey-predator system with cooperation among prey species, Discret. Contin. Dyn. Syst. Ser. S, 15 (2020), 245. https://doi.org/10.3934/dcdss.2020468 doi: 10.3934/dcdss.2020468
[36]	H. J. Alsakaji, S. Kundu, F. A. Rihan, Delay differential model of one-predator two-prey system with Monod-Haldane and holling type Ⅱ functional responses, Appl. Math. Comput., 397 (2021), 125919. https://doi.org/10.1016/j.amc.2020.125919 doi: 10.1016/j.amc.2020.125919
[37]	L. Wang, R. Zhang, Y. Wang, Global exponential stability of reaction-diffusion cellular neural networks with S-type distributed time delays, Nonlinear Anal., 10 (2009), 1101–1113. https://doi.org/10.1016/j.nonrwa.2007.12.002 doi: 10.1016/j.nonrwa.2007.12.002
[38]	L. Wang, D. Xu, Global asymptotic stability of bidirectional associative memory neural networks with S-type distributed delays, Int. J. Syst. Sci., 338 (2002), 869–877. https://doi.org/10.1080/00207720210161777 doi: 10.1080/00207720210161777
[39]	S. Abbas, D. Bahuguna, M. Banerjee, Effect of stochastic perturbation on a two species competitive model, Nonlinear Anal. Hybrid Syst., 3 (2009), 195–206. https://doi.org/10.1016/j.nahs.2009.01.001 doi: 10.1016/j.nahs.2009.01.001
[40]	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
[41]	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
[42]	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
[43]	M. Liu, X. He, J. Yu, Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays, Nonlinear Anal. Hybrid Syst., 28 (2018), 87–104. https://doi.org/10.1016/j.nahs.2017.10.004 doi: 10.1016/j.nahs.2017.10.004
[44]	M. Liu, C. Bai, Dynamics of a stochastic one-prey two-predator model with Lévy jumps, Appl. Math. Comput., 284 (2016), 308–321. https://doi.org/10.1016/j.amc.2016.02.033 doi: 10.1016/j.amc.2016.02.033
[45]	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
[46]	Q. Liu, D. Jiang, N. Shi, T. Hayat, A. Alsaedi, Stochastic mutualism model with Lévy jumps, Commun. Nonlinear Sci. Numer. Simul., 43 (2017), 78–90. https://doi.org/10.1016/j.cnsns.2016.05.003 doi: 10.1016/j.cnsns.2016.05.003
[47]	H. Qiu, W. Deng, Optimal harvesting of a stochastic delay competitive Lotka-Volterra model with Lévy jumps, Appl. Math. Comput., 317 (2018), 210–222. https://doi.org/10.1016/j.amc.2017.08.044 doi: 10.1016/j.amc.2017.08.044
[48]	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
[49]	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
[50]	S. Wang, L. Wang, T. Wei, Optimal harvesting for a stochastic logistic model with S-type distributed time delay, J. Differ. Equation Appl., 23 (2017), 618–632. https://doi.org/10.1080/10236198.2016.1269761 doi: 10.1080/10236198.2016.1269761
[51]	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
[52]	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
[53]	Q. Liu, Q. Chen, Z. Liu, Analysis on stochastic delay Lotka-Volterra systems driven by Lévy noise, Appl. Math. Comput., 235 (2014), 261–271. https://doi.org/10.1016/j.amc.2014.03.011 doi: 10.1016/j.amc.2014.03.011
[54]	X. Mao, Stochastic Differential Equations and Applications, Horwood Publishing Limited, 2007. https://doi.org/10.1533/9780857099402
[55]	D. Applebaum, Lévy Processes and Stochastic Calculus, Cambridge University Press, 2009. https://doi.org/10.1017/CBO9780511809781
[56]	I. Barbalat, Systems dequations differentielles d'osci d'oscillations, Rev. Roumaine Math. Pures Appl., 4 (1959), 267–270.
[57]	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
[58]	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
[59]	G. Prato, J. Zabczyk, Ergodicity for Infinite Dimensional Systems, Cambridge University Press, 1996.
[60]	M. Liu, Optimal harvesting policy of a stochastic predator-prey model with time delay, Appl. Math. Lett., 48 (2015), 102–108. https://doi.org/10.1016/j.aml.2014.10.007 doi: 10.1016/j.aml.2014.10.007

This article has been cited by:

1.	Hardik Joshi, Mehmet Yavuz, Numerical Analysis of Compound Biochemical Calcium Oscillations Process in Hepatocyte Cells, 2024, 8, 2701-0198, 10.1002/adbi.202300647
2.	Yasir Nadeem Anjam, Rubayyi Turki Alqahtani, Nadiyah Hussain Alharthi, Saira Tabassum, Unveiling the Complexity of HIV Transmission: Integrating Multi-Level Infections via Fractal-Fractional Analysis, 2024, 8, 2504-3110, 299, 10.3390/fractalfract8050299
3.	Yasir Nadeem Anjam, Asma Arshad, Rubayyi T. Alqahtani, Muhammad Arshad, Unveiling the dynamics of drug transmission: A fractal-fractional approach integrating criminal law perspectives, 2024, 9, 2473-6988, 13102, 10.3934/math.2024640
4.	Mohamed EL-GAMEL, Nesreen MOHAMED, Waleed ADEL, Genocchi collocation method for accurate solution of nonlinear fractional differential equations with error analysis, 2023, 3, 2791-8564, 351, 10.53391/mmnsa.1373647
5.	Twinkle R. Singh, Approximate‐analytical iterative approach to time‐fractional Bloch equation with Mittag–Leffler type kernel, 2024, 47, 0170-4214, 7028, 10.1002/mma.9955
6.	M. L. Rupa, K. Aruna, K. Raghavendar, Insights into the time Fractional Belousov-Zhabotinsky System Arises in Thermodynamics, 2024, 63, 1572-9575, 10.1007/s10773-024-05770-0
7.	Pushpendra Kumar, Vedat Suat Erturk, A variable-order fractional mathematical model for the strategy to combat the atmospheric level of carbon dioxide, 2024, 10, 2363-6203, 3529, 10.1007/s40808-024-01962-z
8.	Yasir Nadeem Anjam, Saira Tabassum, Muhammad Arshad, Mati ur Rahman, Muhammad Ateeq Tahir, Mathematical insights of social media addiction: fractal-fractional perspectives, 2024, 99, 0031-8949, 055230, 10.1088/1402-4896/ad348c
9.	Asiyeh Ebrahimzadeh, Amin Jajarmi, Dumitru Baleanu, Enhancing Water Pollution Management Through a Comprehensive Fractional Modeling Framework and Optimal Control Techniques, 2024, 31, 1776-0852, 10.1007/s44198-024-00215-y
10.	Ibtehal Alazman, Manvendra Narayan Mishra, Badr Saad Alkahtani, Pranay Goswami, Computational analysis of rabies and its solution by applying fractional operator, 2024, 32, 2769-0911, 10.1080/27690911.2024.2340607
11.	Kottakkaran Sooppy Nisar, M. Sivashankar, S. Sabarinathan, C. Ravichandran, V. Sivaramakrishnan, Evaluating the stability and efficacy of fractal-fractional models in reproductive cancer apoptosis with ABT-737, 2025, 2025, 1029-242X, 10.1186/s13660-024-03249-4
12.	Pasquini Fotsing Soh, Mathew Kinyanjui, David Malonza, Roy Kiogora, Non-local operator as a mathematical tool to improve the modeling of water pollution phenomena in environmental science: A spatio-temporal approach, 2025, 26668181, 101161, 10.1016/j.padiff.2025.101161
13.	Gnanasakthivel Hamsavarthini, Chokkalingam Ravichandran, Fahad Aljuaydi, Kottakkaran Sooppy Nisar, A fractional-order framework for investigating groundwater depletion under environmental and human pressures, 2025, 1598-5865, 10.1007/s12190-025-02514-z

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

Mathematical Biosciences and Engineering

4.4

Metrics

Article views(1938) PDF downloads(120) Cited by(1)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Mathematical Biosciences and Engineering

Optimal harvesting strategy for stochastic hybrid delay Lotka-Volterra systems with Lévy noise in a polluted environment

Related Papers:

Abstract

1. Introduction

2. Preliminaries and basic definitions

3. General solution of the model (1.2) by (LADM)

4. General solution for model (1.2) with (HPM)

5. Numerical results and simulations

6. Conclusion

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Abstract

1. Introduction

2. Preliminaries and basic definitions

3. General solution of the model (1.2) by (LADM)

4. General solution for model (1.2) with (HPM)

5. Numerical results and simulations

6. Conclusion

References

Mathematical Biosciences and Engineering

Optimal harvesting strategy for stochastic hybrid delay Lotka-Volterra systems with Lévy noise in a polluted environment

Related Papers:

Abstract

1. Introduction

2. Preliminaries and basic definitions

3. General solution of the model (1.2) by (LADM)

4. General solution for model (1.2) with (HPM)

5. Numerical results and simulations

6. Conclusion

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

Abstract

1. Introduction

2. Preliminaries and basic definitions

3. General solution of the model (1.2) by (LADM)

4. General solution for model (1.2) with (HPM)

5. Numerical results and simulations

6. Conclusion

References