Techniklexikon

Fourier-Analyse

Mathematische Methoden und Computereinsatz

Vermutlich werden alle Wissenschaftler und Ingenieure, zu diesem Thema befragt, der Ansicht zustimmen, dass die »Fourier-Analyse« eine Angelegenheit von grosser Bedeutung sei; jedoch fragt man dann, auf welche Weise sie bedeutend sei, so muss man sich auf eine Vielzahl verschiedener Antworten einstellen. Für den reinen Mathematiker war der Zeitpunkt, zu dem Fouriers Arbeiten entstanden, nämlich der Beginn des 19. Jh., ein Moment, in dem sich eine tiefgreifende Änderung im Verständnis der Grundlagen der Mathematik vollzog. Fouriers Herangehen an die Probleme der Wärmeleitung - Probleme, die wir aus heutiger Sicht eher dem Bereich der angewandten Mathematik oder den Ingenieurwissenschaften zuweisen würden - hat Fragen von fundamentaler Bedeutung für die gesamte Mathematik aufgeworfen; gleichzeitig lieferte Fouriers Ansatz ein bleibendes und mächtiges Werkzeug, das noch heute, zwei Jahrhunderte später, den Mathematikern und den Ingenieuren eine grosse Hilfe ist.

 

Wärmeleitung in einem dünnen Stab

Wir wollen mit einer Demonstration der Nützlichkeit dieses Werkzeugs beginnen. Die Wärmeleitung in einem Festkörper wird durch die Diffusionsgleichung

,                  (1)

beschrieben, in der k die thermische Diffusivität des Festkörpers und q eine Funktion ist, welche die Temperatur an einem beliebigen Punkt zu einer beliebigen Zeit t angibt. Der Einfachheit halber betrachten wir die Wärmeleitung in nur einer Raumdimension. Man nehme zum Beispiel einen zylindrischen Metallstab, der überall bis auf seine Endflächen isoliert ist, und der so dünn ist, dass kein merklicher Temperaturgradient quer zur Längsausdehnung auftreten kann. Wir wollen annehmen, dass sich der Stab von x = 0 zu x = p erstreckt; es mag unsinnig erscheinen, die Länge des Stabs als p zu wählen, aber es geht hier nur darum, eine passende Längeneinheit für unser Problem festzulegen.

Man sieht leicht (beispielsweise durch Einsetzen in die Gleichung), dass die folgenden Ausdrücke

Lösungen von Gleichung (1) sind, in denen m eine freie Konstante ist. Aus der Linearität der Diffusionsgleichung folgt, dass wir weitere Lösungen durch Linearkombination obiger Lösungen gewinnen können. Diese einfache Beobachtung führt uns direkt zur Fourierschen Methode.

Wir wollen davon ausgehen, dass wir unseren Stab dazu gebracht haben, eine einheitliche Gleichgewichtstemperatur q₀ anzunehmen; dann, zur Zeit t = 0, tauchen wir ihn plötzlich in ein kaltes Bad ein, das in geeigneten Einheiten eine Temperatur von null aufweisen soll. Dies kann mathematisch durch

                   (4)

dargestellt werden. Dieser Ausdruck gibt zusammen mit den Randbedingungen

                   (5)

die anfängliche Temperaturverteilung wieder. Die Randbedingungen legen für alle Zeiten fest, dass die Temperatur an den Enden des Stabes gleich der Temperatur des Bades ist; dies entspricht der Annahme, dass das Bad gut durchmischt ist und eine sehr grosse Wärmekapazität aufweist.

 

Die Darstellung von Funktionen durch trigonometrische Reihen

Wir suchen nach einer Linearkombination von Lösungen der Form (2) und (3), welche die durch Gleichung (4) und (5) formulierten Randbedingungen erfüllt. Fouriers kühne Behauptung war nun, dass beliebige Anfangsbedingungen mit solchen Summen erfüllt werden können, wenn man m über alle positiven ganzen Zahlen variieren lässt. Daher schreiben wir, indem wir die unbekannten Parameter  und  einführen:

                   (6)

und

                                                                               (7)

Der kritische Leser mag sich an diesem Punkt fragen, inwieweit es sinnvoll ist, mit unendlichen Summen von Funktionen zu hantieren, und mag der Behauptung, beliebige Funktionen könnten in dieser Weise dargestellt werden, skeptisch gegenüberstehen; wir werden uns diesen berechtigten Fragen später wieder zuwenden. Im Moment jedoch wollen wir untersuchen, inwieweit eine solche Darstellung bei der Lösung unseres Problems hilfreich sein könnte. Wenn wir die Bildung unendlicher Summen von Funktionen als möglichen Lösungsansatz akzeptieren, und wenn es uns gleichfalls sinnvoll erscheint, diese Summen gliedweise abzuleiten, so lässt sich durch Substitution von Gleichung (7) in Gleichung (1) zeigen, dass es sich tatsächlich um eine Lösung handelt. Darüber hinaus können wir die Einhaltung der Randbedingungen garantieren, indem wir alle Koeffizienten  zu null setzen, so dass

                               (8)

