In der Mathematik ist die diskrete Fourier-Transformation ( DTFT ) eine Form der Fourier-Analyse , die auf eine Folge diskreter Werte anwendbar ist.

Die DTFT wird häufig zur Analyse von Abtastwerten einer stetigen Funktion verwendet . Der Begriff „diskrete Zeit“ rührt daher, dass die Transformation auf diskreten Daten operiert, oft auf Abtastwerten mit einem Zeitintervall. Aus gleichmäßig verteilten Abtastwerten erzeugt sie eine Frequenzfunktion, die die periodische Summe der kontinuierlichen Fouriertransformierten der ursprünglichen stetigen Funktion darstellt. Vereinfacht ausgedrückt: Wendet man die DTFT auf regelmäßig verteilte Abtastwerte eines kontinuierlichen Signals an, erhält man sich wiederholende (und möglicherweise überlappende) Kopien des Frequenzspektrums des Signals, deren Abstände der Abtastfrequenz entsprechen. Unter bestimmten theoretischen Bedingungen, die im Abtasttheorem beschrieben werden , lässt sich die ursprüngliche stetige Funktion perfekt aus der DTFT und somit aus den ursprünglichen diskreten Abtastwerten rekonstruieren. Die DTFT selbst ist eine stetige Frequenzfunktion, deren diskrete Abtastwerte jedoch mithilfe der diskreten Fouriertransformation (DFT) leicht berechnet werden können (siehe § Abtastung der DTFT ). Dies ist die mit Abstand gebräuchlichste Methode der modernen Fourier-Analyse.

Beide Transformationen sind umkehrbar. Die inverse DTFT rekonstruiert die ursprüngliche Abtastdatenfolge, während die inverse DFT eine periodische Summation der ursprünglichen Folge erzeugt. Die schnelle Fourier-Transformation (FFT) ist ein Algorithmus zur Berechnung eines Zyklus der DFT, und ihre Umkehrung erzeugt einen Zyklus der inversen DFT.

Bezug zur Fourier-Transformation

Sei eine stetige Funktion im Zeitbereich . Wir beginnen mit einer üblichen Definition der stetigen Fouriertransformation , wobei die Frequenz in Hertz und die Zeit in Sekunden angibt:

Wir können das Integral durch Abtastung in Sekundenintervallen in eine Summe umwandeln (siehe Fourier-Transformation § Numerische Integration einer Reihe geordneter Paare ). Genauer gesagt können wir durch eine diskrete Folge seiner Abtastwerte, , für ganzzahlige Werte von ersetzen und das Differentialelement durch die Abtastperiode ersetzen . Somit erhalten wir eine Formulierung für die diskrete Fourier-Transformation (DTFT):

Diese Fourierreihe (in der Frequenz) ist eine stetige periodische Funktion , deren Periodizität der Abtastfrequenz entspricht . Der Index unterscheidet sie von der stetigen Fouriertransformation und von der Kreisfrequenzform der DTFT. Letztere erhält man durch die Definition einer Kreisfrequenzvariablen (mit der normierten Einheit Radiant /Abtastwert ), wodurch sich eine periodische Funktion der Kreisfrequenz mit der Periodizität ergibt : [ a ]

Abb. 1. Darstellung einer Fourier-Transformation (oben links) und ihrer periodischen Summation (DTFT) in der unteren linken Ecke. Die untere rechte Ecke zeigt Abtastwerte der DTFT, die mittels diskreter Fourier-Transformation (DFT) berechnet wurden.

Die Nützlichkeit der DTFT beruht auf der Poisson-Summationsformel , die besagt, dass die durch die Fourier-Reihe dargestellte periodische Funktion eine periodische Summe der kontinuierlichen Fourier-Transformation ist : [ b ]

Poisson-Summation

Die Komponenten der periodischen Summation sind um ganzzahlige Werte (bezeichnet mit ) einer normierten Frequenz (Zyklen pro Abtastwert) zentriert. Die physikalische Frequenz (Zyklen pro Sekunde) ist das Produkt aus und der Abtastrate.   Für ausreichend große Werte von ist der Term im Bereich mit geringer oder keiner Verzerrung ( Aliasing ) durch die anderen Terme zu beobachten. Abb. 1 zeigt ein Beispiel, bei dem nicht groß genug ist, um Aliasing zu verhindern.

