Einheit 4 / 9

Signalverarbeitung und Datenanalyse

Gewinne:

  • Fähigkeit, Signalverarbeitungsparameter wie Abtastfrequenz, Filtertyp und Fensterung mit KI genau festzulegen
  • Möglichkeit, Python-Code für FFT, Filterdesign und Rauschanalyse zu generieren sowie Nyquist und Aliasing zu überprüfen
  • Möglichkeit zur Überprüfung der KI-generierten Spektruminterpretation durch Vergleich mit Messungen und physikalischen Erwartungen

Rohdaten von einem Sensor, einem Mikrofon oder einer Stromleitung allein sind oft bedeutungslos; Es ist notwendig, es auf der Achse von Frequenz, Amplitude und Zeit zu verarbeiten und zu interpretieren. Spektrum mit FFT extrahieren, Rauschen filtern, eine Vibrationssignatur erkennen … Das alles ist Signalverarbeitung und KI generiert in diesem Bereich schnellen Python-Code, erklärt Konzepte und berechnet Filterkoeffizienten. Aber ein kleiner Parameterfehler in der Signalverarbeitung (falsche Abtastfrequenz, falsches Fenster, übersprungene Nyquist-Steuerung) macht die gesamte Analyse stillschweigend zunichte; Der Code gibt keine Fehler, sondern nur falsche Ergebnisse. In dieser Einheit behandeln wir die Einstellung der richtigen Parameter mit KI, die Überprüfung des generierten Codes im Hinblick auf Nyquist und Aliasing sowie die Überprüfung der Spektruminterpretation anhand der physikalischen Realität.

Grundparameter: Abtastung, Nyquist, Aliasing

Die erste Regel der digitalen Signalverarbeitung ist das Nyquist-Shannon-Theorem: Um ein Signal genau darzustellen, muss die Abtastfrequenz mindestens doppelt so hoch sein wie die höchste Frequenzkomponente im Signal. Andernfalls erscheinen hohe Frequenzen als niedrige Frequenzen (Aliasing) und dies kann nicht rückgängig gemacht werden.

Höchste Frequenz im Signal: f_max = 2 kHzMinimale Abtastfrequenz: f_s ≥ 2 · f_max = 4 kHzUm in der Praxis Spielraum zu lassen: f_s ≈ 5 · f_max = 10 kHz wird ausgewähltAußerdem ist ein analoger Anti-Aliasing-Filter ein MUSS VOR der Abtastung.

Achtung: Wenn die KI eine Abtastfrequenz vorschlägt, prüfen Sie unbedingt, ob sie Nyquist basierend auf der höchsten Komponente Ihres Signals liefert. Auch bei geringer Abtastfrequenz läuft der Code fehlerfrei; aber das Ergebnis ist aufgrund von Aliasing falsch. Darüber hinaus kann keine Software diesen Fehler ohne einen Anti-Aliasing-Filter in der Hardware rückgängig machen.

Auflösung und Fensterung für FFT

Die Frequenzauflösung der FFT hängt von der Abtastfrequenz und der Anzahl der Abtastwerte ab: Δf = f_s / N. Für eine feinere Auflösung sammeln Sie entweder mehr Abtastwerte oder wählen niedrigere f_s (ohne Nyquist zu brechen). Darüber hinaus führt die Diskontinuität zwischen dem Anfang und dem Ende des Signals zu einer „spektralen Leckage“; Um dies zu reduzieren, werden Fensterfunktionen wie Hann und Hamming angewendet.

f_s = 10 kHz, N = 1024 SamplesFrequenzauflösung: Δf = 10000 / 1024 ≈ 9,77 HzGesamtaufnahmezeit: T = N / f_s = 1024 / 10000 ≈ 102,4 ms

Generieren und Überprüfen von FFT- und Filtercode

Geben Sie die Abtastfrequenz, den Fenstertyp und die Achsenskala explizit an, wenn Sie den FFT-Code von der KI anfordern. Der folgende Code gibt das einseitige Amplitudenspektrum eines Signals aus:

numpy importieren als npfs = 10000 # Abtastfrequenz (Hz) – Nyquist überprüfen! np.hanning(N) # Fensterung: Reduziert spektrales Leck normalisiert durch Verstärkung# Spitzenfrequenzen überprüfenpeaks = f[np.argsort(mag)[-2:]]print("Dominante Frequenzen (Hz):", np.sort(peaks))

Bei der Überprüfung dieses Codes werden Sie gefragt: Ist die Frequenzachse wirklich Hz oder wird der Sample-Index (rfftfreq) verwendet? Ist die Amplitude auf die Fensterverstärkung normiert? Treten an diesen Stellen tatsächlich die erwarteten 500- und 1500-Hz-Spitzen auf? Wenn Sie den Code mit einem bekannten Testsignal (Einzelfrequenz-Sinus) überprüfen, können Sie die Achse und die Normalisierung sicher bestätigen.

Schwache Eingabeaufforderung / Starke Eingabeaufforderung

SCHWACH: „Wenden Sie FFT auf diese Daten an.“ (Ergebnis: Abtastfrequenz, Achsenskalierung und Fensterung sind unklar; das Diagramm wird höchstwahrscheinlich auf der falschen Achse dargestellt.) STARK: „Wenden Sie FFT auf ein abgetastetes 10-kHz-Signal an. Stellen Sie Folgendes bereit: – Zeichnen Sie die Frequenzachse in Hz auf (verwenden Sie rfftfreq). – Wenden Sie das Hann-Fenster an und normalisieren Sie die Amplitude durch die Fensterverstärkung. – Drucken Sie die beiden dominanten Frequenzen digital aus. – Markieren Sie die Nyquist-Grenze (5 kHz). Beachten Sie, dass sich die Achse verschiebt, wenn die Abtastfrequenz falsch ist.

