15. Szeregi Fouriera#
15.1. Harmoniki#
Narration
Niech \(\omega > 0\) oraz \(a, b\in\rr\), przy czym \(a^2 + b^2 \neq 0\).
Definicja 15.1
GEOGEBRA Harmoniką nazywamy funkcję \(f:\rr\to\rr\) o okresie \(\dfrac{2\pi} \omega\) zadaną wzorem
Liczbę \(\omega\) nazywamy pulsacją tej harmoniki.
Uwaga 15.1
W zastosowaniach ważna jest również inna postać harmoniki
gdzie \( r > 0\) oraz \(\varphi \in \rr\).
W celu znalezienia zależności między parametrami zauważmy, że
więc
Odwołując się do liczb zespolonych możemy zapisać równość \( \quad a - bj = r\cdot\left(\cos \varphi + j \sin \varphi\right), \quad \) z której otrzymamy równoważne zależności:
Definicja 15.2
Liczbę \(r > 0\) nazywamy amplitudą, zaś \(\varphi\) – fazą początkową tej harmoniki.
15.2. Szeregi trygonometryczne#
Narration
Przypomnijmy dwa ciągi funkcyjne:
funkcji trygonometrycznych cosinus, gdzie \(\ell\in\rr_+\) jest ustalone i \(f_n:\rr \to \rr\),
\[ \displaystyle f_n(x) = \cos \frac {n\pi x}{\ell}, \qquad n\in \nn_0\]Wspólny okres dla tych funkcji wynosi \(T=2\ell\).
funkcji trygonometrycznych sinus, gdzie \(\ell\in\rr_+\) jest ustalone i \(f_n:\rr \to \rr\),
\[ \displaystyle f_n(x) = \sin \frac {n\pi x}{\ell}, \qquad n\in \nn . \]Wspólny okres dla tych funkcji wynosi \(T=2\ell\).
Ustalmy \(\ell\in\rr_+\) i niech \((a_n)_{n\in\nn_0}\) i \((b_n)_{n\in\nn}\) będą ciągami liczb rzeczywistych.
Rozważmy ciąg funkcyjny \((f_n: \rr \to \rr)_{n\in\nn_0}\), harmonik o pulsacji \(\dfrac {n\pi}{\ell}\) (czyli okresie \(\dfrac {2\ell} n\) dla \(n\in \nn\)):
\[\begin{split} f_0(x) & = \frac{a_0} 2\ifcdot\cos \frac {0\pi x}{\ell} = \frac{a_0} 2, \\ f_n(x) & = a_n\ifcdot\cos \frac {n\pi x}{\ell} + b_n\ifcdot\sin \frac {n\pi x}{\ell} \qquad \text{dla }n\in\nn. \end{split}\]Funkcje \(f_n\) mają wspólny okres \(2\ell\), więc są jednoznacznie wyznaczone przez swoje wartości na przedziale \([-\ell, \ell]\).
Zatem w zależności od potrzeby możemy zajmować się:
funkcjami obciętymi \(f_n : [-\ell, \ell] \to \rr\)
funkcjami okresowymi \(f_n : \rr \to \rr\) o wspólnym okresie \(2\ell\).
Definicja 15.3
GEOGEBRA Wielomianem trygonometrycznym na przedziale \([-\ell, \ell]\) stopnia \(n\) nazywamy \(n\)-tą sumę częściową tego ciągu
Szeregiem trygonometrycznym na przedziale \([-\ell, \ell]\) nazywamy ciąg sum częściowych \((S_n)_{n\in\nn}\) oraz jego granicę
gdzie liczby \(a_n\) dla \(n\in\nn_0\), oraz \(b_n\) dla \(n\in\nn\) nazywamy współczynnikami szeregu trygonometrycznego.
Uwaga 15.2
W powyższych definicjach możemy przyjąć jedną z dwóch interpretacji:
\(S_n, S : [-\ell, \ell] \to \rr\)
\(S_n, S : \rr \to \rr\) są funkcjami okresowymi o wspólnym okresie \(T = 2\ell\).
W teorii sygnałów:
\(f_0 :\rr \to \rr\) to sygnał stały
\(f_n : \rr \to \rr\) to sygnał harmoniczny o okresie podstawowym \(\displaystyle\frac {2\ell} n\).
W szczególności, \(f_1\) jest sygnałem harmonicznym podstawowym o okresie \(2\ell\).
suma częściowa \(S_n : \rr \to \rr\) odpowiada wypadkowej tych sygnałów z odpowiednimi amplitudami.
Z rozważań o harmonikach wiemy, że dla \(n\in\nn\) mamy
\[ \begin{aligned} f_n(x) & = a_n\ifcdot\cos \frac {n\pi x}{\ell} + b_n\ifcdot\sin \frac {n\pi x}{\ell} = r_n \ifcdot\cos\left(\frac {n\pi x}{\ell} + \varphi_n\right), \end{aligned} \]gdzie \(r_n = \Absolute{a_n - b_n j} \geq 0\), oraz \(\varphi_n = \ArgPrincipal\Paren{a_n - b_n j} \in \left(-\pi, \pi \right] \) jeśli \(a_n\cdot b_n\neq 0\).
Stąd szereg trygonometryczny można inaczej zapisać w postaci
\[ S(x) = \frac {a_0} 2 + \Series r_n \ifcdot\cos\left(\frac {n\pi x}{\ell} + \varphi_n\right). \]
15.3. Warunki Dirichleta#
Narration
Niech \(X\subseteq \rr\) oraz \(a,b \in \rr\), \(a< b\).
Definicja 15.4
Funkcję \(f:X\to\rr\), ograniczoną w przedziale \((a, b)\), nazywamy przedziałami monotoniczną w tym przedziale, jeżeli można go podzielić na skończoną liczbę podprzedziałów, wewnątrz których funkcja \(f\) jest monotoniczna.
Przykład 15.1
Funkcja \(f:\rr\to \rr\),
jest przedziałami monotoniczna na każdym przedziale \((a, b)\).
Przykład 15.2
Funkcja \(f:(0, + \infty)\to \rr\),
dla każdego \(b>0\), nie jest przedziałami monotoniczna w przedziale \((0, b)\)
dla każdych \(b > a >0\), jest przedziałami monotoniczna na każdym przedziale \((a, b)\).
Zauważmy, że \(f\) jest złożeniem funkcji \(\sin x\) oraz \(\frac 1 x\), które są przedziałami monotoniczne na każdym podprzedziale \(\rr\setminus\Set{0}\).
Zatem złożenie funkcji przedziałami monotonicznych nie musi mieć tej własności.
Przykład 15.3
Podobnie suma funkcji przedziałami monotonicznych nie musi mieć tej własności.
W tym celu można pokazać, że \(f_1(x) = x^2\sin\frac 1 x + 3x\) oraz \(f_2(x) = -3x\) są monotoniczne (więc również na przedziałami) w przedziale \((0,+\infty)\), ale ich suma \(f_1(x) + f_2(x) = x^2\sin\frac 1 x\) nie jest nawet przedziałami monotoniczna.
Stwierdzenie 15.1
Jeśli ograniczona funkcja \(f:(a,b)\to\rr\) ma skończoną liczbę ekstremów lokalnych, to jest przedziałami monotoniczna.
Definicja 15.5
Funkcja \(f:X\to \rr \), ograniczona w przedziale \((a,b)\), spełnia w przedziale \([a, b]\) warunki Dirichleta , jeżeli:
\(f\) jest przedziałami monotoniczna w \((a, b)\)
\(f\) jest ciągła w \((a, b)\), z wyjątkiem co najwyżej skończonej liczby punktów nieciągłości I-go rodzaju, przy czym w każdym punkcie nieciągłości \(x_0\) spełniony jest warunek
\[ f(x_0) = \frac 1 2 \left(\LimIn{x_0^-} f(x) + \LimIn{x_0^+} f(x) \right) \]na końcach \([a, b]\) spełnione są równości
\[ f(a) = f(b) = \frac 1 2 \left( \LimIn{a^+} f(x) + \LimIn{b^-} f(x) \right). \]Warunki te nazywamy odpowiednio: pierwszym, drugim i trzecim warunkiem Dirichleta.
Ćwiczenie 15.1
Funkcja \(f:[-1, 1]\to \rr\), \(f(x) = x - [x],\) nie spełnia warunków Dirichleta na przedziale \([-1,\ 1]\).
Jak należy zmienić jej wartości w punktach nieciągłości i na brzegu dziedziny, by warunki były spełnione?
Rozwiązanie
Funkcja jest rosnąca na \((-1, 0)\) i na \((0, 1)\) więc jest przedziałami monotoniczna.
Funkcja jest ciągła w przedziale \((-1, 1)\) poza punktem \(x_0 = 0\).
Drugi warunek Dirichleta będzie spełniony jeśli
\[ \begin{aligned} f(0) & = \frac 1 2 \left(\LimIn{0^-} f(x) + \LimIn{0^+} f(x) \right) = \frac 1 2 \Paren{\LimIn{0^-} \Paren{x+1} + \LimIn{0^+} x } = \frac 1 2 \Paren{1 + 0} = \frac 1 2 . \end{aligned} \]Trzeci warunek Dirichleta będzie spełniony jeśli
\[ f(-1) = f(1) = \frac 1 2 \Paren{ \LimIn{-1^+} f(x) + \LimIn{1^-} f(x) } = \frac 1 2 \Paren{ \LimIn{-1^+} \Paren{x+1} + \LimIn{1^-} x } = \frac 1 2 \Paren{ 0 + 1} = \frac 1 2 . \]
15.4. Szereg trygonometryczny Fouriera#
Narration
Niech \(\ell >0\) i niech \(f: [-\ell, \ell] \to \rr\) będzie funkcją całkowalną.
Definicja 15.6
Przyjmijmy:
w szczególności
Szereg trygonometryczny
nazywamy szeregiem trygonometrycznym Fouriera funkcji \(f\) w przedziale \([-\ell, \ell]\).
Twierdzenie 15.2 (Dirichleta)
Niech \(\ell >0\) i niech \(f: [-\ell, \ell] \to \rr\) będzie funkcją całkowalną spełniającą warunki Dirichleta na przedziale \([-\ell,\ell]\).
Wówczas
dla każdego \(x\in [-\ell, \ell]\) zachodzi równość
\[ f(x) = S_f(x) = \frac{a_0} 2 + \Series \left[a_n\ifcdot\cos \frac {n\pi x}{\ell} + b_n\ifcdot\sin \frac {n\pi x}{\ell} \right], \]gdzie współczynniki \(a_n, b_n\) dane są wzorami (15.1)
jeśli dla każdego \(x\in [-\ell, \ell],\)
\[ f(x) = \frac{a_0} 2 + \Series \left[a_n\ifcdot\cos \frac {n\pi x}{\ell} + b_n\ifcdot\sin \frac {n\pi x}{\ell} \right], \]to współczynniki \(a_n, b_n\) dane są wzorami (15.1).
Wniosek 15.3
Jeśli \(f:\rr \to \rr\) jest funkcją całkowalną, okresową o okresie podstawowym \(2\ell\), spełniającą warunki Dirichleta na przedziale \([-\ell,\ell]\), to powyższa równość jest prawdziwa dla \(x\in \rr\) i
dla \(n\in \nn_0\), gdzie \(x_0\in\rr\) można wybrać dowolnie (w szczególności można przyjąć \(x_0 = 0\)).
Uwaga 15.3
Jeśli dla funkcji \(f\) zachodzi równość z powyższego twierdzenia, to mówimy, że funkcja \(f\) rozwija się w szereg Fouriera w przedziale \([-\ell, \ell]\).
Punkt 2. z powyższego twierdzenia mówi, że rozwinięcie w szereg Fouriera funkcji \(f\) jest jednoznaczne jeśli tylko istnieje.
Jeśli \(f:[-\ell, \ell]\to \rr\) jest funkcją przedziałami monotoniczną i nieciągłą w skończonej liczbie punktów, przy czym są to nieciągłości I-go rodzaju, to można „poprawić” jej wartości w tych punktach, tak by warunki Dirichleta były spełnione.
Następnie można rozwinąć „poprawioną” funkcję w szereg Fouriera i równość z twierdzenia Dirichleta, będzie zachodzić we wszystkich punktach ciągłości \(f\) z przedziału \((-\ell, \ell)\).
Jeśli \(f:[-\ell, \ell] \to \rr\) jest nieparzysta i rozwija się w szereg Fouriera, to:
funkcja \(\displaystyle x \mapsto f(x)\cdot \cos \frac {n\pi x}{\ell} \) jest nieparzysta jako iloczyn funkcji nieparzystej i parzystej; to implikuje, że
\[ a_n = \frac 1 \ell \IntDef{-\ell}{\ell}{f(x)\cdot \cos \frac {n\pi x}{\ell} } = 0, \qquad \text{dla } n\in \nn_0 \]funkcja \(\displaystyle x \mapsto f(x)\cdot \sin \frac {n\pi }{\ell} x \) jest parzysta, jako iloczyn funkcji nieparzystych; to implikuje, że
\[ b_n = \frac 2 \ell \IntDef{0}{\ell}{f(x)\cdot \sin \frac {n\pi x}{\ell} }, \quad \text{dla } n\in \nn. \]
Mówimy, że funkcja \(f\) rozwija się w szereg sinusów:
\[ \quad f(x) = \Series b_n\ifcdot\sin \frac {n\pi x}{\ell}. \]Analogicznie, jeśli \(f: [-\ell, \ell] \to \rr\) jest parzysta i rozwija się w szereg Fouriera, to:
funkcja \(\displaystyle x \mapsto f(x)\cdot \cos \frac {n\pi x}{\ell} \) jest parzysta jako iloczyn funkcji parzystych; to implikuje, że
\[ a_n = \frac 2 \ell \IntDef{0}{\ell}{f(x)\cdot \cos \frac {n\pi x}{\ell} }, \quad \text{dla } n\in \nn_0 \]funkcja \(\displaystyle x \mapsto f(x)\cdot \sin \frac {n\pi }{\ell} x \) jest nieparzysta, jako iloczyn funkcji parzystej i nieparzystej; to implikuje, że
\[ b_n = \frac 1 \ell \IntDef{-\ell}{\ell}{f(x)\cdot \sin \frac {n\pi x}{\ell} } = 0, \quad \text{dla } n\in \nn. \]
Mówimy, że funkcja \(f\) rozwija się w szereg cosinusów:
\[ f(x) = \frac{a_0} 2 + \Series a_n\ifcdot\cos \frac {n\pi x}{\ell}. \]
Ćwiczenie 15.2
Znajdźmy rozwinięcie Fouriera (w szczególności współczynniki \(a_n\) i \(b_n\)) funkcji \(f(x) = \sin x\) w przedziale \([-\ell, \ell]\) dla:
\(\ell = 2\pi\)
\(\ell = \pi\)
\(\displaystyle \ell = \frac \pi 2\).
Informacja
Pomoże rozwiązać Ćw. 16.1.
Rozwiązanie
W każdym z trzech przypadków sprawdzamy warunki Dirichleta:
Funkcja jest przedziałami monotoniczna, bo ma skończoną liczbę ekstremów lokalnych.
Funkcja jest ciągła.
Trzeci warunek Dirichleta będzie spełniony jeśli
\[ f(-\ell) = f(\ell) = \frac 1 2 \Paren{ \LimIn{-\ell^+} \sin x + \LimIn{\ell^-} \sin x } = \frac 1 2 \Paren{ \sin(-\ell) + \sin(\ell)} = 0 . \]
W pierwszych dwóch przypadkach warunek jest spełniony, zaś w trzecim należy przyjąć \(f\Paren{-\frac \pi 2} = f\Paren{\frac \pi 2} = 0.\)
Ponadto zauważmy, że funkcja jest nieparzysta, więc \( \quad \Forall{n \in \nn_0} a_n = 0. \)
Jeśli \( \quad \ell = 2\pi\), \(\quad\) to \( \quad \dfrac {n\pi x}{\ell} = \dfrac {n\pi x}{2\pi} = n\cdot \dfrac {x}{2}. \quad \) Zapiszmy funkcję w postaci
\[ f(x) = 1\cdot \sin\Paren{2\cdot \frac x 2 } = 0\cdot \sin\Paren{1\cdot \frac x 2 } + 1\cdot \sin\Paren{2\cdot \frac x 2 } + 0\cdot \sin\Paren{3\cdot \frac x 2 } + \ldots . \]To jest rozwinięcie w szereg Fouriera!
Stąd wniosek, że:
\[ b_1 = 0, \quad b_2 = 1, \quad b_3 = 0, \ldots \qlq b_2 = 1 \ \wedge \ \Forall{n\in \nn\setminus\Set{2}} b_n = 0 \, . \]Jeśli \( \quad \ell = \pi\), \(\quad\) to \( \quad \dfrac {n\pi x}{\ell} = \dfrac {n\pi x}{\pi} = n\cdot x . \quad \) Zapiszmy funkcję w postaci
\[ f(x) = 1\cdot \sin\Paren{1\cdot x } = 1\cdot \sin\Paren{1\cdot x } + 0\cdot \sin\Paren{2\cdot x } + 0\cdot \sin\Paren{3\cdot x } + \ldots . \]To jest rozwinięcie w szereg Fouriera!
Stąd wniosek, że:
\[ b_1 = 1, \quad b_2 = 0, \quad b_3 = 0, \ldots \qlq b_1 = 1 \ \wedge \ \Forall{n\in \nn_2} b_n = 0 \, . \]Jeśli \( \quad \ell = \dfrac \pi 2\), \(\quad\) to \( \quad \dfrac {n\pi x}{\ell} = \dfrac {n\pi x}{\frac \pi 2} = n\cdot 2x . \quad \) Tym razem musimy policzyć całki:
Zatem
Ćwiczenie 15.3
Rozwińmy w szereg Fouriera funkcję
w przedziale \([-\pi, \pi]\) i podajmy wartości całek:
Rozwiązanie
Pozostawiam czytelnikowi sprawdzenie warunków Dirichleta.
Jeśli \( \quad \ell = \pi\), \(\quad\) to \( \quad \dfrac {n\pi x}{\ell} = \dfrac {n\pi x}{\pi} = n\cdot x . \quad \) Korzystając ze wzorów redukcyjnych zapiszmy funkcję w postaci
To jest rozwinięcie w szereg Fouriera.
Stąd wniosek, że:
W celu obliczenia całek zauważamy, że:
Ćwiczenie 15.4
Niech \(A\in \rr \).
Rozwińmy w szereg Fouriera funkcję \(f :[-\pi,\pi] \to \rr\),
Rozwiązanie
Pozostawiam czytelnikowi sprawdzenie warunków Dirichleta.
Ponieważ funkcja jest nieparzysta, to \( \quad \Forall{n \in \nn_0} a_n = 0. \)
Pozostaje obliczenie współczynników (dla \(n\in \nn\))
Zauważmy, że jedynie nieparzyste wyrazy są niezerowe:
Zatem
Przyjmując \(A = \tfrac \pi 2\) otrzymamy następujące wykresy:
15.5. Rozwinięcie Fouriera rozszerzenia do funkcji okresowej#
Narration
Jeżeli funkcja \(f:(a, a+2\ell) \to \rr\) spełnia 1. i 2. warunek Dirichleta w przedziale \((a, a+2\ell)\), to można ją przedłużyć do funkcji okresowej \(\tilde f: \rr \to \rr\) o okresie \(2\ell\), tak aby:
dla każdego \(x \in (a, a+2\ell)\) zachodziło \(\tilde f(x) = f(x)\) (tzn. \(\tilde f|_{(a, a+2\ell)} = f\))
funkcja \(\tilde f\) spełniała wszystkie trzy warunki Dirichleta w przedziale \([-\ell, \ell]\).
Funkcję \(\tilde f\) można rozwinąć w szereg Fouriera \(S_{\tilde f}\) w przedziale \([-\ell, \ell]\).
Z twierdzenia Dirichleta i warunku 1. wynika, że \(f(x) = \tilde f(x) = S_{\tilde f}(x)\) dla każdego \(x \in (a, a+2\ell)\).
W ten sposób rozwinęliśmy funkcję \(f\) w szereg Fouriera w przedziale \((a, a+2\ell)\).
Ćwiczenie 15.5
Zgodnie z powyższą procedurą rozwińmy w szereg Fouriera funkcję
Rozwiązanie
Dla danej funkcji mamy:
Przedłużamy ją do funkcji okresowej jak pokazano na rysunku.
Funkcja nie jest ani parzysta, ani nieparzysta, więc musimy obliczyć wszystkie współczynniki.
Będziemy potrzebować dwóch całek nieoznaczonych:
Korzystając z tych obliczeń otrzymujemy:
Podsumowując,