Wir stellen außerdem fest, dass es sich um die Fouriertransformierte von handelt . Daher lautet eine alternative Definition der DTFT : [ A ]

Die modulierte Dirac- Kammfunktion ist eine mathematische Abstraktion, die manchmal auch als Impulsabtastung bezeichnet wird . [ 3 ]

Inverse Transformation

Eine Operation, die die diskrete Datensequenz aus der DTFT-Funktion wiederherstellt, wird als inverse DTFT bezeichnet . Beispielsweise liefert die inverse kontinuierliche Fourier-Transformation beider Seiten von Gl. 3 die Sequenz in Form einer modulierten Dirac-Kammfunktion :

Da jedoch periodisch ist, sind alle notwendigen Informationen in jedem Intervall der Länge enthalten.   Sowohl in Gl. 1 als auch in Gl. 2 sind die Summen über eine Fourier-Reihe mit Koeffizienten.   Die Standardformeln für die Fourier-Koeffizienten sind auch die inversen Transformationen :

Periodische Daten

Wenn die Eingangsdatensequenz periodisch ist , kann Gl. 2 rechnerisch auf eine diskrete Fourier-Transformation (DFT) reduziert werden, da :

  • Alle verfügbaren Informationen sind in den Proben enthalten.
  • Die Funktion konvergiert überall gegen Null, außer bei ganzzahligen Vielfachen der Frequenz, den sogenannten harmonischen Frequenzen. Bei diesen Frequenzen divergiert die DTFT mit unterschiedlichen frequenzabhängigen Raten. Diese Raten ergeben sich aus der DFT eines Zyklus der Sequenz.
  • Die DTFT ist periodisch, daher ist die maximale Anzahl eindeutiger harmonischer Amplituden

Die DFT eines Zyklus der Sequenz lautet :

Und kann mithilfe der inversen Transformation ausgedrückt werden, die manchmal als diskrete Fourierreihe (DFS) bezeichnet wird : [ 1 ] : S. 542

Mit diesen Definitionen können wir den Zusammenhang zwischen der DTFT und der DFT aufzeigen :

      [ c ] [ B ]

Aufgrund der -Periodizität beider Funktionen lässt sich dies vereinfachen zu :

die die Anforderung der inversen Transformation erfüllt :

Abtastung der DTFT

Wenn die DTFT stetig ist, ist es üblich, eine beliebige Anzahl von Abtastwerten eines Zyklus der periodischen Funktion zu berechnen : [ 1 ] : S. 557–559 & 703 [ 2 ] : S. 76   

wobei es sich um eine periodische Summation handelt :

    (siehe Diskrete Fourier-Reihen )

Die Sequenz ist die inverse DFT. Durch unsere Abtastung der DTFT wird die inverse Transformation periodisch. Das Array der Werte wird als Periodogramm bezeichnet , und der Parameter heißt NFFT in der gleichnamigen Matlab -Funktion. [ 4 ]

Um einen Zyklus numerisch auszuwerten, benötigen wir eine Folge endlicher Länge . Beispielsweise kann eine lange Folge durch eine Fensterfunktion der Länge abgeschnitten werden , was zu drei Fällen führt, die besondere Erwähnung verdienen. Der Einfachheit halber betrachten wir die folgenden Werte als die durch die Fensterfunktion veränderten Werte.

Fallbeispiel: Frequenzdezimierung. für eine ganze Zahl (typischerweise 6 oder 8)

Ein Zyklus reduziert sich auf eine Summe von Segmenten der Länge .   Die DFT hat dann verschiedene Namen, wie zum Beispiel :

Die Dezimierung von Abtastdaten in einem Bereich (Zeit oder Frequenz) führt zu Überlappungen (auch Aliasing genannt ) im anderen Bereich und umgekehrt. Im Vergleich zu einer DFT der Länge λ bewirkt die Summation/Überlappung eine Dezimierung im Frequenzbereich [ 1 ] : S. 558, sodass nur die DTFT-Abtastwerte übrig bleiben, die am wenigsten von spektraler Leckage betroffen sind . Dies ist üblicherweise eine Priorität bei der Implementierung einer FFT -Filterbank ( Kanalisierer ). Mit einer herkömmlichen Fensterfunktion der Länge λ wäre der Scalloping-Verlust inakzeptabel. Daher werden Mehrblock-Fenster mithilfe von FIR-Filterdesign -Tools erstellt [ 14 ] [ 15 ] .   Ihr Frequenzprofil ist am höchsten Punkt flach und fällt in der Mitte zwischen den verbleibenden DTFT-Abtastwerten schnell ab. Je größer der Wert des Parameters λ ist, desto besser ist die potenzielle Leistung.

