Funkcje jako Wektory
Szereg Fouriera to nie tylko narzędzie obliczeniowe — to operacja geometryczna: rzut ortogonalny funkcji na bazę.
Przestrzeń Funkcji
Zbiór funkcji całkowalnych z kwadratem na [0,1] tworzy przestrzeń wektorową L²[0,1]. Dodawanie & mnożenie przez skalar pracują punktowo. Iloczyn skalarny dwóch funkcji f, g:
⟨f, g⟩ = ∫₀¹ f(t) · g*(t) dt
gdzie g* jest sprzężeniem zespolonym g. To spełnia wszystkie aksjomaty iloczynu skalarnego.
Ortogonalność Bazy Fouriera
Funkcje φ_k(t) = e^{i2πkt} tworzą bazę ortonormalną dla L²[0,1]:
⟨φ_k, φ_m⟩ = ∫₀¹ e^{i2πkt} · e^{−i2πmt} dt = ∫₀¹ e^{i2π(k−m)t} dt = δ_{km}
(To równa się 1, jeśli k = m, 0 w pozostałych przypadkach — poprzez całkowanie czystych oscylacji przez pełny okres.)
Współczynnik Fouriera jako Iloczyn Skalarny
k-ty współczynnik Fouriera x(t):
c_k = ⟨x, φ_k⟩ = ∫₀¹ x(t) · e^{−i2πkt} dt
To jest rzut x na wektor bazowy φ_k. Współczynnik mierzy, ile φ_k jest zawarte w x.
Rzut na Podprzestrzeń
Obcięcie szeregu Fouriera do 2N+1 członów rzutuje x na podprzestrzeń rozpinającą {φ_{−N}, …, φ_N}. Obcięty szereg jest rzutem ortogonalnym x na tę skończeniewymiarową podprzestrzeń.
Zgodnie z nierównością Bessela, rzut minimalizuje błąd L²:
‖x − Σ_{k=−N}^{N} c_k φ_k‖² ≤ ‖x − Σ_{k=−N}^{N} a_k φ_k‖² dla dowolnego wyboru a_k
Obcięcie Fouriera jest najlepszym przybliżeniem w L² z tej podprzestrzeni. Minimalizuje średni błąd kwadratowy (kwadrat normy L² różnicy).
Okno Prostokątne → Jądro Sinc
Okno prostokątne w dziedzinie czasu (zachowując tylko współczynniki dla |k| ≤ N) odpowiada mnożeniu przez funkcję rect w indeksie współczynnika.
Mnożenie w jednej dziedzinie odpowiada konwolucji w drugiej dziedzinie.
Transformata Fouriera okna prostokątnego (w dyskretnej przestrzeni współczynników) jest jądrem Dirichleta — funkcją okresową podobną do sinc:
D_N(f) = Σ_{k=−N}^{N} e^{i2πfk} = sin(π(2N+1)f) / sin(πf)
Gdy obcinamy szereg Fouriera, dokonujemy konwolucji idealnej odpowiedzi H_ideal(f) z D_N(f).
Dlaczego Pojawia się Gibbs
Jądro Dirichleta ma duże listki boczne, które maleją powoli. W pobliżu nieciągłości skokowej w H_ideal(f), te listki boczne dzwonią — dodają się spójnie po jednej stronie skoku, dając przebieg ≈9%.
Stała matematyczna: ∫₀^π sin(t)/t dt = Si(π) ≈ 1.8519. Wysokość przebiegu Gibbsa = (2/π)·Si(π) − 1 ≈ 0.0895 = 8.95%. To jest niezależne od N.
Geometria Okna
Gładkie okno (Hamminga, Hanna, Kaisera) ma transformatę Fouriera z mniejszymi liskami bocznymi. Konwolucja H_ideal(f) z jądrem o mniejszych listkach bocznych daje mniej dzwonienia. Kompromis: mniejsze listki boczne zawsze wiążą się z szerszym płatem głównym, poszerzając pasmo przejściowe.
Stała Gibbsa
Przebieg Gibbsa to całka oznaczona, a nie funkcja N.
Pierwsze maksimum N-członowej sumy cząstkowej Fouriera skoku jednostkowego występuje przy f ≈ 1/(2N) od nieciągłości. Gdy N → ∞, to maksimum zbliża się do 1/(2)·(2/π)·Si(π) ≈ 1.0895.
Przebieg: 0.0895 lub w przybliżeniu 8.95% wysokości skoku.
Okna jako Jądra Domeny Częstotliwości
Każda funkcja okna ma transformatę Fouriera, która opisuje jądro używane do wygładzenia idealnej odpowiedzi częstotliwości.
Kluczowe parametry geometryczne jądra:
1. Szerokość płata głównego: określa szerokość pasma przejściowego (szerszy płat główny → szersze przejście).
2. Poziom szczytu listek bocznego: określa falowanie pasma przepustowego & pasma zatrzymania (niższe listki boczne → mniej falowania).
Te dwa parametry nie są niezależne. Dla danej długości okna 2N+1, zmniejszenie wysokości listek bocznego wymaga poszerzenia płata głównego — zawsze.
Okno Kaisera daje użytkownikowi jeden pokręcik (α) do ciągłego kompromisu wysokości listek bocznego vs szerokości płata głównego, zamiast przeskakiwać między stałymi typami okna.
Wgląd Projektowy
Szerokość pasma przejściowego ΔF ≈ szerokość płata głównego / N. Falowanie δ ≈ poziom listek bocznego. Obie formuły są przybliżone; równania Kaisera czynią je dokładnymi.