Teoria de Los Wavelets
Teoria de Los Wavelets
Teoria de Los Wavelets
TEORA DE WAVELETS
Para nuestro estudio diremos que la Teora de Wavelets es el conjunto formado por: las
wavelets, la Transformada Continua de Wavelets y la Transformada Discreta de
Wavelets. De este modo nos referiremos de aqu en adelante a lo que es y lo que
comprende sta Teora la cual es una herramienta matemtica reciente, basada en la
Teora de Conjuntos y representaciones cuadradas integrables, que permiten presentar
una seal o un campo tanto en escala, espacio y posibles direcciones [JAL00]. En este
Captulo se explicar lo que son las wavelets dando algunos ejemplos de ellas, se
definirn las Transformadas Continua y Discreta de Wavelets y se expondrn las
principales ventajas que tiene sobre las Transformadas de Fourier.
2.1
Introduccin
(5)
18
(6)
2.2
2 B,
+
+ , + !
l+l
(7)
y son generadas a partir de funciones madre 2B. A esa funcin madre se le agregan un
par de variables que son la escala (+) que permite hacer dilataciones y contracciones de la
seal y la variable de traslacin (,), que nos permite mover a la seal en el tiempo. Estas
variables son nmeros reales y obviamente para una escala de 0 la funcin se
indetermina. Las wavelets se usan en el procesamiento de seales por sus caractersticas,
Jalali las describe a continuacin:
"They have the best possible simultaneous concentration in time and in
frequency, the set of their finite linear combinations is closed under Fourier
Transforms, pointwise multiplication and convolution. The scalar product of
any two members of this set is given by an explicit formula. And they are among
the very few classes of functions where the transition from one to more
dimensions is immediate." [JAL00].
19
Existen diferentes wavelets que ya son utilizadas de forma constante y que tienen
definiciones establecidas. Sin embargo, la eleccin de un tipo de wavelet depende de la
aplicacin especfica que se le vaya a dar. Actualmente existen muchas aplicaciones en
las que las wavelets actan de manera directa, de acuerdo con Kobayashi [KOB99] las
patentes en Estados Unidos durante el periodo de 1991 a 1998 bajo la palabra wavelet se
incrementaron de 3 a 92, lo que nos permite tener una idea de su utilidad y lo que en un
futuro se puede esperar de ellas. Por otro lado la investigacin en este campo se ha
incrementado notablemente como lo vemos en la figura 4.
2.3
Momentos de desvanecimiento
Este concepto es muy importante, ya que nos permite conocer la forma de la wavelet y es
un parmetro para saber que tan hbil es la wavelet para suprimir un polinomio dado. La
suavidad de la wavelet esta limitada por el nmero de momentos de desvanecimiento que
tenga. El 3-simo momento de la wavelet se calculan con la siguiente integral [POL96]:
(
2BB3 .B !
(8)
7!
(9)
20
2.4
Ejemplos de wavelets
2B "
!
"
#
"
#
B "
otro valor.
(10)
21
#" B# / #
7/B2B
1 $
(11)
22
23
. B#
/ 8
.B
(12)
donde G8 es una constante determinada por # 89<7 de 1+?=B 8 ". Con esta
wavelet slo es posible realizar la Transformada Continua de wavelets y puede ser
simtrica o asimtrica segn el valor de 8. La figura 10 muestra la wavelet Gaussiana de
orden 1.
79<6B / # -9=&B
(13)
24
En la figura 11 podemos ver que la wavelet Morlet es simtrica. Por otro lado
debemos comentar que no posee caractersticas de ortogonalidad ni biortogonalidad,
adems de que slo es til para realizar la Transformada Continua de wavelets. Existen
muchas otras wavelets en la literatura entre ellas podemos encontrar la Meyer, Spline y
Biorthogonal, para mayores detalles se puede consultar [JAL00], [KOB99], [MAT00] y
[TRE00].
2.5
Transformadas de wavelets
25
"... the wavelet transform is a tool that cuts up data or functions or operators
into different frequency components, and then studies each component with a
resolution matched to its scale" [JAL00].
"
B,
2
(
0 aBb.B + , + !
l+l
+
(14)
0
0
( 2 B , 0 aBb.B + , + !
09
0!
(15)
26
27
(a)
(b)
Figura 13. Proceso de anlisis en diferente escala y traslacin.
(a) escala fija a 1 y traslacin variante, (b) escala fija a 20 y traslacin variante.
Ahora veremos un ejemplo de lo que es la Transformada de Continua de
Wavelets, Polikar en [POL96] nos sugiere una seal mostrada en la figura 14. Podemos
ver que las componentes de frecuencia que tiene son de 30,20,10 y 5Hz. La transformada
de la seal se muestra en la figura 15. Los ejes que se manejan son traslacin y escala.
Es importante aclarar que el proceso de transformacin de una seal es reversible,
a esto es a lo que en la literatura se le conoce como la Trasformada Inversa Continua de
Wavelets (ICWT). Esta transformacin inversa esta dada por la siguiente ecuacin:
28
G[ X + ,2 B,
"
+
0 B # ( (
., .+
+#
GG
(16)
donde GG# es una constante que se determina por la wavelet que se ha utilizado en la
transformacin. Esta constante se conoce como constante de admisibilidad y est definida
como:
GG #1 (
lL=l#
.= _
l=l
_
(17)
29
Nyquist, los coeficientes en escala +" en 009" 3/=37+ a la razn de los coeficientes en la
escala +! . Para wavelets discretas los parmetros de escala y traslacin son elegidos tal
que en el nivel 4 la wavelet +!4 2+!4 B es +!4 veces el ancho de 2B Esto significa que el
parmetro de escala es + +!4 4 y el parmetro de traslacin , 5,! +!4 4 5
[ALA03] De este modo la familia de wavelets esta dada por:
4
(18)
(19)
#
""l.45 l# Fll0 Bll#
EF 4 5
(20)
#
"".45 245 B
EF 4 5
(21)
30
(22)
7_
31
Gracias a eso se obtiene la informacin de la seal, tal y como se hizo con las
seales continuas. Este proceso se puede explicar mejor viendo la figura 17. En ella
podemos apreciar como ingresa la seal a analizar a un par de filtros, uno de ellos es pasa
bajas y el otro es pasa altas, los cuales se identifican por la especificacin de ancho de
banda de la seal a la salida de cada uno, la cual va de la mitad del ancho de banda de la
seal de entrada a uno de los extremos. Por ejemplo, tenemos que la seal que entra va de
! a 1. El filtro pasa altas es representado por 18, mientras que el pasa bajas es 28.
La salida de ambos filtros se vuelve a enviar a otro par de filtros de las mismas
caractersticas. De ste modo se va reduciendo el ancho de banda de la seal y eso se
traduce en la reduccin a la mitad de la resolucin. Lo que significa que a mayor nmero
de etapas de filtrado se tendr una mayor resolucin.
La salida de los filtros pasa altas se van eliminando por la regla de Nyquist que
dice que para poder reconstruir una seal a partir de sus muestras es necesarios
muestrearla al menos con el doble de la frecuencia de esa seal, por lo que la informacin
en la salida de los filtros pasa altas no es necesaria. Conforme se van agregando las
etapas de filtrado se va aumentando el nivel N de descomposicin de la seal.
32
245 B # # 2#4 B 5
(23)
_
_
" si 4 4w y 8 8w
! cualquier otro caso
(24)
33
333
a4 B8 V4 B#8 V4"
3@
a5 B8 V! B8 5 V!
33
(25)
4N 5
B8245 #4 8 5.8
(26)
B8945 #4 8 5.8
(27)
34
(28)
(29)
C69A 5 "B82 8 #5
(30)
donde C2312 5 y C69A 5 son las salidas de los filtros pasa altas y pasa bajas
respectivamente. De este modo y cumpliendo las caractersticas requeridas se puede
analizar una seal continua que previamente se muestre para hacerla discreta. Y por
supuesto que el proceso es reversible mediante la Transformada Inversa Discreta de
Wavelets de la ecuacin (25) y donde haciendo algunos despejes se obtiene:
B8 " C2312 51 8 #5 C69A 52 8 #5
_
(31)
5_
la cual se realiza por medio de las salidas de los filtros pasa altas y pasabajas multiplicada
por su respuesta al impulso considerando el proceso de decimado y undecimado. Estos
resultados son sumados desde la primer muestra hasta la ltima y ya con eso se ha
reconstruido la seal discreta en el tiempo. Este proceso se muestra en la figura 19.
2.6
35
Aplicaciones
Los procesos que se han mencionado anteriormente son las principales herramientas que
se utilizan en el procesamiento y anlisis de seales, pero quizs se est preguntando para
que se utilizan estas Transformaciones en la vida real. La respuesta a esta pregunta es lo
que se explicar en sta seccin.
La Teora de Wavelets tiene muchas aplicaciones reales que comprenden la
deteccin de discontinuidades y puntos de ruptura en las seales, la identificacin de
frecuencias puras, la reduccin de ruido en seales, la compresin de seales,
aproximacin de funciones, mtodos espectrales para resolver ecuaciones diferenciales,
anlisis de fluidos turbulentos, entre otros [JAL00]. La que ms nos interesa estudiar es,
lgicamente, aquella aplicacin con la que es posible reducir el ruido de las seales de
inters. Sin embargo, mencionaremos brevemente cmo es que se consiguen algunas
otras aplicaciones de las que se han mencionado anteriormente.
La que resulta muy fcil de entender a simple vista es la identificacin de
frecuencias puras, pues como hemos visto la Transformada de wavelets (en sus 2
versiones) nos permite conocer las frecuencias que tienen la seal analizada, de este
modo y con un tratamiento de extraccin de la informacin adecuado para evitar que sta
informacin se corrompa es muy fcil identificar estas frecuencias puras.
La siguiente es la deteccin de discontinuidades y puntos de ruptura en las
seales, es posible utilizar las wavelets para estos propsitos pues sus procesos de
anlisis de multiresolucin nos permite conocer en que momento se presentan los
cambios de frecuencias en la seal analizada, de ste modo podemos detectar algn punto
en el que la seal tenga un cambio que puede ser un punto de ruptura o una
discontinuidad.
La compresin de seales es, quizs, la aplicacin ms famosa de la Teora de
wavelets, dado que el FBI (Federal Bureau of Investigation) la ha utilizado para guardar
las imgenes de las huellas digitales en sus archivos[JAL00]. Esta aplicacin es posible
gracias a los coeficientes de wavelet que es informacin que sin ser la imagen misma
posee los recursos suficientes, como hemos visto antes, para sintetizar la seal original,
de este modo el total de informacin a guardar se reduce considerablemente, ya que no se
36
guarda la imagen como tal, sino los coeficientes de wavelet que en comparacin son
notablemente menos informacin.
Despus continan otro tipo de aplicaciones que son propiamente los objetivos de
estas aplicaciones y donde no se usan directamente los mtodos aqu estudiados, pero si
cooperan con la obtencin del resultado final. Entre estas se encuentran las mencionadas
anteriormente, mtodos espectrales para resolver ecuaciones diferenciales y anlisis de
fluidos turbulentos. Para mayor detalle se pueden consultar las siguientes fuentes
[KOB99] y [BGG98]. Es importante aclarar que no son todas las aplicaciones que existen
para las wavelets, esto es slo una prueba muy pequea de lo que se puede conseguir
aplicando las Transformadas o las Wavelets por si solas. Tal y como lo marca la historia,
el lmite es nuestra imaginacin.