Fall:

Wird eine symmetrische Fensterfunktion der Länge n ( ) um einen Koeffizienten abgeschnitten, spricht man von einer periodischen oder DFT-geraden Funktion . Dies ist gängige Praxis, jedoch beeinflusst die Abschneidung die DTFT (Spektralleckage) geringfügig. Die Charakterisierung dieses Effekts ist zumindest aus akademischem Interesse. Eine DFT der abgeschnittenen Fensterfunktion der Länge n liefert Frequenzabtastwerte in Intervallen von n statt n.   Die Abtastwerte sind reellwertig ( [ 16 ] : S. 52)   , stimmen aber nicht exakt mit der DTFT der symmetrischen Fensterfunktion überein. Die periodische Summation kann zusammen mit einer DFT der Länge n verwendet werden, um die DTFT in Intervallen von n abzutasten.   Diese Abtastwerte sind ebenfalls reellwertig und stimmen exakt mit der DTFT überein (Beispiel: Datei:Sampling the Discrete-time Fourier transform.svg ). Um das volle symmetrische Fenster für die Spektralanalyse bei dem Abstand zu nutzen , würde man die und Datenproben kombinieren (durch Addition, da das symmetrische Fenster sie gleich gewichtet) und dann das abgeschnittene symmetrische Fenster und die -Länge DFT anwenden.

Abb. 2. DFT von e i2πn/8 für L = 64 und N = 256
Abb. 3. DFT von e i2πn/8 für L = 64 und N = 64

Fallbeispiel: Frequenzinterpolation.

In diesem Fall vereinfacht sich die DFT zu einer bekannteren Form :

Um die Vorteile eines schnellen Fourier-Transformationsalgorithmus zur Berechnung der DFT nutzen zu können, wird üblicherweise über alle Terme summiert, auch wenn einige davon null sind. Daher wird dieser Fall oft als Zero-Padding bezeichnet .

Spektrale Leckage, die mit abnehmendem Wert zunimmt, beeinträchtigt bestimmte wichtige Leistungskennzahlen, wie die Auflösung mehrerer Frequenzkomponenten und das Rauschen pro DTFT-Abtastwert. Diese Aspekte sind jedoch nicht immer relevant, beispielsweise wenn die Sequenz eine rauschfreie Sinuswelle (oder eine Konstante) ist, die durch eine Fensterfunktion geformt wird. In solchen Fällen ist es üblich, Nullen aufzufüllen, um die detaillierten Leckagemuster von Fensterfunktionen grafisch darzustellen und zu vergleichen. Betrachten wir zur Veranschaulichung die folgende Sequenz für ein Rechteckfenster:

Und

Abbildungen 2 und 3 zeigen die Amplituden zweier DFTs unterschiedlicher Größe, wie in den Beschriftungen angegeben. In beiden Fällen liegt die dominante Komponente bei der Signalfrequenz. In Abb. 2 ist außerdem das spektrale Leckagemuster des Rechteckfensters sichtbar. Die Illusion in Abb. 3 entsteht durch die Abtastung der DTFT an ihren Nulldurchgängen. Anstatt der DTFT einer Sequenz endlicher Länge erweckt sie den Eindruck einer unendlich langen Sinusfolge. Zu dieser Illusion tragen die Verwendung eines Rechteckfensters und die Wahl einer Frequenz (1/8 = 8/64) mit genau 8 (einer ganzen Zahl) Zyklen pro 64 Abtastwerten bei. Ein Hann-Fenster würde ein ähnliches Ergebnis liefern, allerdings würde der Peak auf 3 Abtastwerte verbreitert (siehe DFT-gerades Hann-Fenster ).

Faltung

Der Faltungssatz für Folgen lautet :

[ 17 ] : S. 297 [ d ]

