Como se pega uma onda?

September 10, 2017 | Autor: H.m. De Oliveira | Categoria: Higher Education, History of Electrical Engineering
Share Embed


Descrição do Produto

DEPARTAMENTO DE ELETRÔNICA E SISTEMAS

Entrando na Onda… (ou) Como se pega uma “onda”? Prof. Hélio Magalhães de Oliveira

CODEC Grupo de Pesquisas em Comunicações, Departamento de Eletrônica e Sistemas E-mail: [email protected] URL: http://www2.ee.ufpe.br/codec/WEBLET.html AULA DE ENCERRAMENTO DE ENGª ELÉTRICA / ELETRONICA

27/04/05

Uma teoria mais geral introduzida nos anos 80:

A Transformada de Wavelet ORIGEM- Escola Francesa (Morlet, Grossmann, Meyer, Battle, Lemarié, Cohen, Mallat, Coifman, Rioul, etc.) ... Pacotes de ondas acústicas sísmicas.

Gama de aplicações: geologia sísmica visão computacional e humana radar e sonar computação gráfica predição de terremotos e maremotos turbulência, fractais bancos de filtros espectrometria distinção celular (normais vs patológicas) modelos para trato auditivo

compressão de imagens (JPEG 2000) descontaminação de sinais (denoising) detecção de rupturas e bordas análise de seqüências de DNA análise de tons musicais neurofisiologia detecção de curtos eventos patológicos análise de sinais biomédicos (eletrocardiogramas, mamografias, eletroencefalogramas etc.)

espalhamento em banda larga modelagem de sistemas lineares óptica e eletromagnetismo modelagem geométrica

caracterização de sinais acústicos reconhecimento de alvos transitório e falhas em linhas de potência Metalurgia (rugosidade de superfícies) visualização volumétrica Telecomunicações previsão em mercados financeiros Estatística solução de equações diferenciais ordinárias & parciais ...................

Jean-Baptiste-Joseph Fourier 1822:

"La Théorie Analytique de la Chaleur":

Ondas senoidais como elementos de vibrações e ondas periódicas  verdadeiros átomos das flutuações e do fluxo. Os sinais passaram a ser analisados no domínio de Fourier, i.e., no domínio da freqüência. A decomposição evoluiu para a representação via transformada de Fourier.

Por que wavelets? Em que esta ferramenta pode ser mais potente que a análise espectral clássica de Fourier? Conceito de estacionaridade Uma das deficiências da análise de Fourier é que ela não apresenta um caráter local. Todo o sinal, desde o começo dos tempos (-∞) até o fim dos tempos (+∞) é levado em consideração.

(a)

(b)

WAVELETS Década de 80, Alex Grossmann (Université de Marseille) e Jean P. Morlet (Elf Acquitaine) introduziram o conceito de wavelets

Jean Morlet.

(ecce homo!)

Em 1985, Stéphane Mallat (França) estabeleceu a ligação desta teoria com o processamento digital de sinais. Yves Meyer (França) construiu uma das primeiras Wavelets não triviais, continuamente diferenciáveis.

Ingrid Daubechies (Bélgica) construiu o mais usado conjunto de wavelets ortogonais de suporte compacto (tempo-limitada).

Wavelets Contínuas ψ(t) wavelet-mãe

ψ(t)∈L (ℜ). 2

Operações: 1

t

a) escalonamento ψ a (t ) = | a | ψ  a  , a≠0.

b) deslocamento ψ b (t ) = ψ (t − b) .

{ψ (t )} → {ψ a,b (t )} (∀a, a≠0) (∀b∈ℜ).

versão: Comprimida da wavelet mãe, se a1. Define-se +∞

* CWT(a,b):= ∫− ∞ f (t )ψ a ,b (t ) dt =< f(t),ψa,b>.

Analogia com Fourier: F(w)=



+∞

−∞

f (t )e − jwt dt =< f(t),ejwt >.

Esta nova "base" é composta por sinais oscilatórios de curta duração - e não sinais ab aeterno.

