The spectra and periodograms of anti-correlated discrete fractional Gaussian noise
Introduction
Recently, there has been considerable interest in the use of stochastic fractal models to help interpret physiological data [1]. Two models that have been investigated extensively in this context are fractional Brownian motion (FBM) and fractional Gaussian noise (FGN) [2]. These two models are related to one another, but the connection between them and how they are used in practical applications brings up some subtle issues that have not been fully appreciated amongst practitioners and that are the focus of the present manuscript.
Our starting point is FBM, which we denote here as BH(t), 0⩽t⩽∞. FBM is a continuous parameter stochastic process that depends upon a parameter H, the Hurst coefficient, where 0<H<1. The qualifier ‘continuous parameter’ refers to the fact that the independent variable, t, ranges over all non-negative real values. In practical applications, we must deal with sampled data, which leads us to consider BH(t) at just the integers This sampling leads to a discrete parameter version of FBM, which we refer to as dFBM. For clarity, we henceforth refer to the continuous parameter version of FBM as cFBM.
Flandrin [3] shows that, even though cFBM is a non-stationary process, it has a well-defined power spectrum S(f) that obeys a power law exactly over all frequencies; i.e., S(f) is proportional to f−1−2H for −∞<f<∞. Flandrin also considers a form of a derivative of cFBM, namely, a process defined as . He shows that this derivative process also has a spectrum that obeys a power law exactly over all frequencies and is given by S(f)∼f1−2H. The concept of FGN as formulated by Mandelbrot and Van Ness is related to this derivative process in that, rather than letting δ decrease to zero, we fix it at unity. This leads us to what we will refer to as continuous FGN (cFGN) and discrete FGN (dFGN). cFGN is defined as X(t)=BH(t+1)−BH(t), where 0⩽t⩽∞, resulting in a spectrum given by the product of the squared gain function for a first difference filter and the spectrum of cFBM, i.e., , which approaches a power law only as f→0. dFGN is obtained by restricting t to the non-negative integers. Note that, in addition to being the sampled version of cFGN, the process dFGN can be regarded as the first difference of dFBM. We also note that, due to the sampling, the spectra for dFBM and dFGN are even periodic functions with a period of unity, so the frequencies of interest satisfy −1/2⩽f⩽1/2, and their respective spectra also approach power laws only as f→0.
Both dFGN and dFBM have found application in many areas of science and have been applied to discretely sampled time series of many natural processes. The fact that cFBM and its derivative have spectra that exactly obey a power law over −∞<f<∞ has sometimes mistakenly been taken to hold both for the spectra of dFBM and dFGN over −1/2⩽f⩽1/2 and for estimates of these spectra given by the periodogram. Sometimes the power-law expression is cited directly, and it is implied or assumed that the expected estimated spectra as computed by the periodogram for dFGN are a power law everywhere [4], [5], [6]. The idea that the spectra of cFGN are a power law everywhere has been transferred to power-law behavior of dFGN spectral estimates, and the caveat that the actual spectra for dFGN are only a power law in the limit as f→0 appears to be unheeded. Churilla et al. [4] posit separate power-law behavior for the low frequencies and the high frequencies of the periodogram and derive separate Hurst coefficients for both by fitting each band of frequencies with a different power law based on the power-law formulation. The spectral synthesis method of Peitgen and Saupe [7] and that modified by Bassingthwaighte and Raymond [8] generated dFGN by an inverse Fourier transform of the power-law spectral coefficients after phase randomization and multiplying the amplitudes by random Gaussian numbers. The expected value of periodograms from this kind of process will be a power law, and does not have the same spectrum as a dFGN. Pilgram and Kaplan [5] used a spectral synthesis method whose basis is the power-law spectrum, but studied only spectra for H⩾0.5 and therefore did not observe the incorrect spectral representation at the low frequencies though their method does also represent the highest frequencies incorrectly.
If we plot S(f)=f1−2H versus f on a log–log scale, we see a line with a slope of 1−2H. To determine frequencies over which S(f) for dFGN varies as f1−2H, we define a local spectral exponentwhere 0<δ<f. As δ→0, . We also define a local Hurst coefficientBy plotting H(f,δ) as a function of f with δ set to a small number, it is possible to see at what frequencies this function is approximately equal to the Hurst coefficient, showing the region where S(f) varies as f1−2H to a good approximation. It is not widely appreciated how very low the frequency must become when H is also small for the spectra to exhibit power-law behavior.
Marked departure from nominal 1−2H behavior of dFGN spectra occurs under several conditions. At high frequencies near f=0.5 (the Nyquist frequency for the data set sampled from cFGN), the dFGN spectra flatten out and the local spectral exponent is zero. When H is near zero, a significant portion of the frequency band beginning at f=0.5 and extending through mid-range frequencies also deviates from nominal 1−2H power-law behavior. The only spectrum of dFGN that obeys the nominal power law everywhere is for H=0.5, i.e., the ‘white’ noise spectrum with zero exponent everywhere. For H>0.5, the spectra of correlated or persistent dFGN exhibit a significant discrepancy from nominal power-law behavior only near the highest frequencies. The flatness of these spectra near f=0.5 is less remarkable than the overall impression of power-law behavior [9, p. 54]. For H<0.5, the spectra of anti-correlated or anti-persistent dFGN have more extensive regions of non-power-law behavior than what is observed for positively correlated dFGN spectra. More attention has been paid to positively correlated than to anti-correlated series, and consequently the impression of power-law behavior for all spectra of dFGN has been inadvertently reinforced.
The next three sections of this paper are devoted to (1) calculating the local Hurst coefficients, H(f,δ) from Eq. (2) for the spectra of dFGN; (2) comparing the local H(f,δ) from spectra and periodograms and showing the effects of series length, N; and (3) evaluating the utility of the local spectral exponents, e(f,δ), obtained at the lowest available frequencies of the periodograms. The results emphasize that the 1−2H power-law behavior of the spectra applies only to a very small range of frequencies when the Hurst coefficient is small; the expected value of the periodogram exhibits power-law behavior only over a small frequency band, if at all; the expected values of the periodograms at the lowest available frequencies for anti-correlated dFGN do not converge to the actual spectra as N increases; and for very small H, the local spectral exponent of the smallest available frequencies from the periodogram of dFGN is zero.
Section snippets
Calculating the local Hurst coefficients, H(f,δ), for the spectra of dFGN
The spectral density function (SDF) for dFGN is given by Sinai [10] and Beran [9] aswhere σ2 is the variance of the process, 0<H<1 is the Hurst coefficient, andEq. (3) is not useful in calculating the SDF when H<0.5 because of the slow convergence of the infinite summation. The SDF can be approximated and computational efficiency gained by using an Euler–Maclaurin summation [11, p. 280]:
Comparing the local H(f,δ) from dFGN spectra and periodograms with H=0.001 to 0.1
Given a time series X0,X1,…,XN−1 of length N, a basic estimate of the SDF is the periodogram:The expected value of the periodogram is given by the convolution of Fejér's kernel,with the SDF:where S(f) is given by Eq. (3). The contribution to this expectation at a particular frequency from distant frequencies is called spectral leakage and occurs because of an interplay between the shape of Fejér's kernel and
The local spectral exponents, e(f,δ), for periodograms of varied H at the lowest available frequencies
In Fig. 4, for 10−4<H<1, the local spectral exponent, e(f,δ) is shown as the thick line for a set of series, all of length N=215, where f=(1+1/2)/N and δ=(1/2)/N. This choice of f and δ means that e(f,δ) depends on the expected value of the periodogram at its two lowest non-zero Fourier frequencies, namely 1/N and 2/N. The local spectral exponent goes to zero as H goes to zero. On a log–log plot of the periodogram versus frequency, the lowest frequency end behaves as f0 instead of f1−2H as H
Summary and conclusions
Although the SDF for dFGN approaches a power law when the frequency approaches zero, it is not everywhere a power law, except when the Hurst coefficient equals 0.5. For anti-correlated dFGN, the frequency where the SDF becomes a power law to a decent approximation depends on the Hurst coefficient. As H→0, the frequency at which the power-law region becomes discernible approaches zero as well. There is an obvious implication for methods generating dFGN. Attempting to generate dFGN by the Fourier
Acknowledgements
The research was supported by the National Simulation Resource for Circulatory Mass-Transport and Exchange via grant RR-1243 from the National Center for Research Resources of the National Institutes of Health.
References (13)
- J.B. Bassingthwaighte, L.S. Liebovitch, B.J. West, Fractal Physiology, Oxford University Press, New York, London, 1994,...
- et al.
Fractional Brownian motions, fractional noises and applications
SIAM Rev.
(1968) On the spectrum of fractional Brownian motions
IEEE Trans. Inform. Theory
(1989)- et al.
Membrane potential fluctuations of human T-lymphocytes have the fractal characteristics of fractional Brownian motion
Ann. Biomed. Eng.
(1996) - et al.
A comparison of estimators for 1/f noise
Physica D
(1998) - et al.
Establishing the relation between detrended fluctuation analysis and power spectral density analysis for stochastic processes
Phys. Rev. E
(2000)
Cited by (12)
Record length requirement of long-range dependent teletraffic
2017, Physica A: Statistical Mechanics and its ApplicationsCitation Excerpt :Lobato and Robinson [164], Sela and Clifford Hurvich [165], Reisen et al. [166], Sergides and Paparoditis [167], Velasco [168], Raymond et al. [169], Henry [170], Robinson and Marinucci [171], Moulines and Soulier [172], Kokoszka and Mikosch [173], Lobato [174], Bhansali [175], Reisen [176], Hurvich and Beltrao [177], Li and Zhao [178].
Tsallis non-extensive statistical mechanics of El Niño southern oscillation index
2007, Physica A: Statistical Mechanics and its ApplicationsMulti-Fractal Traffic and Anomaly Detection in Computer Communications
2023, Multi-Fractal Traffic and Anomaly Detection in Computer CommunicationsMandelbrot's Stochastic Time Series Models
2019, Earth and Space SciencePreliminary assessment of the scaling relationships of in-situ stress orientation variations indicated by wellbore failure data
2014, Rock Engineering and Rock Mechanics: Structures in and on Rock MassesPreliminary assessment of the scaling relationships of in-situ stress orientation variations indicated by wellbore failure data
2014, Rock Engineering and Rock Mechanics: Structures in and on Rock Masses - Proceedings of EUROCK 2014, ISRM European Regional Symposium