Hauptmenü öffnen

Goertzel-Algorithmus

Methode der Spektralanalyse

Der Goertzel-Algorithmus ist ein Verfahren aus der digitalen Signalverarbeitung und stellt eine besondere Form der diskreten Fourier-Transformation (DFT) dar. Im Gegensatz zu den verschiedenen schnellen Berechnungsmethoden bei der diskreten Fouriertransformation (FFT), welche immer alle diskreten Spektralkomponenten in einem Block berechnen, ist es mit dem Goertzel-Algorithmus möglich, nur einzelne diskrete Spektralanteile zu berechnen. Entwickelt wurde der Algorithmus 1958 von Gerald Goertzel (1919–2002).

FunktionBearbeiten

Der Algorithmus basiert auf einer Struktur bestehend aus einem digitalen Filter, das um eine Zustandssteuerung erweitert ist. Die Zustände unterteilen die Berechnung in den Rückwärtszweig, in dem die im Zeitbereich abgetasteten Eingangswerte geladen werden, und in einen Vorwärtszweig, der das Ausgangssignal liefert. Die Rückwärtsschleife wird bei jedem digitalen Abtastwert (engl. sample) durchlaufen und ist als ein rekursives digitales Filter mit zwei Zustandsspeichern und einem Akkumulator aufgebaut. Der Vorwärtszweig wird erst nach N Abtastwerten einmalig durchlaufen und liefert aus den Zustandsspeichern den berechneten komplexen Ausgangswert, die spektrale Komponente nach Betrag und Phase.

Durch die Wahl der dabei eingesetzten Filterkoeffizienten lässt sich die Frequenzselektivität einstellen. Durch die Wahl der Anzahl der Abtastwerte N lässt sich der Gütefaktor beeinflussen. N kann beliebige natürliche Werte annehmen.

Pro Spektralkomponente ist allerdings eine eigenständige Goertzel-Struktur notwendig. Daher ist dieser Algorithmus vor allem dann vorteilhaft und mit geringerem Rechenaufwand anzuwenden, wenn nicht das komplette Spektrum berechnet werden soll, sondern nur einzelne Spektralkomponenten daraus.

Ausführliche mathematische Herleitungen des Algorithmus finden sich in den unten angegebenen Literaturquellen.

AlgorithmusBearbeiten

Von einem diskreten Signal s wird der Realteil C und der Imaginärteil S einer Spektralkomponente   über einen rekursiven Algorithmus bestimmt. Die Startbedingungen sind  . Über

  für  

ergibt sich

 
 

Auf die Winkelfunktion   und   muss dabei nur je einmal zugegriffen werden.

AufwandsabschätzungBearbeiten

Pro Berechnung einer Spektralkomponente sind bei dem Goertzel-Algorithmus 2N+2 Additionen/Subtraktion und N+2 Multiplikationen notwendig.[1] Vergleicht man diesen Aufwand mit dem Berechnungsaufwand bei der schnellen Fourier-Transformation ist der Goertzel-Algorithmus immer dann effizienter, wenn weniger als 5/6 log2(N) Spektralkomponenten berechnet werden sollen. Denn pro Spektralkomponente (bin) ist eine weitere Goertzelstruktur notwendig, während bei der schnellen Fourier-Transformation der Berechnungsaufwand nur mit Nlog2N ansteigt.

Der Algorithmus kann effizient in digitalen Signalprozessoren implementiert werden.

AnwendungenBearbeiten

Die Anwendungen liegen in der Erkennung einzelner Frequenzen (Tonerkennung) in einem Signal wie beispielsweise bei der Erkennung der Signalisierungsfrequenzen bei dem im Telefonbereich eingesetzten Mehrfrequenzwahlverfahren. In diesem Fall muss nur der Betrag der Spektralkomponente ausgewertet werden, was weitere Vereinfachungen in der Berechnung gestattet.

LiteraturBearbeiten

WeblinksBearbeiten

EinzelnachweiseBearbeiten

  1. Karl Dirk Kammeyer, Kristian Kroschel: Digitale Signalverarbeitung: Filterung und Spektralanalyse mit MATLAB-Übungen. Springer-Verlag, 2009, ISBN 978-3-8348-0610-9, S. 274 (google.de [abgerufen am 24. Juni 2017]).