|
[01 - 02 - 03 - 04 - 05 - 06 - 07 - 08 - 09 - 10 - 11 - 12 - 13] Aleatorización de objetos fractales deterministas - Percolación: ejemplos, renormalización - El modelo de Eden - Agregación limitada por difusión (DLA) - Caracterización fractal de series temporales - Movimiento browniano: autosemejanza estadística, autoafinidad, dimensión box-counting, categorías, simulación - Densidad espectral de Fourier - Procesos 1/f - Dimensión fractal basada en la densidad espectral - Síntesis de fBm basada en la tranformada de Fourier - Extensión del fBm a dimensiones superiores - Método de los desplazamientos en 2d - Costas fractales a partir de paisajes fractales - Fourier filtering en d>1 - Mapas de color - Desplazamiento aleatorio del punto medio no-gaussiano - Series temporales en mercados financieros - Integración de paisajes - Análisis de texturas |
||
|
En 1900 Louis Bachelier (1870-1946) finalizó su tesis doctoral titulada Théorie de la spéculation, donde estableció el primer modelo y las primeras aplicaciones para el movimiento browniano en el contexto de las fluctuaciones de precios en los mercados financieros. Para representar procesos con dependencia a largo plazo (LRD), con aplicación en series temporales de contenido financiero, los modelos basados en fBm no resuelven todos los problemas del modelo, debido, en particular, a que se trata de modelos gaussianos. En 1963 B. Mandelbrot publicó The variation of certain speculatice prices, The Journal of Bussines, 36, donde retomando el trabajo de Bachelier, señalaba, por ejemplo, que los cambios de precio de magnitud elevada, ocurren mucho mas frecuentemente que los que predice una ley gaussiana. Los procesos estables, de los que nos vamos a ocupar seguidamente, se han propuesto para la descripción de porcentajes en cambios de precios en mercados financieros. Estos modelos, resuelven, en particular, el problema apuntado de los procesos gaussianos aunque, como veremo mas adelante, presentan otras carencias. |
||
| PROCESOS ESTABLES |
Consideremos dos variables aleatorias X1, X2, con distribuciones N(m1,s1) y N(m2,s2), respectivamente. Es bien conocido que la variable aleatoria X1+X2, tiene distribución N(m1+m2, (s1)2+(s2)2). Esta propiedad que consiste en la invariancia dentro de una clase de distribuciones, frente al operador suma, se denomina estabilidad. Adoptando otro enfoque, dadas dos variables aleatorias X1, X2, la ecuación funcional s1X1+s2X2=sX, junto con la relación subsidiaria (s1)a+(s2)a=sa, admite la gaussiana como solución para a=2, pero admite también otras soluciones, todas las distribuciones estables en el sentido de Lévy. Una distribución es estable y simétrica de clase a, si su función característica se puede escribir en la forma |
|
|
|
||
|
donde 0<a£2, c³0 y m es la media. Si a=1, se tiene la distribución de Cauchy. Si a=2, se tiene la normal. Si a<<2, los momentos de orden >a no existen. Más generalmente, un parámetro adicional, b, -1£b£1, permite prescindir de la simetría: |
||
![]() |
||
|
Como consecuencia de lo dicho anteriormente, una densidad de probabilidad se denomina estable, si es invariante bajo una operación de convolución consigo misma. La transformada de Fourier de un producto de convolución de dos N(0,1) es una N(1/sqrt(2),0), de forma que recuperamos la estabilidad de Lévy de la distribución normal. En general, una distribución Lévy a-Estable, se caracteriza por Xa(s,b,m), donde 0<a£2, s es un parámetro de escala, b un parámetro que controla la simetría (1£b£1, b=0, distribución simétrica, b=1 o 1, fuerte asimetría) y m el valor medio. Como ya hemos apuntado, el valor a=2 se da para las gaussianas. |
||
| ALGORITMO DE SIMULACIÓN | Presentamos ahora un
algoritmo bien conocido para simular una distribución Xa(1,0,0).
En primer lugar, producimos muestras de una variable aleatoria V
uniformemente distribuida en (-p/2,p/2)
y de una variable aleatoria exponencial W, independiente de la primera,
con media 1. Seguidamente, calculamos |
|
![]() |
||
|
Distribución a-estable (distaest.m) |
||
| INDEPEN- DENCIA DE LA ESCALA |
Como ya hemos indicado, si a<2, los momentos de orden >a no existen. En particular los momentos de orden dos divergen. Para el fBm, el momento de orden uno determina la escala de tiempo, y el de orden dos la escala física de la distribución de las longitudes de los saltos. La divergencia del momento de orden dos en ciertas distribuciones de Lévy, no es una contradicción para nosotros: esta ausencia de escala conduce a comportamientos independientes de la escala, fractales. Para valores grandes de x, la función de distribución de una variable aleatoria estable, satisface la relación de escala |
|
|
|
||
| Aunque cabe la tentación y, de hecho se hace así en ocasiones, el exponente a no tiene el significado de una dimensión fractal. | ||
| EL VUELO DE LÉVY TRUNCADO |
Para ciertos modelos de dinámica de precios en mercados financieros, sin embargo, este comportamiento en cola gruesa, fat tail no se compadece con los resultados empíricos, que presentan varianza finita. Así, Mantegna y Stanley propusieron un modelo, el vuelo de Lévy truncado, que reúne comportamientos mas acordes con aquellos. El TLF se define como un proceso estocástico X(t) caracterizado por la siguiente función de densidad: |
|
![]() |
||
|
donde L es la distribución estable Lévy simétrica de parámetro a y factor de escala s>0, N una constante de normalización y xcut la abscisa de corte. Este modelo explica las siguientes observaciones empíricas:
|
||