Construtor de séries de Fourier

componha ondas complexas a partir de senos simples

A carregar a simulação interativa...

Lição

A teoria — Construtor de séries de Fourier

Uma série de Fourier afirma que qualquer onda periódica, por mais recortada que seja, é uma soma de senos puros cujas frequências são múltiplos inteiros de uma fundamental. A soma é construída termo a termo e desenhada por cima do alvo. A palavra interessante é igual: a soma alcança o alvo onde o alvo é contínuo, e a linha Peak é onde se descobre o que acontece onde não é.

O que significa cada símbolo

N
quantos termos estão a ser somados. O cursor termina em 20, e a linha da série nomeia cada coeficiente que usa.
k
o índice do harmónico — que múltiplo da fundamental um termo carrega. A quadrada e a triangular usam só k ímpares; a dente de serra usa todos.
aₖ
o coeficiente do harmónico k, impresso na linha da série. O seu tamanho diz quanto essa frequência importa, e a rapidez com que os coeficientes encolhem decide a rapidez da convergência.
Peak
max|fN(x)|, o ponto mais alto da soma parcial. O alvo tem amplitude 1, portanto esta linha pergunta na verdade se a soma passou.

De onde vem a fórmula

  1. Lê os coeficientes na página. Para a quadrada, a linha da série imprime 1.273 sin(x) + 0.4244 sin(3x) + 0.2546 sin(5x) + …, ou seja 4/(kπ) em k = 1, 3, 5. Nada é ajustado — os coeficientes vêm de um integral, e cada forma de onda tem os seus.
  2. Repara em que harmónicos faltam. Quadrada e triangular usam só k ímpares porque ambas cumprem f(x + π) = −f(x): deslocadas meio período trocam de sinal, e um harmónico par não consegue fazer isso. A dente de serra não tem essa simetria e por isso precisa de todos os harmónicos — o que a linha «Harmonics used» diz sem rodeios.
  3. O primeiro coeficiente já denuncia a velocidade de convergência. A quadrada dá 4/π = 1.273 e decai como 1/k; a triangular 8/π² = 0.8106 e decai como 1/k²; a dente de serra 2/π = 0.6366 e decai como 1/k. Decair mais depressa significa menos termos para a mesma fidelidade — por isso a triangular parece acabada muito antes das outras duas.
  4. Olha agora para a linha Peak em vez do desenho. Com um alvo sem saltos, o pico devia subir até 1 e ficar. Com um alvo com salto faz algo muito diferente, e é disso que trata o resto desta lição.

Como ler o que vê

Põe a forma em triangular e sobe o N desde 1: o pico marca 0.8106, 0.9006, 0.9331, 0.9596 … 0.9899 em N = 20. Sobe para 1 por baixo e nunca o ultrapassa. Muda para quadrada: 1.273, 1.2, 1.188, 1.182 … e depois simplesmente pára, em 1.179, e nenhum harmónico adicional o move. Por fim a dente de serra: 0.6366, 0.827, 0.9185, depois 1.008 em N = 5 — atravessa o 1 — e em N = 20 vai em 1.131 e continua a subir. Três formas e três histórias completamente diferentes numa só linha. A dente de serra é a que engana: olhando só para ela concluirias que o excesso cresce sem limite. Não cresce. Vai para o mesmo 1.179 da quadrada, apenas por baixo.

Pressupõe
Que vinte termos chegam para ver a tendência, e que o pico é encontrado por amostragem — a linha RMS é avaliada em 801 pontos do período, e o pico é o maior dessas amostras, não um máximo resolvido. Ambas as coisas servem para observar a convergência e nenhuma serviria para uma demonstração.
Falha quando
O excesso nunca desaparece. Aquele 1.179 estável não é a série a ficar sem fôlego: o limite é (2/π)·Si(π) = 1.17898…, e acrescentar harmónicos eternamente estreita o pico sem o baixar. Portanto a série converge para a quadrada em cada ponto isolado e nunca uniformemente — exatamente a distinção que este fenómeno existe para ensinar. Explica também uma contradição aparente nesta página: a ideia acima fala de cerca de 8.95% e esta linha de 17.9%. Ambas estão certas — 0.179 medido contra a amplitude 1 é 0.0895 medido contra o salto de 2. Michelson deparou-se com o pico em 1898 no seu próprio analisador harmónico e suspeitou de uma avaria na máquina; Gibbs escreveu no ano seguinte à Nature a explicar que era a matemática a fazê-lo, e não o latão.

