Next Article in Journal
The Probability Flow in the Stock Market and Spontaneous Symmetry Breaking in Quantum Finance
Previous Article in Journal
Perturbation Theory for Property (VE) and Tensor Product
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Asymptotic and Transient Behaviors of Stochastic Ratio-Dependent Predator–Prey Model

School of Mathematical Sciences, Tiangong University, Tianjin 300387, China
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(21), 2776; https://doi.org/10.3390/math9212776
Submission received: 29 September 2021 / Revised: 20 October 2021 / Accepted: 26 October 2021 / Published: 2 November 2021
(This article belongs to the Topic Dynamical Systems: Theory and Applications)

Abstract

:
In this paper, we focus on the asymptotic and transient dynamics of the studied ecosystem and measure the response to perturbation of the stochastic ratio-dependent predator–prey model. The method we use is mainly based on the Kronecker product and numerical simulation. Firstly, the mean-square stability matrix can be calculated from the Kronecker product, so as to compute three indicators (root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope) of the response to perturbation for the studied ecosystem. Since the above-measured amounts cannot be obtained explicitly, we use numerical simulation to draw the changing figures within the appropriate parameter range. Then we obtain some conclusions by comparing the numerical results. When perturbing any populations, increasing the disturbance intensity will reduce the mean-square stable area of the system. Ecologists can manage the ecosystem, reduce losses and maximize benefits according to the numerical results of the root-mean-square amplification envelope.

1. Introduction

