introducción a wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm analisis...

84
Introducción a Wavelets (ondeletas) Introducción a Wavelets (ondeletas) http://www.jhu.edu/~signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Upload: ascencion-barquero

Post on 23-Jan-2016

213 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Introducción a Wavelets (ondeletas)Introducción a Wavelets (ondeletas)

http://www.jhu.edu/~signals/phasorlecture2/indexphasorlect2.htm

ANALISIS MULTIRESOLUCION

 

  

Page 2: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Transformadas espectrales:Transformadas espectrales:

F uN

f x ejux

N

x

N

( ) ( )

12

0

1

u = 0,1,2, ..., N-1

nucleo

xfuxKuF

)(),()(

Transformada de Fourier:

Page 3: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

xuNj

euxK2

),(

PARTE REAL PARTE IMAGINARIA

Page 4: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Transformada Hartley:Transformada Hartley:

Nucleo o KernelNucleo o Kernel

)(22

cos1

)(1

0

xfxuN

senxuNN

uFN

x

Page 5: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Transformada Discreta Coseno Transformada Discreta Coseno (DCT)(DCT)

)(2

)12(cos

1)(

1

0

xfN

ux

NuF

N

x

Page 6: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

nucleo

xfuxKuF

)(),()(

=

Page 7: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

nucleo

xfuxKuF

)(),()(

=

Page 8: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 9: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Espectro de FourierEspectro de Fourier

f(x)

F(u)

Page 10: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 11: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 12: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 13: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 14: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

STFT (Short time Fourier transform)STFT (Short time Fourier transform)

Or windowed Fourier Or windowed Fourier transformtransform

dtertgtfwrS jwtf

)()(),(

g(t-r) : Window

Page 15: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 16: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

SpectrogramSpectrogram

The square modulus of the windowed Fourier The square modulus of the windowed Fourier transform is the transform is the spectrogramspectrogram of a signal: of a signal:

22

)()(),(),( dtertgtfwrSwrSpec jwtff

Page 17: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 18: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Ventanas:Ventanas:

Ventana de Hamming:Ventana rectangular:

Page 19: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

A segment of a vowel extracted with a rectangular window The amplitude spectrum using a

rectangular window Calculated using Matlab: abs(fft(sig))

Page 20: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

  

  The amplitude spectrum using a hamming window. Calculated using Matlab: abs(fft(hamming(512) .* sig))

A segment of a vowel extracted with a hamming window. Calculated using Matlab: hamming(512) .* sig

Page 21: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

This is the basis for most computer generated spectrograms (display pixel intensity on a log scale by limiting the dynamic range to about 60-80 dB).

Page 22: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Ejemplos de espectrogramasEjemplos de espectrogramas

(Ver MATLAB)(Ver MATLAB)

Page 23: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Ejemplos de espectrogramasEjemplos de espectrogramas

Here is the sum of two parallel linear chirps with its spectrogram.

Page 24: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Here is the sum of two hyperbolic chirps and its spectrogram.

Page 25: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Introducción a WaveletsIntroducción a Wavelets

Page 26: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

four frequency components at different times. The interval 0 to 250 ms is a sinusoid of 300 Hz, and the other 250 ms intervals are sinusoids of 200 Hz, 100 Hz, and 50 Hz

w(t)=exp(-a*(t^2)/2)

Page 27: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 28: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 29: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 30: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 31: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 32: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 33: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 34: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 35: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 36: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 37: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 38: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 39: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 40: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Wavelet de MorletWavelet de Morlet

Page 41: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

The mexican hat The mexican hat waveletwavelet

224

1 2

)1(3

2)(

x

exx

Page 42: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

                                 

            

  Gráficos de varios tipos distintos de wavelets. (a) Wavelet de Haar, (b) Wavelet de Daubechies, (c) Wavelet de Morlet. (Cortesía de Ofer Levi, Universidad de Stanford)

Page 43: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

EscalaEscala

fs

1

Page 44: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

waveletmothert)(

Page 45: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 46: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 47: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 48: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 49: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Ejemplos de escalogramas (CWT)Ejemplos de escalogramas (CWT)Continuous Wavelet TransformContinuous Wavelet Transform

These signals are drawn from a database signals that includes event related potentials of normal people, and patients with Alzheimer's disease.

Page 50: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 51: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

En un espectrograma:En un espectrograma:

Page 52: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

En un escalograma:En un escalograma:

Page 53: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Suma de dos señales CHIRP hiperbólicas

Windowed fourier transform

(Espectrograma)

Continuous Wavelet Transform CWT

(Escalograma)

Page 54: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

AplicacionesAplicaciones

Page 55: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

AnAnálisis de señalesálisis de señales

OxOxímetro de Pulsoímetro de Pulso

Page 56: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 57: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 58: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Representación frecuencia-tiempo para:

(a)Datos muestreados

(b) FT

(c) WFT

(d) CWT

Page 59: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 60: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 61: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 62: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 63: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

DWT (Discrete Wavelet DWT (Discrete Wavelet Transform)Transform)

Page 64: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

A friendly guide to wavelets:A friendly guide to wavelets:

http://perso.orange.fr/polyvalens/clemens/wavelets/wavelets.html#section7

Page 65: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Ahora dejamos fija la Ondeleta y lo que vamos comprimiendo por etapas es la señal

Page 66: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

El análisis multiresolución se El análisis multiresolución se consigue a través de filtrado y consigue a través de filtrado y submuestreo de la señal original.submuestreo de la señal original.

La exploración en tiempo se La exploración en tiempo se consigue a través de operaciones consigue a través de operaciones de convolución (filtrado digital).de convolución (filtrado digital).

Page 67: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

)(nxW

)(nx

Sub-band coding

Page 68: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Sub-band coding algorithmSub-band coding algorithm

Page 69: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 70: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Transformada inversaTransformada inversa

Page 71: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 72: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

2-D Discrete Wavelet Transform2-D Discrete Wavelet Transform

Page 73: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 74: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 75: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 76: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 77: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 78: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 79: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

Wavelet PacketWavelet Packet

Page 80: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION
Page 81: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

En resumen:En resumen:

CWT

DWT

2D - DWT

Page 82: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

... En resumen:... En resumen:

CWT

DWT

2D - DWT

Page 83: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

http://www.gisdevelopment.net/technology/ic/techip0003a.htm

http://www.acm.org/crossroads/xrds6-3/sahaimgcoding.html#Fig6

COMPRESION DE LA DCT A WAVELETS

Page 84: Introducción a Wavelets (ondeletas) signals/phasorlecture2/indexphasorlect2.htm ANALISIS MULTIRESOLUCION

EspectrogramaEspectrograma