und

                   (9)

Da sin(0) = 0 und sin(kp) = 0 gilt, ist Gleichung (5) automatisch erfüllt.

Aber was ist mit der Anfangsbedingung an q(x, 0)? Fouriers Idee war es, die unbekannten Koeffizienten  so zu wählen, dass Gleichung (8) die Form von Gleichung (4) annimmt. Dies ist der Kerngedanke der Fourier-Analyse. Es geht an dieser Stelle ja nicht um die Lösung von Differentialgleichungen, sondern lediglich um die Darstellung einer Funktion als trigonometrische Reihe. In unserem Beispiel handelt es sich um eine Funktion, die auf dem offenen Intervall 0 < x < p einen von null verschiedenen konstanten Wert annimmt, jedoch an den Punkten x = 0 und x = p verschwindet. Der Gedanke, ein solches Objekt als eine einzelne Funktion zu behandeln oder gar als eine Summe von Sinuswellen darzustellen, erschien zu Fouriers Zeit als abwegig. Damals war die Erörterung der Frage, was die Eigenschaften einer Funktion ausmacht, untrennbar mit den Begriffen der Stetigkeit und der Differenzierbarkeit verknüpft, die jedoch ihrerseits noch nicht sauber definiert waren. Es dauerte noch über zwanzig Jahre, bis der strenge Beweis eines Theorems, das Fouriers Behauptung enthielt, veröffentlicht wurde.

 

Der Dirichletsche Satz

Dieses von Dirichlet gefundene Theorem ist, obgleich seither allgemeinere Ergebnisse bewiesen werden konnten, für viele Anwendungen ein nützliches Werkzeug. Dirichlet zeigte, dass unter der Voraussetzung, dass eine auf dem Intervall [-p, p] definierte Funktion F stückweise stetig und stückweise monoton ist, eine eindeutig festgelegte Fourier-Reihe existiert, die an jedem Punkt x, an dem F stetig ist, gegen F(x) konvergiert. An jeder Unstetigkeitsstelle konvergiert die Fourier-Reihe gegen den Mittelwert des rechtsseitigen und des linksseitigen Grenzwertes von F an dieser Stelle.

 

Wie man eine Fourier-Reihe konstruiert

Der direkteste und eleganteste Weg, die passenden Fourier-Koeffizienten zu einer gegebenen Funktion F zu finden, besteht darin, die Orthogonalitätsrelationen der trigonometrischen Funktionen auszunutzen, die für Sinus-Funktionen von der folgenden Form sind:

                   (10)

Analoge Relationen gelten für die Kosinus-Funktion. Darüber hinaus lässt sich auch

                   (11)

zeigen. Zu Demonstrationszwecken wollen wir in einem kleinen Exkurs folgende auf dem Intervall -p < x < p definierte Funktion betrachten:

Diese Funktion, die aus einsichtigen Gründen als Stufenfunktion bezeichnet wird, ist offenbar stückweise stetig und monoton. Sie weist genau eine Unstetigkeitsstelle, nämlich bei x = 0, auf, und sie ist so definiert, dass ihr Wert an dieser Stelle gleich dem Mittelwert aus rechts- und linksseitigem Grenzwert ist. Daher können wir die Stufenfunktion unter Ausnutzung des Dirichletschen Satzes als Fourier-Reihe schreiben:

       (12)

Wenn wir beide Seiten dieser Gleichung mit  multiplizieren und anschliessend über x von -p bis p integrieren, wobei aus der Orthogonalitätsrelation folgt, dass sämtliche der unendlich vielen Integrale auf der rechten Seite bis auf eines verschwinden, so ergibt sich:

                               (13)

In analoger Weise, nämlich durch Multiplikation mit  und anschliessender Integration, bekommt man einen Ausdruck für :

                               (14)

Das Schöne an diesem Verfahren ist, dass es explizite Ausdrücke für die Entwicklungskoeffizienten in Form von Integralen über die bekannte Funktion f liefert.

Für die Stufenfunktion sind diese Integrale besonders leicht auszuführen, da es sich bei ihr um eine ungerade Funktion handelt, die dementsprechend die Symmetriebeziehung f(x) = -f(-x) erfüllt. Damit erhält man:

,                                          (15)

woraus folgt, dass alle Koeffizienten  verschwinden, so dass die in Gleichung (12) dargestellte Reihe nur aus Sinus-Termen besteht. Auf ähnliche Weise findet man

                   (16)

Unsere Fourier-Reihe für f(x) reduziert sich also auf eine Summe von Sinus-Funktionen, deren Argumente ungeradzahlige Vielfache von x sind:

                   (17)