The stability of an ecosystem has always been the focus of biologists and it involves many fields, such as the protection of the ecosystem, economic development, political management and so on. Traditionally, scholars are interested in the asymptotic stability of ecological models and have conducted several qualitative analyses and quantitative calculations. In previous studies, it is known that asymptotic stability can be characterized by resilience. In 1977, Pimm and Lawton [1] proposed the definition of resilience in a biological sense for the first time, which means the asymptotic decay rate of perturbations. Given a linearized (or linear) system
d x d t = A x
The resilience is calculated as v = ( λ m a x ( A ) ) , λ m a x ( A ) is the eigenvalue of A with the largest real part. The greater the resilience, the faster the decay rate of the disturbance in the long run and the more stable the system. Although the long-term asymptotic dynamic behavior is significant, the short-term transient dynamic behavior also cannot be ignored. Some systems are ultimately stable; however, the trajectory of these systems may rapidly move away from the equilibrium in the early stage after being disturbed. For example, the pests that will eventually become extinct can surge in the beginning. The infectious diseases, which will be solved, can cause a large number of deaths at first. Therefore, if we only consider the final state and ignore the temporary behavior of the system, we will suffer serious losses.
In the past few years, scholars have made remarkable achievements in transient dynamics for the ecosystem. In 1997, Neubert and Caswell [2] introduced reactivity as a supplement to resilience, which means the maximum initial rate that a small perturbation can grow. It characterizes the instantaneous responses ( t 0 ). The reactivity is denoted by v 0 = λ m a x ( H ( A ) ) , H ( A ) = A + A T 2 is the Hermitian part of the Jacobian matrix A. The equilibria is called reactive if v 0 > 0 . The larger the reactivity, the more sensitive the system and the faster a perturbation can grow. When λ m a x ( A ) < 0 and λ m a x ( H ( A ) ) > 0 , the system is eventually stable but reactive. This situation is possible and not accidental.
Resilience and reactivity describe two limits of solution behaviors ( t and t 0 ) , but neither describes the transient behaviors between zero and infinity. Therefore, Neubert and Caswell [2] defined the amplification envelope ρ ( t ) . It is the maximum possible amplification that any initial-value perturbation could induce in a solution at time t, given by ρ ( t ) = e A t , where e A t is the Euclidean matrix norm of e x p ( A t ) . In addition, the maximum value of ρ ( t ) is denoted by ρ m a x = max t 0 ρ ( t ) . The time at which ρ m a x occurs is denoted as t m a x .
Based on these conclusions, scholars further studied the resilience, reactivity and envelope amplification of ecosystems. Regarding resilience, Folke, Carpenter and Walker [3] had a heated discussion around resilience and put forward several significant concepts related to ecological resilience, such as adaptability, transformability, etc. In the same year, Hughes and Graham et al. [4] found that resilience played an important role in the management of marine ecosystems (such as coral reef systems), and proposed that we should learn how to avoid unnecessary phase transformations of ecosystems. In 2019, Wright [5] proposed a model of social resilience and sustainable development in the 21st century from the lessons of historical ecology. For reactivity, Snyder [6] analyzed what makes the ecosystem reactive. Then, Vesipa and Ridolfi [7] also focused on the short-term dynamics of long-term stable and reactive ecosystems in 2017. Taking a seasonal forced prey–predator model as an example, it is proved that seasonality can greatly affect the transient dynamics of the system.
In 2014, Evelyn Buckwar and Conall Kelly [8] extend the above measure indicators to stochastic cases for a perturbed predator–prey model, which are called root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope. Since ecosystems exist in a real world and are affected by the external environment continuously, it is more reasonable to consider a random biological model. Compared with the traditional predator–prey model, the ratio-dependent predator–prey model proposed by Arditi and Ginzburg [9] in 1989 avoids the limitation that the trophic function of the traditional predator–prey model ignores the influence of predators. Therefore, it is more meaningful to consider that the stochastic ratio-dependent predator–prey model, whose characteristic is that a functional response depends on the ratio of prey to predator density.
To the best of our knowledge, Berezovskaya, Karev and Arditi [10] made a complete parameter analysis of the stability properties and dynamics of the deterministic ratio-dependent predator–prey model. Recently, some scholars have made careful studies on the dynamic behaviors of the deterministic predator–predator models, such as Chen and Zhang [11] in 2021 and Molla, Sahabuddin and Sajid [12], in 2022. Furthermore, Bandyopadhyay and Chattopadhyay [13] discussed the deterministic and stochastic behaviors and stochastic stability near the internal equilibrium of the ratio-dependent predator–predator model. They also gave the conditions for the deterministic model to enter the Hopf bifurcation. Wu, Huang and Wang [14] researched the dynamical behavior of a stochastic ratio-dependent predator–prey system. They obtained the conclusions of stochastically ultimate boundedness and stochastic permanence, among which persistence is one of the most important characteristics of the biological model. In addition, in 2016, Fan, Teng and Muhammadhaji [15] analyzed the global dynamics of a stochastic ratio-dependent predator–prey system. They gave results for the existence, uniqueness and global attractivity of the solution for the stochastic model. Ji, Jiang and Li [16] described the connection and difference between a deterministic and random ratio-dependent predator–predator model. As mentioned in the article, the previous literature had shown that the solution of the deterministic model is uniformly bounded. In addition, the author had obtained that the solution of the random model is uniformly bounded in mean, which is similar to the deterministic system. They also studied the persistence and extinction of model populations. The conditions under which the stochastic model species is not permanent but the deterministic model species may be permanent are given.
In this paper, we will study the root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope of stochastic ratio-dependent predator–prey model. The layout of the paper is as follows. In Section 2, we briefly introduce the studied model and some basic concepts. Then, in Section 3, we compute the mean-square stability matrix S of stochastic ratio-dependent prey–predator model by the Kronecker product. Meanwhile, we use the numerical method to simulate parameter graphs of the three measure indicators: root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope. Finally, the analysis conclusions are obtained by comparing the numerical figures in Section 4.

2. Preliminaries

2.1. Ratio-Dependent Predator–Prey Model