Ein wichtiger Spezialfall ist die zirkuläre Faltung von Folgen s und y, definiert durch , wobei eine periodische Summation ist. Die diskrete Frequenzcharakteristik von bedeutet, dass auch das Produkt mit der stetigen Funktion diskret ist, was zu einer erheblichen Vereinfachung der inversen Transformation führt :

[ 18 ] [ 1 ] : S. 548

Für s- und y- Folgen, deren von Null verschiedene Dauer kleiner oder gleich N ist , ergibt sich folgende abschließende Vereinfachung :

Die Bedeutung dieses Ergebnisses wird unter zirkuläre Faltung und schnelle Faltungsalgorithmen erläutert .

Beziehung zur Z-Transformation

ist eine Fourierreihe , die auch mithilfe der bilateralen Z-Transformation ausgedrückt werden kann . D.h .:

Die Notation unterscheidet die Z-Transformation von der Fourier-Transformation. Daher können wir auch einen Teil der Z-Transformation mithilfe der Fourier-Transformation ausdrücken :

Beachten Sie, dass sich bei Änderung des Parameters T der Abstand der Terme von S 1/ T ( f ) konstant hält und sich ihre Breite 1/ T ändert.

Tabelle der diskreten Fourier-Transformationen

Einige gängige Transformationspaare sind in der folgenden Tabelle aufgeführt. Es gilt folgende Notation :

  • ist eine reelle Zahl, die die kontinuierliche Kreisfrequenz (in Radiant pro Abtastwert) darstellt. ( wird in Zyklen/Sekunde angegeben, und wird in Sekunden/Abtastwert angegeben.) In allen Fällen in der Tabelle ist die DTFT 2π-periodisch (in ).
  • bezeichnet eine Funktion, die auf definiert ist .
  • bezeichnet eine Funktion, die auf definiert ist und außerhalb von null. Dann:
  • ist die Dirac-Delta-Funktion
  • ist die normalisierte Sinc-Funktion
  • ist die Dreiecksfunktion
  • n ist eine ganze Zahl, die den diskreten Zeitbereich (in Abtastwerten) repräsentiert.
  • ist die diskrete Einheitssprungfunktion
  • ist das Kronecker-Delta
Zeitbereich s [ n ] Frequenzbereich S 2 π ( ω ) Anmerkungen Referenz
[ 17 ] : S. 305
ganze Zahl

    ungerades M     gerades M

ganze Zahl

Der Term muss als Verteilung im Sinne eines Cauchy-Hauptwertes um seine Pole bei interpretiert werden .
[ 17 ] : S. 305
    -π < a < π

reelle Zahl

reelle Zahl mit
reelle Zahl mit
ganze und ungerade ganze Zahl
reelle Zahlen mit
reelle Zahl ,
Es funktioniert als Differenzierungsfilter
reelle Zahlen mit
Hilbert-Transformation
reelle Zahlen komplex

Eigenschaften

Diese Tabelle zeigt einige mathematische Operationen im Zeitbereich und die entsprechenden Auswirkungen im Frequenzbereich.

Eigentum Zeitbereich s [ n ]FrequenzbereichAnmerkungen Referenz
Linearität komplexe Zahlen[ 17 ] : S. 294
Zeitumkehr / Frequenzumkehr [ 17 ] : S. 297
Zeitkonjugation [ 17 ] : S. 291
Zeitumkehr und Konjugation [ 17 ] : S. 291
Tatsächlicher Teil der Zeit [ 17 ] : S. 291
Imaginärer Zeitanteil [ 17 ] : S. 291
Realteil der Frequenz [ 17 ] : S. 291
Imaginärteil der Frequenz [ 17 ] : S. 291
Zeitverschiebung / Frequenzmodulation ganze Zahl k[ 17 ] : S. 296
Frequenzverschiebung / Zeitmodulation reelle Zahl[ 17 ] : S. 300
Dezimierung   [ E ]ganze Zahl
Zeitausdehnung ganze Zahl[ 1 ] : S. 172
Ableitung in der Frequenz [ 17 ] : S. 303
Integration in Frequenz
Zeitliche Unterschiede
Summation in der Zeit
Faltung in der Zeit / Multiplikation in der Frequenz [ 17 ] : S. 297
Multiplikation in der Zeit / Faltung in der Frequenz Periodische Faltung[ 17 ] : S. 302
Kreuzkorrelation
Satz von Parseval[ 17 ] : S. 302

