Fourierreihe

aus Wikipedia, der freien Enzyklopädie
Wechseln zu: Navigation, Suche

Als Fourierreihe (nach Joseph Fourier) bezeichnet man die Reihenentwicklung einer periodischen, abschnittsweise stetigen Funktion in eine Funktionenreihe aus Sinus- und Kosinusfunktionen.

Die Basisfunktionen der Fourierreihe bilden ein bekanntes Beispiel für eine Orthonormalbasis. Im Rahmen der Theorie der Hilberträume werden auch Entwicklungen nach einem beliebigen vollständigen Orthonormalsystem als Fourierreihe bezeichnet.

Geschichte[Bearbeiten]

Joseph Fourier

Bereits im 18. Jahrhundert kannten Mathematiker wie Euler, Lagrange oder die Bernoullis Fourierreihen für einige Funktionen. Zu Beginn des 19. Jahrhunderts behauptete nun Fourier in seinem Werk Théorie analytique de la chaleur, dass es für alle Funktionen solche Reihenentwicklungen gäbe. Diese Behauptung stieß zunächst bei führenden Mathematikern wie Cauchy und Abel auf Ablehnung.

Dirichlet konnte 1829 beweisen, dass Fouriers Behauptung zumindest für Lipschitz-stetige Funktionen zutrifft. Du Bois-Reymond fand 1876 eine stetige Funktion, deren Fourierreihe divergiert. Im 20. Jahrhundert gelangte man schließlich zur Erkenntnis, dass es auch für stetige oder stückweise stetige Funktionen konvergente Fourierreihen gibt, wenn der Konvergenzbegriff geeignet abgeschwächt wird (Lennart Carleson).

Als eine frühe geometrische Vorform der Approximation durch eine Fourierreihe kann die Epizykeltheorie betrachtet werden.

Darstellungsformen[Bearbeiten]

Das dick gezeichnete Signal wird durch die Fourier-Analyse in die beiden dünn gezeichneten Signale zerlegt. Die positive und die negative Halbschwingung sehen gleich aus, weil ein Anteil mit dreifacher Frequenz addiert wurde.
Hier sind positive und die negative Halbschwingung unterschiedlich, weil ein Anteil mit doppelter Frequenz addiert wurde.

Die Partialsummen einer Fourierreihe sind trigonometrische Polynome. Wie diese können Fourierreihen in drei gleichwertigen Formen dargestellt werden. Zu jeder dieser Darstellungen gibt es zugehörige Formeln zum Bestimmen der Koeffizienten bzw. Parameter der Fourierreihenentwicklung einer periodischen Funktion.

Eine Fourierreihenentwicklung einer periodischen Funktion f mit Periode T>0 ist in den folgenden, schrittweise allgemeiner werdenden Fällen möglich:

  1. wenn f stetig und abschnittsweise stetig differenzierbar ist; die Fourierreihe konvergiert dabei punktweise und gleichmäßig.
  2. wenn f eine beschränkte totale Variation über einer Periode hat, und die Funktionswerte von f mit dem Mittel aus den links- und rechtsseitigen Grenzwerten übereinstimmen, 2\,f(x)=f(x+0)+f(x-0) für alle x \in \R; die Fourierreihe konvergiert dann nur punktweise.
  3. wenn f, auf eine Periode [c,c+T] eingeschränkt, dem Funktionenraum L^2([c,c+T]) angehört; mit Konvergenz im Sinne der L²-Norm.

Allgemeine Form[Bearbeiten]

Eine periodische Funktion f mit Periode T>0, die einer der angegebenen Klassen angehört, lässt sich durch eine Reihe von Sinus- und Kosinusfunktionen darstellen, deren Frequenzen ganzzahlige Vielfache der Grundfrequenz \omega=2\pi / T sind,


  f(t)=\frac{a_0}{2} + \sum_{k=1}^\infty (a_k \cdot \cos(k \omega t) + b_k \cdot \sin(k \omega t)).