First, we briefly introduce the ratio-dependent predator–prey model:
d N d t = f ( N ) N g N P P d P d t = e g N P P μ P
where N , P represents the prey abundance and predator abundance, correspondingly. f ( N ) is the natural growth rate of prey, μ and e are the natural death rate and conversion efficiency of predators, respectively. The functional response g ( N P ) is a function of ratio of prey to predator. Normally, we consider f ( N ) as the logistic form: f ( N ) = r 1 N K , and intrinsic growth rate is r, carrying capacity is K. The functional response g ( N P ) will be considered as g ( N P ) = α N P P + α h N with maximum prey consumption rate α and predator handling time h. Therefore, Equation (2) can be written as the following form:
d N d t = r N 1 N K α N P P + α h N d P d t = e α N P P + α h N μ P

2.2. Kronecker Product and Mean-Square Stability Matrix

Definition 1
(Kronecker product (see the literature [17])). Given an m × n matrix A = [ a i j ] and a p × q matrix B = [ b i j ] , the Kronecker product of A and B is defined as follows:
A B = a 11 B a 1 n B a m 1 B a m n B .
Definition 2.
For a linear (or linearized) stochastic system
d X ( t ) = A X ( t ) d t + i = 1 m B i X ( t ) d W i ,
the mean-square stability matrix S of (5) is defined as
S = I d A + A I d + i = 1 m B i B i .
We denote the spectral abscissa α ( A ) of a matrix A by α ( A ) = m a x ( λ i ) , where λ i are the eigenvalues of A.
Lemma 1
(See Reference [8]).
  • The root-mean-square resilience of Equation (2) is denoted by
    1 2 α ( S ) .
  • The root-mean-square reactivity of Equation (2) is denoted by
    1 2 λ m a x ( S + S T 2 ) .
  • The root-mean-square amplification envelope of Equation (2) is denoted by
    e S t .

3. Stochastic Ratio-Dependent Predator–Prey Model

3.1. Nondimensionalized Ratio-Dependent Predator–Prey Model

We use a nondimensionalized method to simplify the deterministic ratio-dependent model. Setting
N = K x P = K α h y t = τ r a = α r b = e r h c = μ r
In this way, we can obtain the new system
d x d τ = x ( 1 x ) a x y x + y d y d τ = b x y x + y c y
For the convenience of description, τ is denoted by t next, that is, the new system is as follows:
d x d t = x ( 1 x ) a x y x + y d y d t = b x y x + y c y
where d is the mortality rate of redefined predator y. Parameter a and b characterize trophic function: a is the maximum prey mortality, b is the maximum predator growth rate. The system has three equilibrium points: O ( x = 0 , y = 0 ) and E 1 ( x = 1 , y = 0 ) are trivial equilibriums, E 2 ( x = x * , y = y * ) is non-trivial equilibrium. Here, we only consider the non-trivial equilibrium E 2 ( x = x * , y = y * ) in biological sense. After calculation, we can obtain
x * = 1 a ( b c ) b y * = b c c x *
Linearize the system (8) at E 2 ( x * , y * ) and obtain the Jacobian matrix:
A = a a c 2 b 2 1 a c 2 b 2 ( b c ) 2 b c + c 2 b

3.2. Generalization to Stochastic Systems

We introduce the stochastic perturbation term into Equation (8) to consider the influence of random fluctuation environment. Refer to [18,19,20] for the method of adding random parameters. This is equivalent to conducting random sensitivity analysis on deterministic systems. Let w 1 and w 2 be independent Wiener processes, and consider the stochastic system
d x = x ( 1 x ) a x y x + y d t + σ 1 ( x x * ) d w 1 d y = [ b x y x + y c y ] d t + σ 2 ( y y * ) d w 2
Next, linearize the system (11) at E 2 ( x * , y * ) and replace ( x x * , y y * ) with ( X , Y ) to rewrite the stochastic system in the matrix form as follows:
d X Y = a a c 2 b 2 1 a c 2 b 2 ( b c ) 2 b c + c 2 b X Y d t
+ σ 1 0 0 0 X Y d w 1 + 0 0 0 σ 2 X Y d w 2
According to Definition 2, we can calculate the mean-square stability matrix S:
σ 1 2 + 2 a 2 a c 2 b 2 2 a c 2 b 2 a c 2 b 2 0 ( b c ) 2 b a + c 2 b a c 2 b 2 c 1 0 a c 2 b 2 ( b c ) 2 b 0 a + c 2 b a c 2 b 2 c 1 a c 2 b 2 0 ( b c ) 2 b ( b c ) 2 b σ 2 2 + 2 c 2 b 2 c
Since it is impossible to explicitly write the eigenvalues in closed form, within the selected parameters range, we use S to simulate the parameter figures of root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope. In the next section, we will display these figures and make a detailed analysis.

