Construction of positivity-preserving numerical method for stochastic SIVS epidemic model

In this paper we propose the balanced implicit numerical techniques for maintaining the nonnegative path of the solution in stochastic susceptible–infected–vaccinated–susceptible (SIVS) epidemic model. We can hardly acquire the explicit solution for the SIVS model, so we often use the numerical scheme to produce approximate solutions. The Euler–Maruyama (EM) method is a useful and effective means in producing numerical solutions of SIVS model. The EM method to simulate the stochastic SIVS model often results in the problem that the numerical solution is not positive. In order to eliminate the negative path of the solution in a stochastic SIVS epidemic model, we construct a numerical method preserving positivity for the SIVS model. It is proved that the balanced implicit method (BIM) can preserve positivity and we show the convergence of the BIM numerical approximate solution to the exact solution. Finally, a numerical example is offered to support the theoretical results and verify the availability of the approach.

[1]  Yanan Zhao,et al.  The threshold of a stochastic SIRS epidemic model with saturated incidence , 2014, Appl. Math. Lett..

[2]  Yaozhong Hu.  Semi-Implicit Euler-Maruyama Scheme for Stiff Stochastic Equations , 1996 .

[3]  Henri Schurz,et al.  Numerical Regularization for SDEs: Construction of Nonnegative Solutions , 1995 .

[4]  E. Platen,et al.  Balanced Implicit Methods for Stiff Stochastic Systems , 1998 .

[5]  Xuerong Mao,et al.  The truncated Euler-Maruyama method for stochastic differential equations , 2015, J. Comput. Appl. Math..

[6]  Chuncheng Wang,et al.  Modeling the effects of vaccination on rotavirus infection , 2015, Advances in Difference Equations.

[7]  Donal O'Regan,et al.  The extinction and persistence of the stochastic SIS epidemic model with vaccination , 2013 .

[8]  Xianning Liu,et al.  Global stability of an SEIR epidemic model with age-dependent latency and relapse ☆ , 2015 .

[9]  Qun Liu,et al.  The threshold of a stochastic SIS epidemic model with imperfect vaccination , 2018, Math. Comput. Simul..

[10]  Lansun Chen,et al.  The effect of impulsive vaccination on an SIR epidemic model , 2009, Appl. Math. Comput..

[11]  Pasquale Vetro,et al.  SIVR epidemic model with stochastic perturbation , 2012, Neural Computing and Applications.

[12]  Fengying Wei,et al.  Persistence and extinction for an age-structured stochastic SVIR epidemic model with generalized nonlinear incidence rate , 2019, Physica A: Statistical Mechanics and its Applications.

[13]  Marija Milosevic,et al.  The Euler-Maruyama approximation of solutions to stochastic differential equations with piecewise constant arguments , 2016, J. Comput. Appl. Math..

[15]  Daqing Jiang,et al.  Stationary distribution of a stochastic SIS epidemic model with vaccination , 2014 .

[16]  M. Bayram,et al.  Numerical methods for simulation of stochastic differential equations , 2018 .

[17]  Michael Günther,et al.  Structure preserving stochastic integration schemes in interest rate derivative modeling , 2008 .

[18]  Lin-Fei Nie,et al.  Dynamic behavior analysis of SIVS epidemic models with state-dependent pulse vaccination , 2018 .

[19]  Yongfeng Guo,et al.  Construction of positivity preserving numerical method for stochastic age-dependent population equations , 2017, Appl. Math. Comput..

[20]  M. Martcheva,et al.  Global threshold dynamics of an SIVS model with waning vaccine-induced immunity and nonlinear incidence. , 2015, Mathematical biosciences.

[21]  L. Szpruch,et al.  Convergence, Non-negativity and Stability of a New Milstein Scheme with Applications to Finance , 2012, 1204.1647.