Transformadas de Fourier: mudanças entre as edições

De MediaWiki do Campus São José
Ir para navegação Ir para pesquisar
Moecke (discussão | contribs)
Moecke (discussão | contribs)
Linha 187: Linha 187:
=== Transformada Rápida de Fourier (FFT) ===
=== Transformada Rápida de Fourier (FFT) ===
As FFTs são algoritmos que reduzem a complexidade do cálculo da DFT, obtendo ordens de <math> \mathrm{ O(log N \times N)} \ </math>.
As FFTs são algoritmos que reduzem a complexidade do cálculo da DFT, obtendo ordens de <math> \mathrm{ O(log N \times N)} \ </math>.
*Para relembrar os conceitos vistos acima, recomendo assistir os vídeos sobre DFT do prof. Steve Brunton da University of Washington
{| class="wikitable"
|-
! The Discrete Fourier Transform (DFT)
! Computing the DFT Matrix
! The Fast Fourier Transform (FFT)
! The Fast Fourier Transform Algorithm
|-
|{{#ev:youtube|nl9TZanwbBk|300}}
|{{#ev:youtube|Xw4voABxU5c|300}}
|{{#ev:youtube|E8HeD-MUrjY|300}}
|{{#ev:youtube|toj_IoCQE-4|300}}
|}


== Soma dos termos de uma Progressão Geométrica ==
== Soma dos termos de uma Progressão Geométrica ==

Edição das 23h38min de 18 de novembro de 2020

1 Transformada de Fourier no tempo contínuo (TFTC)

O sinal x(t) é contínuo no tempo, e o sinal X(jΩ) é contínuo na frequência.

A equação de análise
É uma transformação de um domínio de uma variável real x(t) de tempo continuo em uma variável complexa X(jΩ) de frequência contínua.
 DTDF.
x:.
X(jΩ){x(t)} =defx(t) ejΩtdt
A equação de síntese
É uma transformação de um domínio de uma variável complexa X(jΩ) de frequência contínua em uma variável real x(t) de tempo continuo.
 DFDT.
X:.
x(t)1{X(jΩ)} =def12πX(jΩ) ejΩtdΩ
A transformada de Fourier X(jΩ), por simplificação é muitas vezes representada apenas por X(Ω) ou X(ω), mas é importante lembrar que tratasse de uma frequência complexa, representada pelo eixo imaginário do plano s.

Figura 1 - Sinal x(t) e sua transformada de Fourier X(Ω)
Fonte: Elaborado pelo autor.
  • Para relembrar os conceitos vistos em Sinais e Sistemas, recomendo assitir os vídeos sobre Transformada de Fourier do prof. Luis Antonio Aguirre da UFMG
Introdução a Transformada de Fourier Existência da Transformada Fourier Transformada de Fourier: Propriedades 1 Transformada de Fourier: Propriedades 2 Transformada de Fourier para Sinais Periodicos

2 Transformada Z

A transformada Z de um sinal de tempo discreto x(n) é a função X(z) definida como

X(z)=Z{x(n)}=n=x(n)zn 
onde n é um inteiro; z=Aejϕ  é um número complexo, com A sendo sua magnitude, e ϕ sua frequência angular (em radianos por amostra).

A transformada Z inversa é

x(n)=Z1{X(z)}=12πjCX(z)zn1dz 

onde C  é o caminho antihorário envolvendo a origem, dentro da região de convergencia (ROC) de X(z).

3 Transformada de Fourier de tempo discreto (TFTD)

O sinal x(n) é discreto no tempo, e o sinal X(Ω) é contínuo e periódico em 2π.

A equação de análise
É uma transformação de um domínio de uma variável real x(n) de tempo discreto em uma variável complexa X(ejω) frequência contínua periódica.
 DTDF.
x:.
X(ejω){x(n)} =defn=x(n) ejωn
A equação de síntese
É uma transformação de um domínio de uma variável complexa X(ejω) de frequência contínua periódica em uma variável real x(n) continua.
 DFDT.
X:.
x(n)1{X(ejω)} =def12πX(ejω) ejωndω

Como a transformada de Fourier X(ejω) é periódica com período 2π, pois X(ejω)=X(ej(ω+2πk)), para k, ela pode ser calculada em qualquer faixa de 2π, por exemplo ω[π,π).

x(n)1{X(ejω)} =def12πππX(ejω) ejωndω

Figura 2 - Sinal discreto x(n) e sua TFDT X(ejω)
Fonte: Elaborado pelo autor.

4 Série de Fourier de tempo contínuo (SF)

O sinal x(t) é contínuo e periódico no tempo (com período T) , e o sinal X(k) é discreto na frequência.

A equação de análise
É uma transformação de um domínio de uma variável real x(t) de tempo continuo em uma variável complexa X(k) de frequência discreta.
 DTDF.
x:.
X(k)=1T0Tx(t)ej(2π/T)kt)dt
A equação de síntese
É uma transformação de um domínio de uma variável complexa X(k) de frequência contínua em uma variável real x(t) de tempo continuo.
 DFDT.
X:.
x(t)=k=X(k)ej(2π/T)kt)
A série de Fourier X(k), indicada acima é a série exponencial onde as funções de base são exponenciais complexas ejωt, onde ω=(2π/T)k. Também existem as séries de Fourier usando funções de base senoidais e cossenoidais, as quais podem ser derivadas da série exponencial.

Figura 3 - Sinal x(t) periódico e sua série de Fourier X(Ω)
Arquivo:SerieFourierPlot.png
Fonte: Elaborado pelo autor.

5 Transformada de Discreta de Fourier (TDF)

5.1 Obtenção da TDF a partir da amostragem da TFTD

Sinais discretos no tempo podem ser representados pela sua TFTD, que é uma função continua periódica em 2π de ω :

X(ejω)=n=x(n)ejωn

Para que a mesma possa ser utilizada no processamento de sinais digitais é necessário que a variável frequência seja também discreta. Se amostrarmos uniformemente a frequência ω em N amostras entre 0 e 2π é possível obter a TDF (ou DFT - Discrete Fourier Transform). Assim se tomarmos N frequências ωk=(2π/N)k  com k , and 0kN1 , obtemos o espectro amostrado uniformemente:

X(ejω)=X(ejω)k=δ(ω2πNk).
X(ejω)=n=x(n)ej2πNkn

O sinal equivalente no tempo pode ser obtido aplicando a transformada inversa e a convolução:

x(n)=1{X(ejω)}=x(n)*(N2πp=δ(nNp))=N2πp=x(nNp).

O que mostra que o sinal Esse sinal x(n)  são repetições periódicas (com período N) do sinal discreto x(n)  original.

  • Note que N o período de repetição do sinal x(n)  é o mesmo período de repetição das N amostras da TFTD X(ejω)  original.
  • Se o comprimento L o sinal do x(n)  for maior que N o período de repetição do sinal x(n) , haverá sobreposição das amostras no tempo (time aliasing), e não será possível recuperar o sinal original.
  • Por outro lado, se LN  então x(n)  é a repetição periódica exata de x(n) .
x(n)=2πNx(n) , para 0nN1  ou
x(n)=2πNx(n)[u(n)u(nN) .

Portanto, é possível recuperar as amostras do sinal digital no tempo x(n)  a partir das suas amostras digitais na frequência, desde que o período de repetição das N amostras de frequência seja maior ou igual ao comprimento L do sinal no tempo.

5.2 DFT e IDFT

  • O sinal x(n)  é discreto no tempo pode ser representado pelo o sinal X(ej(2π/N)k)  discreto e periódico em 2π .

Para obter a equação de análise (DFT) pode ser feito o cálculo das amostras do espectro de frequências em ωk=(2π/N)k  em:

X(ejω)=n=x(n)ejωn
X(ej(2π/N)k)=n=x(n)ej(2π/N)kn

Conforme mostrado, o espectro é periódico em N, e portanto é suficiente calcular apenas os valores para 0kN1 . Assim obtém-se

A equação de análise (DFT)
É uma transformação de um domínio de uma variável real x(n)  de tempo discreto em uma variável complexa X(ej(2π/N)k)  frequência discreta periódica.
 DTDF.
x:.
X(ej(2π/N)k){x(n)} =defn=0N1x(n) ej(2π/N)kn , para 0kN1 
A equação de síntese (IDFT)
É uma transformação de um domínio de uma variável complexa X(ej(2π/N)k)  de frequência discreta periódica em uma variável real x(n)  discreta.
 DFDT.
X:.
x(n)1{X(ej(2π/N)k)} =def1Nk=0N1X(ej(2π/N)k) ej(2π/N)kn , para 0nN1 

Ao usar a equação de análise, se o comprimento L de x(n)  for menor que o período de repetição N, é necessário que x(n)  seja preenchido com amostras nulas até atingir o comprimento N (zero-padding).

Simplificação da notação

Para simplificar a notação é podemos utilizar:

X(k)=X(ej(2π/N)k)  para a frequência discreta periódica.

E ainda definir:

WN=ej2π/N  como a frequência fundamental

Essa frequência pode ser representada como uma posição no circulo unitário do plano complexo :z .

Alguns valores de Falhou ao verificar gramática (erro de sintaxe): {\displaystyle \mathrm{W_N \ } que ajudam a lembrar as simplificações:

W1=ej2π=+1 
W1=ej2π=+1 
W2=ejπ=1 
W2=ejπ=1 
W4=ejπ/2=j 
W4=ejπ/2=+j 


Também é importante lembrar que se N  é múltiplo de k  então:

WNk=ej(2π/N)k=ej2π/(N/k)=WN/k 

Dessa forma as equações da DFT e IDFT passam a ser escritas de forma simplificada como:

X(k)=n=0N1x(n)WNkn , para 0kN1 
x(n)=1Nk=0N1X(k)WNkn , para 0nN1 

Dessas equações é possível perceber que o cálculo da DFT e da IDFT requerem a N2  multiplicações e N(N1)  somas, sendo portanto um algoritmo de O(N2) .

5.3 Transformada Rápida de Fourier (FFT)

As FFTs são algoritmos que reduzem a complexidade do cálculo da DFT, obtendo ordens de O(logN×N) .


  • Para relembrar os conceitos vistos acima, recomendo assistir os vídeos sobre DFT do prof. Steve Brunton da University of Washington
The Discrete Fourier Transform (DFT) Computing the DFT Matrix The Fast Fourier Transform (FFT) The Fast Fourier Transform Algorithm

6 Soma dos termos de uma Progressão Geométrica

A soma aritmética dos termos de uma P.G. a partir do primeiro termo, é definida por:

Sn==a1+a1q+a1q2++a1qn1=i=1na1qi1

Caso q1 , a soma pode ser descrita por:

Sn=a1(1qn)1q