4. Numerical Simulation and Analysis

First, choose the appropriate parameter values a = 2.0 , b = 0.7808 , c = 0.5 for the drawing of root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope numerical figures. (The selection of parameters refers to aticle [13]). When a = 2 , b = 0.7808 and c = 0.5 , the model system has a stable Hopf-bifurcating periodic solution and all other trajectories around the limit cycle eventually approach it.

4.1. Root-Mean-Square Resilience and Root-Mean-Square Reactivity

Figure 1 shows the variation of the root-mean-square resilience and root-mean-square reactivity of the studied system within the parameter range when there is no random disturbance. It can be used as a comparison figure for analyzing the influence of random environmental fluctuations on the system. In Figure 1, the parameter ranges of the mean-square stable stage for this system are a [ 0 , 2 ] , b [ 0.5 , 0.7808 ] , c = [ 0.5 , 0.7808 ] , respectively. Resilience shows a trend of monotonous increase first and then decrease. If b and c are changing parameters, the interior equilibria of the system are both reactive (reactivity > 0 ).
The first column of Figure 2 shows the perturbation only applied to the prey, the disturbance intensity is 0.2 and 0.4, correspondingly. The variable parameter is a. It can be seen that the mean-square stable area of the system is reduced from a [ 0 , 2 ] ( σ 1 = 0 , σ 2 = 0 ) to a [ 0 , 1.96 ] ( σ 1 = 0.2 , σ 2 = 0 ) and a [ 0 , 1.8 ] ( σ 1 = 0.4 , σ 2 = 0 ) . The maximum value of resilience (representing the system’s self-recovery ability; the greater the resilience, the more stable the system) is reduced from 0.357 to 0.313 and 0.275. The second column is the disturbance to predator and, similarly, the mean-square stable area is reduced from a [ 0 , 2 ] ( σ 1 = 0 , σ 2 = 0 ) to a [ 0 , 1.96 ] ( σ 1 = 0 , σ 2 = 0.2 ) and a [ 0 , 1.85 ] ( σ 1 = 0 , σ 2 = 0.4 ) . The maximum value of resilience is reduced from 0.357 to 0.298 and 0.210. Regarding the reactivity, there is no obvious change when the disturbance is applied to the prey, but when it acts on the predator, the initial value of the reactivity becomes larger as the intensity of the disturbance increases. Summarizing the above conclusions, the system is more sensitive to predator disturbance and reacts more strongly.
In Figure 3, the range of the change parameter b is [ 0 , 0.7808 ] . When the perturbation is only on the prey, the mean-square stable regions of the system are b [ 0.5 , 0.763 ] ( σ 1 = 0.2 , σ 2 = 0 ) and b [ 0.5 , 0.722 ] ( σ 1 = 0.4 , σ 2 = 0 ) , respectively. The maximum values of resilience are 0.198 and 0.163. When b [ 0.519 , 0.765 ] ( σ 1 = 0 , σ 2 = 0.2 ) and b [ 0.568 , 0.734 ] ( σ 1 = 0 , σ 2 = 0.4 ) , the model is mean-square stable if the predator is disturbed. Comparing the four sub-figures in Figure 3, it is found that increasing the intensity of the disturbance increases the reactivity, the resilience decreases and the mean-square stable area shrinks. Similar to Figure 2, it is still the predator disturbance that has a stronger impact on the system.
Table 1 is a brief description of Figure 4.
The trends of resilience and reactivity are opposite, which also shows that the more stable the system, the stronger the resilience; the more sensitive the system, the stronger the reactivity.
Compare Figure 2, Figure 3 and Figure 4 with Figure 1, we obtain some conclusions as follows:
  • Increasing the random perturbation to any species (predator or prey) will reduce the range of variable parameters of the system in the mean-square asymptotical stable state (where the values of root-mean-square resilience are positive).
  • In Figure 2, Figure 3 and Figure 4, we can clearly see that when perturbation is added to any population, the maximum root-mean-square resilience is visibly reduced. In contrast, when the random perturbation is increased for any species, the value of root-mean-square reactivity increases but the increase is not obvious.
  • The same random disturbance intensity acts on predator and prey, respectively, the system will be more strongly affected when the predator is disturbed.

