DFT – Diskrete Fourier-Transformation
Im folgenden Artikel wird die DFT (Diskrete Fourier Transformation) erklärt und deren Verortung in der Fourier-Analysis dargelegt. Des Weiteren werden wichtige Eigenschaften der DFT gezeigt und außerdem die IDFT (Inverse DFT) erläutert.
Falls dir das zu ausführlich ausfällt, haben wir das Wichtigste zum Thema Diskrete Fourier Transformation für dich in einem anschaulichen Video zusammengefasst.
Inhaltsübersicht
DFT, Fourier Transformation und Fourier Reihen
Die DFT (Diskrete Fourier Transformation) ist eine für die Praxis relevante Transformation aus dem Bereich der Fourier-Analysis.
Die DFT ist mit der kontinuierlichen Fourier-Transformation verwandt und lässt sich aus Überlegungen zu Fourier-Reihen ableiten.
Fourier-Transformation
Mithilfe der Fourier Transformation lassen sich nichtperiodische Funktionen durch eine Linearkombination aus Sinus- und Cosinus-Funktionen verschiedener Frequenzen darstellen. Hierfür ordnet die Fourier-Transformierte jeder Frequenz diejenige Amplitude zu, mit der die Schwingung dieser Frequenz in der Linearkombination auftritt:
Die nichtperiodischen Funktionen können dann durch das Fourier-Integral – die sogenannte inverse Fourier-Transformierte – genähert werden:
Hierbei fasst die Exponentialfunktion die Sinus- und Cosinus-Schwingungen zusammen.
Fourier-Reihen
Für periodische Funktionen f mit der Periode T sind die in der Linearkombination beteiligten Frequenzen diskret. Sie sind ein n-faches der Grundfrequenz und werden mit bezeichnet. Aus dem Integral wird damit die Fourier-Reihe ,
wobei die Amplituden das Analogon zur Fourier-Transformierten darstellen. Sie lassen sich wie folgt berechnen:
,
Durch die Bestimmung dieser Amplituden bzw. der Fourier-Transformierten werden die betrachteten Funktionen in ihr sogenanntes Frequenzspektrum zerlegt.
Definition: Diskrete Fourier Transformierte
Mithilfe der Fourier-Transformation und der Fourier-Reihe lassen sich also die Frequenzspektren von nichtperiodischen und periodischen Funktionen ermitteln. In der Praxis liegen die zu untersuchenden Signalen allerdings nicht als Funktion vor, sondern als Menge diskreter Werte. Die Diskrete Fourier Transformation befasst sich mit der Spektralanalyse solcher diskreten endlichen Signale, welche zur Untersuchung periodisch fortgesetzt werden.
Untersucht werden sollen die komplexen Werte , die die Signalwerte darstellen und in einem Vektor zusammengefasst werden können.
Die Diskrete Fourier Transformierte des Vektors besitzt die Koeffizienten
Zur Effizienten Berechnung der Diskreten Fourier Transformation gibt die FFT (Fast Fourier Transformation) einen Algorithmus an.
Verwandtschaft zwischen diskreter Fourier-Transformation und Fourier-Reihen
Die DFT ist eng verwandt mit der Theorie der Fourier-Reihen bzw. sie lässt sich aus dieser ableiten. Dazu werden die einzelnen Werte ,, des zu untersuchenden diskreten und endlichen Signals als äquidistante Funktionswerte einer -periodischen Funktion angesehen:
mit
Für eine solche -periodischen Funktion f lässt sich eine Fourier-Reihe entwickeln und die Fourier-Koeffizienten lassen sich wie bereits erwähnt durch folgendes Integral berechnen:
,
Dabei können die Grenzen des Integrals auch verschoben werden. Wichtig ist nur, dass über eine komplette Periodendauer integriert wird. Die Koeffizienten können also auch folgendermaßen berechnet werden:
,
Hierbei sind die Frequenzen ein Vielfaches der Grundfrequenz:
Im Fall der DFT sind allerdings nicht alle kontinuierlichen Funktionswerte bekannt, sondern nur die diskreten Messwerte . Das bedeutet, dass sich dieses Integral nicht berechnen lässt. Daher wird das Integral durch eine Riemann-Summe genähert: Die Funktion wird an den Stützstellen ausgewertet, mit der Differenz zwischen zwei Stützstellen multipliziert und dieses Produkt wird über den Laufindex k aufsummiert. Es ergibt sich die folgende Riemann-Summe:
Dies entspricht gerade der Berechnungsformel des n-ten Koeffizienten der Diskreten Fourier Transformierten des Vektors .
Betrachtet man die Berechnungsformel des Koeffizienten , so fällt auf, dass dieser dem Koeffizienten entspricht, weshalb die diskrete Fourier-Transformierte nur zwischen und definiert ist:
Von der kontinuierlichen zur Diskreten Fourier Transformation
Auch zur kontinuierlichen Fourier-Transformation besitzt die DFT große Ähnlichkeiten. Um das deutlich zu machen wird eine nichtperiodische Funktion betrachtet, die außerhalb des Intervalls verschwindet und auf diesem im Abstand jeweils als Funktionswert einen Eintrag aus annimmt:
, für
Sei nun die Zahl der untersuchten Werte N deutlich größer als die Länge des gewählten Intervalls , so lässt sich dort das Integral der Fourier-Transformierten der Funktion f sinnvoll durch eine diskrete Summe nähern:
Wird diese Summe nur für berechnet, so ergibt sich:
Das entspricht bis auf den konstanten Faktor T der Berechnungsformel des n-ten Eintrags der Diskreten Fourier Transformierten .
IDFT: Inverse Diskrete Fourier Transformation
Um zur Formel der Inversen DFT zu gelangen, kann die Fourier-Reihe derselben periodischen Funktion betrachtet werden, welche oben beschrieben wurde.
Im Fall der Diskreten Fourier Transformation werden die Koeffizienten durch die Einträge der Diskreten Fourier Transformierten repräsentiert. Da diese nur für definiert sind, läuft der Index in der Summe auch nur zwischen diesen beiden Grenzen. Berechnet man dementsprechend diese Reihe für die Werte und für die Frequenzen , so erhält man folgendes Ergebnis:
Definition: IDFT
Somit berechnen sich die Koeffizienten der Inversen Diskreten Fourier Transformierten von nach Umbenennung der Indizes auf folgende Art:
Zur Zusammenfassung: Die Koeffizienten der Diskreten Fourier Transformierten von lauten:
Für den Vorfaktor bei der Diskreten Fourier Transformation gibt es unterschiedliche Konventionen. Manchmal wird er auch bei der Inversen DFT statt bei der DFT eingefügt.
Überblick Fourier Analysis
Da nun die Diskrete Fourier Transforamtion und auch ihre Inverse bekannt sind, sollen beide noch einmal in den Kontext der Fourier-Analysis eingeordnet werden und ein kurzer Überblick über die drei wichtigsten Transformationen aus diesem Bereich gegeben werden. Dabei werden nun die in der Signalverarbeitung üblichen Bezeichnungen angewandt. Mit x soll ein Signal im Zeitbereich bezeichnet werden und mit X wird das entsprechende Signal im Frequenzbereich symbolisiert.
DFT Matrix
Wird die Abkürzung für die N-te Einheitswurzel eingeführt, so lauten die Formeln der Diskreten Fourier Transformation:
DFT
IDFT
Im Folgenden werden die einzelnen Komponenten der DFT einmal ausgeschrieben:
Dies lässt sich auch in Matrix-Schreibweise notieren:
Die dort auftretende Matrix
wird als DFT-Matrix oder Fourier-Matrix bezeichnet. Es lässt sich zeigen, dass für ihre Inverse folgender Zusammenhang gilt:
,
wobei die komplex konjugierte Matrix von Z beschreibt.
Somit lauten die Formeln der DFT und der IDFT:
DFT
IDFT
Eigenschaften der DFT
Folgende Eigenschaft der DFT wurde bereits gezeigt. Würde der Koeffizient berechnet werden, entspräche dieser dem Koeffizienten der Diskreten Fourier Transformierten:
Des Weiteren gewinnt man mit der DFT nur Informationen für die Frequenzen mit . Es gilt nämlich:
Werden demnach Werte einer Signalfunktion mit im Abstand , also mit der Abtastfrequenz ermittelt, so liefert die Diskrete Fourier Transformation nur Informationen zum Frequenzspektrum der Funktion für Frequenzen, die unter liegen.
Abtasttheorem
Diese Erkenntnis wird im Abtasttheorem von Shannon zusammengefasst:
Sei ein Signal, in dessen Frequenzspektrum die höchste auftretende Frequenz ist. Dann ist dieses Signal vollständig durch seine Abtastwerte bestimmt, falls die Abtastfrequenz größer als ist. Die Frequenz wird auch als Nyquist-Frequenz bezeichnet.
Linearität der DFT
Eine weitere leicht zu zeigende Eigenschaft der Diskreten Fourier Transformation ist die Linearität:
Diskrete Fourier Transformation Beispiel
In diesem Beispiel soll das Signal betrachtet werden und zu den Zeiten und abgetastet werden.
Die DFT ist also für den Vektor mit und durchzuführen. Der Vektor lautet somit:
Mit der Formel für die DFT berechnet sich die Diskrete Fourier Transformierte mit wie folgt:
Die Diskrete Fourier Transformierte des Vektors lautet also:
Das Ergebnis lässt sich auch mithilfe der Fourier-Matrix berechnen. Das sieht dann folgendermaßen aus:
Werden die passenden Werte eingesetzt, erhält man: