Source author record

Samuel Herrmann

Samuel Herrmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
1topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2023arXiv

Exact simulation of the first passage time through a given level for jump diffusions

Continuous-time stochastic processes play an important role in the description of random phenomena, it is therefore of prime interest to study particular variables depending on their paths, like stopping time for example. One approach consists in pointing out explicit expressions of the probability distributions, an other approach is rather based on the numerical generation of the random variables. We propose an algorithm in order to generate the first passage time through a given level of a one-dimensional jump diffusion. This process satisfies a stochastic differential equation driven by a Brownian motion and subject to random shocks characterized by an independent Poisson process. Our algorithm belongs to the family of rejection sampling procedures, also called exact simulation in this context: the outcome of the algorithm and the stopping time under consideration are identically distributed. It is based on both the exact simulation of the diffusion at a given time and on the exact simulation of first passage time for continuous diffusions. It is therefore based on an extension of the algorithm introduced by Herrmann and Zucca [16] in the continuous framework. The challenge here is to generate the exact position of a continuous diffusion conditionally to the fact that the given level has not been reached before. We present the construction of the algorithm and give numerical illustrations, conditions on the recurrence of jump diffusions are also discussed.

preprint2020arXiv

Exact simulation of diffusion first exit times: algorithm acceleration

In order to describe or estimate different quantities related to a specific random variable, it is of prime interest to numerically generate such a variate. In specific situations, the exact generation of random variables might be either momentarily unavailable or too expensive in terms of computation time. It therefore needs to be replaced by an approximation procedure. As was previously the case, the ambitious exact simulation of exit times for diffusion processes was unreachable though it concerns many applications in different fields like mathematical finance, neuroscience or reliability. The usual way to describe exit times was to use discretization schemes, that are of course approximation procedures. Recently, Herrmann and Zucca \cite{Herrmann-Zucca-2} proposed a new algorithm, the so-called GDET-algorithm (General Diffusion Exit Time), which permits to simulate exactly the exit time for one-dimensional diffusions. The only drawback of exact simulation methods using an acceptance-rejection sampling is their time consumption. In this paper the authors highlight an acceleration procedure for the GDET-algorithm based on a multi-armed bandit model. The efficiency of this acceleration is pointed out through numerical examples.

preprint2016arXiv

Initial-Boundary Value Problem for the heat equation - A stochastic algorithm

The Initial-Boundary Value Problem for the heat equation is solved by using a new algorithm based on a random walk on heat balls. Even if it represents a sophisticated generalization of the Walk on Spheres (WOS) algorithm introduced to solve the Dirich-let problem for Laplace's equation, its implementation is rather easy. The definition of the random walk is based on a new mean value formula for the heat equation. The convergence results and numerical examples permit to emphasize the efficiency and accuracy of the algorithm.

preprint2015arXiv

The first-passage time of the Brownian motion to a curved boundary: an algorithmic approach

Under some weak conditions, the first-passage time of the Brownian motion to a continuous curved boundary is an almost surely finite stopping time. Its probability density function (pdf) is explicitly known only in few particular cases. Several mathematical studies proposed to approximate the pdf in a quite general framework or even to simulate this hitting time using a discrete time approximation of the Brownian motion. The authors study a new algorithm which permits to simulate the first-passage time using an iterating procedure. The convergence rate presented in this paper suggests that the method is very efficient.

preprint2015arXiv

The walk on moving spheres: a new tool for simulating Brownian motion's exit time from a domain

In this paper we introduce a new method for the simulation of the exit time and position of a $δ$-dimensional Brownian motion from a domain. The main interest of our method is that it avoids splitting time schemes as well as inversion of complicated series. The idea is to use the connexion between the $δ$-dimensional Bessel process and the $δ$-dimensional Brownian motion thanks to an explicit Bessel hitting time distribution associated with a particular curved boundary. This allows to build a fast and accurate numerical scheme for approximating the hitting time. Numerical comparisons with existing methods are performed.

preprint2014arXiv

Mean-field limit versus small-noise limit for some interacting particle systems