4.2. Root-Mean-Square Amplification Envelope (Hopf-Bifurcating Periodic Solution)

Figure 5 shows the root-mean-square amplification envelope of the positive equilibrium of (11) with a = 2.0 , b = 0.7808 , c = 0.5 . The positive equilibrium is a stable Hopf-bifurcating periodic solution. Figure 5 illustrates the type of this solution because the change of disturbance is also periodic.
It can be found that adding perturbation to any population at the Hopf-bifurcating equilibrium (even if the perturbation is small) will make the system unstable. In comparison, perturbation to the predator makes the system more sensitive.
Conclusions can be drawn from all the above analyses, and some appropriate suggestions can be made to managers for this ecosystem. For example, when regulating the ecosystem, pay attention to controlling the parameters in the mean-square asymptotically stable state of system (the value of the root-mean-square resilience is positive). Second, according to the numerical graph of root-mean-square amplification envelope, managers can clearly know the maximum value of disturbance magnification: ρ m a x and the time when the maximum disturbance appears: t m a x . If the manager sets a small disturbance to the ecosystem, it needs to consider that the perturbation amplification ρ m a x should be controlled within the range that the system can bear, that is, the system will be mean-square asymptotically stable. t m a x gives critical information, and decision makers should give orders and actively implement them before t m a x , so as to avoid crucial losses and save money and manpower.
In order to make the reader better understand the significance of root-mean-square amplification envelope, the following numerical example is listed.

4.3. Root-Mean-Square Amplification Envelope (Locally Stable Positive Equilibrium)

In this example, select the parameter values a = 1.8 , b = 0.7808 , c = 0.5 . According to literature [13], the interior equilibrium is locally stable in this case. From Figure 6, we find that when the predator is perturbed, and the perturbation intensity reaches 0.4, the system can still restore a stable state through its own regulation. However, the perturbation can no longer disappear and even shows a constant increasing trend when the prey is perturbed with the perturbation intensity 0.4. At the same time, we can see that no matter which species is disturbed, the greater the perturbation intensity, the larger the amplification envelope.

5. Conclusions

By comparing the root-mean-square resilience and root-mean-square reactivity figures of the stochastic model under different disturbance intensities, we reach the following conclusions: As the perturbation intensity increases, the root-mean-square resilience decreases and, on the contrary, the root-mean-square reactivity increases. As a result, the mean square stability area of the random system reduces, and the mean square stability becomes weaker. In addition, numerical results show that the system is more sensitive when the predator is disturbed. Observing and comparing the figures of the root-mean-square envelope amplification, the following points are found: (1) Adding a disturbance to the Hopf-bifurcation stable periodic solution of the deterministic system, no matter how small the disturbance is, the stability of the solution will be destroyed. (2) The greater the disturbance intensity, the greater the root-mean-square envelope amplification. (3) Applying perturbation to the locally stable positive solution, when the perturbation intensity is small enough, the solution will return to the initial state. However, if the perturbation exceeds critical value, the system will be destroyed and the solution cannot be restored.

Author Contributions

Software and writing—original draft, W.L.; project administration and funding acquisition, J.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China OF FUNDER grant number 2018YFC1406403, National Natural Science Foundation of China OF FUNDER grant number Nos. 11501410 and 11471243.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Acknowledgments

