12.07.2015 Views

Cálculo del Exponente de Hurst mediante Wavelets

Cálculo del Exponente de Hurst mediante Wavelets

Cálculo del Exponente de Hurst mediante Wavelets

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Estimación <strong><strong>de</strong>l</strong> Exponent <strong>de</strong> <strong>Hurst</strong> usando <strong>de</strong>nsidad espectral WaveletComo he señalado anteriormente, la estimación wavelet <strong><strong>de</strong>l</strong> exponente <strong>de</strong> <strong>Hurst</strong> es lo que meinició en mi aventura exponente <strong>de</strong> <strong>Hurst</strong>. Mis primeros intentos <strong>de</strong> utilizar la transformadawavelet para calcular el exponente <strong>de</strong> <strong>Hurst</strong> fracasaron, por una variedad <strong>de</strong> razones. Ahora queha sido cubierto el método clásico R/S , es el momento <strong>de</strong> discutir los métodos wavelet.El exponente <strong>de</strong> <strong>Hurst</strong> para un conjunto <strong>de</strong> datos se calcula a partir <strong>de</strong> la <strong>de</strong>nsidad espectralwavelet, que se refiere a veces como escalograma. He cubierto técnicas espectrales wavelet enuna página web Spectral Analysis and Filtering with the Wavelet Transform. Esta página web<strong>de</strong>scribe la estructura wavelet "octava", que se aplica en el cálculo <strong><strong>de</strong>l</strong> exponente <strong>de</strong> <strong>Hurst</strong>. Elplot <strong>de</strong> la <strong>de</strong>nsidad espectral wavelet se genera a partir <strong><strong>de</strong>l</strong> espectro <strong>de</strong> potencia wavelet. Laecuación para calcular la potencia normalizada para la j-ésima octava se muestra en la Ecuación14.Aquí, la potencia se calcula a partir <strong>de</strong> la suma <strong>de</strong> los cuadrados <strong>de</strong> los coeficientes wavelet (elresultado <strong>de</strong> la transformada <strong>de</strong> wavelet hacia a<strong><strong>de</strong>l</strong>ante) para la j-ésima octava. Una octavawavelet contiene 2 j coeficientes wavelet. La suma <strong>de</strong> los cuadrados se normaliza dividiendopor 2 j , dando la potencia normalizada. En el análisis espectral, no siempre es necesario utilizarla potencia normalizada. En la práctica, el cálculo <strong><strong>de</strong>l</strong> exponente <strong>de</strong> <strong>Hurst</strong> requiere potencianormalizada.La Figura 17 muestra la gráfica <strong>de</strong> la <strong>de</strong>nsidad <strong>de</strong> potencia normalizada para la transformadawavelet <strong>de</strong> 1024 puntos <strong>de</strong> datos con un exponente <strong>de</strong> <strong>Hurst</strong> <strong>de</strong> 0.72. En este caso se utilizó elDaubechies D4 wavelet. Puesto que este número <strong>de</strong> octava wavelet es el log 2 <strong><strong>de</strong>l</strong> número <strong>de</strong>coeficientes wavelet, el eje x está en una escala logarítmica.Figura 17El exponente <strong>de</strong> <strong>Hurst</strong> se calcula a partir <strong>de</strong> la <strong>de</strong>nsidad espectral wavelet <strong>mediante</strong> el cálculo <strong>de</strong>una línea <strong>de</strong> regresión lineal a través <strong><strong>de</strong>l</strong> conjunto <strong>de</strong> puntos {x j , y j }, don<strong>de</strong> x j es la octava e y j esel log 2 <strong>de</strong> la potencia normalizada. La pendiente <strong>de</strong> esta línea <strong>de</strong> regresión es proporcional a laestimación para el exponente <strong>de</strong> <strong>Hurst</strong>. Un gráfico <strong>de</strong> regresión para el espectro <strong>de</strong> potenciaDaubechies D4 en la Figura 17 se muestra en la Figura 18.


Figura 18Si <strong>de</strong>sea una respuesta real, en lugar <strong>de</strong> una i<strong>de</strong>a <strong>de</strong> cómo el cambio <strong>de</strong> una cosa causa el cambio<strong>de</strong> otra cosa, la afirmación <strong>de</strong> que algo es proporcional a otra cosa no es muy útil. La ecuaciónpara calcular el exponente <strong>de</strong> <strong>Hurst</strong> partir <strong>de</strong> la pendiente <strong>de</strong> la línea <strong>de</strong> regresión a través <strong>de</strong> la<strong>de</strong>nsidad espectral normalizada se muestra en la Ecuación 15.La Tabla 3 compara la estimación exponente <strong>Hurst</strong> con tres transformada wavelet para elexponente <strong>de</strong> <strong>Hurst</strong> estimado <strong>mediante</strong> el método R/S. La Daubechies D4 es una transformadawavelet "energía normalizada". Se utilizan aquí también formas <strong>de</strong> energía normalizadas <strong><strong>de</strong>l</strong>wavelet <strong>de</strong> Haar e interpolación lineal. Un wavelet <strong>de</strong> energía normalizada parece ser requeridoen la estimación <strong><strong>de</strong>l</strong> exponente <strong>de</strong> <strong>Hurst</strong>.Table 3 Estimación <strong><strong>de</strong>l</strong> exponente <strong>de</strong> <strong>Hurst</strong> para 1024 puntos <strong>de</strong> conjuntos <strong>de</strong> datos sintéticosFunción Wavelet H=0.5 Error H=0.72 Error H=0.8 ErrorHaar 0.5602 0.0339 0.6961 0.0650 0.6959 0.1079Linear Interp. 0.5319 0.1396 0.8170 0.0449 0.9203 0.0587Daubechies D4 0.5006 0.0510 0.7431 0.0379 0.8331 0.0745R/S 0.5791 0.0193 0.7246 0.0149 0.5973 0.0170El error <strong>de</strong> la regresión lineal varía mucho para las diversas funciones wavelet. Aunque elwavelet interpolación lineal parece ser un buen filtro para las formas <strong>de</strong> onda <strong>de</strong> diente <strong>de</strong> sierra,como series <strong>de</strong> tiempo financieras, que no parece funcionar bien para estimar el exponente <strong>de</strong><strong>Hurst</strong>.No he <strong>de</strong>scubierto una manera <strong>de</strong> <strong>de</strong>terminar a priori si una función wavelet dada será un buenestimador <strong><strong>de</strong>l</strong> exponente <strong>de</strong> <strong>Hurst</strong>, aunque experimentalmente se pue<strong>de</strong> <strong>de</strong>terminar a partir <strong><strong>de</strong>l</strong>error <strong>de</strong> regresión. En algunos casos también <strong>de</strong>pen<strong>de</strong> <strong><strong>de</strong>l</strong> conjunto <strong>de</strong> datos. En la Tabla 4 elexponente <strong>de</strong> <strong>Hurst</strong> se estima a partir <strong>de</strong> la gráfica <strong>de</strong> <strong>de</strong>nsidad espectral <strong>de</strong> Haar y a través <strong>de</strong> latécnica <strong>de</strong> R/S. En este caso, el wavelet <strong>de</strong> Haar tenía un error <strong>de</strong> regresión más bajo que elwavelet Daubechies D4.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!