Espectro multifractal
El espectro de singularidad o espectro multifractal es una función usada en análisis multifractal para describir objetos que no tienen una dimensión fractal de Hausdorff bien definida, sino propiedades fractales que diferen de una escala a otra. Es decir, el espectro multifractal describe, para una medida o señal con heterogeneidad a múltiples escalas, la dimensión fractal del conjunto de puntos que comparten el mismo exponente de Hölder local . Intuitivamente, mide la “intensidad” o singularidad local, es decir, cómo se concentra la medida al acercarnos a un punto; mientras que indica “cuántos puntos” presentan dicha singularidad. La formulación moderna se consolidó en física estadística y sistemas dinámicos a mediados de los años 1980 y ha sido desarrollada rigurosamente en matemáticas.[1][2]
Matemáticamente es una función no nula sobre un intervalo abierto dentro del cual es diferenciable y cóncava.[2] Intuitivamente los valores para los cuales la función no se anula representan posibles valores de la dimensión fractal de subconjuntos del objeto multifractal (es decir, si bien no existe un único exponente o dimensión fractal bien definido para todo el conjunto, sí existen subconjuntos para los cuales la dimensión fractal o el exponente de Hölder sí están bien definidos).
Más formalmente, dado un objeto multifractal definido por una medida de Borel , el espectro multifractal es una función real definida como:
donde es la dimensión fractal local de la medida que describe el objeto multifractal y es la dimensión de Hausdorff de un conjunto de puntos. El concepto de espectro multifractal fue introducido por primera vez en 1986, en el trabajo pionero de Halsey et al.[1]
Definición formal
Si es un espacio métrico se denota como el conjunto de todas las medidas de probabilidad sobre . Nótese aquí un objeto multifractal se va a representar como el conjunto del espacio donde una medida de probabilidad no se anula, fuera del conjunto multifractal la medida de probabilidad se tomará como nula. En esas condiciones es cuando se puede definir rigurosamente el espectro multifractal asociado al objeto definido por la medida . Dado un punto se definen los exponentes fractales locales superior e inferior mediante:[2]
(1)
donde es la bola cerrada de radio centrada en .Cuando y coinciden, el valor común se denomina simplemente dimensión fractal local y se designa como . Ahora para cualquier se definen los conjuntos:
(2)
donde es el soporte topológico de la medida . En estas condiciones se puede pensar en la colección de conjuntos de la forma:
(3)
como una descomposición multifractal del objeto multifractal caracterizado por la medida de probabilidad y, entonces, el problema fundamental es caracterizar la dimensión de dichos conjuntos para cada valor de . Usualmente se usa la dimensión de Hausdorff, aunque se pueden usar otras, en estos términos el espectro multifractal viene dado por la función:
(4)
Formalismo práctico: funciones de partición y Legendre
Para estimar el espectro multifractal, se usa el formalismo de partición. Se recubre el soporte del objeto multifractal con cajas de tamaño y pesos y, entonces, se define:[1]
(5)
donde de los exponentes de escala generan las dimensiones :
(6)
y el espectro, mediant la transformada de Legendre:
(7)
En muchos casos relevantes, esta construcción coincide con la definición geométrica presentda más arriba.
Ejemplo: cascada binomial
Para la medida binomial en , con probabilidades y , el exponente de masas es:
(8)
donde y se obtienen analíticamente, produciendo un espectro cóncavo típico. Estas cascadas se usan como modelos de referencia y extensión en geociencias y campos análogos.[3] Este modelo produce un cóncavo; su anchura cuantifica cuánta multifractalidad existe.
Métodos prácticos de estimación
Existen diversos métodos para la estimación del espectro multifractal. Uno de ellos basados en series temporales no estacionarias o MFDFA (Multifractal Detrended Fluctuation Analysis) integra la señal, segmenta en ventanas, elimina tendencias locales y calcula funciones de fluctuación para determinar primero y, a partir de esa función, calcular . Permite distinguir multifractalidad por correlaciones frente a la causada por colas pesadas comparando con series barajadas.[4]
Otro método ampliamente usado en señales e imágenes de naturaleza multifractal es el método basado en ondículas (WTMM: Wavelet Transform Modulus Maxima) que define particiones sobre las crestas de máximos del módulo de ondícula. Es robusto frente a soportes irregulares y se usa para estimar directamente el espectro .[5]
Otro método simple comúnmente usado sn las particiones o box-counting multifractal, que se basa en obtener un ajuste lineal de respecto a para estimar y, vía la transformada Legendre, .[1]
Aplicaciones
El espectro multifractal se emplea para caracterizar intermitencia en turbulencia (por ejemplo, leyes de escala de She-Leveque), texturas e imágenes, finanzas, precipitación y sismología, entre otros.[6][7]
Referencias
- 1 2 3 4 Halsey, T.C.; Jensen, M.H.; Kadanoff, L.P.; Procaccia, I.; Shraiman, B.I. (1986). «Fractal measures and their singularities: The characterization of strange sets». Physical review A 33 (2): 1141. doi:10.1103/PhysRevA.33.1141.
- 1 2 3 Olsen, Lars (1995). «A Multifractal Formalism». Advances in Mathematics 116 (1): 82-196. doi:10.1006/aima.1995.1066.
- ↑ Cheng, Qiuming (2014). «Generalized binomial multiplicative cascade processes and their applications in geosciences». Nonlinear Processes in Geophysics 21: 477-487. doi:10.5194/npg-21-477-2014.
- ↑ Kantelhardt, Jan W.; Zschiegner, Sebastian A.; Koscielny-Bunde, Eva; Havlin, Shlomo; Bunde, Armin; Stanley, H. Eugene (2002). «Multifractal detrended fluctuation analysis of nonstationary time series». Physica A: Statistical Mechanics and its Applications 316 (1-4): 87-114. ISSN 0378-4371. doi:10.1016/S0378-4371(02)01383-3.
- ↑ Muzy, J.-F.; Bacry, E.; Arneodo, A. (1993). «Multifractal formalism for fractal signals: The structure-function approach versus the wavelet-transform modulus-maxima method». Physical Review E 47 (2): 875-884. doi:10.1103/PhysRevE.47.875.
- ↑ She, Z.-S.; Leveque, E. (1994). «Universal scaling laws in fully developed turbulence». Physical Review Letters 72 (3): 336-339. doi:10.1103/PhysRevLett.72.336.
- ↑ Arneodo, A.; Bacry, E.; Muzy, J.-F. (1995). «The thermodynamics of fractals revisited with wavelets». Physica A 213 (1–2): 232-275. doi:10.1016/0378-4371(94)00163-N.