Die Kreisfrequenz \omega skaliert hierbei die Periode 2\pi von Sinus und Kosinus auf die entsprechende Periode T. In der praktischen Anwendung wird man die Reihe häufig nach endlich vielen Reihengliedern abbrechen. Man erhält dann nur eine Approximation von f in Form eines trigonometrischen Polynoms,


  f_n(t)=\frac{a_0}{2} + \sum_{k=1}^n (a_k \cdot \cos(k \omega t) + b_k \cdot \sin(k\omega t)).

Diese endliche Summe wird dann Teilsumme f_n(t) der Fourierreihe genannt. Das so entstehende trigonometrische Polynom ist, unter allen trigonometrischen Polynomen der gleichen Struktur, dasjenige mit minimalem mittlerem quadratischem Fehler zur ursprünglichen Funktion f.

Die Koeffizienten der Entwicklung von f sind:

\begin{align}
  a_k&=\frac{2}{T}\int_{c}^{c+T} f(t) \cdot \cos(k\omega t)\, \mathrm{d}t \\[.7em]
  b_k&=\frac{2}{T}\int_{c}^{c+T} f(t) \cdot \sin(k\omega t)\, \mathrm{d}t
\end{align}

Das c stellt eine Verschiebung des Intervalls dar und kann zur Vereinfachung beliebig gewählt werden.

\displaystyle\frac{a_0}{2}=\frac{1}{T}\int_{c}^{c+T} f(t)\, \mathrm{d}t ist der Gleichanteil (wechsellose Größe)

Einfache Eigenschaften dieser Entwicklung sind, dass

  • b_k=0 für alle k gilt, falls f gerade ist, f(-x)=f(x),
  • a_k=0 für alle k gilt, falls f ungerade ist, f(-x)=-f(x).

Ist f gerade, dann sind also alle b_k=0, und die Koeffizienten a_k können auch über \textstyle a_k=\frac{4}{T} \int_0^{\frac{T}{2}} f(t) \cdot \cos(k\omega t)\, \mathrm{d}t berechnet werden. Dies ist möglich, weil durch die Symmetrie des Kosinus und der Funktion die Werte des Integrals in beiden Halbintervallen [-\tfrac T2,\,0] und [0,\,\tfrac T2] gleich sind. So ergeben sich oft Vereinfachungen. Analog gilt dies auch für ungerades f, das heißt bei a_k=0.

Ist die zugrundeliegende Funktion unbekannt bzw. liegen nur gegebene diskrete Daten (zum Beispiel Messwerte) vor, werden a_k, b_k nur aus den Stützpunkten approximiert (Trigonometrische Interpolation).

Amplituden-Phasen-Notation[Bearbeiten]

In der obigen Darstellung wird das Signal mit Hilfe eines Sinusspektrums und eines Kosinusspektrums dargestellt. Es ist aber auch eine Darstellung mittels Phasenwinkel- und Amplitudenspektrums möglich, da man die additive Überlagerung (Interferenz) einer Sinus- und einer Kosinusschwingung auch als phasenverschobene Kosinusschwingung darstellen kann:

\displaystyle
  f(t)=\frac{a_0}{2} + \sum_{n=1}^\infty (A_n \cos(n\omega t - \varphi_n))

Dabei errechnet sich A_n durch:

\displaystyle
  A_n=\sqrt{a_n^2+b_n^2}
.

Die Berechnung von \varphi_n ist aufwändiger. Wählt man Repräsentanten aus [-\tfrac{1}{2}\pi,\tfrac{3}{2}\pi], so gilt mit der Signumfunktion \sgn:


\varphi_n=\begin{cases} \arctan{\frac {b_n} {a_n}}, & \mbox{wenn }a_n > 0 \\ \arctan{\frac {b_n} {a_n}}+\pi, & \mbox{wenn }a_n < 0 \\ \sgn(b_n)\cdot \frac {\pi} {2}, &\mbox{wenn }a_n=0 \end{cases}
.

