Entendendo a Transformada de Fourier Intuitivamente (Sinais em Frequências)

Pressione três teclas de um piano ao mesmo tempo e uma única onda de pressão emaranhada sai do instrumento. Ela atinge seu tímpano como uma linha oscilante só, não três. Ainda assim, seu ouvido interno separa essa linha de volta em suas três notas sem nenhum esforço, e um músico treinado consegue nomear o acorde. Em algum ponto entre o ar e a sua consciência, um sinal misturado foi desmisturado em suas frequências puras.

Essa separação tem nome. A transformada de Fourier é a máquina matemática que pega qualquer sinal, um som, uma linha de imagem, o preço de uma ação, a vibração de uma ponte, e o reescreve como uma lista das frequências puras que existem dentro dele. É por isso que o seu telefone consegue comprimir áudio, é por isso que máquinas de ressonância magnética produzem imagens, e é por isso que engenheiros olham para um motor batendo e dizem "isso é a componente de 400 hertz".

A maioria dos estudantes conhece a transformada de Fourier como uma fórmula para decorar: uma integral com uma exponencial complexa dentro e nenhuma razão evidente para nada daquilo. Este guia constrói a imagem desde a base. Cobrimos há pouco a transformada de Laplace, que é bem próxima, e ao final você verá que ambas são a mesma máquina com o seletor em outra posição.

Duas descrições da mesma coisa

Pense na receita de um bolo. Você pode descrever o bolo pronto listando o que observaria em cada ponto dentro dele, camada por camada, migalha por migalha. Ou pode descrevê-lo listando os ingredientes: 200 gramas de farinha, 3 ovos, 100 gramas de açúcar. As duas descrições contêm o mesmo bolo. Uma descreve onde as coisas estão; a outra descreve do que ele é feito.

Um sinal tem essas mesmas duas descrições:

• O domínio do tempo lista o valor do sinal a cada instante: a pressão no seu tímpano em 0,001 segundo, em 0,002 segundo, e assim por diante. É assim que os sinais chegam no mundo real.
• O domínio da frequência lista os ingredientes: quanto de uma onda senoidal de 440 hertz está presente, quanto de 523 hertz, quanto de cada frequência possível. É isso que a transformada de Fourier produz.

A percepção de Fourier, publicada pela primeira vez em 1822 enquanto ele estudava o fluxo de calor, foi que essas listas de ingredientes são completas. Qualquer sinal razoável, por mais irregular que seja, pode ser reconstruído exatamente a partir de uma soma de ondas senoidais suaves. Não há um terceiro ingrediente escondido em algum lugar. Os senos bastam.

Por que as ondas senoidais são o alfabeto

Por que senos, e não quadrados, triângulos ou qualquer outro formato? Porque as ondas senoidais são os únicos formatos que sobrevivem às operações com que a física se importa.

Derive uma onda senoidal e você obtém um cosseno, que é a mesma onda deslocada por um quarto de ciclo. Derive de novo e você obtém um seno negativo: o mesmo formato outra vez, apenas escalado. Como exploramos no artigo sobre trigonometria, uma onda senoidal é um círculo visto de lado, e girar um círculo nunca o converte em outra coisa.

Injete uma onda senoidal em uma mola, em um circuito ou em uma sala de concertos, e o que sai é uma onda senoidal de mesma frequência, possivelmente mais alta ou mais baixa e possivelmente atrasada, mas nunca de formato diferente. Engenheiros chamam sistemas com essa propriedade de lineares e invariantes no tempo, e quase todo sistema físico é um deles, pelo menos aproximadamente.

Essa é a razão inteira de o domínio da frequência ser útil. Se você sabe o que um sistema faz com cada frequência pura, e consegue separar qualquer entrada em frequências puras, você sabe o que o sistema faz com toda entrada. O problema difícil de "o que essa ponte faz quando um caminhão passa por ela" vira uma lista de problemas fáceis: "o que essa ponte faz a 1 hertz, a 2 hertz, a 3 hertz".

A máquina de enrolamento

Agora ao mecanismo. Suponha que você tenha um sinal e queira saber quanto de uma frequência específica, digamos 3 ciclos por segundo, está escondido dentro dele. Eis um procedimento físico que responde à pergunta.

