On Spectral Simulation of Fractional Brownian Motion

Abstract

C e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a PNA Probability, Networks and Algorithms Probability, Networks and AlgorithmsOn spectral simulation of fractional Brownian motion . Despite the availability of several exact simulation methods, attention has been paid to approximate simulation (i.e., the output is approximately fBm), particularly because of possible time savings. In this paper, we study the class of approximate methods that are based on the spectral properties of fBm's stationary incremental process, usually called fractional Gaussian noise (fGn). The main contribution is a proof of asymptotical exactness (in a sense that is made precise) of these spectral methods. Moreover, we establish the connection between the spectral simulation approach and a widely used method, originally proposed by Paxson, that lacked a formal mathematical justification. The insights enable us to evaluate the Paxson method in more detail. It is also shown that spectral simulation is related to the fastest known exact method. Mathematics Subject Classification: 65C05, 62M15

Estimation and pricing under long-memory stochastic volatility

“…We refer to [11] for a detailed exposition of all the techniques. Here, we choose to work with the circulant method, since it is faster than others (of order n log n).…”

Modeling and Computation of Bose-Einstein Condensates: Stationary States, Nucleation, Dynamics, Stochasticity

“…If we assume that we want to construct a process of length N ∈ N and being given the stationary process (δW H j ) j∈{1,...,N} , its spectral density [54,123] is such that By using this expression, the following spectral representation of the process [123] can be obtained …”

Bayesian estimation of self-similarity exponent

“…In general, several methods exist to generate numerically fBm data as required in the system I study here [39,67]. I note that here, as well as in the following chapters of this thesis, the synthetic data I use for model testing is based on the Cholesky decomposition of Σ H = L t L [56].…”