Für die üblichere Wahl \varphi_n\in(-\pi,\pi] bietet sich hingegen die Verwendung des Arkustangens mit zwei Argumenten an:

\displaystyle
\varphi_n=\operatorname{atan2}(b_n,a_n)

Für A_n+a_n\neq0 vereinfacht sich dies zu:

\displaystyle
\varphi_n=2\arctan{\frac {b_n} {A_n+a_n}}

Für a_n=b_n=0 sind die obigen Funktionen nicht definiert – in diesem Fall ist aber ohnehin A_n=0. Also spielt \varphi_n keine Rolle.

\varphi_n zeigt in den Quadranten, in welchem auch der Punkt (a_n,b_n) liegt.

Komplexe Fourierreihe[Bearbeiten]

Man kann nun jedes Paar von Amplitude und Verschiebung als komplexe Zahl in Polarkoordinatendarstellung interpretieren. Damit lassen sich die beiden Spektren in eines überführen. Eine Vereinfachung von geraden bzw. ungeraden Funktionen wie im Reellen ist so jedoch nicht möglich.

\displaystyle f(t)=\sum_{n=-\infty}^{\infty} c_n \mathrm{e}^{\mathrm{i}n\omega t}

Dabei ist

\displaystyle c_n =\frac1T\int_{c}^{c+T} f(t) \mathrm{e}^{-\mathrm{i}n\omega t} dt

Häufig wird c = -T/2 gewählt. Damit ergibt sich

\displaystyle c_n =\frac1T\int_{-\frac{T}{2}}^{\frac{T}{2}} f(t) \mathrm{e}^{-\mathrm{i}n\omega t} dt


Die Berechnung ist oft einfacher, da zum einen die e-Funktion leicht zu integrieren ist und zum anderen nur noch ein Koeffizient, c_n, statt zwei zu berechnen ist. Bei c_n handelt es sich in dieser Darstellungsform um eine sogenannte komplexe Amplitude, welche auch die Phaseninformation enthält. Das diskrete und komplexe Amplitudenspektrum stellt ein Linienspektrum dar, der Abstand zweier Spektrallinien \Delta f entspricht dem reziproken Wert der Periodendauer:

\Delta f = \frac{1}{T}

Zusammenhang zwischen reellen und komplexen Fourierkoeffizienten[Bearbeiten]

Reell zu komplex:

\displaystyle c_n = \frac{(a_{-n} + \mathrm{i} b_{-n})}{2} \mbox{ bei } n<0
\displaystyle c_0 = \frac{a_0}{2}
\displaystyle c_n = \frac{(a_n - \mathrm{i} b_n)}{2} \mbox{ bei }n>0

Komplex zu reell:

\displaystyle a_0 = 2 \cdot c_0
\displaystyle a_n = c_n + c_{-n}\!
\displaystyle b_n = \mathrm{i} (c_n - c_{-n})\!

Zusammenhang mit der Fourier-Transformation[Bearbeiten]

Mit Fourierreihen lassen sich nur periodische Funktionen und ihr Spektrum beschreiben. Um auch nichtperiodische Funktionen spektral beschreiben zu können, führt man einen Grenzübergang der Periode T \to \infty durch. Dadurch wird die Frequenzauflösung beliebig fein, was in einem Verschwinden des komplexen Amplitudenspektrums resultiert. Aus diesem Grund führt man das komplexe Amplitudendichtespektrum F ein, ausgehend von der komplexen Fourierreihe zunächst für die diskreten Argumente \omega=n\tfrac{2\pi}{T}:


F(\omega) = \int_{c-\frac{T}{2}}^{c+\frac{T}{2}} f(t) \mathrm{e}^{-\mathrm{i}\omega t} dt\,.