Siehe auch

Anmerkungen

  1. ^ Tatsächlich wird Gl. 2 oft wie folgt begründet : [ 1 ] : S. 143, Gl. 4.6
  2. ^ Aus § Tabelle der diskreten Fourier-Transformationen ergibt sich:
  3. ^ WOLA sollte nicht mit der Überlappungs-Add-Methode der stückweisen Faltung verwechselt werden.
  4. ^ WOLA-Beispiel: Datei:WOLA channelizer example.png
  5. ^ Dieser Ausdruck wird wie folgt hergeleitet: [ 1 ] : S. 168

Seitenzitate

  1. ^ Oppenheim und Schäfer, [ 1 ] S. 147 (4.17), wo:  daher
  2. ^ Oppenheim und Schafer, [ 1 ] S. 147 (4.20), S. 694 (10.1), und Prandoni und Vetterli, [ 2 ] S. 255, (9.33), wobei:    und  
  3. ^ Oppenheim und Schafer, [ 1 ] S. 551 (8.35), und Prandoni und Vetterli, [ 2 ] S. 82 (4.43). Mit Definitionen :     und  dieser Ausdruck unterscheidet sich von den Referenzen um einen Faktor,da dieser beim Übergang vom 3. zum 4. Schritt verloren ging. Insbesondere enthält die DTFT vonin § Tabelle der diskreten Fourier-Transformationen einenFaktor, den die Referenzen ausgelassen haben.    
  4. ^ Oppenheim und Schäfer, [ 1 ] S. 60, (2.169), und Prandoni und Vetterli, [ 2 ] S. 122, (5.21)