Figura. A Wavelet-mãe Symmlet 8 em diferentes escalas e localizações. Uma Transformada inversa de Wavelet pode ser obtida via f (t ) =

1 Cψ

t − b  dadb  1 CWT ( a , b ) ψ ( ) 2  ∫−∞ ∫−∞ a  a .  a +∞ +∞

A Transformada Inversa (CWT-1) e a Condição de Admissibilidade Sejam f(t)∈ L2(ℜ) e ψa,b(t) ∈ L2(ℜ-{0}× ℜ). Sob que condições é possível "pegar uma onda"? (catch the wave...)

Escolher

uma

wavelet

protótipo

tal

que

ψ(t) ↔ Ψ(w)

| Ψ( ζ ) |2 Cψ = ∫ dζ < +∞ −∞ . |ζ | +∞

Exemplos: Um mar de Wavelets Existe um grande número de funções que podem ser eleitas como wavelets mãe.

e

Nome da família de Wavelets 'haar'

Haar wavelet.

'db'

Daubechies wavelets.

'sym'

Symlets.

'coif'

Coiflets.

'bior'

Biorthogonal wavelets.

'rbio'

Reverse biorthogonal wavelets.

'meyr' Meyer wavelet. 'dmey' Discrete approximation of Meyer wavelet. 'gaus'

Gaussian wavelets.

'mexh' Mexican hat wavelet. 'morl'

Morlet wavelet.

'cgau'

Complex Gaussian wavelets.

'shan'

Shannon wavelets.

'deO'

de Oliveira wavelets.

'legd'

Legendre wavelets.

'mth'

Mathieu wavelets.

'Cheb'

Chebyshev wavelets.

'Gegb'

Gegenbauer wavelets.

'beta'

Beta wavelets.

'fbsp'

Frequency B-Spline wavelets.

'cmor' Complex Morlet wavelets. 

Wavelet complexa de Morlet Em sua investigação de sinais geofísicos (exploração de petróleo), Morlet empregou a wavelet complexa:

ψ

( Mor )

(t ) =

Figura. Wavelet complexa de Morlet.

1

π

1/ 4

e

− t 2 / 2 − jw0 t

e

.

Wavelet de Shannon A análise correspondente aos filtros passa-faixa ideais define uma decomposição usando wavelets conhecidas como wavelets de Shannon. No caso da wavelet complexa, pode-se usar ψ

(a)

( CSha )

(t ) = Sinc(t )e − j 2πt .

(b)

Figura. (a) Wavelet de Shannon. (b) A wavelet complexa de Shannon.

Wavelets de Daubechies

Um dos maiores desafios da teoria de wavelets foi a construção de uma família de wavelets ortogonais de suporte compacto.

Figura. Wavelets dbN de Daubechies (N=2,3,4, ...): As Daublets Matlab.

Wavelet de "de Oliveira" Nova família de wavelets ortogonais complexas. Possuem espectro típico passa-faixa ideal (plano), com regiões de "rolamento" assimétricas. 0.435

0.5

0.4

φ( t , 0.1) 0.3 φ( t , 0.2)

 1 φ t ,   3

0.2

0.1

0 − 0.087 0.1

4 −5

alpha=0.1 alpha=0.2 alpha=0.3

2

0 t

2

4 5

Figura. Função escala de "de Oliveira". (esboço para α=0,1, 0,2 e 0,3).

Mostra-se que o módulo wavelet complexa de "de Oliveira" é dado por 0  1 1 (w − π (1 + α ) ) cos  4 α 2 π  1  S ( deO ) ( w) =  2π   1 cos 1 (w − 2π (1 − α ) )  2π 8α 0 

se w < π (1 - α ) se π (1 - α ) < w < π (1 + α ) se π (1 + α ) < w < 2π (1 − α ) se 2π (1 − α ) < w < 2π (1 + α ) se w > 2π (1 + α )

Figura. Módulo da Wavelet de "de Oliveira"

Figura. Wavelet ψ ( deO) (t ) : Parte real

Figura. Wavelet ψ ( deO) (t ) : parte imaginária.