Pegue o gráfico do seu sinal, que é uma linha desenhada ao longo de um eixo de tempo. Agora enrole essa linha em torno de um círculo, a uma taxa de enrolamento de 3 voltas por segundo de tempo do sinal. Onde o sinal é grande, a linha enrolada fica longe do centro; onde o sinal é pequeno, ela se aproxima do centro. Você obtém um rabisco desenhado em volta da origem.

Agora encontre o centro de massa desse rabisco.

Se o sinal não contém nada a 3 hertz, os picos e vales do sinal caem em ângulos aleatórios ao redor do círculo. Eles se cancelam, e o centro de massa fica perto da origem. Mas se o sinal de fato contém uma componente de 3 hertz, seus picos caem todos do mesmo lado do círculo, volta após volta. Eles se acumulam, e o centro de massa se afasta da origem.

Repita isso para cada taxa de enrolamento possível, trace a distância do centro de massa contra a taxa de enrolamento, e você terá desenhado a transformada de Fourier. Os picos do gráfico são as frequências presentes no sinal. É essa a ideia inteira. Todo o resto é notação.

A fórmula, pedaço por pedaço

A máquina de enrolamento, escrita em símbolos, é:

$$\hat{f}(\omega) = \int_{-\infty}^{\infty} f(t)\, e^{-i\omega t}\, dt$$

Cada pedaço corresponde a uma etapa do procedimento.

O enrolamento. O fator $e^{-i\omega t}$ é o giro. Pela fórmula de Euler, $e^{-i\omega t} = \cos(\omega t) - i\sin(\omega t)$, um ponto que se move no sentido horário ao redor do círculo unitário com frequência angular ω. Conforme t avança, esse ponto gira. Esse giro é o enrolamento, e ω é a taxa de enrolamento.

O envolvimento. O produto $f(t)\, e^{-i\omega t}$ é o rabisco. Multiplicar o ponto girante de comprimento unitário pelo valor do sinal f(t) o estica para fora ou o puxa para dentro. Valores grandes do sinal empurram o ponto para longe do centro; valores pequenos o trazem para perto. O produto traça exatamente o rabisco da seção anterior.

O centro de massa é a integral. Somar a posição do rabisco a cada instante, do início ao fim dos tempos, e tratar o resultado como um único número complexo, produz um ponto cuja distância à origem mede com que intensidade o sinal ressoa em ω. A rigor, isso é o total e não a média, então ele escala com a duração do sinal, mas o formato do espectro é o mesmo.

O resultado $\hat{f}(\omega)$ é um número complexo para cada ω, e as duas partes importam. Seu módulo é a amplitude daquela frequência: quanto do ingrediente há na receita. Seu ângulo é a fase: em que ponto do ciclo aquele ingrediente começa. Dois sinais podem ter as mesmas amplitudes em todas as frequências e soar completamente diferentes porque suas fases diferem.

Por que a exponencial em vez de um seno

Os estudantes costumam perguntar por que o núcleo é uma exponencial complexa e não simplesmente uma onda senoidal. Afinal, você poderia multiplicar seu sinal por $\sin(\omega t)$ e integrar.

O problema é a fase. Se o seu sinal contém um cosseno na frequência ω e você o testa contra um seno de mesma frequência, o produto tem média zero, porque seno e cosseno estão a um quarto de ciclo de distância. Você concluiria que a frequência está ausente quando ela está plenamente presente. Para capturar toda fase seriam necessários dois testes, um contra o seno e um contra o cosseno, e depois combinar os resultados.

A exponencial complexa faz os dois testes de uma vez. Sua parte real é o teste do cosseno e sua parte imaginária é o teste do seno, empacotados em um único objeto girante. Com qualquer fase que a frequência chegue, alguma parte da rotação se alinha com ela. O módulo do número complexo resultante é o que você obteria combinando o método de dois testes, e o ângulo lhe dá a fase de graça. No artigo sobre números imaginários descrevemos i como uma rotação; aqui essa rotação faz trabalho honesto de engenharia.

Frequências negativas não são um erro

A integral percorre todo ω, incluindo valores negativos, e a transformada de um sinal real sempre mostra uma imagem espelhada do lado negativo. Isso incomoda as pessoas. O que uma frequência negativa deveria significar fisicamente?

Olhe de novo para a fórmula de Euler. Um cosseno real é a média de duas rotações:

$$\cos(\omega t) = \frac{e^{i\omega t} + e^{-i\omega t}}{2}$$