Durch Bildung des Grenzwertes T \to \infty (wobei gleichzeitig n\to\infty) folgt damit unmittelbar die Fourier-Transformation:

\mathcal{F}(f)(\omega) = \lim_{T \to \infty} F(\omega) = \int_{-\infty}^{\infty} f(t) \mathrm{e}^{-\mathrm{i}\omega t} dt\,.

Beispiele[Bearbeiten]

Dreieckpuls[Bearbeiten]

Verschiedene Näherungen eines Dreieckpulses

Die Dreieckfunktion lässt sich je nach gewünschter Phasenlage mit Sinus- und Kosinustermen approximieren. Mit dem Scheitelwert h lauten die Fourierreihen:

\begin{array}{rl}
f(t)
=& -\frac{8h}{\pi^2}\left[ {\cos{\omega t} + \frac{1}{3^2} \cos{3 \omega t} + \frac{1}{5^2} \cos{5 \omega t} + \cdots}\right] \\[.6em]
=& -\frac{8h}{\pi^2} \sum\limits_{k=1}^\infty \dfrac{ \cos ((2k-1) \omega t)}{(2k-1)^2} 
\end{array}
\begin{array}{rl}
f(t)
=& \frac{8h}{\pi^2}\left[ {\sin {\omega t} - \frac {1}{3^2}\sin{3 \omega t} + \frac {1}{5^2}\sin {5 \omega t} \mp \cdots}\right] \\[.6em]
=& \frac {8h}{\pi^2} \sum\limits_{k=1}^\infty (-1)^{k-1} \dfrac{ \sin((2k-1) \omega t)}{(2k-1)^2} 
\end{array}

Rechteckpuls[Bearbeiten]

Verschiedene Näherungen eines Rechteckpulses

Die Rechteckschwingung ist definiert durch

f(t) = \begin{cases} h, & \mbox{ wenn } 0\leq  t <T/2 \\ -h, & \mbox{ wenn } T/2 \leq t < T \end{cases}  \qquad f(t + T) = f(t)

Die Fourierreihe dazu lautet

\begin{align}
f(t)=& \tfrac{4h}{\pi}\left[\sin \omega t+\tfrac{1}{3}\sin3 \omega t+\tfrac {1}{5}\sin5 \omega t+\tfrac{1}{7}\sin7 \omega t+ \cdots\right]\\
=& \tfrac{4h}{\pi} \sum_{k=1}^\infty \dfrac{ \sin\left( (2k-1)\omega t \right) }{2k-1}
\end{align}

Anhand dieser Funktion erkennt man, dass man eine Rechteckschwingung durch unendlich viele Harmonische darstellen kann. Sie enthält jeweils die ungeraden harmonischen Oberschwingungen, wobei die Amplitude mit steigender Frequenz abnimmt. Aufgrund dessen wird ein Rechtecksignal auch häufig zum Testen elektronischer Schaltungen genommen, da so das Frequenzverhalten dieser Schaltung erkannt wird.

Allgemein enthalten alle periodischen Schwingungen mit der Periodendauer T der Grundschwingung und beliebigem Verlauf innerhalb der Periode nur ungeradzahlige Oberschwingungen, wenn gilt:

f(t + \tfrac T2) = -f(t)
Fourier-Synthese eines Rechtecksignals