This research is supported with funds provided by the National Key Research and Development Program of China (2018YFC1406403) and the National Natural Science Foundation of China (Nos. 11501410 and 11471243). We appreciate the anonymous reviewers for their valuable comments and helpful suggestions which improve this paper significantly.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pimm, S.L.; Lawton, J.H. Number of trophic levels in ecological communities. Nature 1977, 268, 329–331. [Google Scholar] [CrossRef]
  2. Neubert, M.G.; Caswell, H. Alternatives to resilience for measuring the responses of ecological systems to perturbations. Ecology 1997, 78, 653–665. [Google Scholar] [CrossRef]
  3. Folke, C.; Carpenter, S.R.; Walker, B.; Scheffer, M.; Chapin, T.; Rockström, J. Resilience Thinking: Integrating Resilience, Adaptability and Transformability. Ecol. Soc. 2010, 15, 20. [Google Scholar] [CrossRef]
  4. Hughes, T.P.; Graham, N.A.; Jackson, J.B.; Mumby, P.J.; Steneck, R.S. Rising to the challenge of sustaining coral reef resilience. Trends Ecol. Evol. 2010, 25, 633–642. [Google Scholar] [CrossRef] [PubMed]
  5. Wright, D.K. Long-term dynamics of pastoral ecology in northern Kenya: An old model for new resilience. J. Archaeol. Sci. 2019, 55, 101068. [Google Scholar] [CrossRef]
  6. Snyder, R.E. What makes ecological systems reactive? Theor. Popul. Biol. 2010, 77, 243–249. [Google Scholar] [CrossRef] [PubMed]
  7. Vesipa, R.; Ridolfi, L. Impact of seasonal forcing on reactive ecological systems. J. Theoret. Biol. 2017, 419, 23–35. [Google Scholar] [CrossRef] [PubMed]
  8. Buckwar, E.; Kelly, C. Asymptotic and transient mean-square properties of stochastic systems arising in ecology, fluid dynamics, and system control. SIAM J. Appl. Math. 2004, 74, 411–433. [Google Scholar] [CrossRef]
  9. Arditi, R.; Ginzburg, L.R. Coupling in predator-prey dynamics: Ratio-dependence. J. Theoret. Biol. 1989, 139, 311–326. [Google Scholar] [CrossRef]
  10. Berezovskaya, F.; Karev, G.; Arditi, R. Parametric analysis of the ratio-dependent predator-prey model. J. Math. Biol. 2001, 43, 221–246. [Google Scholar] [CrossRef]
  11. Chen, X.; Zhang, X. Dynamics of the predator–prey model with the Sigmoid functional response. Stud. Appl. Math. 2021, 1–19. [Google Scholar] [CrossRef]
  12. Molla, H.; Sarwardi, S.; Sajid, M. Predator-prey dynamics with Allee effect on predator species subject to intra-specific competition and nonlinear prey refuge. J. Math. Comput. Sci. 2022, 25, 150–165. [Google Scholar] [CrossRef]
  13. Bandyopadhyay, M.; Chattopadhyay, J. Ratio-dependent predator–prey model: Effect of environmental fluctuation and stability. Nonlinearity 2005, 18, 913–936. [Google Scholar] [CrossRef]
  14. Wu, Z.; Huang, H.; Wang, L. Dynamical behavior of a Stochastic ratio-dependent predator-prey System. J. Appl. Math. 2012, 2012, 1281–1302. [Google Scholar] [CrossRef]
  15. Fan, X.; Teng, Z.; Muhammadhaji, A. Global dynamics of a stochastic ratio-dependent predator–prey system. Asian-Eur. J. Math. 2016. [Google Scholar] [CrossRef]
  16. Ji, C.; Jiang, D.; Li, X. Qualitative analysis of a stochastic ratio-dependent predator–prey system. J. Comput. Appl. Math. 2011, 235, 1326–1341. [Google Scholar] [CrossRef]
  17. Magnus, J.R.; Neudecker, H. Matrix Differential Calculus with Applications in Statistics and Econometrics; John Wiley and Sons: New York, NY, USA, 2019; ISBN 0-471-98632-1. [Google Scholar]
  18. Baishya, M.C.; Chakrabarti, C.G. Non-equilibrium fluctuation in Lotka–Volterra system. Bull. Math. Biol. 1987, 49, 125–131. [Google Scholar] [CrossRef]
  19. Bandyopadhyay, M.; Chakrabarti, C.G. Deterministic and stochastic analysis of a nonlinear prey–predator system. J. Biol. Syst. 2003, 11, 161–172. [Google Scholar] [CrossRef]
  20. Tapaswi, P.K.; Mukhopadhyay, A. Effects of environmental fluctuation on plankton allelopathy. J. Math. Biol. 1999, 39, 39–58. [Google Scholar] [CrossRef]
