In this paper, we primarily focused on the existence and uniqueness of the initial value problem for fractional order fuzzy ordinary differential equations in a fuzzy metric space. First, definitions and relevant properties of the Gamma function and Beta function within a fuzzy metric space were provided. Second, by employing the principle of fuzzy compression mapping and Choquet integral of fuzzy numerical functions, we established the existence and uniqueness of solutions to initial value problems for fuzzy ordinary differential equations. Finally, several examples were presented to demonstrate the validity of our obtained results.
Citation: Li Chen, Suyun Wang, Yongjun Li, Jinying Wei. New results for fractional ordinary differential equations in fuzzy metric space[J]. AIMS Mathematics, 2024, 9(6): 13861-13873. doi: 10.3934/math.2024674
Related Papers:
[1]
Kimun Ryu, Wonlyul Ko .
Stability and bifurcations in a delayed predator-prey system with prey-taxis and hunting cooperation functional response. AIMS Mathematics, 2025, 10(6): 12808-12840.
doi: 10.3934/math.2025576
[2]
Teekam Singh, Ramu Dubey, Vishnu Narayan Mishra .
Spatial dynamics of predator-prey system with hunting cooperation in predators and type I functional response. AIMS Mathematics, 2020, 5(1): 673-684.
doi: 10.3934/math.2020045
[3]
Na Min, Hongyang Zhang, Xiaobin Gao, Pengyu Zeng .
Impacts of hunting cooperation and prey harvesting in a Leslie-Gower prey-predator system with strong Allee effect. AIMS Mathematics, 2024, 9(12): 34618-34646.
doi: 10.3934/math.20241649
[4]
Weili Kong, Yuanfu Shao .
Bifurcations of a Leslie-Gower predator-prey model with fear, strong Allee effect and hunting cooperation. AIMS Mathematics, 2024, 9(11): 31607-31635.
doi: 10.3934/math.20241520
[5]
Ali Yousef, Ashraf Adnan Thirthar, Abdesslem Larmani Alaoui, Prabir Panja, Thabet Abdeljawad .
The hunting cooperation of a predator under two prey's competition and fear-effect in the prey-predator fractional-order model. AIMS Mathematics, 2022, 7(4): 5463-5479.
doi: 10.3934/math.2022303
[6]
Vikas Kumar, Nitu Kumari .
Controlling chaos in three species food chain model with fear effect. AIMS Mathematics, 2020, 5(2): 828-842.
doi: 10.3934/math.2020056
[7]
Yalong Xue, Fengde Chen, Xiangdong Xie, Shengjiang Chen .
An analysis of a predator-prey model in which fear reduces prey birth and death rates. AIMS Mathematics, 2024, 9(5): 12906-12927.
doi: 10.3934/math.2024630
[8]
Xiaoming Su, Jiahui Wang, Adiya Bao .
Stability analysis and chaos control in a discrete predator-prey system with Allee effect, fear effect, and refuge. AIMS Mathematics, 2024, 9(5): 13462-13491.
doi: 10.3934/math.2024656
[9]
Wei Li, Qingkai Xu, Xingjian Wang, Chunrui Zhang .
Dynamics analysis of spatiotemporal discrete predator-prey model based on coupled map lattices. AIMS Mathematics, 2025, 10(1): 1248-1299.
doi: 10.3934/math.2025059
[10]
Weili Kong, Yuanfu Shao .
The effects of fear and delay on a predator-prey model with Crowley-Martin functional response and stage structure for predator. AIMS Mathematics, 2023, 8(12): 29260-29289.
doi: 10.3934/math.20231498
Abstract
In this paper, we primarily focused on the existence and uniqueness of the initial value problem for fractional order fuzzy ordinary differential equations in a fuzzy metric space. First, definitions and relevant properties of the Gamma function and Beta function within a fuzzy metric space were provided. Second, by employing the principle of fuzzy compression mapping and Choquet integral of fuzzy numerical functions, we established the existence and uniqueness of solutions to initial value problems for fuzzy ordinary differential equations. Finally, several examples were presented to demonstrate the validity of our obtained results.
1.
Introduction
In recent ecological research, the intricate dynamics inherent in three-species food chain models have sparked substantial interest, prompting a deeper understanding of complex ecological relationships [2,8,9,12,23,25,33,43,44]. This fascination has driven investigations into omnivory, a phenomenon where species feed across multiple trophic levels, with particular attention directed towards intraguild predation (IGP) [19,48,50,51]. IGP, a specialized form of omnivory, involves intricate interactions between a resource (prey), an intermediate predator preying exclusively on the prey and a top predator (referred to as an intraguild predator or IG predator) that consumes both the prey and the intermediate predators [14,16]. Given its pronounced impact on species coexistence and the emergence of chaos, the dynamics of IGP have become a central topic in ecological studies. Panja [38] analyzed a fractional-order prey-predator model with intraguild predation and observed instability in the system for integer-order derivatives. In contrast, the system remained stable with fractional-order derivatives. In [20], the authors explored intraguild predation dynamics, focusing on the impact of interference and predation strength on the coexistence of intraguild predators and prey. Gómez-Hernández et al. [13] explored an intraguild predation model and found that interspecific competition consistently reduces infection rates in both species.
Previous research has unveiled various mathematical models to capture the complex dynamics of IGP and its intricate interactions within predator-prey systems. In the realm of ecological modeling, diffusion stands as a fundamental factor closely tied to real environmental processes, influencing a myriad of dynamic behaviors. A wealth of research has delved into the dynamics of diffusive versions of ecological models, shedding light on their intricate mechanisms. Notably, Chen and Wu, in their work [4], illuminated the role of attractive predator-taxis in driving Turing bifurcation and the consequent emergence of spatiotemporal patterns within a predator-prey model characterized by nonlinear harvesting. Conversely, they found that the presence of repulsive predator-taxis or the complete absence thereof did not give rise to such phenomena. Moreover, this study provides valuable insights into how diffusion ratios and nonlinear harvesting practices impact the spatial distribution of species and the direction of Turing bifurcation.
In [5], the authors undertook a study focusing on the stability of bifurcating solutions within a predator-prey diffusive model with prey-taxis, where they not only confirmed the occurrence of steady state bifurcations but also explored the influence of repulsive prey-taxis in promoting this phenomenon. This research further delved into the determination of bifurcating solution stability. These collective findings offer valuable insights into the realm of ecological dynamics and diffusion-based ecological modeling. Hastings and Powell [15] highlighted the potential for chaos within three-species food chain models. On the other hand, Holt and Polis [17] introduced a Lotka-Volterra model characterized by linear functional responses, revealing that increased intensity of IGP could induce destabilization within the system [26]. Research in the field of intraguild predation (IP) has led to diverse insights into predator-prey dynamics. Hsu et al. [18] examined IGP scenarios using Beddington-DeAngelis responses, studying coexistence and competitive exclusion. They highlighted the significance of the interference parameter in shaping stability or oscillation. Kang and Wedekin [21] explored type III functional responses for intermediate predators and IG predators, revealing distinctions between generalist and specialist IG predators, and identifying multiple attractors. Freeze et al. [10] introduced ratio-dependent functional responses, investigating persistence, permanence and stability through innovative methods. These studies have laid crucial groundwork for understanding predator-prey interactions and set the stage for our exploration of cooperative hunting effects in the context of IGP dynamics. Building on this foundation, Tanabe and Namba [47] explored the emergence of chaotic dynamics within IGP scenarios, omitting intraspecific competition among IG prey and predators. Recent advancements by Sen et al. [44] and Majumder et al. [25] extended these models by incorporating saturating functional responses and intraspecific competition among IG prey and predator, leading to enhanced insights into the rich dynamics such as chaos of predator-prey interactions. However, Majumder et al. [25] control chaotic behavior of the system through a memory effect.
Furthermore, cooperative hunting behaviors have emerged as crucial factors shaping ecological stability. A range of adaptive characteristics have evolved among predators, including sharp teeth, claws, fangs and other tools, to facilitate efficient hunting. Predators have developed various strategies such as bait pursuit, collective hunting and other approaches. For instance, carnivores [24], lions [36], wolves [42] and African wild dogs [7] frequently collaborate to capture and subdue their prey. Similarly, ants [27] and spiders [49] exhibit cooperative behaviors in hunting, drawing considerable attention. The concept of hunting cooperation was initially introduced by Cosner et al. [6] in a predator-prey model. He proposed a functional response for predators that engage in spatially linear hunting formations and aggregate when encountering significant prey groups. Cooperative foraging has been explored in diverse predator-prey models, revealing scenarios of both destabilization and stabilization [1,3,6,11,31,32,34]. The incorporation of cooperative hunting behavior among predators in the prey-predator model reveals novel insights into the complex dynamics of ecological systems [45,46]. This growing understanding of cooperative behaviors' influence on ecological dynamics provides a strong motivation for our study.
We undertake a comprehensive investigation into the interplay between hunting cooperation among IG predators and IG prey, aiming to elucidate its impact on the intricate dynamics of predator-prey interactions. Our objective is to decipher how this cooperative dimension influences the emergence of chaotic oscillations and its potential to shape the overall stability or instability of ecological communities. By incorporating this cooperative aspect into our model, we aspire to provide novel insights into the intricate interactions between predator and prey, specifically within the context of IGP and the cooperative behaviors that can arise in these ecological systems. The distinctive feature of our study, setting it apart from existing literature [25,41,44,47] addressing intraguild predation (IGP) scenarios, is the inclusion of cooperative hunting behavior within the model. Through this endeavor, we seek to contribute to the deeper comprehension of the complexities inherent in predator-prey interactions and the factors that govern their dynamics.
2.
Model formulation
We focus on a three-species Lotka-Volterra model that incorporates intraguild predation dynamics. The model involves three major species denoted as x, y and z, representing a basal prey, an intermediate consumer (IG prey) and an omnivorous top predator (IG predator), respectively. We adopt a linear functional response between these species, as discussed in [47]. Additionally, we account for intraspecific competition within the intraguild prey and predator populations, as elaborated in [44]. Given the limited availability of resources, competition becomes essential for their sustenance within the community. In our model, we introduce the concept of hunting cooperation, which confers benefits to both intraguild prey and predator. Moreover, the success in capturing a basal prey and intraguild prey increases with the density of intraguild prey and predator, respectively, as highlighted in [37]. The model can be expressed as follows:
with non-negative initial conditions. The model parameters r11,d11,α11,α12,c11,α13,e11,α22,α23,c22,d22,e22,α31,α33 and d33 are all positive constants and their ecological descriptions are provided in Table 1. Our model is distinct from previous works like Kang and Wedekin [21] and Tanabe and Namba [47] due to the incorporation of intraspecific competitions and hunting cooperation behavior among intraguild prey and predator. Furthermore, compared to models studied by Sen et al. [44] and Majumder et al. [25], our model employs different functional response types between intraguild prey and predator, while also involving cooperative hunting.
Table 1.
Details of parameters with their values and sources.
The positivity of state variables is a fundamental concept in ecological modeling, indicating that the populations of the involved species remain non-negative within the system. Ensuring the positivity of solutions in the proposed model (2.1) is crucial for understanding the persistence of populations.
To establish the positivity of the system (2.1), we integrate its equations and utilize the given initial conditions to obtain the following expressions:
The positivity of these solutions follows from the fact that the right-hand side of each equation is non-negative given the non-negative initial conditions. Consequently, all solutions of the system (2.1) are guaranteed to be non-negative.
While we have established the positivity of solutions, we omit the discussion of the boundedness of solutions, as the boundedness property can be similarly proven and is straightforward in this context.
4.
Equilibria analysis
An equilibrium point in a dynamical system is a state or configuration where the system's variables remain constant over time. In other words, it's a point at which the system's behavior does not change, and the rates of change of the variables are all zero. Equilibrium points are also known as fixed points or steady states. Mathematically, for a dynamical system described by a set of equations, an equilibrium point is found by setting all the derivatives of the variables to zero. Equilibrium points play a crucial role in the analysis of dynamical systems. They serve as important reference states to understand the behavior of the system around those points.
The proposed IGP system has the following feasible equilibria.
(1) Population-free equilibrium E0(0,0,0) exists irrespective of any parametric restriction. From an ecological viewpoint, this equilibrium represents a theoretical scenario where all species involved in the model are entirely absent from the ecosystem. While such a scenario may not directly correspond to a real ecological state, it is a critical reference point.
(2) Basal prey equilibrium E1(r11−d11α11,0,0) exists under the parametric condition r11>d11. From an ecological viewpoint, this equilibrium reflects a scenario where the prey population is the only species present with intermediate and top predator both at zero. Moreover, the existence of E1 suggests that when the prey species experiences a birth rate greater than its death rate, it can maintain its population in the absence of predators.
(3) Basal prey-IG prey equilibrium E2(x2,y2,0), where x2 and y2 are positive solution(s) of following isoclines:
r11−d11−α11x−(α12+c11y)y=0,
(4.1)
e11(α12+c11y)x−α22y−d22=0.
(4.2)
We can obtain x from first equation and then putting the value of x in the second equation we get a cubic equation for y. This is little bit difficult to find the root (i.e., the component of equilibrium point). Therefore, we numerically check the existence E2 by plotting the above isoclines in x−y space as depicted in Figure 1. Figure exhibits unique E2 for the parameter values mentioned in the Table 1. The ecological significance of E2 lies in its representation of a complex ecological interaction where both basal prey and intraguild prey species coexist without the presence of intraguild predators. This equilibrium suggests the existence of resource partitioning mechanisms or the presence of ecological niches that allow both the basal prey and intraguild prey to thrive without direct predation pressure from higher trophic levels.
Figure 1.
Isoclines (4.1) and (4.2) are plotted here. Red dot represents the intersection of these lines. Parameters are chosen from Table 1.
(4) Basal prey-IG predator equilibrium E3(α13d33−α33d11+α33r11α11α33+α13α31,0,α31r11−α11d33−α31d11α11α33+α13α31).E3 exists if the parametric conditions α13d33+α33r11>α33d11 and α31r11>α11d33+α31d11 hold. Ecologically, this equilibrium illustrates the mechanisms that allow basal prey and intraguild predator to coexist without one driving the intermediate predator to extinction. By understanding the conditions that lead to the coexistence of these two species, ecologists can make more informed decisions regarding the preservation of species and ecosystem health.
(5) Coexistence equilibrium E4(x4,y4,z4), where the components x4,y4 and z4 can be obtained from following isoclines.
r11−d11−α11x−(α12+c11y)y−α13z=0,
(4.3)
e11(α12+c11y)x−α22y−(α23+c22z)z−d22=0,
(4.4)
e22(α23+c22z)y+α31x−α33z−d33=0.
(4.5)
Due to the difficulties of finding parametric values of the components of E4, we find it numerically by plotting above isoclines in x−y−z space as illustrated in Figure 2. It is apparent from the figure that for table parameter, E4 exists uniquely. The coexistence equilibrium E4 represents a critical ecological scenario within a model where all three interacting species coexist within an ecosystem. This equilibrium reflects the richness and stability of biodiversity, highlighting the ability of different species to thrive without the extinction of any species.
Figure 2.
Equations (4.3)–(4.5) are plotted here. Red dot denotes intersecting point of the planes, which is the E4. Parameters values are same as in Table 1.
Mathematically, local stability analysis involves linearizing the system's equations of motion around the equilibrium point. This is typically done by calculating the system's Jacobian matrix at the equilibrium, which describes how small changes in the variables affect the rates of change of the variables themselves. The eigenvalues of this Jacobian matrix are crucial indicators of the local stability properties of the system.
Local stability analysis provides valuable insights into the qualitative behavior of a dynamical system around an equilibrium point. However, it is important to note that local stability analysis only considers behavior near the equilibrium and may not fully capture more complex behaviors, such as oscillations, bifurcations and global attractors, that can arise in nonlinear systems. Here, we check local stability of the proposed system around the system's equilibria. First, we evaluate the Jacobian matrix of the system (2.1) as illustrated below.
(1) Eigenvalues of the Jacobian matrix J(E0) are −d11+r11,−d22 and −d33. From this eigenvalue analysis, we can say that E0 is stable under the parametric condition d11>r11. This stability condition is opposite to the existence condition of E1. Thus, E0 and E1 are related via a transcritical bifurcation. Further, we verify the stability of E0 numerically as described in Figure 3. In the figure, we choose r11=0.05 and d11=0.1. It is apparent from this figure that E0 is stable as condition d11>r11 holds.
(2) Eigenvalues of the Jacobian matrix J(E1) are d11−r11, −α12d11e11−α12e11r11+α11d22α11 and −α11d33+α31d11−α31r11α11. Hence, E1 is locally asymptotically stable if the parametric conditions d11<r11, α12d11e11+α11d22>α12e11r11 and α11d33+α31d11>α31r11 hold. E1 loses its stability when either of the conditions fails. We verify its local stability numerically also as described in Figure 4. The parameters used in this figure satisfy stability conditions, as a result, it is stable as depicted in the figure.
Figure 4.
Stability of E1. Here, green and red dots denote the stable and unstable equilibrium points, respectively. Parameters are taken from Table 1 except r11=0.5.
Now, eigenvalues of ˜J(E2) have negative real part if tr(˜J(E2))<0 and det(˜J(E2))>0. Thus, E2 is stable under the following conditions
e22α23y2+α31x2<d33,tr(˜J(E2))<0anddet(˜J(E2))>0.
As it is difficult to follow these conditions for stability of E2, we numerically verify its stability, Figure 5. The parameters used for numerical verification satisfy above stability conditions of E2. As a result E2 is locally asymptotically stable as described in the figure. (5) The Jacobian matrix evaluated at E3 is given by
J(E3)=[B11B12B13B21B22B23B31B32B33],
Figure 5.
Stability of E2. Here, green and red dots denote the stable and unstable equilibrium points, respectively. Parameters are taken from Table 1 except r11=0.75.
By Routh-Hurwitz criteria, E3 is stable if Γ2>0, Γ0>0 and Γ2Γ1>Γ0. The expression of Γi, i=0,1,2, are large to say about the stability E3. So, we numerically investigate its stability as depicted in Figure 6. It is apparent from the figure that E3 is stable since the parameter used here holds its stability conditions.
(6) The Jacobian matrix evaluated at a representative coexisting equilibrium point
E4=(x4,y4,z4)
Figure 6.
Stability of E3. Here, green and red dots denote the stable and unstable equilibrium points, respectively. Parameters are taken from Table 1 except e11=0.01.
The stability of E4 can be determined using the Routh-Hurwitz criteria, which states that E4 is stable if Δ2>0, Δ0>0 and Δ2Δ1>Δ0. The explicit expressions for the components of E4 are not obtainable due to their complexity, leading to the inability to derive parametric restrictions explicitly for the stability of E4. To validate the stability of E4, a numerical example is employed. The chosen parameter values in Figure 7 adhere to the Routh-Hurwitz criteria for the stability of E4. Consequently, E4 is demonstrated to be stable, as depicted in Figure 7.
Figure 7.
Stability of E4. Here, green and red dots denote the stable and unstable equilibrium points, respectively. Parameters are taken from Table 1 except r11=2.5.
Bifurcations in a dynamical system are critical points where the qualitative behavior of the system changes as one or more parameters are varied. In other words, they mark a point at which the stability, periodicity, or overall behavior of the system undergoes a significant transformation. Bifurcations can lead to the emergence of new patterns, such as oscillations, chaos or the creation and destruction of equilibrium points. Bifurcations are fundamental to the study of nonlinear dynamical systems and play a pivotal role in understanding the complex behaviors that can arise in such systems. They are often associated with changes in the number or stability of equilibrium points or the appearance of periodic solutions. Bifurcation analysis involves studying how the system's behavior changes as the parameters are continuously varied. In our proposed IGP system, two types of bifurcations occur: Transcritical and Hopf, which are described in this section.
Theorem 1.The system described by (2.1) undergoes a transcritical bifurcation at the parameter value
r11=r[TB]11=d11
concerning the equilibrium point E0(0,0,0).
Proof. Let us consider the Jacobian matrix J(E0)|r[TB]11. One of its eigenvalues is zero, and the other two are negative. Since two eigenvalues are consistently negative and the remaining one is zero, concluding stability-instability is not feasible at this point.
Let
V=(1,0,0)TandW=(1,0,0)T
be the eigenvectors of J(E0) and J(E0)T corresponding to the zero eigenvalue, respectively. Next, we define
F=[F1,F2,F3]T.
The conditions outlined by Sotomayor's theorem [40] for a transcritical bifurcation are as follows:
(1)WTFr11(E0,r[TB]11)=0,
(2)WT[DFr11(E0,r[TB]11)V]=0,
(3)WT[D2F(E0,r[TB]11)(V,V)]=−2α11≠0.
It's important to note that for a non-degenerate transcritical bifurcation, the second transversality condition should not be zero [40]. Hence, the system defined by (2.1) experiences a degenerate transcritical bifurcation [22] around E0 when
r11=r[TB]11.
Consequently, E0 undergoes a transcritical bifurcation as the parameter r11 crosses the critical value
r11=r[TB]11=d11,
and the conditions outlined in the theorem are satisfied. □ Next, we will establish the presence of a Hopf bifurcation concerning the model parameter r11. We have designated the critical value of r11 as r[HB]11.
Theorem 2.The system defined by (2.1) exhibits a Hopf bifurcation if the condition
(Δ1−3ω2)(dΔ0dr11−ω2dΔ2dr11)+2Δ2ω2dΔ1dr11≠0
is met when the parameter r11 crosses the critical value r11=r[HB]11.
Proof. Considering that the equilibrium point E4 is dependent on the parameter r11, we apply the Hopf bifurcation theorem [29]. The Hopf bifurcation occurs if the variation matrix J(E4) has one eigenvalue with a negative real part, where
Re(dλdr11)r11=r[HB]11≠0,
and the other two eigenvalues are purely imaginary. For purely imaginary eigenvalues, the characteristic Eq (5.2) should satisfy Δ2Δ1−Δ0=0. We assume that iω represents the purely imaginary eigenvalue corresponding to the parameter r[HB]11.
By differentiating both sides of the characteristic Eq (5.2) with respect to r11, we obtain
Therefore, the system described by (2.1) experiences a Hopf bifurcation concerning the parameter r11 if the condition
(Δ1−3ω2)(dΔ0dr11−ω2dΔ2dr11)+2Δ2ω2dΔ1dr11≠0
is satisfied at r11=r[HB]11.□
7.
Numerical simulation
In this section, we will delve into our bifurcation analysis results. We will start by examining the bifurcation structure concerning the parameter r11, as depicted in Figure 8.
Figure 8.
Bifurcation structure of the proposed model. Parameters are adopted from Table 1.
This figure uncovers dynamic behaviors as r11 undergoes variations. For lower values of r11, the system exclusively features the stable equilibrium E0. The presence of any other equilibria is absent within this parameter range. As r11 increases, a critical value is reached, leading to the emergence of the stable equilibrium E1. Simultaneously, the previously stable equilibrium E0 shifts from stability to instability, retaining this behavior for all subsequent higher values of r11. With further increments in r11, another transcritical point comes into view. At this juncture, the stable equilibrium E2 arises, causing the previously stable E1 to transition into an unstable state. This instability of E1 continues for higher values of r11. Finally, the system encounters an additional transcritical point, resulting in the emergence of the stable equilibrium E4, while the previously stable E2 transforms into an unstable state.
Similarly, Figure 9 presents a transcritical bifurcation diagram involving the parameter e11. For lower values of e11, the system displays a stable equilibrium E2. Beyond a critical value of e11, the stable equilibrium E4 comes into play. Consequently, the previously stable E2 shifts into an unstable state. These analyses provide valuable insights into the system's behavior with varying parameters r11 and e11. They reveal critical points where equilibria undergo shifts in stability and new equilibria emerge, contributing to a deeper comprehension of predator-prey interactions within the context of changing parameter values. Transcritical bifurcation is ecologically significant as it marks a critical point in the dynamics of predator-prey relationships and can lead to shifts in population equilibrium, which, in turn, can impact the stability and biodiversity of ecosystems. Understanding these dynamics is crucial for effective ecological management and conservation efforts.
Figure 9.
Transcritical bifurcation. Parameters are adopted from Table 1. Critical value of e11 is e[TC]11=0.083.
Continuing our exploration, we now delve into the system's behavior surrounding the coexistence equilibrium E4. This investigation involves the examination of phase portraits and time series solutions for various values of r11, as showcased in Figure 10. As the parameter r11 undergoes increments, the system's dynamics around the coexistence equilibrium E4 progress through distinct phases, leading to diverse patterns of behavior. For lower values of r11, the system exhibits a stable behavior around the coexistence equilibrium, manifesting as confined trajectories in the phase portrait, Figure 10a. As r11 increases, the system's behavior shifts to one-period oscillations. Trajectories trace out regular, repetitive patterns in the phase portrait, Figure 10b. Further increments in r11 lead to the emergence of two-period oscillations, Figure 10c. Trajectories now complete two cycles before repeating themselves in the phase portrait. The system continues to evolve with higher values of r11, resulting in four-period oscillations, Figure 10d. Trajectories now complete four cycles within the same phase portrait. As r11 further increases, the system transitions into chaotic behavior, Figure 10e. Trajectories become highly sensitive to initial conditions, leading to unpredictable and complex patterns in the phase portrait. Moreover, we observed that the chaotic attractors size increases on increaing r11 (see Figure 11a). Finally, for larger values of r11 both the predator populations extinct making the axial equilibrium point E1 stable (see Figure 11b). In essence, Figure 10 vividly portrays the evolution of the system's dynamics around the coexistence equilibrium E4 in response to changing values of r11. This comprehensive visualization allows us to discern the progression from stable behavior to chaotic complexity and highlights the intricate nature of predator-prey interactions in the context of varying parameter settings. Chaotic dynamics in predator-prey systems can lead to irregular and unpredictable population cycles. This unpredictability can help maintain biodiversity by preventing any one species from dominating the ecosystem for extended periods, fostering coexistence and preventing extinctions.
Figure 10.
Phase portraits and respective time series for the parameter r11=2.58 in (a); r11=2.85 in (b); r11=3.9 in (c); r11=4.35 in (d); r11=5 in (e). Rest are adopted from Table 1.
Moving forward, our analysis takes us to the exploration of species' maximum and minimum densities while systematically varying a key model parameter, as detailed in Figure 12. The insights gleaned from this visualization provide a comprehensive understanding of how the system's behavior evolves in response to specific parameter alterations. Figure 12a provides a revealing depiction of the system's response to changes in the parameter r11. Initially, the system sustains a stable coexistence state. This phase occurs after a transcritical bifurcation point where a switch occurs between equilibrium states. Then the system enters a phase of one-period oscillations, a behavior characterized by regular and predictable cycles. As r11 continues to increase, the system progresses into a phase of many-period oscillations. This pattern is marked by an increasing number of cycles within each oscillatory pattern, indicating a higher level of complexity in the population dynamics. Ultimately, the system experiences a transition into a chaotic state. In this phase, trajectories in the system become highly sensitive to initial conditions, leading to unpredictable and irregular behavior. This chaotic behavior underscores the intricate interactions that can arise within an ecosystem. Similar patterns of behavior are observed for the parameter α23 (Figure 12f), with the primary distinction being the absence of an initial transcritical bifurcation. Furthermore, we observe that parameter c22 (Figure 12b) influences the system's dynamics by inducing two cycles, leading to an increase in the number of cycles and finally chaos. The presence of chaos and transitions between stable states and periodic states can reflect the system's ability to adapt to changing environmental conditions. Ecosystems that exhibit such behavior may be more resilient, can support biodiversity, are capable of adjusting to disturbances and evolving in response to environmental shifts. On the other hand, parameters α11 (Figure 12c), α22 (Figure 12d), d11 (Figure 12e) and α12 (Figure 12g) exert control over the chaotic behavior by modulating the transition from many-period oscillations to one-period oscillations, bringing about a more ordered pattern and ultimately stable coexistence. The ability to control chaotic behavior through specific parameters is ecologically significant because it offers opportunities for managing, restoring and preserving ecosystems. Lastly, Figure 12h demonstrates that the parameter c11 is capable of driving the system towards more orderly behavior from an irregular behavior, as it produces a limit cycle through multiple cycles.
Figure 12.
Bifurcation structure of the proposed system. Parameter values are adopted from Table 1 except r11=4.35 in {b, c, d, e, f, g, h}.
Finally, we present the two-parameter bifurcation structure, varying the IG prey cooperation during hunting basal prey parameter (c11) and Consumption rate of IG predator to IG prey parameter (α23, as shown in Figure 13. This figure reveals the division of the entire parameter space into five subregions, each exhibiting distinct dynamics. The two-parameter bifurcation structure provide a comprehensive view of the model's behavior under varying two parameters, highlighting the occurrence of stable regions, limit cycles, two period, four period and chaotic dynamics. We note that for low values of consumption rate of IG predator, the cooperation of IG predator helps the system to keep in a stable state position. However, increasing in the consumption rate of IG predators push on the system into chaotic oscillations through periodic-doubling oscillations.
Figure 13.
Two parameter bifurcation structure. Table 1 is used for parameter values.
In conclusion, we delve into a three-species Lotka-Volterra model embodying intraguild predation dynamics. The model incorporates complex interactions between a basal prey, an intraguild prey and an omnivorous top predator. These interactions are defined by linear functional responses, alongside considerations for intraspecific competition and cooperative hunting, both of which play essential roles in shaping population dynamics within the ecological community. The proposed model diverges from prior works through its integration of intraspecific competition and cooperative hunting. We analyze equilibria and parametric conditions under which they are stable. Numerical stability of all steady states are justified. Our exploration encompasses a comprehensive bifurcation analysis that reveals intricate behaviors driven by parameter variations.
The examination of bifurcation structures, particularly concerning parameters such as r11 and e11, showcases critical shifts in equilibria stability and the emergence of new equilibrium states. We uncover transitions from stable to unstable states and vice versa, offering insights into the interplay between predator-prey interactions and varying parameter values. Moreover, our investigation into the dynamics around the coexistence equilibrium E4 illuminates a diverse range of behaviors as parameter r11 evolves. These behaviors include stable coexistence, one-period oscillations, multi-period oscillations and ultimately chaotic dynamics. This detailed analysis provides a nuanced understanding of the system's response to changing parameter values.
Furthering our comprehension, the exploration of species' maximum and minimum densities in response to parameter variations demonstrates distinct patterns. These patterns range from stable coexistence to chaotic behavior, each induced by different parameter influences. Parameters such as α12, c11 and d11 exert significant control over the system's chaotic dynamics, promoting stable or orderly behavior. Strong hunting cooperation behavior of intraguild predator produces chaos in the system whereas the same behavior of intraguild prey can control the chaotic nature of the system. In contrast to the findings of Sen et al. [44], our research demonstrates that a chaotic regime can emerge from a state of stable coexistence when considering the consumption rate of intraguild predators towards intraguild prey it may be due to the assumption of hunting cooperation. However, Vandermeer [50] observed chaotic oscillations in an intraguild predation (IGP) system in the absence of self-regulation for the same parameter. A recent study by Namba et al. [35] highlighted the stabilizing influence of strong self-regulation on a similar tri-trophic food web featuring type-I functional responses. Sen et al. [44] reported that robust self-regulation within both IG prey and predator populations can drive the dynamics of the system from a chaotic state to stable steady-state coexistence due to the saturated consumption of IG prey by IG predators. In our research, we extend this concept by revealing that strong self-regulation has the potential to guide the system's dynamics from a chaotic regime to stable coexistence. This stabilization is facilitated by the phenomenon of hunting cooperation. Also for low values of consumption rate of IG predator, the cooperation of IG predator helps the system to keep in a stable state position but increasing in the consumption rate of IG predators push on the system into chaotic oscillations through periodic-doubling oscillations.
In summary, this study offers a comprehensive understanding of complex predator-prey interactions in the context of intraguild predation dynamics. The integrated effects of intraspecific competition and cooperative hunting, coupled with systematic parameter variations, provide valuable insights into the intricate behaviors exhibited by ecological systems. Our findings converge as well as diverge from the existing literature. By shedding light on the intricate interplay between species and parameters, this research contributes to the broader field of ecological modeling and enhances our understanding of the complex dynamics that characterize natural ecosystems.
Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Acknowledgments
This work has been supported by the Deanship of Scientific Research at King Khalid University through the Research Group program with Grant Number RGP2/202/44. The authors are very grateful to the editor and reviewers for their careful reading of the manuscript.
Conflict of interest
The authors declare that they have no conflicts of interest.
References
[1]
K. S. Miller, B. Ross, An introduction to the fractional calculus and fractional differential equations, New York: Wiley, 1993.
[2]
A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Amsterdam, Boston: Elsevier, 2006.
[3]
N. A. Albasheir, A. Alsinai, A. U. K. Niazi, R. Shafqat, Romana, M. Alhagyan, et al., A theoretical investigation of Caputo variable order fractional differential equations: Existence, uniqueness, and stability analysis, Comput. Appl. Math., 42 (2023), 367. https://doi.org/10.1007/S40314-023-02520-6 doi: 10.1007/S40314-023-02520-6
[4]
V. Lakshmikantham, A. S. Vatsala, Basic theory of fractional differential equations, Nonlinear Anal., 69 (2008), 2677–2682. https://doi.org/10.1016/j.na.2007.08.042 doi: 10.1016/j.na.2007.08.042
[5]
M. Al-Refai, D. Baleanu, Comparison principles of fractional differential equations with non-local derivative and their applications, AIMS Mathematics, 6 (2021), 1443–1451. https://doi.org/10.3934/math.2021088 doi: 10.3934/math.2021088
[6]
L. Chen, G. Duan, S. Y. Wang, J. F. Ma, A Choquet integral based fuzzy logic approach to solve uncertain multi-criteria decision making problem, Expert Syst. Appl., 149 (2020), 113303. https://doi.org/10.1016/j.eswa.2020.113303 doi: 10.1016/j.eswa.2020.113303
[7]
M. H. Alshayeji, S. C. Sindhu, S. Abed, Viral genome prediction from raw human DNA sequence samples by combining natural language processing and machine learning techniques, Expert Syst. Appl., 218 (2023), 119641. https://doi.org/10.1016/j.eswa.2023.119641 doi: 10.1016/j.eswa.2023.119641
[8]
K. Tamilselvan, V. Visalakshi, P. Balaji, Applications of picture fuzzy filters: Performance evaluation of an employee using clustering algorithm, AIMS Mathematics, 8 (2023), 21069–21088. https://doi.org/10.3934/math.20231073 doi: 10.3934/math.20231073
[9]
S. Salahshour, A. Ahmadian, B. A. Pansera, M. Ferrara, Uncertain inverse problem for fractional dynamical systems using perturbed collage theorem, Commun. Nonlinear Sci. Numer. Simul., 94 (2021), 105553. https://doi.org/10.1016/j.cnsns.2020.105553 doi: 10.1016/j.cnsns.2020.105553
[10]
M. Z. Ahmad, M. K. Hasan, B. D. Baets, Analytical and numerical solutions of fuzzy differential equations, Inf. Sci., 236 (2013), 156–167. https://doi.org/10.1016/j.ins.2013.02.026 doi: 10.1016/j.ins.2013.02.026
[11]
M. Mazandarani, N. Pariz, A. V. Kamyad, Granular differentiability of fuzzy-number-valued functions, IEEE Trans. Fuzzy Syst., 26 (2018), 310–323. https://doi.org/10.1109/TFUZZ.2017.2659731 doi: 10.1109/TFUZZ.2017.2659731
[12]
Y. B. Shao, Q. Mou, Z. T. Gong, On retarded fuzzy functional differential equations and nonabsolute fuzzy integrals, Fuzzy Sets Syst., 375 (2019), 121–140. https://doi.org/10.1016/j.fss.2019.02.005 doi: 10.1016/j.fss.2019.02.005
[13]
T. Allahviranloo, A. Armand, Z. Gouyandeh, Fuzzy fractional differential equations under generalized fuzzy Caputo derivative, J. Intell. Fuzzy Syst., 26 (2014), 1481–1490.
[14]
N. V. Hoa, H. Vu, T. M. Duc, Fuzzy fractional differential equations under Caputo-Katugampola fractional derivative approach, Fuzzy Sets Syst., 375 (2019), 70–99. https://doi.org/10.1016/j.fss.2018.08.001 doi: 10.1016/j.fss.2018.08.001
[15]
M. Mazandarani, A. V. Kamyad, Modified fractional Euler method for solving fuzzy fractional intial value problem, Commun. Nonlinear Sci. Numer. Simul., 18 (2013), 12–21. https://doi.org/10.1016/j.cnsns.2012.06.008 doi: 10.1016/j.cnsns.2012.06.008
[16]
C. X. Wu, Z. T. Gong, On Henstock integrals of interval-valued and fuzzy-number-valued functions, Fuzzy Sets Syst., 115 (2000), 377–391. https://doi.org/10.1016/S0165-0114(98)00277-2 doi: 10.1016/S0165-0114(98)00277-2
[17]
D. Denneberg, Non-additive measure and integral, Dordrecht: Springer, 1994. https://doi.org/10.1007/978-94-017-2434-0
[18]
Z. T. Gong, L. Chen, G. Duan, Choquet integral of fuzzy-number-valued functions: The differentiability of the primitive with respect to fuzzy measures and choquet integral equations, Abstr. Appl. Anal., 2014 (2014), 953893. https://doi.org/10.1155/2014/953893 doi: 10.1155/2014/953893
[19]
C. V. Negoiţă, D. A. Ralescu, Application of fuzzy sets to system analysis, Basel: Springer, 1975. https://doi.org/10.1007/978-3-0348-5921-9
[20]
L. Stefanini, A generalization of Hukuhara difference and division for interval and fuzzy arithmetic, Fuzzy Sets Syst., 161 (2010), 1564–1584. https://doi.org/10.1016/j.fss.2009.06.009 doi: 10.1016/j.fss.2009.06.009
[21]
P. Diamond, P. Kloeden, Metric topology of fuzzy numbers and fuzzy analysis, In: Fundamentals of fuzzy sets, Boston: Springer, 2000. https://doi.org/10.1007/978-1-4615-4429-6_12
[22]
B. Bede, L. Stefanini, Generalized differentiability of fuzzy-valued functions, Fuzzy Sets Syst., 230 (2013), 119–141. https://doi.org/10.1016/j.fss.2012.10.003 doi: 10.1016/j.fss.2012.10.003
[23]
S. Salahshour, T. Allahviranloo, S. Abbasbandy, D. Baleanu, Existence and uniqueness results for fractional differential equations with uncertainty, Adv. Differ. Equ., 2012 (2012), 112. https://doi.org/10.1186/1687-1847-2012-112 doi: 10.1186/1687-1847-2012-112
This article has been cited by:
1.
Parvaiz Ahmad Naik, Muhammad Amer, Rizwan Ahmed, Sania Qureshi, Zhengxin Huang,
Stability and bifurcation analysis of a discrete predator-prey system of Ricker type with refuge effect,
2024,
21,
1551-0018,
4554,
10.3934/mbe.2024201
2.
Weirun Tao, Zhi-An Wang,
Global well-posedness and Turing–Hopf bifurcation of prey-taxis systems with hunting cooperation,
2025,
0956-7925,
1,
10.1017/S0956792525000026
3.
Alyaa Hussain Naser, Dahlia Khaled Bahlool,
Dynamics of an Intraguild Predation Food Web Cooperation Model Under the Influence of Fear and Hunting,
2025,
13,
2079-3197,
128,
10.3390/computation13060128