Im rechten Bild ist die Fourier-Synthese einer Rechteckschwingung dargestellt. Die Diagramme der ersten Spalte zeigen diejenige Schwingung, die in der jeweiligen Zeile hinzugefügt wird. Die Diagramme in der zweiten Spalte zeigen alle bisher berücksichtigten Schwingungen, die dann in den Diagrammen der dritten Spalte addiert werden, um dem zu erzeugenden Signal möglichst nahe zu kommen. Die Schwingung aus der ersten Zeile nennt sich Fundamentalschwingung, alle weiteren, die hinzugefügt werden, sind Oberschwingungen (Harmonische). Je mehr solcher Vielfache der Grundfrequenz berücksichtigt werden, umso näher kommt man einem idealen Rechtecksignal. An den unstetigen Stellen des Rechtecksignals bildet sich durch die Fourier-Synthese bedingt ein so genannter Überschwinger, der auch bei größerer Approximation nicht verschwindet. Diese Erscheinung wird gibbssches Phänomen genannt, sie weist eine konstante und von der Bandbreite unabhängige Überschwingung von etwa 9 % des vollen Sprungs auf. Die vierte Spalte zeigt das Amplitudenspektrum normiert auf die Grundschwingung.

Sägezahnpuls (steigend)[Bearbeiten]

Verschiedene Näherungen eines Sägezahnpulses

Ebenso lassen sich punktsymmetrische Funktionen aus Sinustermen approximieren. Hier erreicht man eine Phasenverschiebung durch alternierende Vorzeichen:

\begin{array}{rl}
f(t)
=&- \frac{2h}{\pi}\left[ {\sin {\omega t} - \frac {1}{2}\sin{2 \omega t} + \frac {1}{3}\sin {3 \omega t} \mp \cdots}\right] \\[.6em]
=& - \frac {2h}{\pi}\sum\limits_{k=1}^{\infty}(-1)^{k-1} \dfrac {\sin k \omega t}{k} 
\end{array}

Sinuspuls[Bearbeiten]

Verschiedene Näherungen eines Sinuspulses
\begin{array}{rl}
f(t)
=& h\left| \sin {\omega t} \right|\\[.6em]
=& \frac{4h}{\pi}\left[ \frac{1}{2} - \frac { \cos {2 \omega t}}{3}-\frac { \cos {4 \omega t}}{15}-\frac { \cos {6 \omega t}}{35}- \cdots\right] \\[.6em]
=& \frac{2h}{\pi} - \frac{4h}{\pi} \sum\limits_{k=1}^{\infin} \dfrac { \cos {2 k\omega t}}{(2k)^2-1}
\end{array}

Konvergenzaussagen zur Fourierreihe[Bearbeiten]

Man kann zwar bedenkenlos zu einer periodischen Funktion eine Fourierreihe aufstellen, jedoch muss diese Reihe nicht konvergieren. Ist dies der Fall, so erhält man durch diese Transformation auch keine weiteren Informationen. Konvergiert die Reihe, so muss man sich im Klaren sein, in welchem Sinn die Konvergenz vorliegt. Meistens untersucht man Fourierreihen auf punktweise Konvergenz, gleichmäßige Konvergenz oder auf Konvergenz bezüglich der L^2-Norm. Im Folgenden werden einige wichtige Sätze über die Konvergenz von Fourierreihen aufgezählt.

Satz von Dirichlet[Bearbeiten]

Hauptartikel: Dirichlet-Bedingung

Peter Gustav Lejeune Dirichlet bewies, dass die Fourierreihe einer differenzierbaren, 2 \pi periodischen Funktion punktweise gegen die Ausgangsfunktion konvergiert.[1] Unter der Voraussetzung, dass f sogar stetig differenzierbar ist, kann die Aussage noch verbessert werden.

Sei f \colon \R \to \R eine stetig differenzierbare, 2 \pi-periodische Funktion, dann konvergiert die Fourierreihe von f gleichmäßig gegen f.

Satz von Carleson[Bearbeiten]

Der Satz von Carleson ist ein tiefliegendes Resultat zur Konvergenz einer Fourierreihe.

Sei f \in L^2([-\pi,\pi]) eine quadratintegrierbare Funktion, dann konvergiert die Fourierreihe bezüglich f fast überall.