Compressão via wavelets A título de destacar a relevância das wavelets como mecanismos de compressão:

 Inclusão no padrão internacional JPEG 2000.  padrão do FBI (Federal Bureau armazenamento de impressões digitais.

of

Investigations)

para

o

Padrão WSQ, "Padrão de Compressão Wavelet/Quantização escalar" (do Inglês: Wavelet/Scalar Quantization standard).

As imagens são produzidas com qualidade de arquivamento e com taxas de compressão em torno de 20:1.

Figura (a) Imagem original (261 KB); (b) Imagem comprimida por JPEG (28 KB); (c) Imagem comprimida wavelets (5KB).

WAVELET SHIFT-KEYING: A NEW DIGITAL MODULATION H.M. DE OLIVEIRA, H.A.N. SILVA, E.A. BOUTON, Wavelet Shift-Keying: A New Digital Modulation, Simpósio Brasileiro de Telecomunicações, XX SBrT, Setembro 2003, Rio de Janeiro.

Preliminary ideas of a new digital modulation scheme termed wavelet shift keying (WSK), which is based on discrete wavelet transforms.

The WSK Modulation M-ary WSK signalling schemes can be used in data transmission application. • Scaled versions of the mother wavelet are transmitted in each symbol slot. • It is assumed that scale factor of every slot depends on the input binary data.

LOCALIZAÇÃO DE FALTAS EM LINHAS DE TRANSMISSÃO UTILIZANDO A REPRESENTAÇÃO WAVELET MULTIRESOLUÇÃO L.R. SOARES, M.A.CARVALHO JR, H.M. DE OLIVEIRA, Detecção e Classificação de Faltas Utilizando a Transformada de Wavelets, Seminário Nacional de Produção e Transmissão de Energia Elétrica, XVII SNPTEE, Uberlândia, 2003. L.R. SOARES, M. A. DE CARVALHO JR, H.M. DE OLIVEIRA, Localização de Faltas em Linhas de Transmissão Utilizando a Representação Wavelet Multirresolução, Congresso Brasileiro de Automática-CBA, Natal, Setembro, 2002.

LT: Presidente Dutra/Boa Esperança (PDD/BEA), Norte/Nordeste, 205,6km.

As

faltas

foram

simuladas

no

ATP

(Alternative

Transient

Program),

considerando que o monitoramento da linha é realizado no terminal de PDD.

• Localização da falta: método da impedância aparente - versões aproximadas dos sinais de tensão e corrente à terceira escala. • Simulações resultaram erros inferiores a 3,7%: satisfatório para reduzir a área de inspeção por parte das equipes de manutenção.

COMPACTLY SUPPORTED WAVELETS DERIVED FROM LEGENDRE POLYNOMIALS: SPHERICAL HARMONIC WAVELETS M.M.S. LIRA, H.M. DE OLIVEIRA, M.A. CARVALHO Jr, R.M.C. DE SOUZA, Compactly Supported Wavelets Derived from Legendre Polynomials: Spherical Harmonic Wavelets, WSEAS Int. Conf. on Systems, Greece, 2003.

2-dimensional Legendre transform (2D-legdN) for image analysis: Interesting effects in multiresolution decomposition of hand-drawn images.

Fig. - Reconstruction of a hand-drawn draft using 2D-legd2. (a) Original (b) Reconstructed image from 1-level wavelet decomposition.

ELLIPTIC-CYLINDER WAVELETS: THE MATHIEU WAVELETS M.M.S. LIRA, H.M. de OLIVEIRA, R.J. de SOBRAL CINTRA, Elliptic-Cylindrical Wavelets: The Mathieu Wavelets, IEEE Signal Processing Letters, , vol.11, n.1, January, pp. 52 - 55, 2004.. M.M.S. Lira, H.M. de Oliveira, R.J. de Sobral Cintra, R.M.C. de Souza, Wavelets for Elliptical Waveguide Problems, 2002 WSEAS International Conference on Wavelet Analysis and Multirate Systems, Vouliagmeni, Athens, Greece, December 29-31, 2002.