O primeiro termo gira no sentido anti-horário à taxa ω; o segundo gira no sentido horário à mesma taxa. A rotação horária em ω é o que a transformada registra como frequência menos ω. Um sinal de valores reais não consegue distinguir os dois sentidos, então ele sempre contém ambos em igual medida, e o espectro é simétrico.

A metade negativa não carrega informação nova para sinais reais, e na prática os engenheiros muitas vezes traçam apenas o lado positivo. Mas manter as duas metades deixa as fórmulas simétricas e permite que a mesma transformada trate sinais de valores complexos, como os do rádio e da ressonância magnética, onde os dois sentidos de rotação de fato diferem.

A viagem de volta

Converter uma receita de volta em um bolo é a transformada inversa de Fourier:

$$f(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} \hat{f}(\omega)\, e^{i\omega t}\, d\omega$$

Leia isso como uma instrução: para cada frequência ω, pegue a quantidade $\hat{f}(\omega)$ daquele ingrediente, gire-a no sentido anti-horário àquela taxa e some tudo. O sinal no expoente passou de menos para mais porque agora você está desenrolando em vez de enrolando. O fator 2π é contabilidade, corrigindo o fato de que uma volta completa tem 2π radianos; alguns livros o dividem simetricamente como $1/\sqrt{2\pi}$ nas duas transformadas, e outros o escondem por completo medindo a frequência em hertz em vez de radianos por segundo. A física não muda, apenas a constante.

O par de fórmulas é quase perfeitamente simétrico. Essa simetria é uma pista profunda: tempo e frequência são duas vistas de um mesmo objeto, e nenhuma é mais fundamental que a outra.

Derivadas viram multiplicação de novo

A transformada de Laplace ganhou seu espaço ao converter derivadas em multiplicação por s. A transformada de Fourier faz o mesmo truque com iω.

Derive a fórmula da transformada inversa em relação a t. O único lugar em que t aparece é dentro de $e^{i\omega t}$, e sua derivada puxa para baixo um fator iω:

$$\frac{d}{dt} f(t) = \frac{1}{2\pi} \int_{-\infty}^{\infty} \bigl(i\omega\, \hat{f}(\omega)\bigr)\, e^{i\omega t}\, d\omega$$

Então a transformada de Fourier de $f'(t)$ é simplesmente $i\omega\, \hat{f}(\omega)$. Derivar no tempo multiplica cada componente de frequência por iω: frequências altas são amplificadas, frequências baixas são suprimidas, e toda componente gira um quarto de volta. É exatamente por isso que a derivação deixa sinais ruidosos ainda mais ruidosos. O ruído vive nas frequências altas, e a derivada as multiplica por números grandes.

A mesma lógica converte uma equação diferencial de coeficientes constantes em uma equação algébrica em ω, que é como os físicos resolvem as equações da onda e do calor. Fourier inventou a ferramenta com exatamente esse propósito.

Fourier e Laplace: uma máquina, dois ajustes do seletor

Coloque as duas transformadas lado a lado:

$$F(s) = \int_0^{\infty} f(t)\, e^{-st}\, dt, \qquad \hat{f}(\omega) = \int_{-\infty}^{\infty} f(t)\, e^{-i\omega t}\, dt$$

No artigo sobre Laplace descrevemos s como um seletor com dois botões, $s = \sigma + i\omega$: σ para o decaimento exponencial, ω para a oscilação. Coloque o botão do decaimento em zero, de modo que σ = 0 e s = iω, e o núcleo de Laplace $e^{-st}$ vira o núcleo de Fourier $e^{-i\omega t}$. A transformada de Fourier é a transformada de Laplace avaliada ao longo do eixo imaginário do plano s.

A outra diferença é o limite inferior. Laplace começa em t = 0 porque foi feita para sistemas de causa e efeito, com um instante de acionamento e condições iniciais. Fourier percorre todo o tempo porque foi feita para sinais estacionários que sempre estiveram lá e sempre estarão: uma nota sustentada, uma onda portadora, um zumbido constante.

Qual usar decorre da pergunta. Transitórios, chaveamentos e condições iniciais: Laplace. Espectros, filtros e oscilação estacionária: Fourier. O botão do decaimento é o que permite a Laplace lidar com sinais que crescem sem limite, que a integral de Fourier não consegue absorver; a integral bilateral é o que permite a Fourier descrever sinais sem começo.

O compromisso do qual você não escapa

Uma última propriedade explica uma quantidade surpreendente de física. Um sinal estreito no tempo é largo em frequência, e um sinal estreito em frequência é largo no tempo. Você não pode ter os dois.