In the nonlinear diffusion framework, stochastic processes of McKean-Vlasov type play an important role. In some cases they correspond to processes attracted by their own probability distribution: the so-called self-stabilizing processes. Such diffusions can be obtained by taking the hydrodymamic limit in a huge system of linear diffusions in interaction. In both cases, for the linear and the nonlinear processes, small-noise asymptotics have been emphasized by specific large deviation phenomenons. The natural question, therefore, is: is it possible to interchange the mean-field limit with the small-noise limit? The aim here is to consider this question by proving that the rate function of the first particle in a mean-field system converges to the rate function of the hydrodynamic limit as the number of particles becomes large.

preprint2014arXiv

Simulation of hitting times for Bessel processes with non integer dimension

In this paper we pursue and complete the study of the simulation of the hitting time of some given boundaries for Bessel processes. These problems are of great interest in many application fields as finance and neurosciences. In a previous work, the authors introduced a new method for the simulation of hitting times for Bessel processes with integer dimension. The method was based mainly on the explicit formula for the distribution of the hitting time and on the connexion between the Bessel process and the Euclidean norm of the Brownian motion. This method does not apply anymore for a non integer dimension. In this paper we consider the simulation of the hitting time of Bessel processes with non integer dimension and provide a new algorithm by using the additivity property of the laws of squared Bessel processes. We split each simulation step in two parts: one is using the integer dimension case and the other one considers hitting time of a Bessel process starting from zero.

preprint2013arXiv

Hitting time for Bessel processes - walk on moving spheres algorithm (WoMS)

In this article we investigate the hitting time of some given boundaries for Bessel processes. The main motivation comes from mathematical finance when dealing with volatility models, but the results can also be used in optimal control problems. The aim here is to construct a new and efficient algorithm in order to approach this hitting time. As an application we will consider the hitting time of a given level for the Cox-Ingersoll-Ross process. The main tools we use are on one side, an adaptation of the method of images to this particular situation and on the other side, the connection that exists between Cox-Ingersoll-Ross processes and Bessel processes.

preprint2013arXiv

Statistics of transitions for Markov chains with periodic forcing

The influence of a time-periodic forcing on stochastic processes can essentially be emphasized in the large time behaviour of their paths. The statistics of transition in a simple Markov chain model permits to quantify this influence. In particular the first Floquet multiplier of the associated generating function can be explicitly computed and related to the equilibrium probability measure of an associated process in higher dimension. An application to the stochastic resonance is presented.

preprint2012arXiv

Persistent random walks, variable length Markov chains and piecewise deterministic Markov processes

A classical random walk $(S_t, t\in\mathbb{N})$ is defined by $S_t:=\displaystyle\sum_{n=0}^t X_n$, where $(X_n)$ are i.i.d. When the increments $(X_n)_{n\in\mathbb{N}}$ are a one-order Markov chain, a short memory is introduced in the dynamics of $(S_t)$. This so-called "persistent" random walk is nolonger Markovian and, under suitable conditions, the rescaled process converges towards the integrated telegraph noise (ITN) as the time-scale and space-scale parameters tend to zero (see Herrmann and Vallois, 2010; Tapiero-Vallois, Tapiero-Vallois2}). The ITN process is effectively non-Markovian too. The aim is to consider persistent random walks $(S_t)$ whose increments are Markov chains with variable order which can be infinite. This variable memory is enlighted by a one-to-one correspondence between $(X_n)$ and a suitable Variable Length Markov Chain (VLMC), since for a VLMC the dependency from the past can be unbounded. The key fact is to consider the non Markovian letter process $(X_n)$ as the margin of a couple $(X_n,M_n)_{n\ge 0}$ where $(M_n)_{n\ge 0}$ stands for the memory of the process $(X_n)$. We prove that, under a suitable rescaling, $(S_n,X_n,M_n)$ converges in distribution towards a time continuous process $(S^0(t),X(t),M(t))$. The process $(S^0(t))$ is a semi-Markov and Piecewise Deterministic Markov Process whose paths are piecewise linear.