como vetores rotativos constroem uma forma de onda 🖖

Cada vetor rotativo contribui com uma componente senoidal no índice harmônico k. Somando os vetores ponta a ponta forma-se a soma parcial fN(x). À medida que N cresce, o conteúdo de alta frequência afia os cantos, enquanto perto das descontinuidades de salto permanece um sobressinal de Gibbs limitado.

Toda onda periódica é um acorde de senoides 🖖

A ideia de Fourier: qualquer sinal periódico, por mais irregular que seja, é uma soma de ondas senoidais puras cujas frequências são múltiplos inteiros de uma fundamental. Os controles apenas escolhem quantos desses tons ingredientes empilhar. Numa onda quadrada surgem só harmônicos ímpares, cada um mais fraco que o anterior: o 3.º carrega 1/3 da amplitude da fundamental, e o 5.º apenas 1/5. Mais termos afinam a soma rumo à forma alvo.

O sobressinal de 9% que nunca desaparece 🖖

Perto de um salto, a soma parcial ultrapassa o alvo em cerca de 8.95% da altura do salto — e esse pico nunca encolhe, por mais harmônicos que você some; apenas fica mais estreito. O físico Albert Michelson deparou com ele em 1898, quando seu analisador harmônico mecânico insistia em desenhar picos nas bordas de uma onda quadrada, e ele suspeitou de uma falha do aparelho. J. Willard Gibbs provou que a máquina estava certa: o sobressinal é matemática real. O mesmo 'ringing' assombra imagens JPEG perto de bordas nítidas.

TEORIA DE FOURIER — POR QUE ALGUMAS ONDAS PRECISAM DE MAIS HARMÔNICOS

Ler uma forma de onda: com que rapidez seus harmônicos decaem?

Toda onda periódica é uma soma de senos. Quantos você precisa até a soma ficar certa não é questão de gosto: é fixado por quão suave é a forma de onda alvo. Um salto custa muito mais harmônicos que um canto, a simetria pode cortar o trabalho pela metade de graça, e há um erro específico que não some por mais termos que você acrescente.

Quadrada — um salto, então decaimento 1/n cn ∝ 1/n, n = 2k+1
Dente de serra — um salto sem simetria, então todos os harmônicos cn ∝ 1/n, ∀n
Triangular — um canto, não um salto, então decaimento 1/n² cn ∝ 1/n2, n = 2k+1
Um harmônico — o preço de parar cedo ε ≈ ∑n>N cn

01

Quadrada — um salto, então decaimento 1/n

A forma de onda tem: Um salto vertical entre dois níveis, mais simetria de meia onda.

Decaimento dos coeficientes: cn ∝ 1/n, n = 2k+1

O que se vê: 10 harmônicos dão bordas nítidas que ainda assim ultrapassam cerca de 9% da altura do salto.

Abrir este caso: quadrada, 10 harmônicos
Quadrada — um salto, então decaimento 1/n. Dez harmônicos: as bordas ficam nítidas, mas a ultrapassagem continua ali. Um salto vertical entre dois níveis, mais simetria de meia onda.
Dez harmônicos: as bordas ficam nítidas, mas a ultrapassagem continua ali.

02

Dente de serra — um salto sem simetria, então todos os harmônicos

A forma de onda tem: Um salto por período e nenhuma simetria de meia onda para cancelar nada.

Decaimento dos coeficientes: cn ∝ 1/n, ∀n

