An algorithm for simulating Brownian increments on a sphere

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Mijatović, Aleksandar, Mramor, Veno, Bravo, Gerónimo Uribe
Format: Preprint
Veröffentlicht: 2020
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866909556352745472
author Mijatović, Aleksandar
Mramor, Veno
Bravo, Gerónimo Uribe
author_facet Mijatović, Aleksandar
Mramor, Veno
Bravo, Gerónimo Uribe
contents This paper presents a novel formula for the transition density of the Brownian motion on a sphere of any dimension and discusses an algorithm for the simulation of the increments of the spherical Brownian motion based on this formula. The formula for the density is derived from an observation that a suitably transformed radial process (with respect to the geodesic distance) can be identified as a Wright-Fisher diffusion process. Such processes satisfy a duality (a kind of symmetry) with a certain coalescent processes and this in turn yields a spectral representation of the transition density, which can be used for exact simulation of their increments using the results of Jenkins and Spanò (2017). The symmetry then yields the algorithm for the simulation of the increments of the Brownian motion on a sphere. We analyse the algorithm numerically and show that it remains stable when the time-step parameter is not too small.
format Preprint
id arxiv_https___arxiv_org_abs_2012_12018
institution arXiv
publishDate 2020
record_format arxiv
spellingShingle An algorithm for simulating Brownian increments on a sphere
Mijatović, Aleksandar
Mramor, Veno
Bravo, Gerónimo Uribe
Statistical Mechanics
Numerical Analysis
Probability
65C05
This paper presents a novel formula for the transition density of the Brownian motion on a sphere of any dimension and discusses an algorithm for the simulation of the increments of the spherical Brownian motion based on this formula. The formula for the density is derived from an observation that a suitably transformed radial process (with respect to the geodesic distance) can be identified as a Wright-Fisher diffusion process. Such processes satisfy a duality (a kind of symmetry) with a certain coalescent processes and this in turn yields a spectral representation of the transition density, which can be used for exact simulation of their increments using the results of Jenkins and Spanò (2017). The symmetry then yields the algorithm for the simulation of the increments of the Brownian motion on a sphere. We analyse the algorithm numerically and show that it remains stable when the time-step parameter is not too small.
title An algorithm for simulating Brownian increments on a sphere
topic Statistical Mechanics
Numerical Analysis
Probability
65C05
url https://arxiv.org/abs/2012.12018