| Fourier-Transformationen |
|---|
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 ]
| Gl. 1 |

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 ]
| Gl. 2 |
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 ]
| Gl. 3 |
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 :
| Gl. 4 |
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 :
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 :
- Fenster-Vorsummen-FFT [ 5 ]
- Gewichtung, Überlappung, Addition (WOLA) [ 6 ] [ 7 ] [ 8 ] [ 9 ] [ 10 ] [ 11 ] [ C ] [ D ]
- Polyphasen-DFT [ 9 ] [ 10 ]
- Polyphasen-Filterbank [ 12 ]
- Mehrfachblock-Fensterung und Zeitaliasing . [ 13 ]
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.


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 :
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 :
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.
- ist die diskrete Faltung zweier Sequenzen
- ist die komplex Konjugierte von
| Eigentum | Zeitbereich s [ n ] | Frequenzbereich | Anmerkungen | 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
- Spektralanalyse nach der Methode der kleinsten Quadrate
- Multidimensionale Transformation
- Zak verwandelt sich
Anmerkungen
- ^ Tatsächlich wird Gl. 2 oft wie folgt begründet : [ 1 ] : S. 143, Gl. 4.6
- ^ Aus § Tabelle der diskreten Fourier-Transformationen ergibt sich:
- ^ WOLA sollte nicht mit der Überlappungs-Add-Methode der stückweisen Faltung verwechselt werden.
- ^ WOLA-Beispiel: Datei:WOLA channelizer example.png
- ^ Dieser Ausdruck wird wie folgt hergeleitet: [ 1 ] : S. 168
Seitenzitate
- ^ Oppenheim und Schäfer, [ 1 ] S. 147 (4.17), wo: daher
- ^ Oppenheim und Schafer, [ 1 ] S. 147 (4.20), S. 694 (10.1), und Prandoni und Vetterli, [ 2 ] S. 255, (9.33), wobei: und
- ^ 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.
- ^ Oppenheim und Schäfer, [ 1 ] S. 60, (2.169), und Prandoni und Vetterli, [ 2 ] S. 122, (5.21)
Referenzen
- ^ 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-2
Abtastwerte 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
. - ^ 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
. - ^ Rao, R. (2008). Signale und Systeme . Prentice-Hall of India Pvt. Limited. ISBN 9788120338593Die
- ^ "Periodogramm - Schätzung der spektralen Leistungsdichte eines Periodogramms - MATLAB" . www.mathworks.com .
- ^ 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) - ^ Crochiere, RE; Rabiner, LR (1983). „7.2“. Multirate Digital Signal Processing . Englewood Cliffs, NJ: Prentice-Hall. S. 313–326 . ISBN 0136051626Die
- ^ 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 .
- ^ 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 .
- ^ 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.
- ^ 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).
- ^ 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.
- ^ Chennamangalam, Jayanth (18.10.2016). „Die Polyphasen-Filterbanktechnik“ . CASPER Group . Abgerufen am 30.10.2016 .
- ^ Dahl, Jason F. (06.02.2003). Time Aliasing Methods of Spectrum Estimation (Dissertation). Brigham Young University . Abgerufen am 31.10.2016 .
- ^ 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 .
- ^ Harris, Frederic J. (24.05.2004). „9“. Multiraten-Signalverarbeitung für Kommunikationssysteme . Upper Saddle River, NJ: Prentice Hall PTR. S. 226–253 . ISBN 0131465112Die
- ^ 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 .
- ^ 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.
- ^ 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