Mikrocontroller.net - Benutzerbeiträge [de]

Entwurf: Spektralanalyse mit der FFT

2011-10-08T22:01:17Z

89.204.153.229: /* Was macht eine FFT? */

{{Baustelle}}

Oft findet man im [http://www.mikrocontroller.net/forum/dsp DSP-Forum] Fragen wie: "Ich habe 1024 Signalpunkte gespeichert und darauf eine FFT ausgeführt, was mache ich mit dem Ergebnis?"

Da muss man zuerst mal die Gegenfrage stellen: was willst du überhaupt mit dem Ergebnis machen?

In der Regel geht es dem Frager darum, ein Leistungsdichtespektrum darzustellen, also einen Plot von der Signalleistung über der Frequenz, ähnlich wie es eine Audio-Spektrumanalysator anzeigt.

== Was macht eine FFT? ==

'''FFT''' (Fast Fourier Transform) ist eigentlich nur eine Bezeichnung für einen Algorithmus zur schnellen Berechnung der '''DFT''' (Discrete Fourier Transform) mit bestimmten Randbedingungen. Da die DFT sehr oft mit dem FFT-Algorithmus berechnet wird, verwenden viele beide Begriffe gleichwertig.

Die DFT/FFT berechnet, vereinfacht ausgedrückt, welche Frequenzen im Originalsignal enthalten sind, welche Sinusschwingungen also mit welcher Phase addiert werden müssten, um auf das ursprüngliche Signal zu kommen. Konkret bekommt man für jeden Frequenzpunkt (0 bis Abtastfrequenz/2) eine komplexe Zahl, die Amplitude und Phase der entsprechenden Sinusschwingung repräsentiert.

Grafisch lässt sich die FFT so interpretieren, dass ein orthogonales Koordinatensystem bestehend aus Sin/Cos auf das Signal angewendet wird, um dessen Vektoranteile in diesem Koordinatensystem zu bestimmen, die dann zu Betrag (Länge) und Phase (Winkel) umgerechnet werden.

Eine FFT arbeitet daher grundsätzlich mit komplexen Signalen, das heißt man gibt ein Signal mit einem Real- und Imaginärteil herein, und heraus kommt ein Spektrum mit einem Real- und Imaginärteil. Arbeitet man mit rein reellen Signalen (der Normalfall), dann setzt man den Imaginärteil des Eingangssignals einfach auf 0. Ansonsten gibt es auch spezielle FFT-Implementierungen für rein reelle Signale, und auch ein paar Tricks um die Transformation von rein reellen Signalen mit komplexen FFTs beschleunigen kann. Darauf soll hier aber erst mal nicht eingegangen werden.

== Leistungsdichte ==

Oft interessiert man sich nicht für die Phasen der Sinuskomponenten des Signals, sondern nur für die Beträge, bzw. nach Quadrierung die Leistung. Wenn man die Beträge der komplexen Zahlen quadriert erhält man das sogenannte Leistungsdichtespektrum des Signals, kurz '''PSD''' (Power Spectral Density).

(Plot Zeitsignal -> PSD)

== Fensterung ==

(Plot 1. abgeschnittener Sinus, 2. mit Fensterung)

== Realität: stochastische Signale ==

In realen Anwendungen hat man meist kein Signal mit einer definierten Form wie einen Sinus, sondern Signale von denen man nicht wissen wie sie aussehen, die sich ständig ändern, und von denen man nur Erwartungswerte kennt, z. B. den Mittelwert oder die Leistung. Solche Signale nennt man Zufallssignale oder '''stochastische Signale'''. Wenn wir ein Stück aus so einem Signal, hier als Beispiel ein Sprachsignal, ausschneiden, fenstern, FFT-transformieren, den Betrag bilden und quadrieren, dann sieht das Ergebnis zum Beispiel so aus:

(Plot gefenstertes Sprachsignal, Spektrum)

== Mittelung: Die Welch-Methode ==

Diese Varianz zu reduzieren geht glücklicherweise sehr einfach: man bildet einen Mittelwert über die FFTs von mehreren Blöcken. Wie viele Blöcke man verwendet bestimmt die Reaktionszeit des Spektrums auf Änderungen. Ein weiterer Faktor ist die Überlappung zwischen den Blöcken, hierfür wird meistens 50% gewählt.

(Plot gemitteltes Spektrum)

== Mittelung: Rekursiv ==

Ein Nachteil bei der Welch-Methode ist der hohe Speicherbedarf, da die Ergebnisse aller Blöcke über die gemittelt werden soll ja irgendwo gespeichert werden müssen. Alternativ kann man daher eine rekursive Mittelung verwenden.

<math>S_{neu} = S_{alt} \cdot \lambda + \mathrm{FFT}(aktueller Block) \cdot (1 - \lambda)</math>

Beispiel:

<math>S_{neu} = S_{alt} \cdot 0.99 + \mathrm{FFT}(aktueller Block) \cdot 0.01</math>

Jetzt ist neben dem Speicher für das FFT-Ergebnis und das Spektrum kein zusätzlicher Speicher mehr nötig.

Auch hier stellt sich wieder die Frage welche Fensterung verwendet wird, und ob aufeinanderfolgende Datenblöcke überlappen sollen oder nicht. Hierfür gelten die selben Überlegungen wie für die Welch-Methode.

Außerdem gilt es natürlich den Parameter <math>\lambda</math> zu wählen. In der Praxis ist i.d.R. Ausprobieren angesagt. Wenn man schon eine Mittelungslänge für die Welch-Methode gefunden hat, dann lässt sich mit folgender Faustregel ein äquivalenter Lambda-Faktor für die rekursive Mittelung berechnen<ref>Porat, B.: Second-order equivalence of rectangular and exponential windows in least-squares estimation of Gaussian autoregressive processes</ref>:

<math>R = (1 + \lambda)/(1 - \lambda)</math>

(Plot gefiltertes Spektrum)

=== Effiziente Implementierung ===

Da Multiplikationen auf vielen Mikrocontrollern aufwändig sind, ist es sinnvoll sich auf einen Faktor <math>\lambda = 1 - 2^{-x}</math> zu beschränken und die Berechnung mit Schiebeoperationen durchzuführen:

<math>S_{neu} = S_{alt} \cdot (1 - 2^{-x}) + \mathrm{FFT}(aktueller Block) \cdot 2^{-x}</math>

Mit ein bisschen Umstellen und Ersetzen der Multiplikation mit <math>2^{-x}</math> durch x-faches Rechtsschieben kommt man auf:

<math>Differenz = \mathrm{FFT}(aktueller Block) - S_{alt}</math>

<math>Aenderung = Differenz >> x</math>

<math>S_{neu} = S_{alt} + Aenderung</math>

Ein Nachteil der rekursiven Mittelung ist, dass Lambda-Werte sehr nahe bei 1 (bzw. ein großer Wert x) dazu führen können, dass "Aenderung" durch die begrenzte Genauigkeit immer 0 ist, und sich <math>S_{neu}</math> deshalb nicht ändert.

<references />
[[Category:DSP]]

Abtasttheorem

2011-10-08T21:55:59Z

89.204.153.229:

Das sogenannte 'Abtasttheorem' von Nyquist und Shannon beschreibt die Tatsache, dass sich ein beliebig geformtes analoges Signal immer dann exakt wiederherstellen lässt, wenn die Abtastrate beim Digitalisieren und damit die Dichte der Stützpunkte größer ist, als die höchste im Signal enthaltene Oberwelle.

Dies erklart sich dadurch, dass beim späteren Rekonstruieren ein ideales Tiefpassfilter theoretisch wieder einen oberwellenfreien Sinus zwischen den Stützstellen generieren kann, mit dem es möglich ist, die maximale Krümmung im Ursignal, die ein Mass für die maximale Änderungsgeschwindigkeit darstellt, abzubilden.

Höhere Frequenzen, die im Ursignal enthalten waren, können nicht wieder hergestellt werden. Im Gegenteil: Sie verursachen beim Digitalisieren sogenannte aliasing Fehler, die dazu führen, dass falsche Amplitudenwerte erfasst werden, die bei der Rekonstruktion zu ebenfalls falschen Amplitudenverläufen des erzeugten Signals führen.

Ebenso ist festzustellen, dass aufgrund der Unzulänglichkeiten realer Filter die theoretisch exakt mögliche Rekonstruktion des Signals in der Praxis nicht gelingt. Bei der künstlichen Erzeugung von Sinuswellen z.B., wie bei der [[DDS]], liegen mathematisch exakte Amplitudenwerte vor, dennoch kommt man mit den generierten Sinuswellen praktisch nur bis an 70%-80% der halben Abtastfrequenz heran.

Daher muss je nach Qualitätsanspruch sowohl bei der Digitalsierung als auch der Rekonstruktion des Signals mit deutlich höheren Abtastfrequenzen gearbeitet werden.

[[Kategorie:Signalgeneratoren]]
[[Kategorie:Grundlagen]]
[[Kategorie:Datenübertragung]]
[[Kategorie:Signalverarbeitung]]

Abtasttheorem

2011-10-08T21:53:27Z

89.204.153.229: kat

Abtasttheorem

2011-10-08T21:50:27Z

89.204.153.229: neu