29
Transformada de Fourier: fundamentos matemáticos, implementação e aplicações musicais MAC 0337 – Computação Musical Jorge H. Neyra-Araoz IME – USP 22/11/2007

fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

Embed Size (px)

Citation preview

Page 1: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

Transformada de Fourier: fundamentos matemáticos, implementação e aplicações

musicais

MAC 0337 – Computação Musical

Jorge H. Neyra-AraozIME – USP

22/11/2007

Page 2: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

2

Resumo

● Série de Fourier para funções periódicas● Transformada de Fourier: definição e

exemplo● Transformada de Fourier Discreta

– Tranformada Rápida de Fourier

Page 3: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

3

Série de Fourier para funções periódicas

● A Série de Fourier é uma representação de uma função periódica como uma soma de funções periódicas:

Page 4: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

4

Série complexa de Fourier

● Segundo a Fórmula de Euler, as séries podem ser expressadas equivalemente

Page 5: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

5

Equações de Síntese e Análise

● Equação de Síntese:

● Equação de Análise:

● Uma propriedade fundamental dos valores de é a seguinte:

Page 6: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

6

Exemplo ...

Considere a função

Temos . Para

Page 7: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

7

Exemplo ...(2)

Page 8: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

8

Exemplo ...(3)

● O n-ésimo harmônico de f é dado por:

Page 9: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

9

Transformada de Fourier

● Quando uma função não é periódica é impossível escrevê-la como combinação linear de uma família de senos e cossenos harmonicamente relacionados.

● No entanto, muitas vezes é possível escrevê-la como combinação linear de todos os senos e cossenos que existem, utilizando todas as freqüencias disponíveis:

Page 10: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

10

Transformada de Fourier (2)

● Para funções que satisfazem o critério

entre outras, pode-se provar que os valores deque satisfazem a equação acima são

dados por

Page 11: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

11

Transformada Inversa de Fourier

● Se F(ω) e f(t) estão relacionadas pelas equações de análise e síntese

denotamos esta relação por

indicando que F é a Transformada de Fourier de f, e f é a Transformada Inversa de Fourier de F.

Page 12: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

12

Transformada de Fourier: Exemplo

● Considere a função

Page 13: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

13

Transformada de Fourier: Exemplo

● Para , temos que:

● e para , temos que:

Page 14: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

14

Transformada de Fourier: Exemplo

● Assim, pela equação de síntese,

● Freqüentemente, representamos a transforma de Fourier através dos espectro de amplitude e fase.

Page 15: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

15

Transformada de Fourier: Exemplo

● A escolha entre na expressão deembora arbitrária, deve preservar a propriedade de que é uma função ímpar.

Page 16: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

16

Transformada de Fourier Discreta (DFT)

● Dado considere a família de números complexos que divide o círculo unitário em N partes iguais:

Page 17: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

17

DFT: Relações de Ortogonalidade

● A família de número complexos anterior possui uma propriedade de ortogonalidade em relação ao produto interno emdefinido por

● Vamos utilizar o símbolo de tal forma que

Lema (Relações de Ortogonalidade)

Page 18: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

18

DFT: Exemplo

● Considere . Então, para os vetores

temos

Page 19: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

19

Transformada Rápida de Fourier (FFT)

● O método de cálculo da transformada discreta de Fourier a partir da expressão

utiliza produtos entre números complexos e somas, possuindo assim complexidade computacional

● O método FFT permite obter o mesmo resultado em tempo

Page 20: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

20

FFT: Implementação

● Existem muitas implementações:– Mapeamento em Índices Multidimensionais,

transforma uma DFT de uma dimensão em uma DFT de duas ou mais dimensões.

– Algoritmo COOLEY-TURKEY: usa o mapeamento de índices os índices de tempo e de freqüência.

– FFT por Dizimição no Tempo, descompõe a seqüência de tamanho N em seqüências sucessivas menores.

Page 21: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

21

FFT: Implementação

● FFT por dizimição em Freqüência– A mais común e a mais eficiente forma de

FFT.– Usa todas as dimensões do mesmo tamanho.– Calcula uma DFT de tamanho N, em termos

de duas DFTs de tamanho N/2.

Page 22: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

22

FFT: ilustração de três etapas

Page 23: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

23

FFT: Exemplo

● Consider o vetorSua transformada de Fourier discreta poder ser calculada com auxílio do diagrama do FFT do slide anterior.

Page 24: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

24

FFT: Algoritmo

Page 25: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

25

Janelas

● Também chamadas Funções de Ponderação.

● Modificam as características de resposta em freqüência dos algoritmos de DFT.

● Usadas em conjunto com análise espectral via DFT para reduzir o espalhamento em freqüência.

Page 26: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

26

Janela Rectangular e Triangular

Page 27: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

27

Aplicações da FT

● A Transfomada de Fourier (FT) é uma ferramenta largamente empregada em processamento de sinais (sons, imagens).

● A teoria de Fourier demonstrou que um som periódico (como os sons instrumentais utilizados na música) podia decompor-se numa soma de sons puros (sinusoidais) de frequência 1f, 2f, 3f,...

Page 28: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

28

Aplicações da FT

● Usado amplamente no análise, síntese, e processamento de voz e sons musicais.

● Pode ser usada para obter o Espectro a partir da Forma de Onda.

Page 29: fundamentos matemáticos, implementação e aplicações ...kon/MAC5900/seminarios/seminario_Jorge.pdf · 2 Resumo Série de Fourier para funções periódicas Transformada de Fourier:

29

Obrigado