Lição
A teoria — Transformada Discreta de Fourier
A transformada discreta de Fourier de N amostras x0 … xN−1 é constituída por N números Xk = Σn xn e−2πikn/N, um por banda k. Cada número mede quanto do sinal é composto por uma onda que completa k ciclos nas N amostras. A transformada trata as amostras como um período de um sinal que se repete. Deste modo, é exata para tons que completam ciclos inteiros e espalha qualquer outro tom por todas as bandas: a fuga espectral. Uma janela wn multiplica primeiro as amostras para suavizar a junção.
O que significa cada símbolo
A_k- a amplitude que a banda k implica, o que as barras mostram: 2|Xk|/(N·G) para 0 < k < N/2, e |Xk|/(N·G) em k = 0 e k = N/2, os quais não possuem banda espelho para partilhar. Um tom de amplitude 1 numa banda regista 1
G- o ganho coerente da janela, a média de wn: 1 sem janela, ½ para a de Hann
- Pressupõe
- Cada tom começa no seu pico, como um cosseno. Os valores da cauda distante variam com a fase inicial. A banda 31 do tom de 10,5 ciclos encontra-se 26,3 dB abaixo da barra mais alta para um cosseno e 31,2 dB abaixo para um seno, ao passo que os 81% nas duas bandas mais próximas quase não se movem. A janela de Hann utiliza a forma periódica, ½ − ½cos(2πn/N), cuja transformada ocupa exatamente três bandas.
Problema resolvido na íntegra
-
Corrigir uma leitura e eliminar a necessidade de o fazer 4 passos
Um analisador recolhe amostras 8000 vezes por segundo e transforma 64 amostras de cada vez, sem usar qualquer janela. Se lhe fornecer um tom de teste de 1312,5 Hz, a barra mais alta registará 0,32. Qual é a verdadeira amplitude do tom e quantas amostras por bloco colocariam o tom exatamente numa banda? Prima As setas curvam-se: a página reproduz um tom com essa frequência.
Abrir este problema na ferramenta
-
As bandas estão separadas por 8000/64, logo o tom completa 10,5 ciclos por bloco: a meio caminho entre as bandas 10 e 11.
-
A meio caminho, a barra mais alta da página regista 64,61% do tom. Basta dividir:
-
O tom cai numa banda quando um bloco contém um número inteiro dos seus ciclos. Com 128 amostras as bandas ficam separadas por 62,5 Hz:
-
Se o analisador tivesse aplicado uma janela de Hann às 64 amostras, a barra mais alta teria registado 84,89% do tom. Assim, o mesmo tom mostraria:
Resposta
Cerca de 0,495 e 128 amostras por bloco. Qualquer múltiplo de 128 também funciona, visto que 1312,5 Hz necessitam de 128 amostras para completar um número inteiro de ciclos (21). A página está fixada em 64, por isso o caso de 128 amostras ultrapassa o que consegue mostrar. Com 128 amostras a barra lê a totalidade dos 0,495, e a janela deixaria de ser necessária para ler a amplitude, servindo apenas para afastar a cauda de um tom forte de um outro mais fraco.
-
Percurso de aprendizagem
A transformada de Fourier na prática
Referências (2)
- The standard reference on windows. Leakage defined as the page shows it: of all frequencies “only those which coincide with the basis will project onto a single basis vector; all other frequencies will exhibit non zero projections on the entire basis set”, because “the periodic extension of a signal not commensurate with its natural period exhibits discontinuities at the boundaries of the observation” (p. 52). Table I (p. 55) gives coherent gain 1.00 and 0.50 and scalloping loss 3.92 dB and 1.42 dB for the rectangle and Hann windows: F. J. Harris, “On the use of windows for harmonic analysis with the discrete Fourier transform”, Proceedings of the IEEE, vol. 66, no. 1, pp. 51–83 (1978).
- The DFT's definition, and what happens to a frequency that is not a bin: the sum “is N if ωₖ = ωₓ and zero at ωₗ, for l ≠ k. However, the sum is nonzero at all other frequencies”: J. O. Smith III, Mathematics of the Discrete Fourier Transform (DFT), with Audio Applications, 2nd ed., W3K Publishing (2007), “DFT Definition” and “Frequencies in the ‘Cracks’”.