Referenzen

  1. ^ a b c d e f g h i j k Oppenheim, Alan V. ; Schafer, Ronald W. ; Buck, John R. (1999). „4.2, 8.4“. Diskrete Signalverarbeitung (2. Aufl.). Upper Saddle River, NJ: Prentice Hall. ISBN 0-13-754920-2Abtastwerte der Fouriertransformierten einer aperiodischen Folge x[n] können als DFS-Koeffizienten einer periodischen Folge betrachtet werden, die durch Summieren periodischer Kopien von x[n] erhalten wird . 
  2. ^ a b c d Prandoni, Paolo; Vetterli, Martin (2008). Signalverarbeitung für die Kommunikation (PDF) (1. Aufl.). Boca Raton, FL: CRC Press. S. 72, 76. ISBN 978-1-4200-7046-0Abgerufen am 4. Oktober 2020. Die DFS-Koeffizienten für das periodisierte Signal bilden eine diskrete Wertemenge für seine DTFT .
  3. ^ Rao, R. (2008). Signale und Systeme . Prentice-Hall of India Pvt. Limited. ISBN 9788120338593Die
  4. ^ "Periodogramm - Schätzung der spektralen Leistungsdichte eines Periodogramms - MATLAB" . www.mathworks.com .
  5. ^ Gumas, Charles Constantine (Juli 1997). „Window-Presum FFT erreicht hohen Dynamikbereich und hohe Auflösung“ . Personal Engineering & Instrumentation News : 58–64 . Archiviert vom Original am 10.02.2001.{{cite journal}}: CS1 maint: bot: original URL status unknown (link)
  6. ^ Crochiere, RE; Rabiner, LR (1983). „7.2“. Multirate Digital Signal Processing . Englewood Cliffs, NJ: Prentice-Hall. S.  313–326 . ISBN 0136051626Die
  7. ^ Wang, Hong; Lu, Youxin; Wang, Xuegang (16. Oktober 2006). „Kanalisierter Empfänger mit WOLA-Filterbank“. 2006 CIE International Conference on Radar . Shanghai, China: IEEE. S.  1–3 . doi : 10.1109/ICR.2006.343463 . ISBN 0-7803-9582-4. S2CID  42688070 .
  8. ^ Lyons, Richard G. (Juni 2008). „DSP-Tricks: Bau eines praktischen Spektrumanalysators“ . EE Times . Abgerufen am 19.09.2024 .   Beachten Sie jedoch, dass es einen Link mit der Bezeichnung „weighted overlap-add structure“ enthält , der fälschlicherweise zur Overlap-add-Methode führt .
  9. ^ a b Lillington, John (März 2003). „Vergleich von Breitbandkanalisierungsarchitekturen“ (PDF) . Dallas: International Signal Processing Conference. S. 4 (Abb. 7). S2CID 31525301. Archiviert vom Original (PDF) am 08.03.2019 . Abgerufen am 06.09.2020 . Die „Weight Overlap and Add“-Methode (WOLA) bzw. ihre Teilmethode „Polyphase DFT“ etabliert sich zunehmend und ist besonders effizient, wenn große, hochwertige Filterbänke benötigt werden. 
  10. ^ a b Lillington, John. „A Review of Filter Bank Techniques - RF and Digital“ (PDF) . armms.org . Isle of Wight, UK: Libra Design Associates Ltd. S. 11. Abgerufen am 06.09.2020 . Glücklicherweise gibt es eine wesentlich elegantere Lösung, wie in Abbildung 20 unten dargestellt, bekannt als Polyphase- oder WOLA-FFT (Weight, Overlap and Add).
  11. ^ Hochgürtel, Stefan (2013), „2.5“, Effiziente Implementierungen hochauflösender Breitband-FFT-Spektrometer und ihre Anwendung auf eine APEX-Liniendurchmusterung im galaktischen Zentrum (PDF) , Bonn: Rheinische Friedrich-Wilhelms-Universität Bonn, S.  26–31 , Bibcode : 2013PhDT.......427H , abgerufen am 19.09.2024 . Um eine M-fache WOLA für eine N-Punkt-DFT durchzuführen, werden M·N reelle Eingangswerte a<sub> j</sub> zunächst mit einer Fensterfunktion w<sub> j </sub> gleicher Größe multipliziert.
  12. ^ Chennamangalam, Jayanth (18.10.2016). „Die Polyphasen-Filterbanktechnik“ . CASPER Group . Abgerufen am 30.10.2016 .
  13. ^ Dahl, Jason F. (06.02.2003). Time Aliasing Methods of Spectrum Estimation (Dissertation). Brigham Young University . Abgerufen am 31.10.2016 .
  14. ^ Lin, Yuan-Pei; Vaidyanathan, PP (Juni 1998). „Ein Kaiser-Fenster-Ansatz für den Entwurf von Prototypfiltern für kosinusmodulierte Filterbänke“ (PDF) . IEEE Signal Processing Letters . 5 (6): 132–134 . Bibcode : 1998ISPL....5..132L . doi : 10.1109/97.681427 . S2CID 18159105. Abgerufen am 16.03.2017 . 
  15. ^ Harris, Frederic J. (24.05.2004). „9“. Multiraten-Signalverarbeitung für Kommunikationssysteme . Upper Saddle River, NJ: Prentice Hall PTR. S.  226–253 . ISBN 0131465112Die
  16. ^ Harris, Fredric J. (Jan. 1978). „On the use of Windows for Harmonic Analysis with the Discrete Fourier Transform“ (PDF) . Proceedings of the IEEE . 66 (1): 51–83 . Bibcode : 1978IEEEP..66...51H . CiteSeerX 10.1.1.649.9880 . doi : 10.1109/PROC.1978.10837 . S2CID 426548 .  
  17. ^ a b c d e f g h i j k l m n o p q Proakis, John G.; Manolakis, Dimitri G. (1996). Digitale Signalverarbeitung: Prinzipien, Algorithmen und Anwendungen (3. Aufl.). New Jersey: Prentice-Hall International. Bibcode : 1996dspp.book.....P . ISBN 9780133942897. sAcfAQAAIAAJ.
  18. ^ Rabiner, Lawrence R.; Gold, Bernard (1975). Theorie und Anwendung der digitalen Signalverarbeitung . Englewood Cliffs, NJ: Prentice-Hall, Inc. S. 59 (2.163). ISBN 978-0139141010Die

Weiterführende Literatur

  • Porat, Boaz (1996). Einführung in die digitale Signalverarbeitung . John Wiley and Sons. S. 27–29 und 104–105. ISBN 0-471-14961-6Die
  • Siebert, William M. (1986). Schaltungen, Signale und Systeme . MIT Electrical Engineering and Computer Science Series. Cambridge, MA: MIT Press. ISBN 0262690950Die
  • Lyons, Richard G. (2010). Understanding Digital Signal Processing (3. Aufl.). Prentice Hall. ISBN 978-0137027415Die