{"title":"Positivity-Preserving Numerical Method for a Stochastic Multi-Group SIR Epidemic Model","authors":"Han Ma, Qimin Zhang, X. Xu","doi":"10.1515/cmam-2022-0143","DOIUrl":null,"url":null,"abstract":"Abstract The stochastic multi-group susceptible–infected–recovered (SIR) epidemic model is nonlinear, and so analytical solutions are generally difficult to obtain. Hence, it is often necessary to find numerical solutions, but most existing numerical methods fail to preserve the nonnegativity or positivity of solutions. Therefore, an appropriate numerical method for studying the dynamic behavior of epidemic diseases through SIR models is urgently required. In this paper, based on the Euler–Maruyama scheme and a logarithmic transformation, we propose a novel explicit positivity-preserving numerical scheme for a stochastic multi-group SIR epidemic model whose coefficients violate the global monotonicity condition. This scheme not only results in numerical solutions that preserve the domain of the stochastic multi-group SIR epidemic model, but also achieves the “ order - 1 2 {\\mathrm{order}-\\frac{1}{2}} ” strong convergence rate. Taking a two-group SIR epidemic model as an example, some numerical simulations are performed to illustrate the performance of the proposed scheme.","PeriodicalId":48751,"journal":{"name":"Computational Methods in Applied Mathematics","volume":null,"pages":null},"PeriodicalIF":1.0000,"publicationDate":"2022-12-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computational Methods in Applied Mathematics","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1515/cmam-2022-0143","RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
Abstract The stochastic multi-group susceptible–infected–recovered (SIR) epidemic model is nonlinear, and so analytical solutions are generally difficult to obtain. Hence, it is often necessary to find numerical solutions, but most existing numerical methods fail to preserve the nonnegativity or positivity of solutions. Therefore, an appropriate numerical method for studying the dynamic behavior of epidemic diseases through SIR models is urgently required. In this paper, based on the Euler–Maruyama scheme and a logarithmic transformation, we propose a novel explicit positivity-preserving numerical scheme for a stochastic multi-group SIR epidemic model whose coefficients violate the global monotonicity condition. This scheme not only results in numerical solutions that preserve the domain of the stochastic multi-group SIR epidemic model, but also achieves the “ order - 1 2 {\mathrm{order}-\frac{1}{2}} ” strong convergence rate. Taking a two-group SIR epidemic model as an example, some numerical simulations are performed to illustrate the performance of the proposed scheme.
期刊介绍:
The highly selective international mathematical journal Computational Methods in Applied Mathematics (CMAM) considers original mathematical contributions to computational methods and numerical analysis with applications mainly related to PDEs.
CMAM seeks to be interdisciplinary while retaining the common thread of numerical analysis, it is intended to be readily readable and meant for a wide circle of researchers in applied mathematics.
The journal is published by De Gruyter on behalf of the Institute of Mathematics of the National Academy of Science of Belarus.