Diese Aussage ist sogar für alle L^p-Räume mit p \in ]1,\infty[ richtig und heißt in dieser allgemeinen Form Satz von Carleson–Hunt. Dass die Aussage für p=1 falsch ist, konnte Kolmogorov 1923 durch ein Gegenbeispiel zeigen. Nikolai Nikolajewitsch Lusin vermutete schon 1915 die Richtigkeit des Satzes von Carleson, konnte sie jedoch nicht beweisen. Der Beweis gelang erst Lennart Carleson im Jahr 1964.[2]

Satz von Fejér[Bearbeiten]

Hauptartikel: Satz von Fejér

Leopold Fejér bewies, dass das arithmetische Mittel der Partialsummen der Fourierreihe einer stetigen, 2 \pi-periodischen Funktion gleichmäßig gegen die Funktion konvergieren.

Sei f \colon \R \to \R eine stetige, 2 \pi-periodische Funktion und \textstyle S(f) = \sum_{k = - \infty}^\infty c_k e^{ikx} die Fourierreihe von f. Mit \textstyle s_n(f) = \sum_{k = - n}^n c_k e^{ikx} wird die n-te Partialsumme dieser Reihe beschrieben. Dann besagt der Satz von Fejér, dass die Partialsummen \textstyle \frac{1}{m+1}\sum_{n = 1}^m s_n(f) gleichmäßig gegen f konvergieren. Es gilt also

\lim_{m \to \infty} \frac{1}{m+1} \sum_{n = 1}^m s_n(f) = \lim_{m \to \infty} \frac{1}{m+1}\sum_{n = 1}^m \sum_{k = - n}^n c_k e^{ikx} = f(x),

wobei die Konvergenz gleichmäßig ist.

Gibbssches Phänomen[Bearbeiten]

Gibbssches Phänomen bei einer Rechteckkurve
Hauptartikel: Gibbssches Phänomen

In der Umgebung von Sprungstellen konvergiert die Fourierreihe nicht mehr gleichmäßig, sondern nur noch punktweise. Es entstehen dort in den Partialsummen der Reihe typische Über- und Unterschwinger von etwa 9 % der Sprunghöhe. Dieser Effekt hat weitreichende Auswirkungen in der Signalverarbeitung.

Verallgemeinerte Fourier-Reihe[Bearbeiten]

Sei (H, \langle ., .\rangle_H) ein Hilbertraum mit einer Orthonormalbasis S. Dann kann man jedes Element f \in H des Hilbertraums durch

f=\sum_{u \in S} \langle f, u \rangle_H\, u

darstellen. Diese Reihendarstellung wird auch (verallgemeinerte) Fourier-Reihe genannt.

Verallgemeinerungen[Bearbeiten]

Verallgemeinerungen der Fourier-Reihe, die sich zwar auch als Darstellungen in Orthonormalbasen beschreiben lassen, aber zusätzlich ähnlich der Fourier-Reihe bestimmte Eigenschaften in Bezug zu Symmetrien aufweisen, untersucht die harmonische Analyse. Die Pontrjagin-Dualität verallgemeinert dabei die Fourier-Reihe auf Funktionen auf beliebigen abelschen lokalkompakten topologischen Gruppen, der Satz von Peter-Weyl auf kompakten topologischen Gruppen.

Weblinks[Bearbeiten]

Einzelnachweise[Bearbeiten]

  1. St. Goebbels, St. Ritter: Mathematik verstehen und anwenden - von den Grundlagen bis zu Fourier-Reihen und Laplace-Transformation. Spektrum, Heidelberg 2011, ISBN 978-3-8274-2761-8, S. 696, 704 - 706
  2.  S. A. Telyakovskii: Carleson theorem. In: Michiel Hazewinkel (Hrsg.): Encyclopaedia of Mathematics. Springer-Verlag, Berlin 2002, ISBN 1-4020-0609-8 (Online).

Quelle[Bearbeiten]

  •  Konrad Königsberger: Analysis 1. Springer, Berlin 2004, ISBN 3-540-41282-4 (Kapitel 16).