Samuel Livingstone, Nikolas Nüsken, Giorgos Vasdekis, Rui-Yang Zhang
{"title":"Skew-symmetric schemes for stochastic differential equations with non-Lipschitz drift: an unadjusted Barker algorithm","authors":"Samuel Livingstone, Nikolas Nüsken, Giorgos Vasdekis, Rui-Yang Zhang","doi":"arxiv-2405.14373","DOIUrl":null,"url":null,"abstract":"We propose a new simple and explicit numerical scheme for time-homogeneous\nstochastic differential equations. The scheme is based on sampling increments\nat each time step from a skew-symmetric probability distribution, with the\nlevel of skewness determined by the drift and volatility of the underlying\nprocess. We show that as the step-size decreases the scheme converges weakly to\nthe diffusion of interest. We then consider the problem of simulating from the\nlimiting distribution of an ergodic diffusion process using the numerical\nscheme with a fixed step-size. We establish conditions under which the\nnumerical scheme converges to equilibrium at a geometric rate, and quantify the\nbias between the equilibrium distributions of the scheme and of the true\ndiffusion process. Notably, our results do not require a global Lipschitz\nassumption on the drift, in contrast to those required for the Euler--Maruyama\nscheme for long-time simulation at fixed step-sizes. Our weak convergence\nresult relies on an extension of the theory of Milstein \\& Tretyakov to\nstochastic differential equations with non-Lipschitz drift, which could also be\nof independent interest. We support our theoretical results with numerical\nsimulations.","PeriodicalId":501215,"journal":{"name":"arXiv - STAT - Computation","volume":null,"pages":null},"PeriodicalIF":0.0000,"publicationDate":"2024-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - STAT - Computation","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2405.14373","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
We propose a new simple and explicit numerical scheme for time-homogeneous
stochastic differential equations. The scheme is based on sampling increments
at each time step from a skew-symmetric probability distribution, with the
level of skewness determined by the drift and volatility of the underlying
process. We show that as the step-size decreases the scheme converges weakly to
the diffusion of interest. We then consider the problem of simulating from the
limiting distribution of an ergodic diffusion process using the numerical
scheme with a fixed step-size. We establish conditions under which the
numerical scheme converges to equilibrium at a geometric rate, and quantify the
bias between the equilibrium distributions of the scheme and of the true
diffusion process. Notably, our results do not require a global Lipschitz
assumption on the drift, in contrast to those required for the Euler--Maruyama
scheme for long-time simulation at fixed step-sizes. Our weak convergence
result relies on an extension of the theory of Milstein \& Tretyakov to
stochastic differential equations with non-Lipschitz drift, which could also be
of independent interest. We support our theoretical results with numerical
simulations.