Main Content


Parameter estimation of fractional Brownian motion



    hest = wfbmesti(X) returns estimates of the fractal index H of the input signal X.


    collapse all

    This example shows how to estimate the Hurst index of a fractional Brownian motion. The example simulates 1,000 realizations of fractional Brownian motion with H=0.6. Each realization consists of 10,000 samples. At the end of the simulation, the three estimates of the Hurst index are compared.

    Initialize the random number generator for repeatable results. Set the Hurst index equal to 0.6 and the length of the realizations to be 10,000.

    rng default
    H = 0.6;
    len = 10000;

    Generate 1,000 realizations of fractional Brownian motion and compute the estimates of the Hurst parameter.

    n = 1000; 
    Hest = zeros(n,3);
    for ii = 1:n
    	fBm06 = wfbm(H,len);
    	Hest(ii,:) = wfbmesti(fBm06);

    Compare the estimates.

    title("Discrete Second Derivative Estimator (DSOD)")
    title("Wavelet Version of DSOD") 
    title("Wavelet Details Regression Estimator")
    xlabel("True value of the parameter H = 0.6")

    Input Arguments

    collapse all

    Input signal, specified as a vector. The signal X is assumed to be a realization of fractional Brownian motion with Hurst index H.

    Data Types: double

    Output Arguments

    collapse all

    Fractal index estimates, returned as a 1-by-3 vector. hest contains three estimates of the fractal index H.

    • The first two elements of hest are estimates based on the second derivative with the second computed in the wavelet domain.

    • The third estimate is based on the linear regression in loglog plot, of the variance of detail versus level.

    More About

    collapse all

    Fractional Brownian Motion

    A fractional Brownian motion (fBm) is a continuous-time Gaussian process depending on the Hurst parameter 0 < H < 1. It generalizes the ordinary Brownian motion corresponding to H = 0.5 and whose derivative is the white noise. The fBm is self-similar in distribution and the variance of the increments is given by

    Var(fBm(t)-fBm(s)) = v |t-s|^(2H),

    where v is a positive constant.

    This special form of the variance of the increments suggests various ways to estimate the parameter H. One can find in Bardet et al [2] a survey of such methods. The wfbmesti function provides three different estimates. The first one, due to Istas and Lang [4], is based on the discrete second-order derivative. The second one is a wavelet-based adaptation and has similar properties. The third one, proposed by Flandrin [3], estimates H using the slope of the loglog plot of the detail variance versus the level. A more recent extension can be found in Abry et al [1].


    [1] Abry, Patrice, Patrick Flandrin, Murad S. Taqqu, and Darryl Veitch. “Self-Similarity and Long-Range Dependence Through the Wavelet Lens.” In Theory and Applications of Long-Range Dependence, edited by P. Doukhan, G. Oppenheim, and M. Taqqu, 527–56. Birkhäuser, 2003.

    [2] Bardet, Jean-Marc, Gabriel Lang, Georges Oppenheim, Anne Philippe, Stilian Stoev, and Murad S. Taqqu. “Semi-Parametric Estimation of the Long-Range Dependence Parameter: A Survey.” In Theory and Applications of Long-Range Dependence, edited by P. Doukhan, G. Oppenheim, and M. Taqqu, 557–77. Birkhäuser, 2003.

    [3] Flandrin, Patrick. “Wavelet Analysis and Synthesis of Fractional Brownian Motion.” IEEE Transactions on Information Theory 38, no. 2 (March 1992): 910–17.

    [4] Istas, Jacques, and Gabriel Lang. “Quadratic Variations and Estimation of the Local Hölder Index of a Gaussian Process.” Annales de l’Institut Henri Poincare (B) Probability and Statistics 33, no. 4 (1997): 407–36.

    Version History

    Introduced before R2006a

    See Also