Figure 1. Resilience and reactivity (in the mean-square sense) of the positive equilibrium of (11), σ 1 = σ 2 = 0 . In (a), b = 0.7808 , c = 0.5 . In (b), a = 2.0 , c = 0.5 . In (c), a = 2.0 , b = 0.7808 .
Figure 1. Resilience and reactivity (in the mean-square sense) of the positive equilibrium of (11), σ 1 = σ 2 = 0 . In (a), b = 0.7808 , c = 0.5 . In (b), a = 2.0 , c = 0.5 . In (c), a = 2.0 , b = 0.7808 .
Mathematics 09 02776 g001
Figure 2. Resilience and reactivity (in the mean-square sense) of the positive equilibrium of (11). The variable parameter is a, fixed parameters b = 0.7808 , c = 0.5 . In (a), σ 1 = 0.2 , σ 2 = 0 . In (b), σ 1 = 0.4 , σ 2 = 0 . In (c), σ 1 = 0 , σ 2 = 0.2 . In (d), σ 1 = 0 , σ 2 = 0.4 .
Figure 2. Resilience and reactivity (in the mean-square sense) of the positive equilibrium of (11). The variable parameter is a, fixed parameters b = 0.7808 , c = 0.5 . In (a), σ 1 = 0.2 , σ 2 = 0 . In (b), σ 1 = 0.4 , σ 2 = 0 . In (c), σ 1 = 0 , σ 2 = 0.2 . In (d), σ 1 = 0 , σ 2 = 0.4 .
Mathematics 09 02776 g002
Figure 3. Resilience and reactivity (in the mean-square sense) of the positive equilibrium of (11). The variable parameter is b, fixed parameters a = 2.0 , c = 0.5 . In (a), σ 1 = 0.2 , σ 2 = 0 . In (b), σ 1 = 0.4 , σ 2 = 0 . In (c), σ 1 = 0 , σ 2 = 0.2 . In (d), σ 1 = 0 , σ 2 = 0.4 .
Figure 3. Resilience and reactivity (in the mean-square sense) of the positive equilibrium of (11). The variable parameter is b, fixed parameters a = 2.0 , c = 0.5 . In (a), σ 1 = 0.2 , σ 2 = 0 . In (b), σ 1 = 0.4 , σ 2 = 0 . In (c), σ 1 = 0 , σ 2 = 0.2 . In (d), σ 1 = 0 , σ 2 = 0.4 .
Mathematics 09 02776 g003
Figure 4. Resilience and reactivity (in the mean-square sense) of the positive equilibrium of (11). The variable parameter is c, fixed parameters a = 2.0 , b = 0.7808 . In (a), σ 1 = 0.2 , σ 2 = 0 . In (b), σ 1 = 0.4 , σ 2 = 0 . In (c), σ 1 = 0 , σ 2 = 0.2 . In (d), σ 1 = 0 , σ 2 = 0.4 .
Figure 4. Resilience and reactivity (in the mean-square sense) of the positive equilibrium of (11). The variable parameter is c, fixed parameters a = 2.0 , b = 0.7808 . In (a), σ 1 = 0.2 , σ 2 = 0 . In (b), σ 1 = 0.4 , σ 2 = 0 . In (c), σ 1 = 0 , σ 2 = 0.2 . In (d), σ 1 = 0 , σ 2 = 0.4 .
Mathematics 09 02776 g004
Figure 5. Root-mean-square amplification envelope of the positive equilibrium of (11). In (a), blue solid line: σ 1 = 0 , σ 2 = 0 ; red dash line: σ 1 = 0.1 , σ 2 = 0 ; green solid line: σ 1 = 0.15 , σ 2 = 0 . In (b), blue solid line: σ 1 = 0 , σ 2 = 0 ; red dash line: σ 1 = 0 , σ 2 = 0.1 ; green solid line: σ 1 = 0 , σ 2 = 0.15 .
Figure 5. Root-mean-square amplification envelope of the positive equilibrium of (11). In (a), blue solid line: σ 1 = 0 , σ 2 = 0 ; red dash line: σ 1 = 0.1 , σ 2 = 0 ; green solid line: σ 1 = 0.15 , σ 2 = 0 . In (b), blue solid line: σ 1 = 0 , σ 2 = 0 ; red dash line: σ 1 = 0 , σ 2 = 0.1 ; green solid line: σ 1 = 0 , σ 2 = 0.15 .
Mathematics 09 02776 g005
Figure 6. Amplification envelope(in the mean-square sense) of the positive equilibrium of (11). In (a), blue solid line: σ 1 = 0 , σ 2 = 0 ; red dash line: σ 1 = 0.2 , σ 2 = 0 ; black dash line: σ 1 = 0.4 , σ 2 = 0 . In (b), blue solid line: σ 1 = 0 , σ 2 = 0 ; red dash line: σ 1 = 0 , σ 2 = 0.2 ; black dash line: σ 1 = 0 , σ 2 = 0.4 .
Figure 6. Amplification envelope(in the mean-square sense) of the positive equilibrium of (11). In (a), blue solid line: σ 1 = 0 , σ 2 = 0 ; red dash line: σ 1 = 0.2 , σ 2 = 0 ; black dash line: σ 1 = 0.4 , σ 2 = 0 . In (b), blue solid line: σ 1 = 0 , σ 2 = 0 ; red dash line: σ 1 = 0 , σ 2 = 0.2 ; black dash line: σ 1 = 0 , σ 2 = 0.4 .
Mathematics 09 02776 g006
Table 1. Information summary of Figure 4.
Table 1. Information summary of Figure 4.
Disturbance IntensityParameter Value Range of Mean Square Stable StageMaximum Value of Resilience
σ 1 = 0.2 , σ 2 = 0 d [ 0.52 , 0.78 ] 0.226
σ 1 = 0.4 , σ 2 = 0 d [ 0.51 , 0.76 ] 0.213
σ 1 = 0 , σ 2 = 0.2 d [ 0.54 , 0.78 ] 0.190
σ 1 = 0 , σ 2 = 0.4 d [ 0.53 , 0.70 ] 0.142
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, W.; Feng, J. Analysis of Asymptotic and Transient Behaviors of Stochastic Ratio-Dependent Predator–Prey Model. Mathematics 2021, 9, 2776. https://doi.org/10.3390/math9212776

AMA Style

Liu W, Feng J. Analysis of Asymptotic and Transient Behaviors of Stochastic Ratio-Dependent Predator–Prey Model. Mathematics. 2021; 9(21):2776. https://doi.org/10.3390/math9212776

Chicago/Turabian Style

Liu, Wen, and Jianfeng Feng. 2021. "Analysis of Asymptotic and Transient Behaviors of Stochastic Ratio-Dependent Predator–Prey Model" Mathematics 9, no. 21: 2776. https://doi.org/10.3390/math9212776

APA Style

Liu, W., & Feng, J. (2021). Analysis of Asymptotic and Transient Behaviors of Stochastic Ratio-Dependent Predator–Prey Model. Mathematics, 9(21), 2776. https://doi.org/10.3390/math9212776

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

Article Metrics

Back to TopTop