Broad scope of wave-guide problems involving elliptical geometry, including: i) analysis for weak guiding for step index elliptical core optical fibres ii) power transport of elliptical wave-guides iii) evaluating radiated waves of elliptical horn antennas

Figure - Mathieu Wavelet with ν=1 and q=5 for 2,4 and 6 iteration.

ON WAVELET DECOMPOSITION OVER FINITE FIELDS DE OLIVEIRA, H. M., FALK, T. H. On Wavelet Decomposition over Finite Fields In: 19º Simpósio Brasileiro de Telecomunicações, 2001, Fortaleza CE. Anais do 19º Simpósio Brasileiro de Telecomunicações, 2001. (artigo selecionado para publicação do n. especial da Rev. da Soc. Bras. de Telecom.) TAVORA, R. G. F., SOUZA, R. M. C., DE OLIVEIRA, H. M. Um Algoritmo Rápido para a Transformada Wavelet em Corpos Finitos In: 19 Simpósio Brasileiro de Telecomunicações, 2001, Fortaleza. Anais do 19 Simpósio Brasileiro de Telecomunicações , 2001. DE OLIVEIRA, H.M., FALK T.H., TÁVORA R.F.G., Decomposição de Wavelets sobre Corpos Finitos, Revista da Sociedade Brasileira de Telecomunicações. Campinas, SP: v.17, n.1, p.38 - 47, 2002.

•Foundations of wavelets over Galois fields. •Standard orthogonal finite-field wavelets including FF-Haar and FF-Daubechies. •FF-wavelets to design spread-spectrum sequences.

APLICAÇÃO: PROJETO DE SEQÜÊNCIAS DE ESPALHAMENTO ESPECTRAL As portadoras wavelet digitais têm a mesma duração T de um símbolo de entrada a ser modulado, tendo N chips por símbolo de dados. Nova classe de esquemas CDM/CDMA baseado em transformadas de wavelet sobre corpos finitos:

Figura: Sistema Multiplex utilizando wavelets FF- Haar.

Análise wavelet para Seqüências de DNA

As FF-Wavelets podem ser mais adequadas! (pesquisa)

CHEBYSHEV WAVELETS R.J.S. CINTRA, L.R. SOARES, H.M. DE OLIVEIRA, On Filter Banks and Wavelets Based on Chebyshev Polynomials, submitted to WSEAS Int. Conf. on Systems, Greece, 2003.

New family of wavelets, named Chebyshev wavelets, which are derived from conventional Chebyshev polynomials. The convergence of the cascade algorithm of 2nd kind Chebyshev wavelets is proved by using properties of Markov chains.

Aplicações: Wavelets Polinomiais em Eletrogastrografia Acometimentos patológicos do sistema digestivo, ligados ao estômago (e.g. gastroparesia, gastrites, úlceras pépticas, constipações idiopáticas, refluxo gástrico-esofágico simples, neoplasias gástricas, taquigastrias, doenças gástrico-congênitas).

Figura. Trecho de 20 segundos de um sinal de EGG amostrado a uma freqüência de 2 Hz Sinal Típico de EGG e a função wavelet de Chebyshev.

CONCLUSÕES • RESUMO DOS PRINCIPAIS DESENVOLVIMENTOS NA UFPE

Novas wavelets Wavelets em corpos finitos Wavelet de "de Oliveira" Wavelets de Chebyshev wavelets de Gegenbauer Wavelets de Legendre Wavelets de Mathieu Wavelets Beta

Potenciais aplicações Multiplex, Acesso múltiplo, Análise de seqüências de DNA Modulação digital, Modems ADSL Análise em biomédica (eletrogastrografia, ECG...) Desenhos, Tomografia, RMN, Imagens médicas Óptica, eletromagnetismo (antenas...) Sinais com formato tipo quase senoidal

• IMPORTANTE E PODEROSA FERRAMENTA NA ANÁLISE DE SINAIS

• NUMEROSAS APLICAÇÕES E PERSPECTIVAS Sucesso profissional para todos os caros novos colegas!

Lihat lebih banyak...

Comentários

Copyright © 2017 DADOSPDF Inc.