DE102013217181A1 - Verfahren und System zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit - Google Patents
Verfahren und System zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit Download PDFInfo
- Publication number
- DE102013217181A1 DE102013217181A1 DE201310217181 DE102013217181A DE102013217181A1 DE 102013217181 A1 DE102013217181 A1 DE 102013217181A1 DE 201310217181 DE201310217181 DE 201310217181 DE 102013217181 A DE102013217181 A DE 102013217181A DE 102013217181 A1 DE102013217181 A1 DE 102013217181A1
- Authority
- DE
- Germany
- Prior art keywords
- sampled data
- signal
- filter
- sets
- fft
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/175—Spectrum analysis; Fourier analysis by delay means, e.g. tapped delay lines
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/165—Spectrum analysis; Fourier analysis using filters
- G01R23/167—Spectrum analysis; Fourier analysis using filters with digital filters
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
- Radar Systems Or Details Thereof (AREA)
- Image Processing (AREA)
Abstract
Description
- HINTERGRUND
- Herkömmliche Verfahren zum Durchführen einer Spektralanalyse von Signalen setzen voraus, dass die Eingangssignale beim Abtastvorgang während eines Beobachtungsfensters einen konstanten Oberschwingungsanteil besitzen. In vielen Anwendungen sind Signale zeitabhängig (nicht stationär oder transient) und können für Beobachtungszwecke nur dann als stationär angesehen werden, wenn die Dauer oder Länge des Beobachtungsfensters angemessen kurz ist. Die Länge des Beobachtungsfensters sollte beispielsweise kleiner als die Unterbrechungsdauer oder die Abklingdauer sein, so dass die Signalstatistik während des gesamten Beobachtungsfensters stationär bleibt. Demgegenüber wird die Anzahl der Samples erhöht, was die Verwendung eines längeren Beobachtungsfensters erfordert, um eine gute Spektralauflösung des Signals zu erzielen.
- Im Allgemeinen ist eine Fast Fourier Transformation (FFT) nicht geeignet, wenn ein schmales Frequenzband eines nicht stationären Signals analysiert werden soll, da eine gute Spektralauflösung mit einem kurzen Beobachtungsfenster nicht möglich ist. Um sowohl die Spektralauflösung als auch die Auflösung bei Verwendung der FFT zu verbessern, wird ein längeres Beobachtungsfenster, d. h. eine größere Anzahl an Punkten, benötigt. Ein nicht stationäres Signal kann jedoch nicht über die gesamte Dauer eines langen Beobachtungsfensters als stationär angesehen werden, da nicht stationäre Signale nicht konstant sind. Daher ist es schwierig, eine gute Spektralauflösung bei gleichzeitiger Reduzierung der Länge des Beobachtungsfensters in Anpassung an ein nicht stationäres Signal zu erzielen. Anders ausgedrückt ist es eventuell nicht möglich, die Länge des Beobachtungsfensters zu vergrößern, da das nicht stationäre Signal nur für ein kurzes Zeitintervall als stationär angenommen werden kann.
- Eine weitere Beschränkung der FFT ist dadurch gegeben, dass sie die Analyse des gesamten Spektrums erfordert, selbst wenn der spektrale Gehalt des Signals nur in einem begrenzten Frequenzband ausgewertet wird. Dementsprechend muss, da eine bestimmte Auflösung für das betrachtete Frequenzband benötigt wird, das gesamte Spektrum mit derselben Auflösung analysiert werden. Daraus kann sich ein langes Beobachtungsfenster ergeben, welches die Annahme, dass das Signal stationär ist, wie oben beschrieben, verhindert.
- Verglichen damit erfordert ein Algorithmus mit einer Chirp-Z-Transformation (CZT) nicht die Analyse des gesamten Spektrums mit derselben Auflösung, so wie bei der FFT. Auf diese Weise liefert die CZT eine beschränkte Spektralanalyse mit besserer Spektralauflösung, da die CZT nicht wie die FFT (und andere klassische Verfahren, wie das nach Welch) nur von der Länge des Beobachtungsfensters abhängt. Die CZT hängt vielmehr von dem Verhältnis zwischen analysiertem Spektralband und Abtastfrequenz ab. In diesem Fall ist es nach Wahl eines Frequenzbandes möglich, auch mit einer erheblich niedrigeren Anzahl an Samples und einem entsprechend kürzeren Beobachtungsfenster eine gute Spektralauflösung zu erreichen. Auf diese Weise erlaubt die CZT eine Optimierung aller Leistungsparameter, selbst wenn ein kürzeres Beobachtungsfenster benötigt wird.
- Die CZT ermöglicht, nur durch Verwendung der abgetasteten Daten eine gezoomte Signaldarstellung zu erzielen, ohne irgendwelche weiteren Operationen, um Leck- oder Amplitudenfehler zu korrigieren. Ebenso ist es möglich, die CZT zu verwenden, um das vollständige Signalspektrum, selbst für nicht stationäre Signale, wiederherzustellen, indem verschiedene CZTs an begrenzten Spektrumsabschnitten unter Verwendung derselben gepufferten, abgetasteten Daten ausgewertet werden. CZT allein liefert keine verbesserte Spektralauflösung für die gezoomte oder die gesamte Bandbreite des Signals. Daher besteht ein Bedarf, eine genauere Signalanalyse durch Verbessern der Spektralauflösung zu erreichen, was nicht durch Verwendung von CZT allein geschehen kann.
- KURZBESCHREIBUNG
- In einer charakteristischen Weiterbildung wird ein Verfahren zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit zur Verfügung gestellt. Das Verfahren umfasst: Abtasten des nicht stationären Signals unter Verwendung eines Beobachtungsfensters mit einer Länge, die kurz genug ist, um ein stationäres Signal zu approximieren, zum Bereitstellen eines ersten Satzes abgetasteter Daten; Pufferung des ersten Satzes abgetasteter Daten, um mehrere gepufferte Sätze abgetasteter Daten zu erhalten; Filtern des ersten Satzes abgetasteter Daten und der gepufferten Sätze abgetasteter Daten unter Verwendung entsprechender Filterantworten, um mehrere gefilterte Sätze abgetasteter Daten zu erhalten; und Durchführen einer Chirp-Z-Transformation (CZT) der gefilterten Sätze abgetasteter Daten, um einen Satz von Koeffizienten von Diskreten Fouriertransformationen (DFT) bereitzustellen. Ein vollständiges Signalspektrum des nicht stationären Signals wird unter Verwendung des Satzes von DFT Koeffizienten wiederhergestellt.
- In einer weiteren charakteristischen Weiterbildung wird ein Verfahren zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit zur Verfügung gestellt. Das Verfahren umfasst: Abtasten des nicht stationären Signals unter Verwendung eines kurzen Beobachtungsfensters zum Approximieren eines stationären Signals, um abgetastete Signale bereitzustellen; Filtern der abgetasteten Signale unter Verwendung eines ersten Mehrphasenfilters, um erste gefilterte, abgetastete Signale bereitzustellen; Verzögern der abgetasteten Signale, um verzögerte, abgetastete Signale bereitzustellen; Filtern der verzögerten, abgetasteten Signale unter Verwendung eines zweiten Mehrphasenfilters, um zweite gefilterte, abgetastete Signale bereitzustellen; Durchführen von CZTs der ersten und zweiten gefilterten, abgetasteten Signale, um DFT Koeffizienten bereitzustellen; und Wiederherstellen eines vollständigen Signalspektrums des nicht stationären Signals unter Verwendung der DFT Koeffizienten.
- In einer weiteren charakteristischen Weiterbildung wird ein System zum Durchführen einer Spektralanalyse eines Eingangssignals von einem Messobjekt (DUT = Device Under Test) in Echtzeit zur Verfügung gestellt. Das System umfasst einen Abtaster, der dazu ausgebildet ist, das Eingangssignal abzutasten, um erste abgetastete Daten bereitzustellen, einen ersten Mehrphasenfilter, der dazu ausgebildet ist, die ersten abgetasteten Daten zu filtern, um erste gefilterte, abgetastete Daten bereitzustellen, einen ersten Puffer, der dazu ausgebildet ist, die ersten abgetasteten Daten um eine erste Verzögerung zu verzögern, um zweite abgetastete Daten bereitzustellen, einen zweiten Mehrphasenfilter, der dazu ausgebildet ist, die zweiten abgetasteten Daten zu filtern, um zweite gefilterte, abgetastete Daten bereitzustellen, einen zweiten Puffer, der dazu ausgebildet ist, die zweiten abgetasteten Daten um eine zweite Verzögerung zu verzögern, um dritte abgetastete Daten bereitzustellen, und einen dritten Mehrphasenfilter, der dazu ausgebildet ist, die dritten abgetasteten Daten zu filtern, um dritte gefilterte, abgetastete Daten bereitzustellen. Das System umfasst weiterhin einen Addierer, der dazu ausgebildet ist, die ersten, zweiten und dritten gefilterten, abgetasteten Daten zu addieren, um addierte, gefilterte Daten bereitzustellen, und eine CZT-Schaltung, die dazu ausgebildet ist, eine CZT an den addierten, gefilterten Daten durchzuführen, um Frequenzdomänensamples bereitzustellen, die dem Eingangssignal entsprechen. Die Frequenzdomänensamples ermöglichen ein Wiederherstellen eines vollständigen Signalspektrums des Eingangssignals.
- KURZBESCHREIBUNG DER ZEICHNUNGEN
- Die veranschaulichenden Weiterbildungen können am besten verstanden werden, wenn die folgende detaillierte Beschreibung zusammen mit den begleitenden Zeichnungen gelesen wird. Es wird hervorgehoben, dass die einzelnen Elemente nicht notwendigerweise maßstabsgetreu gezeichnet sind. Vielmehr können die Abmessungen zur klareren Darstellung beliebig vergrößert oder reduziert sein. Überall, wo es möglich und praktisch ist, beziehen sich gleiche Bezugszeichen auf gleiche Elemente.
-
1 ist ein Flussdiagram eines Verfahrens zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit unter Verwendung eines auf einer kontinuierlichen Filterbank beruhenden CZT-Algorithmus gemäß einer charakteristischen Weiterbildung. -
2 ist ein funktionales Blockdiagramm, welches einen Chirp-Z-Transformationsalgorithmus (CZT) darstellt. -
3 ist ein funktionales Blockdiagramm, welches einen CZT-Algorithmus gemäß einer Implementierung in der Frequenzdomäne darstellt. -
4 ist ein funktionales Blockdiagramm, welches einen auf einer fensterbewerteten FFT beruhenden CZT-Algorithmus darstellt, der nicht stationäre Signale abbildet. -
5 ist ein funktionales Blockdiagramm, welches den auf der fensterbewerteten FFT beruhenden CZT-Algorithmus, der in4 gezeigt ist, mit kombinierten Vormultiplizierern darstellt. -
6 ist ein funktionales Blockdiagramm, welches eine Filterbank eines Spektralanalysators darstellt. -
7 ist ein funktionales Blockdiagramm, welches eine kontinuierliche Filterbank eines Spektralanalysators darstellt. -
8 ist ein funktionales Blockdiagramm, welches einen auf einer kontinuierlichen Filterbank beruhenden CZT-Algorithmus gemäß einer charakteristischen Weiterbildung darstellt. -
9 ist ein funktionales Blockdiagramm, welches einen auf einer kontinuierlichen Filterbank beruhenden CZT-Algorithmus gemäß einer charakteristischen Weiterbildung darstellt. -
10 ist ein funktionales Blockdiagramm, welches einen auf einer kontinuierlichen Filterbank beruhenden CZT-Algorithmus gemäß einer charakteristischen Weiterbildung darstellt. -
11 ist ein Diagramm, welches eine spiralförmige CZT-Kontur in der z-Ebene zum Berechnen einer CZT zeigt. - DETAILLIERTE BESCHREIBUNG
- In der folgenden detaillierten Beschreibung werden zum Zwecke der Erläuterung und nicht der Beschränkung veranschaulichende Weiterbildungen, die spezifische Details offenbaren, dargelegt, um ein tiefgreifendes Verständnis der Weiterbildungen gemäß der vorliegenden Lehre zu ermöglichen. Es wird jedoch jedem, der in den Vorzug der vorliegende Offenbarung gekommen ist, klar sein, dass weitere Weiterbildungen gemäß der vorliegenden Lehre, welche von den spezifischen, hierin offenbarten Details abweichen, unter den Umfang der beigefügten Ansprüche fallen. Darüber hinaus können Beschreibungen wohlbekannter Vorrichtungen und Verfahren entfallen, um die Beschreibung der beispielhaften Weiterbildungen nicht zu verschleiern. Derartige Verfahren und Vorrichtungen liegen innerhalb des Rahmens der vorliegenden Lehre.
- In der Regel kann davon ausgegangen werden, dass die Zeichnungen und die einzelnen, darin dargestellten Elemente nicht maßstabsgetreu gezeichnet sind. Des Weiteren werden relative Begriffe, wie etwa ”über”, ”unter”, ”oben”, ”unten”, ”oberer”, ”unterer”, ”links”, ”rechts”, ”vertikal” und ”horizontal” verwendet um den Bezug der einzelnen Elemente zueinander, so wie sie in den beigefügten Zeichnungen dargestellt sind, zu beschreiben. Es versteht sich, dass diese relativen Begriffe zusätzlich zu der in den Zeichnungen dargestellten Ausrichtung verschiedenste Ausrichtungen der Vorrichtung und/oder der Elemente umfassen sollen. Falls die Vorrichtung beispielsweise in Bezug auf die Ansicht in den Zeichnungen umgekehrt werden sollte, würde ein Element, das als ”über” einem anderen Element beschrieben ist, dann beispielsweise ”unter” diesem Element sein. Entsprechend würde in dem Fall, dass die Vorrichtung um 90 Grad in Bezug auf die Ansicht in den Zeichnungen gedreht werden sollte, ein Element, dass als ”vertikal” beschrieben ist, dann ”horizontal” sein.
- Eine Reihe von Weiterbildungen stellen Abtastsysteme und -verfahren zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit zur Verfügung. In der Regel wird ein nicht stationäres Signal unter Verwendung eines kurzen Beobachtungsfensters abgetastet, um ein stationäres Signal zu approximieren. Das abgetastete Signal wird mittels einer Reihe von Filterbänken gefiltert, welche eine verbesserte Spektralauflösung für die gesamte Bandbreite des abgetasteten Signals bereitstellen. An dem gefilterten Signal wird ein CZT-Algorithmus mittels einer CZT-Schaltung durchgeführt, die eine Zoom-Darstellung beliebiger Bandbreiten nach Wahl durch den Anwender (beispielsweise mit minimaler und maximaler Frequenz [fmin, fmax]) erlaubt. Die Kombination verschafft die Möglichkeit, Merkmale bei verbesserter Auflösung über die gesamte Bandbreite des abgetasteten Signals zu zoomen und zu identifizieren.
-
1 ist ein Flussdiagramm eines Verfahrens zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit unter Verwendung eines auf einer kontinuierlichen Filterbank beruhenden CZT Algorithmus gemäß einer charakteristischen Weiterbildung. - Mit Berg auf
1 wird ein nicht stationäres Signal in Block S110 empfangen und abgetastet, um einen ersten Satz abgetasteter Daten bereitzustellen, wobei ein Beobachtungsfenster verwendet wird, dessen Länge kurz genug ist, um ein stationäres Signal zu approximieren. Das Signal kann ein diskretes Signal sein, welches von einem Messobjekt (DUT = Device Under Test), beispielsweise in Reaktion auf einen Auslöseimpuls, ausgegeben wird. - In Block S120 wird der erste Satz abgetasteter Daten gepuffert, um mehrere gepufferte Sätze abgetasteter Daten zu erhalten. Jeder von dem ersten Satz abgetasteter Daten und den mehreren gepufferten Sätzen abgetasteter Daten entspricht einem begrenzten Spektrumsabschnitt eines vollständigen Signalspektrums des nicht stationären Signals. Gemeinsam emulieren die gepufferten Sätze abgetasteter Daten in effektiver Weise Daten, die unter Verwendung eines längeren Beobachtungsfensters gesammelt werden, ausgehend von dem nicht stationären Signal, das mit dem kurzen Beobachtungsfenster abgetastet wird. Aus diesem Grund kann eine Spektralauflösung des nicht stationären Signals aufgrund eines im Verhältnis größeren Eingangsdatensatz, der den ersten und die gepufferten Sätze abgetasteter Daten umfasst, verbessert werden, wodurch beispielsweise ein Frequenzauflösungsfehler abgeschwächt wird.
- Beispielsweise kann der erste Satz abgetasteter Daten an einen ersten Puffer (oder eine andere Verzögerungsschaltung) geliefert werden, um um eine vorbestimmte erste Verzögerungszeit verzögert zu werden. Die Ausgabe des ersten Puffers (d. h. erste verzögerte Daten) kann dann an einen zweiten Puffer geliefert werden, um um eine vorbestimmte zweite Verzögerungszeit, welche dieselbe oder verschieden von der ersten Verzögerungszeit sein kann, verzögert zu werden. Die Ausgabe des zweiten Puffers (d. h. zweite verzögerte Daten) kann dann an einen dritten Puffer geliefert werden, um um eine vorbestimmte dritte Verzögerungszeit, die gleich oder verschieden von der ersten und zweiten Verzögerungszeit sein kann, verzögert zu werden. In verschiedenen Weiterbildungen kann die Gesamtzahl der Puffer eine beliebige positive, von dem Anwender wählbare ganze Zahl sein. Die Gesamtzahl der Puffer ist direkt proportional zu der gewünschten (emulierten) Länge des Beobachtungsfensters.
- In Block S130 werden der erste Satz abgetasteter Daten und die mehreren gepufferten Sätze abgetasteter Daten mittels Mehrphasenfilter unter Verwendung entsprechender Filterantworten gefiltert, um mehrere gefilterte Sätze abgetasteter Daten zu erhalten. Jede der Filterantworten kann einen Tiefpassfilter, beispielsweise unter Verwendung eines vom Anwender wählbaren Filterkoeffizienten, bereitstellen. Allgemein ist die Filterantwort (H(ω)) die FFT des Filterkoeffizienten (h(n)) für einen bestimmten Filter, wobei ω der Frequenzindex und n der Koeffizientenindex ist. In vielen Weiterbildungen können die Filterantworten von einem einzigen Filter abgeleitet werden. Ähnlich wie bei den Puffer ist die Gesamtzahl der Mehrphasenfilter in einer bestimmten Implementierung direkt proportional zu der gewünschten (emulierten) Länge des Beobachtungsfensters. Beispielsweise würden zwei Filter effektiv die Länge des Beobachtungsfensters, das zum Sammeln des ersten Satzes abgetasteter Daten verwendet wird, verdoppeln. Die Zahl der Mehrphasenfilter (und der entsprechenden Puffer) kann durch die gesamte Bandbreite des Spektrums des empfangenen nicht stationären Signals bestimmt werden. Die gefilterten Sätze abgetasteter Daten, die von den jeweiligen Mehrphasenfilter ausgegeben werden, werden in Block S140 addiert.
- In Block S150 wird eine Chirp-Z-Transformation (CZT) an den addierten, gefilterten Sätzen abgetasteter Daten unter Verwendung eines CZT-Algorithmus durchgeführt, um einen Satz von Koeffizienten von Diskreten Fouriertransformationen (DFT) bereitzustellen. In verschiedenen Weiterbildungen kann ein Teil des CZT-Algorithmus, wie beispielsweise Vormultiplizierer, mittels der Mehrphasenfilter durchgeführt werden. Das Durchführen des CZT-Algorithmus kann ein Durchführen einer Fast Fourier Transformation (FFT) an den addierten, gefilterten Sätzen abgetasteter Daten, ein Filtern der FFT Sätze abgetasteter Daten unter Verwendung einer FFT einer unten beschriebenen, komplexwertigen Chirp-Funktion, ein Durchführen einer inversen FFT (IFFT) an den gefilterten FFT Sätzen abgetasteter Daten und ein Filtern der IFFT Sätze abgetasteter Daten unter Verwendung einer weiteren komplexwertigen Chirp-Funktion umfassen, um den Satz von DFT Koeffizienten bereitzustellen.
- Das vollständige Signalspektrum des empfangenen nicht stationären Signals kann dann in Block S160 unter Verwendung des Satzes von DFT Koeffizienten wiederhergestellt werden. Das Wiederherstellen des vollständigen Signalspektrums des nicht stationären Signals kann ein Berechnen des Satzes von DFT Koeffizienten beispielsweise in einer z-Ebene umfassen, wobei der Satz von DFT Koeffizienten auf kreisförmigen oder spiralförmigen Konturen liegt, die an einem frei wählbaren Punkt in der z-Ebene wie in
4 angedeutet und unten beschrieben beginnen. - Die Blöcke in
1 können mittels eines Verarbeitungsgeräts, wie etwa eines Prozessors, eines digitalen Signalprozessors (DSP), eines oder mehrerer anwendungsspezifischer integrierter Schaltkreise (ASICs), eines oder mehrerer feldprogrammierbarer Gate-Arrays (FPGAs) oder Kombinationen davon unter Verwendung von Software, Firmware, festverdrahteter logischer Schaltungen oder Kombinationen davon implementiert werden. Bei Verwendung eines Computerprozessors und/oder eines DSP beispielsweise kann ein Speicher zum Speichern von ausführbarer Software/Firmware und/oder ausführbarem Code vorgesehen sein, der es ihm erlaubt, die verschiedenen Funktionen auszuführen. Der Speicher kann ein nicht flüchtiges, computer-lesbares Medium sein und kann eine beliebige Zahl, Typ und Kombination von beispielsweise Direktzugriffsspeicher (RAM) und nicht flüchtigem Speicher (z. B. Festspeicher (ROM)) aufweisen. Eine Filterbank, die vollständig oder teilweise mittels des Verarbeitungsgeräts implementiert werden kann, kann die Funktionalität der Mehrphasenfilter und der entsprechenden Filterantworten bereitstellen. Die Filterbank empfängt das nicht stationäre Signal von einem Analog-Digital-Wandler (ADC) und ermöglicht eine Darstellung des gesamten oder eines ausgewählten Abschnitts des Signalspektrums auf einem Anzeigegerät. - Hier sind die
2 –10 funktionale Blockdiagramme, welche verschiedene CZT-Algorithmen darstellen. Insbesondere sind die8 –10 funktionale Blockdiagramme, die auf einer kontinuierlichen Filterbank beruhende CZT-Algorithmen gemäß charakteristischen Weiterbildungen darstellen. Die2 –7 werden zur Erklärung der auf einer kontinuierlichen Filterbank beruhenden CZT-Algorithmen in den8 –10 herangezogen. - Im Allgemeinen ist ein CZT-Algorithmus eine effektive Methode zum numerischen Berechnen einer Z-Transformation einer Folge von N Samples (wobei N eine positive ganze Zahl ist) von einem empfangenen Eingangssignal. Der CZT-Algorithmus ermöglicht ein Berechnen der Z-Transformation an M Punkten (wobei M eine positive ganze Zahl ist) in der z-Ebene, die auf kreisförmigen oder spiralförmigen Konturen liegen, die beispielsweise wie in
11 gezeigt an einem frei wählbaren Punkt in der z-Ebene beginnen. Der CZT-Algorithmus beruht auf der Tatsache, dass die Werte der Z-Transformation auf den kreisförmigen oder spiralförmigen Konturen als eine diskrete Faltung dargestellt werden können. Die DFT für ein diskretes Signal x(n) wird durch Xk angegeben und durch Gleichung (1) beschrieben, wobei x(n), 0 ≤ n ≤ N – 1 eine vorgegebene N Punkte Folge ist: -
-
-
- Sodann können die Frequenzen θ0 und φ0 gemäß Gleichung (7) bzw. (8) bestimmt werden.
-
- Für eine allgemeine Spektralanalyse ist jede der komplexen Zahlen A0 und W0 gleich eins, wobei die DFT auf dem Einheitskreis berechnet wird. Es ist offensichtlich, dass die Frequenz θ0 die Startfrequenz [fmin] bestimmt und dass die Frequenz φ0 die Endfrequenz [fmax] und die Auflösung bestimmt.
-
-
-
-
2 ist ein funktionales Blockdiagramm, welches den CZT-Algorithmus aus Gleichung (13) darstellt. Jeder der Blöcke in2 (ebenso wie in den funktionalen Blockdiagrammen, die hier in Hinblick auf3 –10 besprochen werden) kann beispielsweise ein Verarbeitungsmodul, wie etwa ein Softwaremodul, darstellen, dass dazu ausgebildet ist, die dargestellte Funktion durchzuführen. Die Verarbeitungsmodule können von einem Verarbeitungsgerät ausgeführt werden, das durch einen Prozessor, einen DSP, einen oder mehrere ASICs, ein oder mehrere FPGAs oder Kombinationen davon unter Verwendung von Software, Firmware, festverdrahteten logischen Schaltungen oder Kombinationen davon wie oben beschrieben implementiert werden kann. - Bezugnehmend auf
2 bestimmt der Multiplizierer215 das Produkt aus abgetastetem, diskretem Signal xn (wobei n = 0, 1, ..., N – 1) und dem Funktionsblock210 , das an den Faltungsblock220 ausgegeben wird. Der Multiplizierer235 bestimmt das Produkt aus der Ausgabe des Faltungsblocks220 und Funktionsblock230 , um Xk (wobei k = 0, 1, ..., M – 1) bereitzustellen, welches die DFT für das diskrete Signal xn ist. Der Faltungsblock220 führt die Faltungsfunktion y(n)·v(n) durch, die in Gleichung (13) gezeigt ist. Der Funktionsblock210 führt die Funktion A–nWn2/2 durch und der Funktionsblock230 führt die Funktion Wk2/2 durch, die ebenso oben mit Bezug auf Gleichung (13) bestimmt ist. -
-
3 ist ein funktionales Blockdiagramm, welches den CZT-Algorithmus aus Gleichung (14) gemäß einer Implementierung in der Frequenzdomäne darstellt.3 zeigt das Produkt aus abgetastetem, diskretem Signal xn und dem Funktionsblock210 , das durch den FFT Block321 von dem Multiplizierer215 empfangen wird. Der Multiplizierer325 bestimmt das Produkt aus der Ausgabe des FFT Blocks321 und dem FFT Funktionsblock326 , der eine FFT an der Funktion Wn2/2 durchführt. Die Ausgabe des Multiplizierers325 wird an den IFFT Block322 geliefert, der eine inverse FFT durchführt. Der Multiplizierer235 bestimmt das Produkt aus der Ausgabe des IFFT Blocks322 und dem Funktionsblock230 , um Xk bereitzustellen, welches die DFT für das diskrete Signal xn ist. - Die DFT in der z-Ebene wird durch einen Einheitskreis dargestellt. Da die CZT ein allgemeinerer Algorithmus ist, ist es offenkundig, dass sie die Z-Transformation an Punkten sowohl innerhalb als auch außerhalb des Einheitskreises berechnen kann, der beispielhaft in
11 gezeigt ist, die eine spiralförmige CZT Kontur in der z-Ebene zum Berechnen der Transformation zeigt. Um zirkuläre Faltungseffekte zu vermeiden, sind die Größen der FFTs auf L ≥ M + N – 1 beschränkt, wobei L die Größe der FFT ist. - Für eine hochaufgelöste Spektralanalyse kann der Anwender sich für die Fähigkeit des CZT-Algorithmus, effizient hochaufgelöste Spektren mit schmalem Frequenzband, gekennzeichnet durch ΔF, zu berechnen, interessieren. Bei Verwendung von üblichen FFT Verfahren werden zum Erreichen einer Frequenzauflösung von ≤ ΔF, bei Abtastfrequenz der Daten von 1/T, N ≥ 1/(T·ΔF) Punkte benötigt. Für sehr kleines ΔF bedeutet dies sehr große Werte für N. Eine hohe Auflösung kann jedoch für einen begrenzten Frequenzbereich benötigt werden, während eine geringe Auflösung für den Rest des Spektrums verwendet werden kann. Ein Beispiel für eine solche Anwendung findet sich in der Breitbandspektralanalyse, wo eine Zoomfunktion verwendet werden kann, um Einzelheiten der Frequenzantwort in einem schmalen Frequenzband von Interesse zu betrachten und der Rest des Spektrums nur grob betrachtet wird. Ein weiteres Beispiel ist der Entwurf von Bandpass- oder Tiefpassfiltern. Üblicherweise ist es wünschenswert, eine mikroskopische Betrachtung der Details der Frequenzantwort in dem Durchlassband und nur eine grobe Betrachtung außerhalb des Durchlassbandes bereitzustellen. Die Verwendbarkeit des CZT-Algorithmus für solche Frequenzspreizungen liefert ein mächtiges Werkzeug für die eingehende Untersuchung von kleinen Frequenzbändern, sowie für die Fehlerbeseitigung bei Implementierungen von Digitalfiltern. Es kann auf einfache Art und Weise geprüft werden, ob ein gewünschter Filter seine Entwurfsspezifikation bzgl. Rauschens im Band, Übergangsverhältnis, usw. erfüllt.
- Eine Situation, in der der CZT-Algorithmus hilfreich ist, ergibt sich beim Berechnen des Spektrums einer extrem langen Folge, für die ein feinkörniges Spektrum über ein schmales Frequenzband erwünscht ist, was normalerweise sehr lange FFTs erfordert. Beispielsweise kann der CZT-Algorithmus für eine Folge von P Samples (wobei P eine positive ganze Zahl ist), wenn die spektralen Samples für M << P gewünscht werden, in r Summen über N zerlegt werden, wobei rN >> P. Jede der r Summen kann unter Verwendung des CZT-Algorithmus berechnet werden, wobei ein Speicher von ungefähr 3(N + M – 1) Speicherstellen benötigt wird. Zusätzlich werden 2M Speicherstellen benötigt, um die M komplexen Werte der Transformation zu speichern. Obwohl zwei FFTs und 2M komplexwertige Multiplikationen für jede der r Transformationen benötigt werden, ist es immer noch möglich, dass eine Einsparung in der Gesamtdauer im Vergleich zur Berechnung einer P Punkte Transformation unter Verwendung eines großen externen Speichers oder Plattenspeichers erzielt werden kann. Der CZT-Algorithmus ist für solche Fälle gut geeignet, da er eine Auswahl von Anfangsfrequenz und Frequenzabstand unabhängig von der Anzahl der Zeitsamples erlaubt, wie es beispielsweise von L. R. Rabiner et al., The Chirp-Z transform Algorithm and Its Application, BELL SYSTEM TECHNICAL JOURNAL 48 (No. 5), S. 1249–1292 (Mai-Juni 1969) beschrieben wurde, was hiermit durch Verweisung aufgenommen wird. Somit können hochaufgelöste Daten kostensparend über einen schmalen Frequenzbereich erzielt werden.
- Eine alternative spektrale Zoomtechnik verwendet einen digitalen Abwärtswandler, um eine Zoomfunktion vor der FFT durch eine klassische Überlagerungstechnik, so wie sie im Stand der Technik bekannt ist, bereitzustellen. Bei dieser Technik wird das Breitbandsignal mit einem Träger bei oder nahe der Trägerfrequenz des interessierenden Bandes gemischt, um diesen Kanal auf Basisband zu legen. Das resultierende Signal wird anschließend gefiltert und dezimiert, um das interessierende Band von dem Breitbandsignal zu isolieren und zu extrahieren, und um die Gesamtabtastrate auf das Minimum zu reduzieren, das erforderlich ist, um das gewünschte Band zu erhalten. In Abhängigkeit von der Bandbreite des interessierenden Kanals können typischerweise verschiedene Filter- und Dezimierungsverfahren eingesetzt werden. Beispielsweise wird für Signale mit breiterer Bandbreite, welche Dezimierungen von acht oder weniger benötigen, ein Breitband-FIR (Finite Impulse Response) Filter direkt auf den Basisband-Mischer folgend verwendet. Die Ausgabe dieses Filters wird dann um einen geeigneten Betrag dezimiert. Für Signale mit schmalerer Bandbreite werden die Filter- und Dezimierungsfunktionen üblicherweise in mehrere Stufen aufgeteilt. Eine erste Stufe wird mittels eines kaskadierten Integrator-Differentiator-Filters (CIC) vorgesehen, der eine annehmbare Kanalisolierung in der ersten Stufe gewährt, während er die Anzahl der komplexen Rechenoperationen, die vor der Dezimierung durchgeführt werden müssen, minimiert. Ein Folgefiltern kann anschließend mittels eines Schmalband-FIR Filters vorgesehen sein, der programmierbare Taps haben kann, um eine Optimierung des Filters für das interessierende Band zu erlauben. Die Berechnungskomplexität dieser Algorithmen hängt von den Einzelheiten der jeweiligen Implementierung ab. Der Breitband DDC benötigt mindestens sechs Rechenoperationen für die Mischfunktion und sechs Rechenoperationen pro Tap (zwei Multiplikationen, zwei Additionen und zwei Verschiebungen) für das FIR Filtern. Für große Dezimierungs- und Zoomanforderungen wird ein sehr großer FIR Filter benötigt, wodurch eine Implementierung in der Frequenzdomäne unter Verwendung von FFTs leichter möglich wird. In diesem Fall kann der Ansatz mit einem CZT-Algorithmus kosteneffizienter sein, da alle Abstimmungs- und Filteroperationen in der CZT Architektur integriert sind, und er muss nicht auf Zweierpotenzen an Frequenzpunkten beschränkt sein.
- Bei der digitalen Spektralanalyse endlicher Folgen wird eine Fenstertechnik auf das Eingangssignal angewendet, um ein Frequenzleck, das aufgrund des Berechnens der DFT für Signale endlicher Länge auftritt, zu reduzieren. Die Auswirkung der Fenstertechnik auf Effizienz und Implementierung des CZT-Algorithmus wird untenstehend analysiert.
-
-
-
4 ist ein funktionales Blockdiagramm, welches den auf einer fensterbewerteten FFT beruhenden CZT-Algorithmus aus Gleichung (16) darstellt.5 ist ein funktionales Blockdiagramm, welches den in4 gezeigten, auf einer fensterbewerteten FFT beruhenden CZT-Algorithmus darstellt, wobei die Vormultiplizierer kombiniert sind. - Bezugnehmend auf
4 bestimmt der Multiplizierer405 das Produkt aus dem abgetasteten, diskreten Signal xn und dem Funktionsblock410 , der eine Fensterfunktion w(n) ausführt. Der Multiplizierer215 bestimmt wie oben beschrieben das Produkt aus der Ausgabe des Multiplizierers405 und dem Funktionsblock210 . Der FFT Block321 führt eine FFT an der Ausgabe des Multiplizierers215 durch und der Multiplizierer325 bestimmt wie oben beschrieben das Produkt aus der Ausgabe des FFT Blocks321 und dem Funktionsblock326 . Die Ausgabe des Multiplizierers325 wird an den IFFT Block322 geliefert, der eine inverse FFT durchführt. Der Multiplizierer235 bestimmt wie oben beschrieben das Produkt aus der Ausgabe des IFFT Blocks322 und dem Funktionsblock230 , um Xk bereitzustellen, welches die DFT für das diskrete Signal xn ist. In5 sind die Funktionen der Funktionsblöcke410 und210 zu einem Funktionsblock510 kombiniert, der somit die Funktion w(n)·A–nWn2/2 ausführt. - Spektralanalyseverfahren setzen voraus, dass das verarbeitete Signal während des Abtastvorgangs (d. h. während der Dauer des Beobachtungsfensters) einen konstanten Oberschwingungsanteil besitzt. Bei vielen üblichen Anwendungen sind die Signale zeitabhängig und können üblicherweise nur dann als stationär angesehen werden, wenn die Länge des Beobachtungsfensters angemessen kurz ist. Um jedoch eine gute Spektralauflösung eines Signals zu erzielen, muss die Anzahl der Samples erhöht werden, was die Länge des Beobachtungsfensters vergrößert.
- Die FFT ist nicht geeignet, wenn ein schmales Frequenzband eines nicht stationären Signals analysiert wird, da eine gute Spektralauflösung nicht mit einem kurzen Beobachtungsfenster möglich ist. Um sowohl die Spektralauflösung als auch die Auflösung bei Verwendung der FFT zu verbessern, ist ein relativ langes Beobachtungsfenster (d. h. eine große Zahl von Samples) erforderlich. Ein transientes Signal kann jedoch nicht in einem langen Fenster als stationär angesehen werden, da zeitabhängige Signale nicht konstant sind. Daher ist es ein Problem, eine gute Spektralauflösung zu erzielen, während die Länge des Beobachtungsfensters reduziert wird.
- Unter gewissen Umständen ist es eventuell nicht möglich, die Länge des Beobachtungsfensters zu erhöhen, da das Signal nur für ein kurzes Zeitintervall als stationär angenommen werden darf. Zusätzlich erfordert eine Anwendung der FFT eine Analyse des gesamten Signalspektrums, selbst wenn eine Auswertung des spektralen Gehalts des Signals nur in einem begrenzten Frequenzband benötigt wird. Da eine festgesetzte Auflösung für das betrachtete Frequenzband benötigt wird, bedeutet dies, dass das gesamte Spektrum mit derselben Auflösung analysiert werden muss, was ein langes Beobachtungsfenster hervorrufen kann. Das lange Beobachtungsfenster widerspricht der Annahme, dass das Eingangssignal als stationär angesehen werden soll.
- Der CZT-Algorithmus benötigt nicht eine Analyse des gesamten Spektrums mit derselben Auflösung, so wie es die FFT tut. Daher bietet der CZT-Algorithmus eine beschränkte Spektralanalyse mit einer erheblich besseren Spektralauflösung, da der CZT-Algorithmus, anders als klassische, konventionelle Verfahren (z. B. FFT, Welch, usw.), nicht nur von der Länge des Beobachtungsfensters abhängt, sondern eher von dem Verhältnis zwischen dem analysierten spektralen Frequenzband und der Abtastfrequenz. In diesem Fall ist es nach Auswahl eines Frequenzbands möglich, eine gute Spektralauflösung selbst mit einer erheblich geringeren Zahl von Samples und einem dementsprechend kürzeren Beobachtungsfenster zu erzielen. Auf diese Weise ermöglicht es der CZT-Algorithmus, alle Leistungsparameter zu optimieren, selbst wenn ein kürzeres Beobachtungsfenster benötigt wird.
- Der CZT-Algorithmus erlaubt es weiterhin dem Anwender, nur die abgetasteten Daten zu verwenden, um eine exaktere Signalanalyse zu erhalten, ohne irgendeine weitere Operation zum Korrigieren von Leck- oder Amplitudenfehlern. Zusätzlich ist es möglich, den CZT-Algorithmus zu verwenden, um das vollständige Signalspektrum selbst für nicht stationäre Signale wiederherzustellen, indem verschiedene CZTs auf begrenzten Spektrumsabschnitten verarbeitet werden, wobei dieselben gepufferten, abgetasteten Daten verwendet werden. Daher eignet sich der CZT-Algorithmus sowohl zur Begrenzung des spektralen Frequenzfehlers aufgrund nicht stationärer Signale als auch zur Verbesserung der Spektralauflösung.
- Insbesondere liegt ein Nachteil der oben beschriebenen konventionellen fensterbewerteten Spektralanalyse darin, dass die Länge des Beobachtungsfensters an die Länge der FFT gekoppelt ist. Daher müssen zum Verbessern der Leistung der Fenstertechnik sowohl die Länge des Beobachtungsfensters als auch der FFT vergrößert werden, was für eine Spektralanalyse in Echtzeit eventuell nicht praktikabel ist.
- Im Allgemeinen erlaubt die Verwendung von Filterbänken eine Vergrößerung der Länge des Beobachtungsfensters unabhängig von der Länge der FFT, wodurch möglicherweise eine erhöhte Leistung der Fenstertechnik mit kürzeren FFTs bereitgestellt wird.
6 ist ein funktionales Blockdiagramm, welches eine Filterbank eines Spektralanalysators darstellt. - In
6 sind Filter H0(z)631 bis Filter HN-1(z)633 Mehrphasenfilter, die eine Fensterfunktion repräsentieren. Die Mehrphasenfilter H0(z)631 bis Filter HN-1(z) können von einem einzigen Prototypfilter H(z) wie in der untenstehenden Gleichung (17) abgeleitet sein. Wenn jeder der Filter H0(z)631 bis Filter HN-1(z)633 beispielsweise die Länge eins hat, vereinfacht sich die Struktur zu einem klassischen Fenstertechnikfall, bei dem der Mehrphasenfilterkoeffizient h gleich der Fensterfunktion w ist. In der allgemeinen Filterbankstruktur kann die Länge jedes Filterkoeffizienten hk eine beliebige ganze Zahl sein. Die Mehrphasenfilter werden wie in Gleichung (17) gezeigt von den Filterkoeffizienten hk abgeleitet:hk(n) = h(k + nN) (17) -
6 zeigt weiterhin ein diskretes Signal x(n), das über den Downsampling-Block623 (IN) auf den Filter HN-1(z)633 , über den Verzögerungsblock611 (Z–1) und den Downsampling-Block622 (IN) auf den Filter H1(z)632 , und über Verzögerungsblöcke611 und612 (Z–1) (sowie etwaige dazwischen liegende Verzögerungsblöcke) und den Downsampling-Block621 (IN) auf den Filter H0(z)631 angewendet wird. Jeder der Downsampling-Blöcke621 bis623 behält eines von N Samples bei und jeder der Verzögerungsblöcke611 bis612 verzögert die Daten x(n) um ein Sample und kann beispielsweise unter Verwendung eines Speicherregisters oder eines Speichers implementiert werden. Die Ausgabe von jedem der Filter H0(z)631 bis Filter HN-1(z)633 wird an einen N-Punkt FFT Block640 geliefert, der entsprechende DFTs X0 bis XN-1 für das diskrete Signal x(n) ausgibt. -
7 ist ein funktionales Blockdiagramm, welches eine kontinuierliche Filterbank eines Spektralanalysators darstellt und somit eine effizientere Umsetzung der6 zeigt, in der eine Filterstruktur verwendet wird, während sich K Filterkoeffizienten pro Multiplizierer abwechseln, was als eine auf einer kontinuierlichen Filterbank beruhende FFT Analyse bezeichnet werden kann. - In
7 wird ein Eingangssignal mit einer Abtastrate FS empfangen. Das Eingangssignal wird fortlaufend einer Reihe von Puffer zugeführt, die entsprechende Verzögerungen aufweisen und durch Puffer731 bis733 angedeutet sind. Der Puffer731 verzögert das Eingangssignal um eine erste Verzögerung, der Puffer732 verzögert das von dem Puffer731 ausgegebene, verzögerte Signal um eine zweite Verzögerung, und der (N – 1)-te Puffer733 verzögert das von dem vorhergehenden Puffer (z. B. Puffer732 ) ausgegebene, verzögerte Signal um eine (N – 1)-te Verzögerung. Zusätzlich multipliziert der Multiplizierer721 das Eingangssignal mit mehreren Taps, die den Filterkoeffizienten h0, h1...hN-1 entsprechen. Der Multiplizierer722 multipliziert das von Puffer731 ausgegebene, verzögerte Signal mit mehreren Taps, die den Filterkoeffizienten hK, hK+1, ... h2N-1 entsprechen. Der Multiplizierer723 multipliziert das von Puffer732 ausgegebene, verzögerte Signal mit mehreren Taps, die den Filterkoeffizienten h2K, h2K+1, ... h3N-1 entsprechen. Der Multiplizierer724 multipliziert das von Puffer733 ausgegebene, verzögerte Signal mit mehreren Taps, die den Filterkoeffizienten h(N-1)K, h(N-1)K+1...hsN-1 entsprechen. Die Ausgaben der Multiplizierer721 bis724 werden durch den Addierer740 addiert (oder summiert), um I-Daten der IQ-Daten des Eingangssignals bereitzustellen. Währenddessen wird ein entsprechender Prozess an den Q-Daten der IQ-Daten aus dem Eingangssignal wie durch Block760 angedeutet durchgeführt. Die I-Daten und die Q-Daten werden an den N-Punkt FFT Block750 geliefert, der IQ-Daten in der Frequenzdomäne aus dem Eingangssignal als Mehrphasen-DFT des diskreten Signals x(n) ausgibt. -
3 und7 können wie in8 gezeigt kombiniert werden, um einen auf einer allgemeinen Filterbank beruhenden CZT-Algorithmus umzusetzen. Insbesondere ist8 ein funktionales Blockdiagramm, welches einen auf einer kontinuierlichen Filterbank beruhenden CZT-Algorithmus mit doppelter Fensterlänge darstellt, der den in5 gezeigten, auf einer fensterbewerteten FFT beruhenden CZT-Algorithmus und die in7 gezeigte kontinuierliche Filterbank kombiniert. Der auf einer kontinuierlichen Filterbank beruhende CZT-Algorithmus besitzt wie unten beschrieben die doppelte Fensterlänge aufgrund einer Implementierung von zwei Mehrphasenfiltern. - Mit Bezug auf
8 umfasst eine Filterschaltung die Multiplizierer815 und825 , zwei Mehrphasenfilterblöcke810 und820 , den Puffer830 und den Addierer840 . Eine CZT-Schaltung umfasst die Multiplizierer215 ,325 und235 , die Funktionsblöcke510 ,326 und230 , den FFT Block321 und den IFFT Block322 , die oben in Zusammenhang mit den4 und5 erläutert wurden. - Der Multiplizierer
815 ist dazu ausgebildet, das Produkt von jedem von 2N Samples des diskreten Signals xn (wobei n = 0, 1, ..., 2N – 1) mit dem Mehrphasenfilterblock810 zu bilden, der die Filterkoeffizienten h(0...N – 1) implementiert. Der Multiplizierer825 ist dazu ausgebildet, das Produkt von N Samples des diskreten, durch den Puffer830 verzögerten Signals xn und dem Mehrphasenfilterblock820 zu bilden, der die Filterkoeffizienten h(N...2N) implementiert. Die von den Multiplizierern815 und825 ausgegebenen Produkte werden mittels des Addierers840 addiert und an den Multiplizierer215 ausgegeben. Der Multiplizierer215 bildet das Produkt aus der Ausgabe des Addierers840 und dem Funktionsblock510 , das an den FFT Block321 weitergegeben wird. Der Rest der CZT-Verarbeitung wird wie oben in Hinblick auf die3 bis5 beschrieben durchgeführt. -
9 ist ein funktionales Blockdiagramm, welches den in8 gezeigten, auf einer kontinuierlichen Filterbank beruhenden CZT-Algorithmus darstellt, wobei der Vormultiplizierer der CZT-Schaltung mit den Filterkoeffizienten in der Filterschaltung kombiniert ist. Anders ausgedrückt zeigt9 eine Kombination der Funktionen des Funktionsblocks510 der CZT-Schaltung und des Filterblocks810 der Filterschaltung in einem kombinierten Funktionsblock910 , welcher dann also die Funktion h(0...N – 1)·(A–nWn)/2 ausführt. Entsprechend sind die Funktionen des Funktionsblocks510 der CZT-Schaltung und des Filterblocks820 der Filterschaltung in einem kombinierten Funktionsblock920 vereint, welcher dann also die Funktion h(N...2N)·(A–nWn)/2 ausführt. Jede der kombinierten Funktionen wird durch zwei dividiert, um die zwei Filterblöcke810 und820 (und die doppelte Fensterlänge) zu berücksichtigen. -
10 ist ein funktionales Blockdiagramm, welches einen auf einer allgemeinen kontinuierlichen Filterbank beruhenden CZT-Algorithmus für einen Filter der Länge sN darstellt, wobei s eine ganze Zahl größer oder gleich eins ist. Im Vergleich hierzu stellt9 den Spezialfall für s = 1 dar. Wenn s = 1, liegt der klassische fensterbewertete Fall vor, und wenn alle Filterkoeffizienten h gleich eins sind, reduziert sich die Struktur auf den grundlegenden CZT-Algorithmus in3 . Ohne Einschränkung der Allgemeinheit kann die in10 dargestellte Struktur verwendet werden, um die DFT für das diskrete Signal xn, die mit Xk bezeichnet sind und die auf kreisförmigen oder spiralförmigen Konturen, die an einem frei wählbaren Punkt in der z-Ebene beginnen, liegen, auszuwerten, indem einfach die A0 und W0 Parameter wie in11 gezeigt modifiziert werden, für welche Anwendungen jenseits der klassischen, auf DFT basierenden Spektralanalyse denkbar sind. -
- Insbesondere ist wie in
10 gezeigt der CZT Vormultiplizierer aus8 in der Filterschaltung kombiniert, die so verallgemeinert wurde, dass sie s + 1 Filterblöcke umfasst, die als stellvertretende, kombinierte Filterblöcke1010 ,1020 und1030 gezeigt sind, welche stellvertretende Mehrphasenfilter sind. Somit führt der kombinierte Filterblock1010 die Funktion h(0...N – 1)·(A–nWn)/(s + 1) aus, der kombinierte Filterblock1020 führt die Funktion h(N...2N)·(A–nWn)/(s + 1) aus, und der kombinierte Filterblock1030 führt die Funktion h(N...sN)·(A–nWn)/(s + 1) aus. Der Multiplizierer1015 empfängt sN Samples des diskreten Signals xn (wobei n = 0, 1, ..., sN – 1) und ist dazu ausgebildet, das Produkt aus dem abgetasteten, diskreten Signal xn und dem kombinierten Filterblock1010 zu bilden. Der Multiplizierer1025 ist dazu ausgebildet, das Produkt von N Samples des diskreten, durch den Puffer1031 verzögerten Signals xn mit dem kombinierten Filterblock1020 zu bilden. Der Multiplizierer1035 ist dazu ausgebildet, das Produkt von N Samples des diskreten, durch alle vorhergehenden Puffer, die durch die stellvertretenden Puffer1031 und1032 angedeutet sind, verzögerten Signals xn mit dem kombinierten Filterblock1030 zu bilden. Die von den Multiplizierern1015 ,1025 und1035 ausgegeben Produkte werden mittels des Addierers1040 addiert und an den FFT Block321 ausgegeben. Der Rest der CZT Verarbeitung wird wie oben erläutert durchgeführt, um Xk (wobei k = 0, 1, ..., M – 1) bereitzustellen, welches die DFT für das diskrete Signal xn ist. - Folglich gestatten die vielfältigen Weiterbildung eine Zoomdarstellung einer beliebigen schmalen Bandbreite des abgetasteten Signals, beispielsweise durch den Anwender ausgewählt, während sie es gleichzeitig ermöglichen, Merkmale bei über die gesamte Bandbreite des abgetasteten Signals erhöhter Auflösung zu identifizieren. Zurückverweisend auf
6 stellen beispielsweise die Ausgabesignale Xk der Filterbank unabhängige, frequenzselektive Kanäle dar, die durch die Form des Mehrphasenfilters bestimmt sind, der die gesamte Bandbreite in N gleichmäßig verteilte Kanäle unterteilt. Mit der Struktur der in10 offenbarten Weiterbildung kann das System auch als ein kanalisierendes System mit beliebiger Bandbreite beschrieben werden, welches N Kanäle zwischen beliebigen Anfangs- und Endfrequenzen, nicht nur die gesamte Bandbreite, erzeugt. - Während die Offenbarung Bezug auf beispielhafte Weiterbildungen nimmt, ist es für den Fachmann offenbar, dass vielfältige Änderungen und Modifikationen durchgeführt werden können, ohne vom Geist und vom Umfang der vorliegenden Lehre abzuweichen. Daher sollte Einverständnis darüber herrschen, dass die oben beschriebenen Weiterbildungen nicht einschränkend, sondern veranschaulichend sind.
- ZITATE ENTHALTEN IN DER BESCHREIBUNG
- Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
- Zitierte Nicht-Patentliteratur
-
- L. R. Rabiner et al., The Chirp-Z transform Algorithm and Its Application, BELL SYSTEM TECHNICAL JOURNAL 48 (No. 5), S. 1249–1292 (Mai-Juni 1969) [0049]
Claims (13)
- Ein Verfahren zum Durchführen von Spektralanalyse eines nicht stationären Signals in Echtzeit, wobei das Verfahren umfasst: Abtasten des nicht stationären Signals unter Verwendung eines Beobachtungsfensters mit einer Länge, die kurz genug ist, um ein stationäres Signal zu approximieren, zum Bereitstellen eines ersten Satzes abgetasteter Daten (S110); Pufferung des ersten Satzes abgetasteter Daten, um eine Vielzahl von gepufferten Sätzen abgetasteter Daten zu erhalten (S120); Filtern des ersten Satzes abgetasteter Daten und der Vielzahl von gepufferten Sätzen abgetasteter Daten unter Verwendung einer entsprechenden Vielzahl von Filterantworten, um eine Vielzahl von gefilterten Sätzen abgetasteter Daten zu erhalten (S310); Durchführen einer Chirp-Z-Transformation (CZT) der Vielzahl von gefilterten Sätzen abgetasteter Daten, um einen Satz von Koeffizienten von Diskreten Fouriertransformationen (DFT) bereitzustellen (S150); und Wiederherstellen eines vollständigen Signalspektrums des nicht stationären Signals unter Verwendung des Satzes von DFT Koeffizienten (S160).
- Das Verfahren nach Anspruch 1, wobei sich die Vielzahl von Filterantworten von einem einzigen Filter ableiten.
- Das Verfahren nach Anspruch 1 oder 2, wobei die Pufferung des ersten Satzes abgetasteter Daten ein Verzögern des ersten Satzes abgetasteter Daten unter Verwendung einer Vielzahl von Puffer umfasst, wobei jeder Puffer einem aus der Vielzahl von gepufferten Sätzen abgetasteter Daten entspricht.
- Das Verfahren nach einem der Ansprüche 1 bis 3, wobei jeder von dem ersten Satz abgetasteter Daten und der Vielzahl von gepufferten Sätzen abgetasteter Daten einem begrenzten Spektrumsabschnitt des vollständigen Signalspektrums des nicht stationären Signals entspricht.
- Das Verfahren nach einem der Ansprüche 1 bis 4, wobei eine Spektralauflösung aufgrund eines im Verhältnis größeren Eingangsdatensatzes, der den ersten Satz abgetasteter Daten und die Vielzahl von gepufferten Sätzen abgetasteter Daten umfasst, verbessert wird, wodurch ein Frequenzauflösungsfehler abgeschwächt wird.
- Das Verfahren nach einem der Ansprüche 1 bis 5, wobei das Durchführen der CZT der Vielzahl von gefilterten Sätzen abgetasteter Daten umfasst: Addieren der Vielzahl von gefilterten Sätzen abgetasteter Daten; Durchführen einer Fast Fourier Transformation (FFT) an den addierten, gefilterten Sätzen abgetasteter Daten; Filtern der FFT Sätze abgetasteter Daten unter Verwendung einer FFT einer komplexwertigen Chirp-Funktion; Durchführen einer inversen FFT (IFFT) an den gefilterten FFT Sätzen abgetasteter Daten; und Filtern der IFFT Sätze abgetasteter Daten unter Verwendung einer weiteren komplexwertigen Chirp-Funktion, um den Satz von DFT Koeffizienten bereitzustellen.
- Das Verfahren nach Anspruch 6, wobei das Filtern des ersten Satzes abgetasteter Daten umfasst: Multiplizieren eines Filterkoeffizienten für jeden abgetasteten Datenpunkt in dem ersten Satz abgetasteter Daten mit einem Quotienten aus einer komplexen Zahl, auf der Grundlage einer Abtastfrequenz des ersten Satzes abgetasteter Daten, und der Anzahl der gepufferten Sätze plus den ersten Satz abgetasteter Daten.
- Das Verfahren nach Anspruch 7, wobei das Filtern jedes der Vielzahl von gepufferten Sätzen abgetasteter Daten umfasst: Multiplizieren eines Filterkoeffizienten für jeden abgetasteten Datenpunkt in jedem der gepufferten Sätze abgetasteter Daten mit dem Quotienten aus der komplexen Zahl und der Anzahl der gepufferten Sätze plus den ersten Satz abgetasteter Daten.
- Das Verfahren nach einem der Ansprüche 1 bis 8, wobei das Wiederherstellen des vollständigen Signalspektrums des nicht stationären Signals ein Berechnen des Satzes von DFT Koeffizienten in einer z-Ebene umfasst, wobei der Satz von DFT Koeffizienten auf kreisförmigen oder spiralförmigen Konturen, die an einem frei wählbaren Punkt in der z-Ebene beginnen, liegt.
- Das Verfahren nach einem der Ansprüche 1 bis 9, wobei eine Dauer des Beobachtungsfensters kleiner als eine Unterbrechungsdauer oder eine Abklingdauer des nicht stationären Signals ist.
- Ein System zum Durchführen einer Spektralanalyse eines Eingangssignals von einem Messobjekt (DUT = Device Under Test) in Echtzeit, wobei das System umfasst: einen Abtaster, der dazu ausgebildet ist, das Eingangssignal abzutasten, um erste abgetastete Daten bereitzustellen; einen ersten Mehrphasenfilter (
1010 ), der dazu ausgebildet ist, die ersten abgetasteten Daten zu filtern, um erste gefilterte, abgetastete Daten bereitzustellen; einen ersten Puffer (1031 ), der dazu ausgebildet ist, die ersten abgetasteten Daten um eine erste Verzögerung zu verzögern, um zweite abgetastete Daten bereitzustellen; einen zweiten Mehrphasenfilter (1020 ), der dazu ausgebildet ist, die zweiten abgetasteten Daten zu filtern, um zweite gefilterte, abgetastete Daten bereitzustellen; einen zweiten Puffer (1032 ), der dazu ausgebildet ist, die zweiten abgetasteten Daten um eine zweite Verzögerung zu verzögern, um dritte abgetastete Daten bereitzustellen; einen dritten Mehrphasenfilter (1030 ), der dazu ausgebildet ist, die dritten abgetasteten Daten zu filtern, um dritte gefilterte, abgetastete Daten bereitzustellen; einen Addierer (1040 ), der dazu ausgebildet ist, die ersten, zweiten und dritten gefilterten, abgetasteten Daten zu addieren, um addierte, gefilterte Daten bereitzustellen; und eine Chirp-Z-Transformations-Schaltung (CZT) (321 ,322 ,325 ,326 ,230 ,235 ), die dazu ausgebildet ist, eine CZT an den addierten, gefilterten Daten durchzuführen, um Frequenzdomänensamples bereitzustellen, die dem Eingangssignal entsprechen, wobei die Frequenzdomänensamples ein Wiederherstellen eines vollständigen Signalspektrums des Eingangssignals ermöglichen. - Das System nach Anspruch 11, wobei der Abtaster einen Analog-digital-Wandler aufweist.
- Das System nach Anspruch 11 oder 12, wobei die CZT-Schaltung umfasst: eine Fast Fourier Transformationsschaltung (FFT), die dazu ausgebildet ist, eine FFT Operation an den addierten, gefilterten Daten durchzuführen, um ein FFT Signal bereitzustellen; einen vierten Filter, der dazu ausgebildet ist, das FFT Signal zu filtern, um ein gefiltertes FFT Signal bereitzustellen; eine Inverse-FFT (IFFT) Schaltung, die dazu ausgebildet ist, eine IFFT Operation an dem gefilterten FFT Signal durchzuführen, um ein IFFT Signal bereitzustellen; und einen fünften Filter, der dazu ausgebildet ist, das IFFT Signal zu filtern, um die Frequenzdomänensamples bereitzustellen.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/661,834 US10371732B2 (en) | 2012-10-26 | 2012-10-26 | Method and system for performing real-time spectral analysis of non-stationary signal |
US13/661,834 | 2012-10-26 |
Publications (1)
Publication Number | Publication Date |
---|---|
DE102013217181A1 true DE102013217181A1 (de) | 2014-04-30 |
Family
ID=50479869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE201310217181 Pending DE102013217181A1 (de) | 2012-10-26 | 2013-08-28 | Verfahren und System zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit |
Country Status (3)
Country | Link |
---|---|
US (1) | US10371732B2 (de) |
CN (1) | CN103792427B (de) |
DE (1) | DE102013217181A1 (de) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106546817A (zh) * | 2016-11-08 | 2017-03-29 | 电子科技大学 | 一种具有反馈功能的频率估计和能量估计方法 |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104007317A (zh) * | 2014-06-04 | 2014-08-27 | 哈尔滨工业大学 | 一种获取伺服***频率特性的方法及装置 |
CN104007318B (zh) * | 2014-06-17 | 2016-11-09 | 中国科学院电子学研究所 | 获取信号时频函数的方法 |
CN105158494A (zh) * | 2015-08-06 | 2015-12-16 | 浙江工业大学 | 一种往复式制冷压缩机转速测量的方法 |
CN105137175B (zh) * | 2015-08-06 | 2017-11-07 | 湖南大学 | 一种加窗短时傅里叶变换三点插值动态频率测量方法 |
JP6634159B2 (ja) * | 2015-12-18 | 2020-01-22 | オリンパス株式会社 | 改善されたスライディングdftを使用する正確かつ効率的なスペクトル推定のための方法および装置 |
CN106018955B (zh) * | 2016-05-07 | 2018-11-13 | 天津大学 | 快速卷积可调滤波器组的低速率端频率估计方法 |
CN106092492B (zh) * | 2016-06-01 | 2018-08-28 | 南京航空航天大学 | 一种滤波去噪方法 |
US10330773B2 (en) * | 2016-06-16 | 2019-06-25 | Texas Instruments Incorporated | Radar hardware accelerator |
JP6928208B2 (ja) * | 2017-07-24 | 2021-09-01 | アイオワ、ステイト、ユニバーシティー、リサーチ、ファウンデーション、インコーポレイテッドIowa State University Research Foundation,Inc. | O(nlog(n))時間及びびO(n)メモリでチャープZ変換を逆変換するためのシステム及び方法 |
DE102018127132A1 (de) | 2018-10-30 | 2020-04-30 | Sma Solar Technology Ag | Wechselrichter mit mindestens zwei Gleichspannungswandlern und Verwendung eines solchen Wechselrichters in einer Photovoltaikanlage |
CN111579868B (zh) * | 2020-05-25 | 2022-06-10 | 云南电网有限责任公司电力科学研究院 | 一种高次谐波的测量方法及装置 |
CN112328963A (zh) * | 2020-09-29 | 2021-02-05 | 国创新能源汽车智慧能源装备创新中心(江苏)有限公司 | 信号的有效值计算方法、装置 |
CN112559439B (zh) * | 2020-12-21 | 2024-04-26 | 中国电子科技集团公司第十四研究所 | 一种实现czt变换的fpga架构 |
CN112664177B (zh) * | 2020-12-29 | 2022-06-21 | 重庆邮电大学 | 一种基于归一化卷积与自适应滤波的油井动液面测量方法 |
CN112883326B (zh) * | 2021-03-19 | 2022-07-08 | 吉林大学 | 一种基于流式算法的自适应时频变换方法 |
US11874308B2 (en) * | 2021-10-29 | 2024-01-16 | Keysight Technologies, Inc. | ASIC implementing real-time spectrum analysis |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO831718L (no) | 1983-05-13 | 1984-11-14 | Vingmed As | Fremgangsmaate og apparat ved blodstroem-hastighetsmaaling med ultralyd for dannelse av todimensjonal avbildning av blodets hastighet |
US5075619A (en) | 1990-04-06 | 1991-12-24 | Tektronix, Inc. | Method and apparatus for measuring the frequency of a spectral line |
MX9702434A (es) * | 1991-03-07 | 1998-05-31 | Masimo Corp | Aparato de procesamiento de señales. |
US5495432A (en) * | 1994-01-03 | 1996-02-27 | Industrial Technology Research Institute | Apparatus and method for sampling rate conversion |
US8019400B2 (en) | 1994-10-07 | 2011-09-13 | Masimo Corporation | Signal processing apparatus |
US6218896B1 (en) * | 1999-08-27 | 2001-04-17 | Tachyon, Inc. | Vectored demodulation and frequency estimation apparatus and method |
DE102004051386A1 (de) * | 2004-09-28 | 2006-04-06 | Rohde & Schwarz Gmbh & Co. Kg | Verfahren und Vorrichtung zur Spektrumanalyse in mehreren Frequenzbändern mit verschiedener Frequenzauflösung |
US7774396B2 (en) * | 2005-11-18 | 2010-08-10 | Dynamic Hearing Pty Ltd | Method and device for low delay processing |
US7612275B2 (en) | 2006-04-18 | 2009-11-03 | Nokia Corporation | Method, apparatus and computer program product for providing rhythm information from an audio signal |
US20080052335A1 (en) * | 2006-08-01 | 2008-02-28 | Gee Edward C | Systems and methods for time domain to frequency domain conversion using frequency shifting |
CN101231315A (zh) | 2007-01-24 | 2008-07-30 | 涂亚庆 | 频率估计的多段采样信号融合处理方法 |
US8249386B2 (en) | 2008-03-28 | 2012-08-21 | Tektronix, Inc. | Video bandwidth resolution in DFT-based spectrum analysis |
JP4512835B2 (ja) * | 2008-05-02 | 2010-07-28 | テクトロニクス・インターナショナル・セールス・ゲーエムベーハー | コンスタレーション表示方法及び装置 |
US8255179B2 (en) * | 2009-02-11 | 2012-08-28 | Tektronix, Inc. | Time qualified frequency mask trigger |
US8626810B2 (en) * | 2009-05-15 | 2014-01-07 | Texas Instruments Incorporated | Method and system for finite impulse response (FIR) digital filtering |
US8179305B2 (en) * | 2009-08-12 | 2012-05-15 | Tektronix, Inc. | Enhanced impulse response measurement of an FM radar transmitter pulse |
WO2011026004A1 (en) * | 2009-08-31 | 2011-03-03 | Bae Systems National Security Solutions Inc. | Digital signal conditioning |
US9031834B2 (en) * | 2009-09-04 | 2015-05-12 | Nuance Communications, Inc. | Speech enhancement techniques on the power spectrum |
US8090337B2 (en) * | 2010-01-25 | 2012-01-03 | Bae Systems Information And Electronic Systems Integration Inc. | Chirp fourier transform method and apparatus for canceling wide band interference |
US8521460B2 (en) * | 2010-09-28 | 2013-08-27 | Tektronix, Inc. | Multi-domain test and measurement instrument |
US9083471B2 (en) * | 2011-03-03 | 2015-07-14 | Technion Research And Development Foundation Ltd. | Coherent and self-coherent signal processing techniques |
CN102650658B (zh) | 2012-03-31 | 2014-03-19 | 机械工业第三设计研究院 | 一种时变非平稳信号时频分析方法 |
US8694306B1 (en) * | 2012-05-04 | 2014-04-08 | Kaonyx Labs LLC | Systems and methods for source signal separation |
US9886419B2 (en) * | 2012-07-26 | 2018-02-06 | Tektronix, Inc. | System for improving probability of transient event detection |
-
2012
- 2012-10-26 US US13/661,834 patent/US10371732B2/en active Active
-
2013
- 2013-08-28 DE DE201310217181 patent/DE102013217181A1/de active Pending
- 2013-09-26 CN CN201310445073.7A patent/CN103792427B/zh active Active
Non-Patent Citations (1)
Title |
---|
L. R. Rabiner et al., The Chirp-Z transform Algorithm and Its Application, BELL SYSTEM TECHNICAL JOURNAL 48 (No. 5), S. 1249-1292 (Mai-Juni 1969) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106546817A (zh) * | 2016-11-08 | 2017-03-29 | 电子科技大学 | 一种具有反馈功能的频率估计和能量估计方法 |
CN106546817B (zh) * | 2016-11-08 | 2019-01-22 | 电子科技大学 | 一种具有反馈功能的频率估计和能量估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103792427B (zh) | 2018-04-06 |
US20140122025A1 (en) | 2014-05-01 |
US10371732B2 (en) | 2019-08-06 |
CN103792427A (zh) | 2014-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102013217181A1 (de) | Verfahren und System zum Durchführen einer Spektralanalyse eines nicht stationären Signals in Echtzeit | |
DE102006047197B3 (de) | Vorrichtung und Verfahren zum Verarbeiten eines reellen Subband-Signals zur Reduktion von Aliasing-Effekten | |
DE112008003098B4 (de) | Verfahren und Vorrichtung zur Berechnung von Interpolationsfaktoren in Abtastratenwandlungssystemen | |
DE60317722T2 (de) | Verfahren zur Reduzierung von Aliasing-Störungen, die durch die Anpassung der spektralen Hüllkurve in Realwertfilterbanken verursacht werden | |
DE69328828T2 (de) | Verfahren und gerät zur analyse von zeitlich veränderlichem frequenzspektrum | |
DE60221927T2 (de) | Vorrichtung und Programm zur Schallcodierung | |
DE10236898A1 (de) | Videofiltern mit verbesserter Spur unter Verwendung von Wavelet-Entrauschungstechniken | |
DE102019119776B4 (de) | Zeitverschachtelte digital-analog-wandler-korrektur | |
EP1794602B1 (de) | Verfahren und vorrichtung zur spektrumanalyse eines nutz- oder rauschsignals | |
EP0696848A1 (de) | Verfahren zur digitalen Interpolation von Signalen | |
DE102005012977B4 (de) | Verfahren und System zur Rauschmessung mit den kombinierbaren Unterverfahren Messen, Identifizieren und Beseitigen von sinusförmigen Störsignalen in einem Rauschsignal | |
DE2608515A1 (de) | Doppel ungerade diskrete fourier- transformationsanordnung | |
EP3459173B1 (de) | Verfahren und vorrichtung zum senden bzw. empfangen von mindestens einem hochfrequenzsignal mit paralleler und unterabgetasteter basisbandsignalverarbeitung | |
DE69028782T2 (de) | Digitales filter und entwicklungsmethode | |
DE102010002111A1 (de) | Verfahren und Anordnung zur Verteilung der Rechenlast in Datenverarbeitungseinrichtungen bei einer Durchführung von blockbasierten Rechenvorschriften sowie ein entsprechendes Computerprogramm und ein entsprechendes computerlesbares Speichermedium | |
DE4317869B4 (de) | System zur Bestimmung einer Instabilität in einem IIR-Filter | |
EP0608281B1 (de) | Verfahren zur reduzierung des frequenzübersprechens bei der übertragung und/ oder speicherung akustischer oder optischer signale | |
EP3332478A1 (de) | Vorrichtung zum verändern einer abtastrate, system umfassend eine vorrichtung zum verändern einer abtastrate und verfahren zum verändern einer abtastrate | |
DE69515509T2 (de) | Sprachverarbeitung | |
DE60210479T2 (de) | Audiokodierer mit unregelmässiger filterbank | |
EP1537750B1 (de) | Verfahren zum bestimmen der hüllkurve eines modulierten signals | |
DE68928228T2 (de) | Verfahren und Gerät, um zwischen Datenproben zu interpolieren | |
DE4022381C2 (de) | Verwendung langer Digitalfilter bei Vorkommnis von Abrundungsfehlern | |
DE102018208465A1 (de) | Verfahren und system zum verhindern von durch spiegelfrequenzen verursachten interferenzen | |
DE69522883T2 (de) | Vorrichtung zur bestimmung des signalspektrums eines digitalen breitbandsignals und zur ableitung einer bitzuweisungsinformation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R163 | Identified publications notified |
Effective date: 20140618 |
|
R082 | Change of representative |
Representative=s name: SCHOPPE, ZIMMERMANN, STOECKELER, ZINKLER & PAR, DE |
|
R082 | Change of representative |
Representative=s name: SCHOPPE, ZIMMERMANN, STOECKELER, ZINKLER & PAR, DE |
|
R081 | Change of applicant/patentee |
Owner name: KEYSIGHT TECHNOLOGIES, INC. (N.D.GES.D.STAATES, US Free format text: FORMER OWNER: AGILENT TECHNOLOGIES, INC., SANTA CLARA, CALIF., US Effective date: 20141111 |
|
R082 | Change of representative |
Representative=s name: SCHOPPE, ZIMMERMANN, STOECKELER, ZINKLER & PAR, DE Effective date: 20141111 Representative=s name: SCHOPPE, ZIMMERMANN, STOECKELER, ZINKLER & PAR, DE Effective date: 20140808 Representative=s name: SCHOPPE, ZIMMERMANN, STOECKELER, ZINKLER, SCHE, DE Effective date: 20140808 Representative=s name: SCHOPPE, ZIMMERMANN, STOECKELER, ZINKLER, SCHE, DE Effective date: 20141111 |
|
R012 | Request for examination validly filed |