|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| [[File:Karperien Strange Attractor 200.gif|thumb|300px|{{anchor|multifractal}}A [[Strange Attractor]] that exhibits [[multifractal]] scaling]]
| | Hi there. My title is Sophia Meagher even though it is not the name on my birth certification. To climb is some thing she would by no means give up. I am an invoicing officer and I'll be promoted soon. North Carolina is where we've been living for years and will never transfer.<br><br>Here is my homepage :: clairvoyants; [http://www.article-galaxy.com/profile.php?a=143251 related web-site], |
| [[File:WF111-Anderson transition-multifractal.jpeg|thumbnail|Example of a multifractal electronic eigenstate at the [[Anderson localization]] transition in a system with 1367631 atoms.]]
| |
| 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.<ref>{{cite book | last = Harte | first = David | title = Multifractals | publisher = Chapman & Hall | location = London | year = 2001 | isbn = 978-1-58488-154-4 }}</ref>
| |
| | |
| Multifractal systems are common in nature, especially [[geophysics]]. They include [[Turbulence|fully developed turbulence]], [[stock market]] time series, real world scenes, the Sun’s magnetic field time series, [[Cardiac cycle|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. The origin of multifractality in sequential (time series) data has been attributed, to mathematical convergence effects related to the [[central limit theorem]] that have as foci of convergence the family of statistical distributions known as the [[Tweedie distributions|'''Tweedie exponential dispersion models''']]<ref name=Kendal2011b>Kendal WS & Jørgensen BR (2011) Tweedie convergence: a mathematical basis for Taylor's power law, ''1/f'' noise and multifractality. ''Phys. Rev E'' 84, 066120</ref> as well as the geometric Tweedie models.<ref name=Jørgensen2011> Jørgensen B, Kokonendji CC (2011) Dispersion models for geometric sums. ''Braz J Probab Stat'' 25, 263-293</ref> The first convergence effect yields monofractal sequences and the second convergence effect is responsible for variation in the fractal dimension of the monofractal sequences.<ref name=Kendal2014>Kendal WS (2014) Multifractality attributed to dual central limit-lie convergence effects. ''Physica A'' 401, 22-33</ref>
| |
| | |
| From a practical perspective, multifractal analysis uses the mathematical basis of multifractal theory to investigate datasets, often in conjunction with other methods of [[fractal analysis]] and [[lacunarity]] analysis. The technique entails distorting datasets extracted from patterns to generate multifractal spectra that illustrate how scaling varies over the dataset. The [[#techniques|techniques of multifractal analysis]] have been applied in a variety of practical situations such as predicting earthquakes and interpreting medical images.<ref>{{cite pmid|19535282}}</ref><ref name="chhabra">{{cite doi|10.1186/1471-2164-12-506}}</ref><ref>{{cite pmid|22101185}}</ref>
| |
| | |
| == Definition ==
| |
| In a multifractal system <math>s</math>, the behavior around any point is described by a local [[power law]]:
| |
| | |
| :<math>s(\vec{x}+\vec{a})-s(\vec{x}) \sim a^{h(\vec{x})}.</math>
| |
| | |
| The exponent <math>h(\vec{x})</math> is called the [[singularity exponent]], as it describes the local degree of [[Mathematical singularity|singularity]] or [[regularity]]{{disambiguation needed|date=May 2012}} around the point <math>\vec{x}</math>.
| |
| | |
| 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 <math>s</math>.
| |
| | |
| In practice, the multifractal behaviour of a physical system <math>X</math> is not directly characterized by its singularity spectrum D(h). Data analysis rather gives access to the ''multiscaling exponents'' <math>\zeta(q),\ q\in{\mathbb R}</math>. Indeed, multifractal signals generally obey a ''scale invariance'' property which yields power law behaviours for multiresolution quantities depending on their scale <math>a</math>. Depending on the object under study, these multiresolution quantities, denoted by <math>T_X(a)</math> in the following, can be local averages in boxes of size <math>a</math>, gradients over distance <math>a</math>, wavelet coefficients at scale <math>a</math>... For multifractal objects, one usually observes a global power law scaling of the form:
| |
| | |
| :<math>\langle T_X(a)^q \rangle \sim a^{\zeta(q)}\ </math>
| |
| | |
| at least in some range of scales and for some range of orders <math>q</math>. When such a behaviour is observed, one talks of scale invariance, self-similarity or multiscaling.<ref>{{cite journal |author=A.J. Roberts and A. Cronin |title=Unbiased estimation of multi-fractal dimensions of finite data sets |journal=Physica A |volume=233 |year=1996 |pages=867–878 |doi=10.1016/S0378-4371(96)00165-3 }}</ref>
| |
| | |
| == Estimation ==
| |
| | |
| Using the so-called ''multifractal formalism'', it can be shown that, under some well-suited assumptions, there exists a correspondence between the singularity spectrum <math>D(h)</math> and the multiscaling exponents <math>\zeta(q)</math> through a [[Legendre transform]]. While the determination of <math>D(h)</math> calls for some exhaustive local analysis of the data, which would result in difficult and numerically unstable calculations, the estimation of the <math>\zeta(q)</math> relies on the use of statistical averages and linear regressions in log-log diagrams. Once the <math>\zeta(q)</math> are known, one can deduce an estimate of <math>D(h)</math> thanks to a simple Legendre transform.
| |
| | |
| Multifractal systems are often modeled by stochastic processes such as [[multiplicative cascade]]s. Interestingly, the <math>\zeta(q)</math> receives some statistical interpretation as they characterize the evolution of the distributions of the <math>T_X(a)</math> as <math>a</math> 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 ({{harvnb|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 [http://www.maths.adelaide.edu.au/anthony.roberts/multifractal.php]).
| |
| | |
| {{anchor|techniques}}
| |
| | |
| == Practical application of multifractal spectra ==
| |
| [[File:Distort.gif|thumb|Multifractal analysis is analogous to viewing a dataset through a series of distorting lenses to home in on differences in scaling. The pattern shown is a [[Hénon map]]]]
| |
| | |
| {{anchor|distort}}Multifractal analysis has been used in several fields in science to characterize various types of datasets.<ref>{{cite doi|10.1364/OE.20.003015}}</ref> In essence, multifractal analysis applies a distorting factor to datasets extracted from patterns, to compare how the data behave at each distortion. This is done using graphs known as '''multifractal spectra''' that illustrate how the distortions affect the data, analogous to viewing the dataset through a "distorting lens" as shown in the [[#distort|illustration]].<ref name="bcmf"/> Several types of multifractal spectra are used in practise.
| |
| | |
| === D<sub>Q</sub> vs Q ===
| |
| {{anchor|dqvsq}}[[File:Dqvsq.gif|thumb|D<sub>Q</sub> vs Q spectra for a non-fractal circle (empirical box counting dimension = 1.0), mono-fractal [[List of fractals by Hausdorff dimension#cross|Quadric Cross]] (empirical box counting dimension = 1.49), and multifractal [[Hénon map]] (empirical box counting dimension = 1.29).]]
| |
| | |
| {{anchor|dimensional ordering}}One practical multifractal spectrum is the graph of D<sub>Q</sub> vs Q, where D<sub>Q</sub> is the '''generalized dimension''' for a dataset and Q is an arbitrary set of exponents. The expression ''generalized dimension'' thus refers to a set of dimensions for a dataset (detailed calculations for determining the generalized dimension using [[box counting]] are described [[#generalized dimension|below]]).
| |
| | |
| ==== Dimensional ordering ====
| |
| The general pattern of the graph of D<sub>Q</sub> vs Q can be used to assess the scaling in a pattern. The graph is generally decreasing, sigmoidal around Q=0, where D<sub>(Q=0)</sub> ≥ D<sub>(Q=1)</sub> ≥ D<sub>(Q=2)</sub>. As illustrated in the [[#dqvsq|figure]], variation in this graphical spectrum can help distinguish patterns. The image shows D<sub>(Q)</sub> spectra from a multifractal analysis of binary images of non-, mono-, and multi-fractal sets. As is the case in the sample images, non- and mono-fractals tend to have flatter D<sub>(Q)</sub> spectra than multifractals.
| |
| | |
| The generalized dimension also offers some important specific information. D<sub>(Q=0)</sub> is equal to the [[Capacity Dimension]], which in the analysis shown in the figures here is the [[box counting dimension]]. D<sub>(Q=1)</sub> is equal to the [[Information Dimension]], and D<sub>(Q=2)</sub> to the [[Correlation Dimension]]. This relates to the "multi" in multifractal whereby multifractals have multiple dimensions in the D<sub>(Q)</sub> vs Q spectra but monofractals stay rather flat in that area.<ref name="bcmf"/><ref name="chaabra"/>
| |
| | |
| === <math>f(\alpha)</math> vs <math>\alpha</math> ===
| |
| Another useful multifractal spectrum is the graph of <math>f(\alpha)</math> vs <math>\alpha</math> (see [[#calculations|calculations]]). These graphs generally rise to a maximum that approximates the [[fractal dimension]] at Q=0, and then fall. Like D<sub>Q</sub> vs Q spectra, they also show typical patterns useful for comparing non-, mono-, and multi-fractal patterns. In particular, for these spectra, non- and mono-fractals converge on certain values, whereas the spectra from multifractal patterns are typically humped over a broader extent.
| |
| | |
| == Estimating multifractal scaling from box counting ==
| |
| {{anchor|calculations}}
| |
| Multifractal spectra can be determined from [[box counting]] on digital images. First, a box counting scan is done to determine how the pixels are distributed; then, this "mass distribution" becomes the basis for a series of calculations.<ref name="bcmf">{{citation | author=Karperien, A |title=What are Multifractals? | publisher=ImageJ | accessdate=2012-02-10|url=http://rsbweb.nih.gov/ij/plugins/fraclac/FLHelp/Multifractals.htm |year=2002 |archive=http://www.webcitation.org/65LENzkV8}}</ref><ref name="chaabra">{{cite doi|10.1103/PhysRevLett.62.1327}}</ref><ref>{{cite doi|10.2136/sssaj2001.6551361x}}</ref> The chief idea is that for multifractals, the probability, <math>P</math>, of a number of pixels, <math>m</math>, appearing in a box, <math>i</math>, varies as box size, <math>\textstyle\epsilon</math>, to some exponent, <math>\textstyle\alpha</math>, which changes over the image, as in {{EquationNote|Eq.0.0}}. ''NB: For mono[[fractals]], in contrast, the exponent does not change meaningfully over the set.'' <math>P</math> is calculated from the box counting pixel distribution as in {{EquationNote|Eq.2.0}}.
| |
|
| |
| {{NumBlk|:|<math>P_{[i,\epsilon]} \varpropto \epsilon^{-\alpha_i} \therefore\alpha_i \varpropto \frac{\log{P_{[i,\epsilon]}}}{\log{\epsilon^{-1}}}</math>|{{EquationRef|Eq.0.0}}}}
| |
| | |
| :<math>\epsilon</math> = an arbitrary scale ([[Box counting|box size]] in box counting) at which the set is examined
| |
| | |
| :<math>i</math> = the index for each box laid over the set for an <math>\epsilon</math>
| |
| | |
| :<math>m_{[i,\epsilon]}</math> = the number of pixels or ''mass'' in any box, <math>i</math>, at size <math>\epsilon</math>
| |
| | |
| :<math>N_\epsilon</math> = the total boxes that contained more than 0 pixels, for each <math>\epsilon</math>
| |
| | |
| {{NumBlk|:|<math>M_\epsilon = \sum_{i=1}^{N_\epsilon}m_{[i,\epsilon]} = </math> the total mass or sum of pixels in all boxes for this <math>\epsilon</math>|{{EquationRef|Eq.1.0}}}}
| |
| | |
| {{NumBlk|:|<math>P_{[i,\epsilon]} = \frac{m_{[i,\epsilon]}}{M_\epsilon} = </math> the probability of this mass at <math>i</math> relative to the total mass for a box size|{{EquationRef|Eq.2.0}}}}
| |
| | |
| <math>P</math> is used to observe how the pixel distribution behaves when distorted in certain ways as in {{EquationNote|Eq.3.0}} and {{EquationNote|Eq.3.1}}:
| |
| | |
| :<math>Q</math> = an arbitrary range of values to use as exponents for distorting the data set
| |
| | |
| {{NumBlk|:|<math>I_{{(Q)}_{[\epsilon]}} = \sum_{i=1}^{N_\epsilon} {P_{[i,\epsilon]}^Q} = </math> the sum of all mass probabilities distorted by being raised to this Q, for this box size|{{EquationRef |Eq.3.0}}}}
| |
| :*When <math>Q=1</math>, {{EquationNote|Eq.3.0}} equals 1, the usual sum of all probabilities, and when <math>Q=0</math>, every term is equal to 1, so the sum is equal to the number of boxes counted, <math>N_\epsilon</math>.
| |
| {{NumBlk|:|<math>\mu_{{(Q)}_{[i,\epsilon]}} = \frac{P_{[i,\epsilon]}^Q}{I_{{(Q)}_{[\epsilon]}}} = </math> how the distorted mass probability at a box compares to the distorted sum over all boxes at this box size|{{EquationRef|Eq.3.1}}}}
| |
| | |
| These distorting equations are further used to address how the set behaves when scaled or resolved or cut up into a series of <math>\epsilon</math>-sized pieces and distorted by Q, to find different values for the dimension of the set, as in the following:
| |
| | |
| :*An important feature of {{EquationNote|Eq.3.0}} is that it can also be seen to vary according to scale raised to the exponent <math>\textstyle\tau</math> in {{EquationNote|Eq.4.0}}:
| |
| {{NumBlk|:|<math>I_{{(Q)}_{[\epsilon]}} \varpropto \epsilon^{\tau_{(Q)}}</math>|{{EquationRef|Eq.4.0}}}}
| |
| | |
| Thus, a series of values for <math>\tau_{(Q)} </math> can be found from the slopes of the regression line for the log of {{EquationNote|Eq.3.0}} vs the log of <math>\epsilon</math> for each <math>Q</math>, based on {{EquationNote|Eq.4.1}}:
| |
| | |
| {{NumBlk|:|<math>\tau_{(Q)} = {\lim_{\epsilon\to0}{\left[ \frac {ln{I_{{(Q)}_{[\epsilon]}}}} {ln{\epsilon}} \right ]}} </math>|{{EquationRef|Eq.4.1}}}}
| |
|
| |
| :*{{anchor|generalized dimension}}For the generalized dimension:
| |
| {{NumBlk|:|<math>D_{(Q)} = {\lim_{\epsilon\to0} { \left [ \frac{ln{I_{{(Q)}_{[\epsilon]}}}}{ln{\epsilon^{-1}}} \right ]}} {(1-Q)^{-1}} </math>|{{EquationRef|Eq.5.0}}}}
| |
| {{NumBlk|:|<math>D_{(Q)} = \frac{\tau_{(Q)}}{Q-1}</math>|{{EquationRef|Eq.5.1}}}}
| |
| {{NumBlk|:|<math>\tau_{{(Q)}_{}} = D_{(Q)}\left(Q-1\right)</math>|{{EquationRef|Eq.5.2}}}}
| |
| {{NumBlk|:|<math>\tau_{(Q)} = \alpha_{(Q)}Q - f_{\left(\alpha_{(Q)}\right)}</math>|{{EquationRef|Eq.5.3}}}}
| |
| | |
| :*<math>\textstyle\alpha_{(Q)}</math> is estimated as the slope of the regression line for {{math|log A<sub><math>\epsilon</math>,Q</sub>}} vs {{math|log <math>\epsilon</math>}} where:
| |
| | |
| {{NumBlk|:|<math>A_{\epsilon,Q} = \sum_{i=1}^{N_\epsilon}{\mu_{{i,\epsilon}_{Q}}{P_{{i,\epsilon}_{Q}}}} </math>|{{EquationRef|Eq.6.0}}}}
| |
| :*Then <math>f_{\left(\alpha_{{(Q)}}\right)}</math> is found from {{EquationNote|Eq.5.3}}.
| |
| | |
| :*The mean <math>\textstyle\tau_{(Q)}</math> is estimated as the slope of the log-log regression line for <math>\textstyle\tau_{{(Q)}_{[\epsilon]}}</math> vs <math>\textstyle\epsilon</math>, where:
| |
| | |
| {{NumBlk|:|<math>\tau_{(Q)_{[\epsilon]}} = \frac{\sum_{i=1}^{N_\epsilon} {P_{[i,\epsilon]}^{Q-1}}} {N_\epsilon} </math>|{{EquationRef|Eq.6.1}}}}
| |
| | |
| In practise, the probability distribution depends on how the dataset is sampled, so optimizing algorithms have been developed to ensure adequate sampling.<ref name="bcmf"/>
| |
| | |
| == See also ==
| |
| * [[Multifractal Model of Asset Returns]] (MMAR)
| |
| * [[Multifractal Random Walk model]] (MRW)
| |
| * [[Fractional Brownian motion]]
| |
| * [[Mandelbrot cascade]], [[continuous cascade]] and [[lognormal cascade]]
| |
| * [[Detrended fluctuation analysis]]
| |
| * [[Tweedie distributions]]
| |
| | |
| ==References==
| |
| <references />
| |
| :11. Veneziano, D. and Essiam, A.K. (2003). Flow through porous media with multifractal hydraulic conductivity.
| |
| :Water Resources Research 39: doi: 10.1029/2001WR001018. issn: 0043-1397.
| |
| | |
| ==External links==
| |
| *{{cite journal |author=Stanley H.E., Meakin P. |title=Multifractal phenomena in physics and chemistry |journal=Nature |volume=335 |year=1988 |pages=405–9 |url=http://polymer.bu.edu/hes/articles/sm88.pdf |format=Review |doi=10.1038/335405a0 |issue=6189}}
| |
| | |
| *{{cite journal |author=Alain Arneodo, ''et al.'' |title=Wavelet-based multifractal analysis |journal=Scholarpedia |volume=3 |issue=3 |pages=4103 |year=2008 |url=http://www.scholarpedia.org/article/Wavelet-based_multifractal_analysis |doi=10.4249/scholarpedia.4103}}
| |
| | |
| * [http://www2.warwick.ac.uk/fac/sci/physics/research/theory/research/disqs/media Movies of visualizations of multifractals]
| |
| | |
| {{Fractals}}
| |
| | |
| {{DEFAULTSORT:Multifractal System}}
| |
| [[Category:Fractals]]
| |
| [[Category:Dimension theory]]
| |