Filterdesign und Geräuschanalyse

KI kann Niedrig-/Hoch-/Bandpassfilterkoeffizienten (FIR/IIR) erzeugen. Allerdings müssen Grenzfrequenz, Ordnung und Stabilität des Filters überprüft werden. Dies kann insbesondere bei IIR-Filtern zu Instabilität hoher Ordnung und Phasenverzerrung führen. Zeichnen Sie nach dem Entwurf des Filters den Frequenzgang (Bode-ähnlich) auf, um sicherzustellen, dass die Grenzfrequenz an der richtigen Stelle liegt und keine unerwünschten Wellen im Durchlassbereich auftreten.

Die typische Aufgabe bei der Rauschanalyse besteht darin, das reale Signal vom Rauschen zu trennen. Bewerten Sie hier den von der KI vorgeschlagenen Schwellenwert oder Filter anhand der physikalischen Bandbreite des Signals: Wenn Ihr reales Signal bei 0–500 Hz liegt, handelt es sich bei einer Komponente bei 2 kHz höchstwahrscheinlich um Rauschen oder Interferenzen. KI liefert Statistiken; Sie machen die physikalische Interpretation.

Validierung der Spektruminterpretation mit der physischen Realität

Bei der Interpretation eines Spektrums zieht die KI vernünftige, aber nicht immer genaue Schlussfolgerungen wie „Es gibt einen Peak bei 50 Hz, das ist eine Netzstörung.“ Vergleichen Sie diese Kommentare mit der physikalischen Erwartung: Verfügt die Anlage wirklich über einen 50-Hz-Netzanschluss? Sind die erwarteten Rotationsfrequenzen und Harmonischen im Schwingungsspektrum eines Motors annähernd vorhanden? Wenn möglich, wiederholen Sie die Messung mit einer zweiten Methode (anderer Sensor, andere Software) und sehen Sie die gleichen Peaks.

Hinweis: Für jede Spektralinterpretation: „Welchem ​​physikalischen Phänomen entspricht dieser Peak?“ Stellen Sie die Frage. 50/100/150-Hz-Netze und ihre Harmonischen, die Rotationsfrequenz und ihre Vielfachen in rotierenden Maschinen sowie die Schaltfrequenz in Schaltnetzteilen sind bekannte Signaturen. Ein unbekannter Peak ist entweder ein reales Ereignis oder ein Mess-/Verarbeitungsfehler; Recherchieren Sie beides.

Mini-Koffer

Ein F&E-Ingenieur misst die Vibration eines Lüfters und lässt die KI das Spektrum interpretieren. „Es gibt eine dominante Spitze bei 1200 Hz, das könnte ein Lagerschaden sein“, sagt AI. Der Techniker steuert die Drehzahl des Lüfters: 1800 U/min = 30 Hz. Es wird erwartet, dass die Lager- und Schaufelübergangsfrequenzen bei bestimmten Werten liegen, die mit diesem Zyklus kompatibel sind; 1200 Hz entsprechen keinem davon. Es steuert die Abtastfrequenz: Das Signal überträgt Inhalte bis zu 2 kHz, aber nur 2 kHz werden abgetastet, also eine Nyquist-Verletzung. Die Spitze bei 1200 Hz ist eine Geisterfrequenz, die durch Aliasing entsteht. Wenn ich die Abtastfrequenz auf 8 kHz erhöhe, verschwindet der falsche Peak. Lektion: Testen Sie die Spektruminterpretation immer mit physikalischen Zyklus-/Frequenzerwartungen und korrekter Abtastung.

Häufige Fehler

  • Auswahl der Abtastfrequenz ohne Überprüfung des Nyquist-Kriteriums; Aliasing übersehen.
  • Belassen der FFT-Frequenzachse beim Sample-Index statt bei Hz.
  • Spektrale Leckage wird als Signal ohne Fensterung verwechselt.
  • Die Amplitude wird nicht anhand der Anzahl der Fenster/Samples normalisiert.
  • IIR-Filterstabilität und Phasenverzerrung werden nicht überprüft.
  • Interpretieren von Spektrumsspitzen, ohne sie in Erwartung der physikalischen Drehzahl/Frequenz zu testen.

Zusammenfassend

  • Die Abtastfrequenz muss mindestens das Doppelte der höchsten Signalkomponente betragen; Andernfalls wird durch Aliasing die gesamte Analyse unterbrochen.
  • FFT-Auflösung Δf = f_s/N; Wählen Sie N und f_s entsprechend Ihrem Bedarf.
  • Fensterung reduziert spektrale Leckage; Normalisieren Sie die Amplitude entsprechend der Fensterverstärkung.
  • Überprüfen Sie den AI-Code mit einem bekannten Testsignal, um Achse und Normalisierung zu bestätigen.
  • Überprüfen Sie die Filterstabilität und den Frequenzgang, indem Sie sie grafisch darstellen.
  • Testen Sie die Spektruminterpretation immer unter Berücksichtigung der physikalischen Frequenz.

Anwendungsaufgabe

Erzeugen Sie ein Testsignal mit bekannter Frequenz (z. B. 500 Hz + 1500 Hz Sinus). Fordern Sie FFT- und Filtercode von AI an. Dann: (1) überprüfen, ob das Nyquist-Kriterium erfüllt ist, (2) überprüfen, ob die FFT-Peaks tatsächlich bei 500 und 1500 Hz auftreten, (3) absichtlich die Abtastfrequenz unter Nyquist reduzieren und beobachten, wie Aliasing einen falschen Peak erzeugt. Notieren Sie Ihre Beobachtungen und den von Ihnen korrigierten Parameter.