🌀 Epicykle Fouriera

Zjawisko Gibbsa to przeregulowanie pojawiające się, gdy szereg Fouriera przybliża sygnał ze skokową nieciągłością (jak fala prostokątna). W pobliżu nieciągłości suma częściowa przekracza prawdziwą wartość o około 9% wysokości skoku, niezależnie od liczby uwzględnionych wyrazów — przeregulowanie nie znika, tylko przesuwa się bliżej nieciągłości. Widać to w animacji epicykli jako „dzwonienie” przy ostrych narożnikach. Można je zredukować, stosując funkcje okna (czynnik sigma Lanczosa) lub sumowanie Fejéra zamiast sum częściowych.

Tak, po modyfikacjach. Dla ścieżki nieokresowej można albo potraktować skończony fragment jako jeden okres (co wprowadza skokową nieciągłość na końcach, powodując dzwonienie Gibbsa), albo zastosować dyskretną transformatę kosinusową (DCT) zamiast DFT — DCT niejawnie odbija ścieżkę, tworząc okresowe rozszerzenie bez skoku. Alternatywnie ścieżkę można sparametryzować tak, by ją zamknąć (np. dodając powrotny odcinek), czyniąc ją prawdziwie okresową. Standard kompresji obrazów JPEG wykorzystuje właśnie z tego powodu bloki DCT o rozmiarze 8×8 pikseli.

Rzeczywisty szereg Fouriera wyraża funkcję jako Σ (aₙ cos(nωt) + bₙ sin(nωt)), używając par wyrazów sinus i kosinus dla każdej częstotliwości. Zespolony szereg Fouriera wykorzystuje zespolone wykładniki Σ cₙ e^(inωt) dla n = −∞ do ∞, gdzie c₋ₙ = cₙ* (sprzężenie zespolone) dla sygnałów rzeczywistych. Forma zespolona jest bardziej zwarta i bezpośrednio odwzorowuje się na epicykle: dodatnie n daje okręgi przeciwne do ruchu wskazówek zegara, ujemne n — zgodne z ruchem wskazówek zegara. Większość wyprowadzeń w przetwarzaniu sygnałów korzysta z formy zespolonej dla wygody algebraicznej.

Symulator sortuje epicykle według promienia |cₙ| w kolejności malejącej — od największego. Jest to najbardziej czytelne wizualnie uporządkowanie, ponieważ największe okręgi (składowe o najniższej amplitudzie i najgrubszej częstotliwości) ustalają ogólny kształt, a mniejsze, szybsze okręgi dodają drobne detale. Matematycznie dowolne uporządkowanie daje ten sam ostateczny kształt; wybór wpływa jedynie na wizualny przebieg animacji i pośrednie pozycje końcówki łańcucha.

Zachowanie tylko nieparzystych harmonicznych (n = 1, 3, 5, 7, …) odpowiada szeregowi Fouriera sygnału o symetrii nieparzystej — takiego, dla którego f(t + T/2) = −f(t). Wizualnie powstała ścieżka ma symetrię obrotową 180°: druga połowa rysunku to pierwsza połowa obrócona o pół obrotu. Fale prostokątne, trójkątne i piłokształtne (w postaci symetrycznej względem połowy fali) składają się wyłącznie z nieparzystych harmonicznych, dlatego ich łańcuchy epicykli wykazują tę elegancką dwukrotną symetrię.

Tak — jest to popularna wizualizacja matematyczna od czasu, gdy spopularyzował ją film 3Blue1Brown z 2019 roku. Ścieżka odręcznej litery lub konturu portretu jest próbkowana, obliczana jest jej DFT, a powstałe epicykle animują odtwarzanie kształtu. Typowy kontur portretu próbkowany w 1000 punktach potrzebuje mniej więcej 200–300 istotnych epicykli, by wyglądać przekonująco. Technika jest zasadniczo taka sama jak kompresja ścieżek SVG: zachowanie tylko dominujących składowych częstotliwościowych odrzuca drobne, niedostrzegalne szczegóły.