In der Praxis wäre es oftmals wünschenswert, diese unendliche Summe durch ihre ersten (möglichst wenigen) Glieder annähern zu können. Wie sinnvoll dies ist, hängt davon ab, wie schnell die Reihe konvergiert - dass die Reihe überhaupt konvergiert, und zwar gegen die richtige Funktion, wissen wir bereits aus dem Dirichletschen Satz. Aus [23]Abb. 1, in der verschiedene Näherungen von Gleichung (17) dargestellt sind, können wir entnehmen, dass die Reihe doch eher langsam konvergiert - selbst wenn man die Glieder bis zu k = 30 berücksichtigt, zeigt die Näherung besonders in der Nähe der Singularität bei x = 0, aber auch noch in der Nähe von x = -p und x = p, starke Oszillationen.

Da trigonometrische Funktionen auf dem gesamten Bereich der reellen Zahlen wohldefiniert sind, kann jede beliebige Fourier-Reihe wie eine Funktion von ganz  behandelt werden; dies sogar dann, wenn man sie damit über den Definitionsbereich der Funktion, die dargestellt werden soll, hinaus erweitert. Dieser Aspekt impliziert, dass jede Fourier-Reihe (und damit auch die zugrundeliegende Funktion f) periodisch sein muss; in unserem Beispiel ergibt sich f(x + 2p) = f(x) aus der Periodizität der Sinus-Funktionen in den Gliedern der Reihe. In diesem Sinne stellt sich die Fourier-Reihe, obwohl sie auf dem beschränkten Definitionsbereich der Originalfunktion f wie eine Stufenfunktion aussieht, auf dem erweiterten Definitionsbereich als eine periodische, stückweise konstante Funktion dar. Letztere weist zusätzlich zu der Unstetigkeitsstelle bei x = 0 eine solche bei x = p und, dies erfordert die Periodizität, auch bei x = -p auf. Man kann aus dem Dirichletschen Satz oder auch schlicht aus der Tatsache, dass sin((2k + 1)p) = 0 gilt, folgern, dass f(p) = 0 gilt.

 

Die Lösung der Wärmeleitungsgleichung

Dass wir diese Überlegungen bei der Lösung des Wärmeleitungsproblems gebrauchen können, erkennt man, wenn man die Funktion f durch ihre Fourier-Reihe darstellt, deren Definitionsbereich man auf das Intervall 0 £ x £ p einschränkt. Dann nämlich kann man die Anfangsverteilung der Temperatur als

                   (18)

schreiben. Daraus ergibt sich, dass die gesamte Zeitentwicklung der Temperaturverteilung des Stabes durch

       (19)

gegeben ist. Diese Lösung ist in [24]Abb. 2 dargestellt; offenbar verliert die Verteilung ihre anfängliche charakteristische rechteckige Form sehr schnell. Dies hat mit der langsamen Konvergenz der Fourier-Reihe für f(x) zu tun: Jede Näherung der vollen Fourier-Reihe durch endlich viele Glieder ist ein Versuch der Darstellung einer unstetigen Funktion mit einer glatten, denn der Gradient einer solchen Näherung kann eine Grenze gemäss

       (20)

nicht überschreiten. Intuitiv erwarten wir, dass die Unstetigkeiten von f(x) durch solche Stellen der endlichen Reihe angenähert werden, die einen sehr steilen Verlauf aufweisen; um dies zu erreichen, brauchen wir grosse K, also schnell oszillierende Sinus-Funktionen. Für die Lösung der Diffusionsgleichung (19) gilt, dass die Sinus-Funktion um so schneller bei der zeitlichen Entwicklung abfällt, je höher ihre Frequenz ist; daher glättet sich die eckige Form der Anfangsverteilung schnell.

 

Weitere Entwicklungen und Anwendungen

Wir haben bereits weiter oben darauf hingewiesen, dass es heikel ist, unendliche Summen wie endliche zu behandeln. Um die Fourier-Analyse zur Lösung von Differentialgleichungen verwenden zu können, müssen wir wissen, ob es erlaubt ist, die Reihe gliedweise abzuleiten. Das ist offenbar nicht selbstverständlich; die Fourier-Reihe der Stufenfunktion beispielsweise ist nicht einmal stetig, geschweige denn stellt sie eine Summe differenzierbarer Funktionen dar. Damit eine Fourier-Reihe stetig ist, reicht es jedoch schon hin, wenn sie gleichmässig konvergent ist - dies ist im wesentlichen eine Bedingung, die sicherstellt, dass die Reihe über den ganzen Definitionsbereich schnell genug konvergiert. Damit eine Fourier-Reihe differenzierbar ist, genügt es, dass die Reihe, die durch gliedweise Ableitung der Originalreihe entsteht, gleichmässig konvergiert. Zu der Zeit, zu der diese Gedanken ausgearbeitet wurden, herrschte eine beträchtliche Verwirrung bezüglich des Zusammenhangs zwischen Stetigkeit und Differenzierbarkeit. Man hatte irgendwie das Gefühl, dass eine stetige Funktion nur an wenigen besonderen Stellen nicht differenziert werden kann. Diese Illusion zerschlug sich, als Weierstrass 1872 auf folgende Funktion hinwies, die explizit als trigonometrische Reihe definiert ist:

                   (21)

