Multifractal system

A multifractal system is a generalization of a fractal system in which a single exponent (the fractal dimension) is not enough to describe its dynamics; instead, a continuous spectrum of exponents (the so-called singularity spectrum) is needed.

Multifractal systems are common in nature, especially geophysics. They include fully developed turbulence, stock market time series, real world scenes, the Sun’s magnetic field time series, heartbeat dynamics, human gait, and natural luminosity time series. Models have been proposed in various contexts ranging from turbulence in fluid dynamics to internet traffic, finance, image modeling, texture synthesis, meteorology, geophysics and more.

In a multifractal system s, the behavior around any point is described by a local power law:

s(\vec{x}%2B\vec{a})-s(\vec{x}) \sim a^{h(\vec{x})}.

The exponent h(\vec{x}) is called the singularity exponent, as it describes the local degree of singularity or regularity around the point \vec{x}.

The ensemble formed by all the points that share the same singularity exponent is called the singularity manifold of exponent h, and is a fractal set of fractal dimension D(h). The curve D(h) versus h is called the singularity spectrum and fully describes the (statistical) distribution of the variable s.

In practice, the multifractal behaviour of a physical system X is not directly characterized by its singularity specrum D(h). Data analysis rather gives access to the multiscaling exponents \zeta(q),\ q\in{\mathbb R}. Indeed, multifractal signals generally obey a scale invariance property which yields power law behaviours for multiresolution quantities depending on their scale a. Depending on the object under study, these multiresolution quantities, denoted by T_X(a) in the following, can be local averages in boxes of size a, gradients over distance a, wavelet coefficients at scale a... For multifractal objects, one usually observes a global power law scaling of the form:

<T_X(a)^q> \sim a^{\zeta(q)}\

at least in some range of scales and for some range of orders q. When such a behaviour is observed, one talks of scale invariance, self-similarity or multiscaling.

Contents

Estimation

Thanks to the so-called multifractal formalism, it can be shown that, under some well-suited assumptions, there exists a correspondence between the singularity spectrum D(h) and the multiscaling exponents \zeta(q) through a Legendre transform. While the determination of D(h) calls for some exhaustive local analysis of the data, which would result difficult and numerically unstable, the estimation of the \zeta(q) relies on the use of statistical averages and linear regressions in log-log diagrams. Once the \zeta(q) are known, one can deduce an estimate of D(h) thanks to a simple Legendre transform.

Multifractal systems are often modeled by stochastic processes such as multiplicative cascades. Interestingly, the \zeta(q) receives some statistical interpretation as they characterize the evolution of the distributions of the T_X(a) as a goes from larger to smaller scales. This evolution is often called statistical intermittency and betrays a departure from Gaussian models.

Modelling as a multiplicative cascade also leads to estimation of multifractal properties for relatively small datasets (Roberts & Cronin 1996). A maximum likelihood fit of a multiplicative cascade to the dataset not only estimates the complete spectrum, but also gives reasonable estimates of the errors (see the web service [1]).

See also

References

External links