O que se vê: 10 harmônicos: oscilação em cada salto, e todo harmônico 1, 2, 3… contribui.

Abrir este caso: dente de serra, 10 harmônicos
Dente de serra — um salto sem simetria, então todos os harmônicos. Todos os harmônicos contribuem, e a oscilação fica em cada salto. Um salto por período e nenhuma simetria de meia onda para cancelar nada.
Todos os harmônicos contribuem, e a oscilação fica em cada salto.

03

Triangular — um canto, não um salto, então decaimento 1/n²

A forma de onda tem: A onda em si é contínua; só a inclinação dela salta.

Decaimento dos coeficientes: cn ∝ 1/n2, n = 2k+1

O que se vê: 8 harmônicos já ficam quase exatamente sobre o alvo, sem ultrapassagem nenhuma.

Abrir este caso: triangular, 8 harmônicos
Triangular — um canto, não um salto, então decaimento 1/n². Oito harmônicos acompanham a triangular quase exatamente, sem oscilação. A onda em si é contínua; só a inclinação dela salta.
Oito harmônicos acompanham a triangular quase exatamente, sem oscilação.

04

Um harmônico — o preço de parar cedo

A forma de onda tem: Você cortou a soma num N pequeno. A pergunta é o que ficou para trás.

Decaimento dos coeficientes: ε ≈ ∑n>N cn

O que se vê: Uma quadrada com n = 1 é um único seno: o período certo, nada da forma.

Abrir este caso: quadrada, 1 harmônico
Um harmônico — o preço de parar cedo. n = 1 contra n = 10: o período acerta muito antes da forma. Você cortou a soma num N pequeno. A pergunta é o que ficou para trás.
n = 1 contra n = 10: o período acerta muito antes da forma.
Referências (4)

Problema resolvido na íntegra

  1. Cinco ondas senoidais a aproximar uma onda quadrada e a atingir 1,182 5 passos

    Cinco ondas sinusoidais a aproximar uma onda quadrada ultrapassam-na, atingindo 1,182 contra o 1 da onda quadrada. Descubra se adicionar mais termos resolve o problema.

    1. Uma onda quadrada é ímpar e tem simetria de meia-onda, o que elimina todos os harmónicos pares antes de qualquer aritmética acontecer. Cinco termos significam, portanto, k = 1, 3, 5, 7, 9 — a lista de harmónicos do painel é uma consequência da simetria, não uma definição.

    2. A soma parcial sofre um sobressinal perto do salto. Para N = 5 o pico é 1,182 e o RMS é 0,979, já perto do RMS exato da onda, que é 1 — portanto, a energia está quase correta enquanto a forma perto da extremidade não está.

    3. Teste agora o remédio óbvio. N = 1 tem um pico em 1,273, N = 5 em 1,182, N = 50 em 1,179, N = 200 em 1,17898.

    4. Está a convergir — mas não para 1. O limite é (2/π)·Si(π) = 1,17898, e o sobressinal nunca desaparece para qualquer número de termos.

    5. O que acontece é que o sobressinal fica mais estreito. Desliza em direção à descontinuidade e ocupa cada vez menos largura, razão pela qual o erro RMS tende para zero, embora o erro de pico não o faça.

    Resposta

    A ferramenta indica N = 5, harmónicos 1, 3, 5, 7, 9, um pico de 1,182 e RMS de 0,979. O resultado a reter é que o pico converge para 1,17898, e não para 1. Trata-se do fenómeno de Gibbs, e citá-lo exige cuidado: o sobressinal é de 8,95% do salto — e esta onda salta 2, de −1 a +1, pelo que o mesmo facto se lê como 17,9% da amplitude. Ambos os números estão corretos e nenhum é o outro. Também não é um artefacto numérico: é o aspeto da convergência pontual quando o limite é descontínuo, e é a razão pela qual as arestas nítidas oscilam em JPEG e em filtros de áudio.

Percurso de aprendizagem

Círculos e ondas são o mesmo movimento

Problemas de exemplo