Ist a eine ungerade ganze Zahl, und gilt 0 < b < 1 und ab > 1 + 3p/2, so ist diese Funktion für alle x stetig, jedoch nirgendwo differenzierbar. Aus der Sicht heutiger Entwicklungen der Geometrie fraktaler Mengen und ihrer Anwendung auf die Beschreibung von Naturphänomenen müssen wir den Graphen dieser Funktion als ein frühes Beispiel eines mathematisch konstruierten Fraktals betrachten.

Die Fourier-Analyse hat aufgrund der Nützlichkeit des Verfahrens beim Lösen linearer Differentialgleichungen einen weiten Anwendungsbereich. Die Diffusionsgleichung (1) und ihr verwandte Gleichungen tauchen in so verschiedenartigen Gebieten wie der elektrischen Leitfähigkeit, der selbständigen Durchmischung von in Flüssigkeiten und Gasen gelösten Materialien, der Wahrscheinlichkeitstheorie und neuerdings sogar bei der Berechnung der zukünftigen Entwicklungen von Finanzmärkten auf. Der Diffusionsgleichung sehr nah verwandt ist die Schrödinger-Gleichung, die von zentraler Bedeutung für die Quantenmechanik ist. Fourier-Reihen können eingesetzt werden, um die Wellengleichung zu lösen, die von der mathematischen Nachbildung schwingender Saiten und Orgelpfeifen bis hin zur Beschreibung der Ausbreitung elektromagnetischer Wellen nahezu überall Anwendung findet. Selbst bei der Behandlung nichtlinearer Gleichungen kann die Fourier-Analyse von grossem Nutzen sein, weil es oftmals möglich ist, die Stabilität einer Lösung der vollen nichtlinearen Gleichung zu überprüfen, indem man lineare Gleichungen ableitet, welche die Entwicklung kleiner Störungen des gegeben Zustands beschreiben. In solchen Analysen kann man auf Lösungen, die der Form von Gleichung (7) ähneln, stossen, wobei jedoch einige der Koeffizienten exponentiell ansteigen. Solche Lösungen bilden eine neue Generation räumlicher Strukturen auf dem Gebiet nichtlinearer Gleichungen.

Schliesslich wollen wir darauf hinweisen, dass es noch weit mehr mathematische Techniken gibt, die sich mit Fouriers Namen verbinden, als hier im zur Verfügung stehenden Raum auch nur erwähnt werden kann. Von diesen vielleicht am weitesten verbreitet ist das Fourier-Integral, das als Erweiterung der Fourier-Reihe verstanden werden kann. Wenn wir den ganzzahligen Index k in Gleichung (12) als reelle Variable uminterpretieren und die Summenzeichen durch Integrale ersetzen, so gelangen wir zu einem Konstruktionsverfahren für Funktionen f(x), die nicht mehr zwingend periodisch sein müssen. Bei dieser Verallgemeinerung sind die Fourier-Koeffizienten  und  Funktionen der reellen Variablen k, die als Fourier-Transformierte von f(x) bezeichnet werden. Diese Verallgemeinerung ist bei der Behandlung von Differentialgleichungen über einem unendlichen Definitionsbereich, bei der Zeitreihenanalyse und in der Spektroskopie von grösster Bedeutung.

Ganz zum Schluss sollten wir uns vielleicht zurücklehnen und einen Moment darüber nachsinnen, dass die Geschichte der Fourier-Reihen ein sehr schönes Beispiel für das ist, was die Mathematik und die Wissenschaften allgemein in den letzten drei Jahrhunderten hat so erfolgreich sein lassen. Es ist eine Geschichte darüber, wie Anwendungen die Entwicklung der reinen Mathematik vorantreiben und dann wiederum die Entwicklungen in der reinen Mathematik zu neuartigen, oftmals unerwarteten Anwendungen führen.

 

Fourier-Analyse 1: Darstellung vier verschiedener Näherungen von Gleichung (17), an denen man die Auswirkung der Hinzunahme weiterer Glieder beobachten kann. Man beachte, dass der Graph um so mehr Oszillationen aufweist, je mehr Glieder berücksichtigt werden.

bottom:.0001pt\'>Fourier-Analyse 2: Eine Darstellung von q(x, t) für den abkühlenden Stab.