File Exchange

image thumbnail

Fractional Brownian motion generator

version 2.1.0.0 (1.73 KB) by Zdravko Botev
Generates fractional Brownian motion with a given Hurst parameter using the FFT.

10 Downloads

Updated 20 Jun 2016

View License

fast one dimensional fractional Brownian motion (FBM) generator
output is 'W_t' with t in [0,T] using 'n' equally spaced grid points;
code uses Fast Fourier Transform (FFT) for speed.
INPUT:
- Hurst parameter 'H' in [0,1]
- number of grid points 'n', where 'n' is a power of 2;
if the 'n' supplied is not a power of two,
then we set n=2^ceil(log2(n)); default is n=2^12;
- final time 'T'; default value is T=1;
OUTPUT:
- Fractional Brownian motion 'W_t' for 't';
- time 't' at which FBM is computed;
If no output it invoked, then function plots the FBM.
Example: plot FBM with hurst parameter 0.95 on the interval [0,10]
[W,t]=fbm1d(0.95,2^12,10); plot(t,W)
Reference:
Kroese, D. P., & Botev, Z. I. (2015). Spatial Process Simulation.
In Stochastic Geometry, Spatial Statistics and Random Fields(pp. 369-404)
Springer International Publishing, DOI: 10.1007/978-3-319-10064-7_12

Comments and Ratings (1)

Sanghai Guan

Very good.

Updates

2.1.0.0

- bug fix (courtesy of Viktor Bezborodov): W=sqrt(T)*W ---> W=T^H*W; (unless final time T=1, this would affect the scaling).

2.0.0.0

- Vectorized for speed
- Written as an 'm' file with error control

MATLAB Release Compatibility
Created with R2015b
Compatible with any release
Platform Compatibility
Windows macOS Linux

Discover Live Editor

Create scripts with code, output, and formatted text in a single executable document.


Learn About Live Editor