Uma onda senoidal pura tem uma única frequência, então seu espectro é um único pico. Mas uma onda senoidal pura dura para sempre. Corte-a em um segundo e o espectro se espalha em torno da frequência original; corte-a em um milissegundo e o espectro se espalha muito mais. No limite, um estalo instantâneo contém todas as frequências em igual medida, e é por isso que uma palma soa como ruído branco enquanto um diapasão soa como um tom.

A imagem do enrolamento deixa isso óbvio. Uma rajada curta de sinal se enrola em torno do círculo apenas algumas vezes em qualquer taxa de enrolamento, então o centro de massa nunca tem chance de se cancelar direito. Muitas taxas de enrolamento vizinhas dão um resultado parecido e moderado. Um sinal longo se enrola milhares de vezes, e só a taxa exatamente correspondente sobrevive ao cancelamento.

Na mecânica quântica, posição e momento são pares de Fourier, e esse mesmo compromisso é o princípio da incerteza de Heisenberg. Não é uma limitação da medição. É uma propriedade das ondas.

De onde vêm os erros

Estudantes que trabalham com transformadas de Fourier tropeçam em três lugares previsíveis.

O primeiro é descartar a fase. Como o espectro de amplitude é o que se costuma traçar, os estudantes passam a pensar na transformada como uma lista apenas de amplitudes. Aí uma transformada inversa sai errada porque os ângulos de $\hat{f}(\omega)$ foram jogados fora. A transformada é uma lista de números complexos, e metade de cada número é o seu ângulo.

O segundo é a confusão de sinais e constantes. Livros diferentes colocam o 2π em lugares diferentes e escolhem sinais diferentes no expoente, e todos estão corretos desde que a transformada direta e a inversa concordem. Escolha uma convenção, escreva as duas fórmulas no topo da página e verifique se transformar e destransformar uma função simples a devolve inalterada.

O terceiro é esperar que um sinal finito tenha um espectro limpo. Medições reais são finitas, então seus espectros sempre mostram o espalhamento da seção anterior, mais ondulações provocadas pelo corte abrupto nas bordas. Isso não é um erro de cálculo. É o compromisso se impondo, e as funções janela existem para administrá-lo.

Onde o Math Zen se encaixa

A transformada de Fourier é construída sobre três habilidades que precisam ser automáticas antes que a intuição consiga levar você adiante: identidades trigonométricas, exponenciais complexas e integrais de produtos.

Se converter entre $\cos(\omega t)$ e exponenciais complexas exige esforço consciente, ou se integrar $e^{-at} e^{-i\omega t}$ à mão parece incerto, a transformada sempre vai parecer uma caixa-preta. O Math Zen treina essas peças em sessões curtas e diárias:

• Trigonometria: Fluência com as identidades de soma e de arco duplo e com a imagem do círculo por trás delas, no módulo de prática de trigonometria.
• Números complexos: Forma polar, fórmula de Euler e rotação por multiplicação, no módulo de prática de números complexos.
• Transformadas: A mesma mecânica de integral com núcleo, incluindo transformadas de derivadas, no módulo de prática de transformadas de Laplace.

Quando isso vira reflexo, a transformada de Fourier deixa de ser uma fórmula e passa a ser o que ela é: uma máquina de enrolamento que lê a receita de um sinal.

A conclusão

A transformada de Fourier reescreve um sinal de "qual é o valor em cada instante" para "quanto de cada frequência está presente". Ela funciona porque as ondas senoidais são os formatos que os sistemas físicos deixam inalterados, então conhecer o conteúdo de frequências de um sinal diz como qualquer sistema linear vai responder a ele.

A fórmula é uma máquina de enrolamento: $e^{-i\omega t}$ enrola o sinal em torno de um círculo à taxa ω, a integral encontra o centro de massa, e um centro de massa longe da origem significa que aquela frequência está presente. A exponencial complexa captura toda fase de uma vez, as frequências negativas são a imagem espelhada que todo sinal real carrega, e a transformada inversa desenrola girando no sentido oposto.

Coloque o botão do decaimento da transformada de Laplace em zero e você tem a transformada de Fourier. Use Laplace para chaveamentos e transitórios, Fourier para espectros e sinais estacionários, e lembre-se de que, seja qual for a escolha, um sinal curto no tempo é sempre largo em frequência.