Cholera, as an endemic disease around the world, has imposed great harmful effects on human health. In addition, from a microscopic viewpoint, the interference of random factors exists in the process of virus replication. However, there are few theoretical studies of viral infection models with biologically reasonable stochastic effects. This paper studied a stochastic cholera model used to describe transmission dynamics in China. In this paper, we adopted a special method to simulate the effect of environmental perturbations to the system instead of using linear functions of white noise, i.e., the transmission rate of environment to human was satisfied Ornstein–Uhlenbeck processes, which is a more practical and interesting. First, it was theoretically proved that the solution to the stochastic model is unique and global, with an ergodic stationary distribution. Moreover, by solving the corresponding Fokker–Planck equation and using our developed algebraic equation theory, we obtain the exact expression of probability density function around the quasi-equilibrium of the stochastic model. Finally, several numerical simulations are provided to confirm our analytical results.
Citation: Ying He, Bo Bi. Threshold dynamics and density function of a stochastic cholera transmission model[J]. AIMS Mathematics, 2024, 9(8): 21918-21939. doi: 10.3934/math.20241065
Cholera, as an endemic disease around the world, has imposed great harmful effects on human health. In addition, from a microscopic viewpoint, the interference of random factors exists in the process of virus replication. However, there are few theoretical studies of viral infection models with biologically reasonable stochastic effects. This paper studied a stochastic cholera model used to describe transmission dynamics in China. In this paper, we adopted a special method to simulate the effect of environmental perturbations to the system instead of using linear functions of white noise, i.e., the transmission rate of environment to human was satisfied Ornstein–Uhlenbeck processes, which is a more practical and interesting. First, it was theoretically proved that the solution to the stochastic model is unique and global, with an ergodic stationary distribution. Moreover, by solving the corresponding Fokker–Planck equation and using our developed algebraic equation theory, we obtain the exact expression of probability density function around the quasi-equilibrium of the stochastic model. Finally, several numerical simulations are provided to confirm our analytical results.
[1] | R. M. Anderson, R. M. May, Infectious diseases of humans: dynamics and control, Oxford: Oxford University Press, 1991. https://doi.org/10.1093/oso/9780198545996.001.0001 |
[2] | M. J. Keeling, P. Rohani, Modeling infectious diseases in humans and animals, Princeton: Princeton University Press, 2008. https://doi.org/10.1515/9781400841035 |
[3] | M. Marathe, A. K. S. Vullikanti, Computational epidemiology, Commun. ACM., 56 (2013), 88–96. https://doi.org/10.1145/2483852.2483871 doi: 10.1145/2483852.2483871 |
[4] | J. P. Tian, J. Wang, Global stability for cholera epidemic models, Math. Biosci., 232 (2011), 31–41. https://doi.org/10.1016/j.mbs.2011.04.001 doi: 10.1016/j.mbs.2011.04.001 |
[5] | G. Q. Sun, J. H. Xie, S. H. Huang, Z. Jin, M. T. Li, L. Liu, Transmission dynamics of cholera: mathematical modeling and control strategies, Commun. Nonlinear Sci. Numer. Simul., 45 (2017), 235–244. https://doi.org/10.1016/j.cnsns.2016.10.007 doi: 10.1016/j.cnsns.2016.10.007 |
[6] | Y. Cai, J. Jiao, Z. Gui, Y. Liu, W. Wang, Environmental variability in a stochastic epidemic model, Appl. Math. Comput., 329 (2018), 210–226. https://doi.org/10.1016/j.amc.2018.02.009 doi: 10.1016/j.amc.2018.02.009 |
[7] | B. Han, B. Zhou, D. Jiang, T. Hayat, A. Alsaedi, Stationary solution, extinction and density function for a high-dimensional stochastic SEI epidemic model with general distributed delay, Appl. Math. Comput., 405 (2021), 126236. https://doi.org/10.1016/j.amc.2021.126236 doi: 10.1016/j.amc.2021.126236 |
[8] | O. M. Otunuga, Estimation of epidemiological parameters for COVID-19 cases using a stochastic SEIRS epidemic model with vital dynamics, Results Phys., 28 (2021), 104664. https://doi.org/10.1016/j.rinp.2021.104664 doi: 10.1016/j.rinp.2021.104664 |
[9] | X. Zhang, R. Yuan, A stochastic chemostat model with mean-reverting Ornstein–Uhlenbeck process and Monod-Haldane response function, Appl. Math. Comput., 394 (2021), 125833. https://doi.org/10.1016/j.amc.2020.125833 doi: 10.1016/j.amc.2020.125833 |
[10] | W. Wang, Y. Cai, Z. Ding, Z. Gui, A stochastic differential equation SIS epidemic model incorporating Ornstein–Uhlenbeck process, Phys. A., 509 (2018), 921–936. https://doi.org/10.1016/j.physa.2018.06.099 doi: 10.1016/j.physa.2018.06.099 |
[11] | D. Jiang, N. Shi, X. Li, Global stability and stochastic permanence of a non-autonomous logistic equation with random perturbation, J. Math. Anal. Appl., 340 (2008), 588–597. https://doi.org/10.1016/j.jmaa.2007.08.014 doi: 10.1016/j.jmaa.2007.08.014 |
[12] | C. Huang, S. Gan, D. Wang, Delay-dependent stability analysis of numerical methods for stochastic delay differential equations, J. Comput. Appl. Math., 236 (2012), 3514–3527. https://doi.org/10.1016/j.cam.2012.03.003 doi: 10.1016/j.cam.2012.03.003 |
[13] | D. Li, The stationary distribution and ergodicity of a stochastic generalized logistic system, Stat. Probabil. Lett., 83 (2013), 580–583. https://doi.org/10.1016/j.spl.2012.11.006 doi: 10.1016/j.spl.2012.11.006 |
[14] | M. Liu, K. Wang, Staionary distribution, ergodicity and extinction of a stochastic generalized logistic system, Appl. Math. Lett., 25 (2012), 1980–1985. https://doi.org/10.1016/j.aml.2012.03.015 doi: 10.1016/j.aml.2012.03.015 |
[15] | S. Zhang, X. Meng, T. Feng, T. Zhang, Dynamics analysis and numerical simulations of a stochastic non-autonomous predator–prey system with impulsive effects, Nonlinear Anal.: Hybrid Syst., 26 (2017), 19–37. https://doi.org/10.1016/j.nahs.2017.04.003 doi: 10.1016/j.nahs.2017.04.003 |
[16] | X. Mao, Stochastic differential equations and their applications, Chichester: Horwood Publishing, 1997. |
[17] | B. Zhou, D. Jiang, Y. Dai, T. Hayat, Stationary distribution and density function expression for a stochastic SIQRS epidemic model with temporary immunity, Nonlinear Dyn., 105 (2021), 931–955. https://doi.org/10.1007/s11071-020-06151-y doi: 10.1007/s11071-020-06151-y |
[18] | B. Zhou, D. Jiang, B. Han, T. Hayat, Threshold dynamics and density function of a stochastic epidemic model with media coverage and mean-reverting Ornstein–Uhlenbeck process, Math. Comput. Simulat., 196 (2022), 15–44. https://doi.org/10.1016/j.matcom.2022.01.014 doi: 10.1016/j.matcom.2022.01.014 |
[19] | N. H. Du, D. H. Nguyen, G. G. Yin, Conditions for permanence and ergodicity of certain stochastic predator–prey models, J. Appl. Probab., 53 (2016), 187–202. https://doi.org/10.1017/jpr.2015.18 doi: 10.1017/jpr.2015.18 |
[20] | Z. Shi, D. Jiang, Dynamical behaviors of a stochastic HTLV-I infection model with general infection form and Ornstein–Uhlenbeck process, Chaos Solitons Fract., 165 (2022), 112789. https://doi.org/10.1016/j.chaos.2022.112789 doi: 10.1016/j.chaos.2022.112789 |
[21] | H. Roozen, An asymptotic solution to a two-dimensional exit problem arising in population dynamics, SIAM. J. Appl. Math., 49 (1989), 1793–1810. https://doi.org/10.1137/0149110 doi: 10.1137/0149110 |