Das vorgestellte erklärbarkeitsbasierte KI-gestützte cyber-physische Kollaborationssystem wurde experimentell anhand von etwa 12.000 Audiodaten der MIMII-Datensammlung getestet, die sowohl normale als auch abnormale Maschinenzustände enthält. Der Trainingsdatensatz bestand aus ungefähr 9.600 Datensätzen, während der Testdatensatz etwa 2.400 Audiosignale umfasste. Basierend auf den Ergebnissen erreichte das vorgeschlagene CNN-Transformer-Modell eine Klassifikationsgenauigkeit von 98,5 % auf dem Testdatensatz, mit einem Bereich von 0,97–0,98. Dies bezieht sich auf die Angaben in Tabelle 3, die Konfusionsmatrix sowie die Analyse der statistischen Stabilität. Die Integration des Transformer-Encoders verbessert die zeitliche Modellierungsfähigkeit der Architektur, während der CNN diskriminierende räumliche Merkmale aus dem Eingangsspektrogramm extrahiert.
Neben der Klassifikationsleistung wurden auch andere systembezogene Eigenschaften bewertet. Erstens überstieg die Inferenzlatenz für eine Edge-basierte Implementierung nicht 20 ms, wodurch eine Echtzeiterkennung von Anomalien ermöglicht wurde. Zweitens trug das vorgestellte System durch die Reduzierung des Kommunikationsbedarfs zur Verbesserung der Energieeffizienz bei. Schließlich bietet der vorgeschlagene Erklärbarer-KI-Ansatz (Explainable AI) SHAP-Interpretierbarkeit, wodurch das Modell aus der Sicht eines fachkundigen Experten leichter analysierbar wird. Die MIMII-Daten wurden mithilfe der Kurzzeit-Fourier-Transformation (STFT) in Spektrogramme umgewandelt, wobei eine Fenstergröße von 1024, eine Hop-Größe von 512 und eine FFT-Größe (Fast Fourier Transform) von 1024 verwendet wurden. Die resultierenden logarithmischen Mel-Spektrogramme wurden auf 224 × 224 Pixel skaliert und anschließend mittels z-Score-Normalisierung vorverarbeitet. Zur Steigerung der Robustheit des trainierten Modells kamen Methoden der Datenvermehrung (Data Augmentation) zum Einsatz, darunter Zeitmaskierung, Frequenzmaskierung und Hinzufügen von Gaußschem Rauschen. Die Architektur des CNN bestand aus vier Faltungsschichten mit jeweils 32, 64, 128 und 256 Filtern sowie Batch-Normalisierung, ReLU- und Max-Pooling-Schichten. Die extrahierten Merkmalskarten wurden in Sequenzen umgewandelt und als Eingaben für den Transformer-Encoder verwendet, der aus vier Encodern, acht Attention-Heads, einem 256-dimensionalen Embedding und einer 1024-dimensionalen Feed-Forward-Dimension bestand. Eine Dropout-Wahrscheinlichkeit von 0,3 wurde angewandt, um Überanpassung (Overfitting) zu minimieren. Zum Training wurde der Adam-Optimierer mit einer Lernrate von 0,0001 und einer Gewichtsreduzierung (Weight Decay) von 1 × 10⁻5 verwendet. Die Anzahl der Epochen und die Batch-Größe beim Training betrugen jeweils 100 und 32. Der Zufallsstartwert (Random Seed) wurde auf 42 festgelegt, um die Reproduzierbarkeit des Experiments sicherzustellen. Zur Überprüfung der statistischen Gültigkeit der Ergebnisse wurden folgende Methoden angewandt: Alle Experimente wurden unabhängig 10-mal mit unterschiedlichen Zufallsstartwerten und Datenmischungen durchgeführt. Die Durchschnittswerte und Standardabweichungen von Genauigkeit (Accuracy), Präzision (Precision), Sensitivität (Recall), F1-Score, Fläche unter der ROC-Kurve (ROC-AUC) und Fläche unter der Präzisions-Recall-Kurve (PR-AUC) wurden berechnet. Außerdem wurden 95 % Konfidenzintervalle berechnet, um die Unsicherheit der Leistung abzuschätzen. Schließlich wurde ein gepaarter t-Test durchgeführt, um die Ergebnisse unseres CNN-Transformer-Ansatzes mit dem leistungsstärksten Baseline-Modell zu vergleichen. Die Methode des frühen Abbruchs (Early Stopping) wurde mit einer Geduld (Patience) von 10 Epochen angewandt. Tabelle 4 zeigt den Umgebungsaufbau und die Hyperparameter-Konfiguration, und Tabelle 5 enthält die Details der Simulationsumgebung.
Das Bewertungsexperiment in der vorliegenden Arbeit wurde mithilfe des MIMII-Datensatzes durchgeführt, der verschiedene Maschinenbetriebszustände und akustische Fehler-Szenarien umfasst. Obwohl die Ergebnisse des vorgeschlagenen Frameworks das Modell bestätigen, ist eine weitere Validierung anhand verschiedener Industriedatensätze und Fertigungsumgebungen erforderlich.
Der vorgeschlagene Rahmen wurde auf einer Edge-Computing-Plattform implementiert, die mit einem Multicore-Prozessor, 16 GB Systemspeicher und einem Grafikprozessor (GPU)-basierten Edge-Beschleuniger ausgestattet ist, um die Echtzeit-Inferenz und Erklärbarkeit zu unterstützen. Industrielle akustische Sensoren übertrugen Audiodatenströme über die IIoT-Kommunikationsschicht an das Edge-Gateway. Edge-Knoten wurden zur Durchführung von Aufgaben wie der Erstellung von Spektrogrammen, der Anomalie-Inferenz und der Generierung von SHAP-Erklärungen genutzt, während der Cloud-Server zur Speicherung langfristiger Daten und zum Aktualisieren des Modells verwendet wurde. Eine solche Konfiguration reduziert nicht nur den erforderlichen Kommunikationsaufwand, sondern ermöglicht auch die Echtzeit-Interaktion zwischen Sensoren, KI-Modulen und menschlichen Bedienern in cyber-physischen Systemen.
Der oben eingeführte CNN-Transformer-Ansatz wird mit mehreren bestehenden Methoden verglichen, einschließlich herkömmlicher ML-Modelle wie Support Vector Machine (SVM) und Random Forest, CNN-Modellen ohne Transformer sowie auf Long Short-Term Memory (LSTM) basierenden Verfahren. Während herkömmliche ML-Modelle aufgrund mangelnder Fähigkeit, räumliche und zeitliche Korrelationen in den Daten zu verarbeiten, keine zufriedenstellenden Ergebnisse erzielen, können CNN-Modelle räumliche Merkmale effizient verarbeiten, während LSTM-Modelle die zeitlichen Aspekte akustischer Signale bewältigen können. Dennoch sind erstere hinsichtlich Rechenaufwand und Fähigkeit zur parallelen Verarbeitung den letzteren unterlegen.
Um einen fairen Vergleich zwischen dem vorgeschlagenen CNN-Transformer-Modell und den Baseline-Modellen zu gewährleisten, wurde ersteres mit derselben Aufteilung des MIMII-Datensatzes (80 % Training und 20 % Test), denselben Vorverarbeitungsschritten, dem gleichen Verfahren zur Erstellung der Spektrogramme und denselben Evaluationsmaßen trainiert und bewertet. Bei der SVM-Methode kam der Radial-Basisfunktion (RBF)-Kernel mit C = 10 und γ = 0,01 zum Einsatz. Die Random-Forest-Methode umfasste 200 Bäume mit einer maximalen Baumtiefe von 20. Was das CNN betrifft, so enthielt es vier Faltungsschichten (32, 64, 128 und 256 Filter), denen volle Verbindungen ohne ein Transformer-Modul folgten. Schließlich wurden für die LSTM-Methode zwei gestapelte LSTM-Schichten mit 128 versteckten Einheiten und einer Dropout-Rate von 0,3 verwendet.
Im Gegenteil ermöglicht die Verwendung des CNN-Transformer-Frameworks die besten Ergebnisse, da sowohl die räumlichen als auch die zeitlichen Eigenschaften akustischer Signale effizient genutzt werden. Außerdem erlaubt der Einsatz von Edge-Computing eine Reduzierung der Rechenkosten, verbesserte Latenzzeiten und Energieeinsparungen. Darüber hinaus kann die Erklärbarkeit von Vorhersagen, die mithilfe von SHAP implementiert wird, als weiterer Vorteil angesehen werden, der die vorgeschlagene Methode von bestehenden Techniken unterscheidet.
Trotz jüngster Fortschritte bei der Entwicklung neuartiger Modelle für die akustische Anomalieerkennung, wie beispielsweise dem Vision Transformer (ViT), dem Audio Spectrogram Transformer (AST), dem Conformer und modellbasierten Ansätzen mit selbstüberwachtem Lernen, wurde im vorliegenden Studie kein experimenteller Vergleich dieser Architekturen durchgeführt. In Zukunft ist geplant, diese Modelle als Referenzwerte zur Bewertung des vorgeschlagenen Frameworks zu verwenden.
Metriken wie Genauigkeit, Präzision, Sensitivität (Recall) und F1-Score dienen zur Bewertung der Leistung des vorgeschlagenen Modells. Während die Präzision das Verhältnis der richtigen Positiven zur Gesamtanzahl der richtigen und falschen Positiven misst, bewertet die Genauigkeit die Gesamtkorrektheit der Vorhersagen. Die Sensitivität liefert eine Schätzung der Fähigkeit des Modells, Anomalien in einem bestimmten System vorherzusagen. Dies ist entscheidend, da das Versäumen einer Fehlerprognose zu Systemausfällen führen könnte. Die folgenden Formeln zeigen, wie sich diese Metriken mathematisch berechnen lassen:
(25)
(26)
(27)
(28)
Neben diesen werden zusätzliche Metriken bezüglich der Systemleistung, wie Latenz, Energieverbrauch und Skalierbarkeit, bei der Bewertung der realen Leistung des vorgeschlagenen CPS-Modells basierend auf Edge-Computing berücksichtigt.
Neben standardmäßigen Leistungsmessverfahren gibt es viele weitere Maßnahmen, die helfen können, die vorgeschlagene Lösung aus einer ganzheitlicheren Perspektive zu bewerten, insbesondere angesichts ihrer Anwendbarkeit auf reale, unausgeglichene und strombasierte Datenfallstudien. Beispielsweise wird die Fähigkeit des Algorithmus, zwischen den beiden Gruppen anhand unterschiedlicher Schwellenwerte zu unterscheiden, durch die Fläche unter der Receiver-Operating-Characteristic-Kurve (ROC-AUC) gemessen. Aufgrund seiner Robustheit und Trennfähigkeit sollte der ROC-AUC-Wert hoch sein. Außerdem ist das Precision-Recall-Flächen-unter-der-Kurve-Maß (PR-AUC) wichtig für die Anomalieerkennungs-Fallstudie, da die abnormalen Klassen selten sind.
Ein wichtiges Maß, das berücksichtigt werden sollte, ist der Matthews-Korrelationskoeffizient (MCC). Er berücksichtigt alle vier Elemente der Konfusionsmatrix und liefert auch bei Klassenungleichgewicht eine genaue Bewertung. Dieses Maß wird wie folgt berechnet:
(29)
Die Werte des MCC liegen zwischen -1 und +1, wobei +1 einer vollständig korrekten Vorhersage entspricht. Eine weitere anwendbare Metrik ist die Spezifität (True Negative Rate), die in einer industriellen Anwendung nützlich wäre, um zu beurteilen, ob falsch-positive Vorhersagen vorliegen:
(30)
Um sicherzustellen, dass keine Inkonsistenz zwischen der Konfusionsmatrix und den Leistungskennzahlen bestand, wurden alle Evaluierungskriterien basierend auf dem endgültigen Testdatensatz berechnet. In den folgenden Gleichungen werden durchgehend die Bezeichnungen für true positive (TP), true negative (TN), false positive (FP) und false negative (FN) verwendet. Genauigkeit (Accuracy), Präzision (Precision), Sensitivität (Recall), Spezifität (Specificity), F1-Score und MCC wurden gemäß den Gleichungen (25)–(30) ermittelt. Zusätzlich wurden ROC-AUC und PR-AUC anhand der Wahrscheinlichkeitsausgaben des CNN-Transformer-Modells durch unabhängige Schwellwertbewertung berechnet.
Zu den wichtigsten Kennzahlen auf Systemebene gehören Latenz, Durchsatz und Rechenaufwand. Die Latenz ist die Zeit, die für eine Vorhersage benötigt wird. Der Durchsatz gibt die Anzahl der pro Sekunde verarbeiteten Datensätze an. Der Rechenaufwand, der die Energieeffizienz widerspiegelt, ist entscheidend, um eine effiziente Berechnung sicherzustellen. Darüber hinaus können auch Erklärbarkeitskennzahlen wie SHAP-Konsistenz und Stabilität der Merkmalswichtigkeit bei der Bewertung des Modells berücksichtigt werden.
Aus den Ergebnissen in Tabelle 3 für den MIMII-Datensatz geht hervor, dass das vorgeschlagene CNN-Transformer-Modell alle anderen maschinellen Lern- und Deep-Learning-Verfahren übertrifft. Die Leistung der SVM- und Random-Forest-Modelle ist akzeptabel, jedoch begrenzt, mit Genauigkeiten von 88,6 % bzw. 91,2 %. Dies lässt sich durch die Unfähigkeit dieser Modelle erklären, komplexe zeitliche und räumliche Merkmale aus akustischen Signalen zu extrahieren. LSTM übertrifft hingegen SVM und Random Forest und erreicht eine Genauigkeit von 94,5 %, indem es Zeitreihendaten modelliert. Der vorgeschlagene CNN-Transformer-Ansatz weist jedoch die beste Leistung auf, mit verbesserter Genauigkeit (98,5 %), Präzision (98,06 %), Sensitivität (99,02 %) und F1-Score (98,54 %). Somit kann der vorgeschlagene Algorithmus normale und abnormale Maschinenzustände präzise klassifizieren. Die verbesserte Spezifität (97,96 %) zeigt die Wirksamkeit des Ansatzes bei der Reduzierung falsch positiver Ergebnisse, was für die Industrie von entscheidender Bedeutung ist. Es ist anzumerken, dass die Leistung des CNN-Transformers durch die Kombination aus CNN und Transformer erklärt werden kann, wobei ersteres diskriminierende Merkmale extrahiert und letzteres zeitliche Beziehungen in langen Sequenzen erfasst.
Die ROC-Kurve in Abbildung 6 zeigt die vergleichende Leistung von SVM, zufälligem Wald, LSTM und dem vorgeschlagenen CNN-Transformer-Ansatz hinsichtlich ihrer Klassifikationsfähigkeit bezüglich des MIMII-Datensatzes. Wie in Abbildung 6 dargestellt, wurde für den vorgeschlagenen CNN-Transformer die höchste Fläche unter der Kurve (AUC) mit einem Wert von 0,994 beobachtet. Offensichtlich bietet die vorgeschlagene Methode eine bessere diskriminative Leistung bei der Unterscheidung zwischen normalen und abnormalen Kategorien. Im Gegensatz dazu erzielt LSTM leicht niedrigere Ergebnisse als der vorgeschlagene CNN-Transformer-Ansatz, sodass seine AUC bei etwa 0,97 liegt. Sowohl der zufällige Wald als auch SVM weisen relativ niedrige AUC-Werte von etwa 0,958 bzw. 0,913 auf. Dies wird auf die Unfähigkeit beider Ansätze zurückgeführt, komplexe akustische Muster aus den vorliegenden Daten zu lernen.
Bei der Anomalieerkennung, bei der abnormale Proben selten sind, ermöglicht die Precision-Recall-(PR)-Kurve eine bessere Analyse. In der PR-Kurve erreicht das CNN-Transformer-Modell die höchste durchschnittliche Präzision (Average Precision, AP) von etwa 0,97–0,98, was auf eine überlegene Fähigkeit hinweist, abnormale Instanzen mit weniger Fehlalarmen zu erkennen, wie in Abbildung 7 dargestellt. Umgekehrt belegt der LSTM-Klassifikator den zweiten Platz mit einer AP zwischen 0,92 und 0,94. Der Zufallswald (Random Forest) und der SVM-Klassifikator belegen hingegen den dritten bzw. vierten Platz mit AP-Werten von 0,88 und 0,84. Der vorgeschlagene hybride Klassifikator erreicht über alle Recall-Stufen hinweg eine höhere Genauigkeit, was für praktische industrielle Anwendungen wichtig ist, da diese sowohl eine hohe Präzision als auch eine geringe Zahl an Fehlalarmen erfordern.
Das CNN-Transformer-Modell wurde anhand des MIMII-Datensatzes mit bestehenden Modellen anhand der in Tabelle 6 gezeigten erweiterten Leistungsmetriken verglichen. Der ROC-AUC-Wert zeigt die Diskriminationsfähigkeit des Modells zwischen den beiden Klassen über verschiedene Schwellenwerte hinweg an. Die SVM- und Random-Forest-Klassifikatoren erreichen ROC-AUC-Werte von 0,913 bzw. 0,958, während das LSTM-Modell einen ROC-AUC-Wert von 0,970 erreicht und das CNN einen Wert von 0,98 erzielt. Der vorgeschlagene CNN-Transformer-Ansatz erreicht den höchsten ROC-AUC-Wert von 0,994 und demonstriert damit eine überlegene Unterscheidung zwischen normalen und abnormalen Maschinenzuständen. Ebenso zeigen die PR-AUC-Werte, dass das vorgeschlagene Modell mit Klassenungleichgewicht effektiv umgehen kann, was ein entscheidender Faktor bei der Anomalieerkennung ist. Der vorgeschlagene CNN-Transformer erreicht den höchsten PR-AUC-Wert (durchschnittliche Präzision) von 0,982, gefolgt von CNN (0,950), LSTM (0,912), Random Forest (0,891) und SVM (0,765), was eine überlegene Präzisions- und Rückrufleistung für die industrielle akustische Anomalieerkennung anzeigt. Darüber hinaus zeigt der Matthews-Korrelationskoeffizient (MCC), der das Gleichgewicht der Modellleistung misst, für das vorgeschlagene Modell einen relativ hohen Wert von 0,97 an. Frühere Methoden wie SVM (0,73) und Random Forest (0,78) wiesen hingegen eine geringere Vorhersagegenauigkeit auf.
Um die Stabilität des vorgestellten Ansatzes zu überprüfen, wurde dasselbe Experiment 10-mal mit unterschiedlichen Initialisierungsseeds und Datenmischungsverfahren durchgeführt. Die Ergebnisse in Tabelle 7 für die CNN-Transformer-Architektur waren wie folgt: Genauigkeit = 98,50 % ± 0,19 %, Präzision = 98,06 % ± 0,23 %, Sensitivität = 99,02 % ± 0,22 % und F1-Score = 98,54 % ± 0,24 %. Die relativ geringe Standardabweichung weist auf die Stabilität des Modells hin. Die 95-%-Konfidenzintervalle wurden für alle Evaluierungsmaße berechnet. Es zeigte sich, dass das Konfidenzintervall für die Genauigkeit [98,1 %, 98,9 %] betrug, was darauf hindeutet, dass das Modell konsistent gut abschneidet. Auch für ROC-AUC, PR-AUC und MCC ergaben sich enge Konfidenzintervalle. Der gepaarte t-Test zum Vergleich des vorgeschlagenen CNN-Transformer-Modells mit dem leistungsstärksten Baseline-Modell (LSTM) ergab einen p-Wert < 0,05, was darauf hinweist, dass der beobachtete Leistungszuwachs signifikant ist und nicht zufällig zustande kam.
Um die Robustheit der vorgeschlagenen CNN-Transformer-Architektur zu überprüfen, wurden die Experimente 10-mal mit unterschiedlichen Zufallsinitialisierungs-Seeds und Mischkonfigurationen wiederholt, wie in Abbildung 8 dargestellt. Die durchschnittlichen Werte für Genauigkeit, Präzision, Sensitivität und F1-Score, die vom Modell erzielt wurden, betrugen jeweils 98,50 % ± 0,19 %, 98,06 % ± 0,23 %, 99,02 % ± 0,22 % und 98,54 % ± 0,24 %. Aus der Boxplot-Analyse geht hervor, dass die Leistungskennzahlen des Modells nur geringfügige Schwankungen aufweisen, was auf Stabilität und Robustheit hinweist. Folglich zeigt der vorgeschlagene Ansatz über mehrere Experimente hinweg eine hohe Stabilität ohne signifikante Variationen, was auf gute Generalisierungsfähigkeiten des Modells schließen lässt.
Zusätzlich wurde ein statistischer Test basierend auf den Ergebnissen mehrerer Durchläufe des Experiments durchgeführt. Die geringen Standardabweichungen bei allen Evaluierungsmaßen deuten auf eine geringe Anfälligkeit gegenüber Zufallsvariationen oder Trainingseinflüssen hin. Dies bestätigt, dass das vorgeschlagene CNN-Transformer-Modell stabile und konsistente Leistungen in realen cyber-physikalischen Fertigungssystemen erbringen kann.
Die SHAP-Heatmap vermittelt einen Eindruck von der Bedeutung der Zeit-Frequenz-Bereiche im Spektrogramm für die Vorhersage, ob sich der Maschinenzustand im Normal- oder im Fehlerzustand befindet. SHAP-Werte zeigen den Beitrag jedes Prädiktors zur endgültigen Vorhersage an.
In der SHAP-Wert-Heatmap in Abbildung 9A stellen positive SHAP-Werte (z. B. zwischen +0,6 und +1,2) einen Beitrag zur Abnormalität dar, während negative SHAP-Werte (d. h. zwischen −0,3 und −0,8) einen Beitrag zur Normalitätskategorie anzeigen. Der intensivste SHAP-Wert wurde im mittleren Frequenzbereich (40–80 Frequenzbänder) in den Zeitfenstern 50–100 gefunden. Diese Bereiche sind die entscheidenden Regionen, in denen kritische Fehler signale in akustischen Messungen erfasst werden. Außerdem zeigt das niederfrequente Band (weniger als 20 Bänder) SHAP-Werte nahezu null (−0,1 bis +0,1), was auf eine vernachlässigbare Auswirkung auf die Entscheidungsfindung des Modells hinweist. Dies zeigt, dass das Modell Hintergrundgeräusche, die für die Fehlervorhersage irrelevant sind, ignoriert. Die zeitliche Konsistenz hoher SHAP-Werte innerhalb bestimmter Zeitbänder unterstreicht zudem die Fähigkeit des Transformer-Moduls, langreichweitige Abhängigkeiten in akustischen Signalen zu erfassen. Zusammenfassend macht die SHAP-Heatmap deutlich, dass das CNN-Transformer-Modell auf spezifische Zeit-Frequenz-Bänder mit hoher Diskriminierungskraft fokussiert. Dies verbessert die Interpretierbarkeit des Modells und hilft dem menschlichen Bediener, Fehler visuell zu erkennen.
Es wird gezeigt, dass die SHAP-Methode ausschließlich als Interpretierbarkeitstechnik eingesetzt wird und nicht Teil des Trainingsprozesses des CNN-Transformer-Modells ist. Die Klassifizierungsgenauigkeit wird ausschließlich durch die Parameter bestimmt, die von den CNN- und Transformer-Komponenten gelernt wurden. SHAP wird nach dem Modell-Inferenzschritt verwendet, um die Vorhersagen zu interpretieren, indem die einflussreichsten Zeit-Frequenz-Fenster identifiziert werden, die die Vorhersagen in Richtung normaler oder abnormer Kategorien lenken. Dadurch verbessert SHAP die Transparenz und das Verständnis für den Menschen, ohne die Klassifizierungsgenauigkeit zu beeinflussen. Die vom vorgeschlagenen Framework bereitgestellten SHAP-Werte ermöglichen eine klare Visualisierung, die Wartungsoperatoren dabei unterstützt, Modellvorhersagen zu überprüfen, indem sie erkannte Bereiche mit dem tatsächlichen Maschinenverhalten und Wartungsdatensätzen vergleichen. Somit wird die Zusammenarbeit zwischen Mensch und Maschine bei vorausschauenden Wartungsmaßnahmen erleichtert.
Zur Bewertung des Erklärbarkeitsgrads wurde eine globale SHAP-Analyse der Merkmalswichtigkeit durchgeführt, wie in Abbildung 9B veranschaulicht. Die Ergebnisse zeigen, dass einige mittlere Frequenzkomponenten eine bedeutende Rolle bei der Anomalieerkennung spielen. Darüber hinaus zeigt eine in Abbildung 9C dargestellte Fallstudie, dass SHAP fehlerrelevante Bereiche in abnormalen Maschinenspektrogrammen korrekt identifiziert. Die Studie umfasste eine qualitative Fallstudie unter Verwendung abnormaler Datensätze von Pumpen und Ventilen aus dem MIMII-Datensatz. Bei defekten Pumpen identifizierte SHAP Frequenzbereiche von 2–4 kHz, die mit Kavitation und Verschleiß assoziiert sind. Im Fall defekter Ventile war die SHAP-Aktivierung in Bereichen mit wiederholten harmonischen Signalen, die auf Leckagen hindeuten, am stärksten ausgeprägt. Bei repräsentativen Testbeispielen identifizierten die SHAP-Visualisierungen zuverlässig diskriminierende Zeit-Frequenz-Zonen, die mit anomalen Maschinenzuständen verbunden sind. Indem sie die Spektralbereiche hervorheben, die am stärksten zu jeder Vorhersage beitragen, liefern diese Erklärungen Einblicke in die Entscheidungsfindung des Modells. Die SHAP-Analyse wird als Interpretationswerkzeug zum Verständnis des Modellverhaltens angeboten, nicht jedoch als durch Experten validierte Fehlerverifizierung, da eine formale Validierung durch Wartungsfachleute außerhalb des Untersuchungsumfangs lag.
Die Konfusionsmatrix verdeutlicht, wie gut das vorgeschlagene CNN-Transformer-Modell die Klassifizierung des in Abbildung 10 gezeigten MIMII-Datensatzes durchführt. Insgesamt wurden 960 normale Proben korrekt als normal identifiziert (True Negatives), während 20 normale Proben fälschlicherweise als abnormal eingestuft wurden (False Positives). Darüber hinaus wurden von allen Proben insgesamt 1010 abnormale Proben korrekt erkannt (True Positives), während 10 abnormale Proben als normal klassifiziert wurden (False Negatives).
Daher wurde eine hohe Erkennungsleistung erzielt, insbesondere bei anomalen Proben, die für die vorausschauende Wartung entscheidend sind. Wie in Abbildung 10 gezeigt, tritt ein Fehlalarm vom Typ „False Negative“ äußerst selten auf (10 Proben) – dies ist für jede Industrie von zentraler Bedeutung, um Sicherheit und Zuverlässigkeit sicherzustellen. Umgekehrt kommt der Fehler vom Typ „False Positive“ (20 Proben) etwas häufiger vor, bleibt jedoch auf einem relativ niedrigen Niveau. Die in Abbildung 10 dargestellte Konfusionsmatrix wurde anhand einer ausgewogenen Teilmenge von 2.000 Testaufnahmen (mit gleicher Anzahl normaler und anormaler Proben) erstellt, um eine klare Visualisierung der Klassifikatorleistung zu ermöglichen. Alle angegebenen Evaluierungsmetriken, einschließlich Genauigkeit, Präzision, Sensitivität, Spezifität und F1-Score, wurden anhand des vollständigen Testdatensatzes (etwa 2.400 Aufnahmen) berechnet.
Unter Berücksichtigung der in Abbildung 10 dargestellten Konfusionsmatrix und von Tabelle 8 (TN = 960, FP = 20, TP = 1010, FN = 10) lassen sich folgende Berechnungen durchführen: Genauigkeit = (TP + TN) / (TP + TN + FP + FN) = 98,50 %, Präzision = TP/(TP + FP) = 98,06 %, Sensitivität = TP/(TP + FN) = 99,02 %, Spezifität = TN / (TN + FP) = 97,96 %, F1-Score = 98,54 % und MCC = 0,97. Außerdem betrugen die ROC-AUC- und PR-AUC-Werte 0,994 bzw. 0,982. Die in Abbildung 10 gezeigte Konfusionsmatrix wurde mithilfe einer ausgewogenen, repräsentativen Teilmenge von 2.000 Aufnahmen aus den Testdaten lediglich zu Visualisierungszwecken erstellt. Alle in dieser Studie angegebenen quantitativen Leistungsmetriken wurden anhand des vollständigen Testdatensatzes (etwa 2.400 Aufnahmen) berechnet. Insgesamt zeigen die Ergebnisse, dass die vorgeschlagene CNN-Transformer-Architektur eine Klassifizierungsgenauigkeit von 98,5 % erreicht, was mit den in Tabelle 8 angegebenen Ergebnissen und der Ablassanalyse übereinstimmt.
Zur Überprüfung der Effektivität der kantenbasierten Implementierung wurden zusätzliche Messungen durchgeführt. Die durchschnittliche Inferenzlatenz des vorgeschlagenen Modells beträgt 18,7 ms/Probe, wodurch etwa 53 Proben pro Sekunde verarbeitet werden können, wie in Tabelle 9 gezeigt. Der Energieverbrauch während der Inferenz betrug 8,9 W und führte zu einem geschätzten Energieverbrauch von 0,167 J/Vorhersage. Zur Bewertung der Nachhaltigkeitsauswirkungen des vorgeschlagenen Ansatzes in Tabelle 10 wurden Kommunikationsaufwand, Energieverbrauch und Ressourcennutzung berücksichtigt. Da der Inferenzprozess lokal auf Edge-Knoten erfolgt, werden nur zusammengefasste Diagnosedaten an die Cloud-Umgebung übertragen. Die experimentellen Ergebnisse zeigten einen Rückgang des Kommunikationsaufwands um 73 % im Vergleich zu herkömmlichen, cloud-basierten CPS-Systemen. Außerdem verringerte die lokale Inferenz den Energieverbrauch um 32 %, während sich die Kommunikationszeit von 75 ms auf 20 ms reduzierte.
Abtragungsstudie
Der MIMII-Datensatz wird einer Ablationsforschung unterzogen, um die Auswirkungen jeder Komponente der auf erklärbarem KI-basierten CPHS-Architektur zu bewerten. Diese Ablationsstudie untersucht die Effekte der Verwendung einer CNN zur Erfassung räumlicher Merkmale, eines Transformer-Encoders zur Erfassung zeitlicher Merkmale sowie des SHAP-Erklärbarkeitsalgorithmus selbst. Mehrere Modifikationen des vorgeschlagenen Modells werden betrachtet, und zwar: (i) eine nur CNN-basierte Version, (ii) eine nur Transformer-basierte Version, (iii) das Modell mit kombinierter CNN- und Transformer-Architektur, jedoch ohne SHAP-Algorithmus, und (iv) das vollständige vorgeschlagene Modell mit kombinierter CNN- und Transformer-Architektur sowie SHAP. Zur Bewertung werden gängige Metriken wie Genauigkeit, Präzision, Sensitivität und F1-Score verwendet.
Diese Ergebnisse der Ablationsstudie verdeutlichen die Bedeutung jedes Elements im Rahmenwerk, das in Tabelle 11 dargestellt ist. Das ausschließlich auf CNN basierende Modell erreicht eine Genauigkeit von 93,8 %, was zeigt, dass die Extraktion räumlicher Merkmale aus Spektrogrammen wirksam ist, jedoch keine zeitlichen Abhängigkeiten in akustischen Signalen erfassen kann. Ähnliches gilt für das ausschließlich auf dem Transformer basierende Modell, das eine leicht niedrigere Genauigkeit (92,6 %) erreicht, da es zwar zeitliche Merkmale extrahiert, dabei aber räumliche Informationen vernachlässigt. Schließlich erzielte das vorgeschlagene CNN-Transformer-Modell eine Genauigkeit von 98,5 %, eine Präzision von 98,06 %, einen Recall von 99,02 % und einen F1-Score von 98,54 %, was mit den Ergebnissen des Hauptexperiments übereinstimmt. Da SHAP nur für die nachträgliche Interpretation und nicht für das Modelltraining verwendet wird, entspricht die Klassifikationsgenauigkeit der des trainierten CNN-Transformer-Modells. Anschließend wird SHAP verwendet, um Erklärungen zur Merkmalzugehörigkeit zu erstellen. Der verbesserte Recall (99,02 %) zeigt, dass das Modell effektiv dabei ist, abnormale Zustände in den Maschinen zu erkennen, wodurch sichergestellt wird, dass keine fehlerhaften Bauteile übersehen werden – eine entscheidende Voraussetzung für die Implementierung eines Systems zur vorausschauenden Wartung.
VERFÜGBARKEIT VON DATEN:
Der in dieser Studie verwendete Datensatz „Malfunctioning Industrial Machine Investigation and Inspection“ (MIMII) ist öffentlich im offiziellen Zenodo-Repository verfügbar. Die Experimente wurden mithilfe der öffentlich zugänglichen Audiodaten und Annotationen durchgeführt, die von den Autoren des Datensatzes bereitgestellt wurden. Auf den Datensatz kann unter https://zenodo.org/records/3384388 zugegriffen werden. Der Implementierungscode, der dieser Studie beigefügt ist, steht öffentlich über das Zenodo-Repository unter https://zenodo.org/records/21405292 zur Verfügung. Das Repository enthält den Implementierungscode für den vorgeschlagenen Ansatz, einschließlich der Daten-Preprocessing-Pipeline, der Modellarchitektur, des Trainingsablaufs, der Evaluierungsskripte, Konfigurationsdateien und unterstützender Hilfsprogramme, um das Verständnis und die unabhängige Reproduktion der vorgeschlagenen Methodik zu erleichtern. Der MIMII-Datensatz wird nicht zusammen mit dem Implementierungscode weiterverbreitet und muss separat aus dem offiziellen Repository bezogen werden.

Abbildung 1: Aufbau der vorgeschlagenen Arbeit. Architektonisches Design eines erklärbarkeitsbasierten, KI-gesteuerten cyber-physischen Mensch-System, das die Vorverarbeitung, räumliches Lernen auf Basis von CNNs, zeitliche Modellierung auf Basis von Transformatoren, SHAP-basierte Erklärbarkeit und Edge-Computing für eine nachhaltige Produktion kombiniert. Klicken Sie bitte hier, um eine vergrößerte Darstellung dieser Abbildung anzusehen.

Abbildung 2: Verarbeitungspipeline für die Vorverarbeitung. Die Vorverarbeitungspipeline für akustische Signale umfasst Rauschunterdrückung, Normalisierung, Segmentierung, STFT-Transformation und die Extraktion des Spektrogramms. Klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.

Abbildung 3: Räumliche Merkmalsextraktion mittels CNN. Architektur zur räumlichen Merkmalsextraktion durch Faltungsschichten (CNNs) unter Verwendung von Spektrogrammen über die Prozesse der Faltung, Aktivierung und des Poolings. Klicken Sie hier, um eine vergrößerte Ansicht dieser Abbildung anzuzeigen.

Abbildung 4: Zeitliche Modellierung mithilfe eines Transformer-Encoders. Architektur eines Transformer-Encoders, um zeitliche Abhängigkeiten und langreichweitige Abhängigkeiten durch die Repräsentation einer Spektrogrammsequenz zu erfassen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 5: Klassifizierungsschicht. Eine Klassifizierungsschicht, die die kodierten zeitlichen Merkmale auf Wahrscheinlichkeitswerte abbildet, um den Maschinenzustand vorherzusagen. Klicken Sie bitte hier, um eine vergrößerte Version dieser Abbildung anzusehen.

Abbildung 6: Ergebnis der ROC-Kurve. Vergleich der ROC-Kurven für SVM-, Random-Forest-, LSTM-, CNN- und vorgeschlagene CNN-Transformer-Algorithmen auf dem MIMII-Datensatz. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 7: Ergebnis der Präzisions-Recall-Kurve. Vergleich der Präzisions-Recall-Kurven, der die Effizienz der Anomalieerkennung verschiedener Methoden des maschinellen Lernens und des Deep Learnings darstellt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 8: Analyse der Stabilität der Leistung bei mehrfacher Durchführung. Boxplot, der die Stabilität der Leistung hinsichtlich Genauigkeit, Präzision, Rückruf und F1-Score-Ergebnisse zeigt, die aus zehn unabhängigen Experimenten mit der vorgeschlagenen CNN-Transformer-Architektur erzielt wurden. Die geringen Schwankungen deuten auf Stabilität und Robustheit des Modells hin. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 9: SHAP-basierte Erklärbarkeitsanalyse des vorgeschlagenen Anomalieerkennungs-Frameworks. (A) SHAP-basierte Erklärbarkeitsanalyse. Heatmap, die die SHAP-Methode verwendet, um die relevantesten Zeit-Frequenz-Bereiche anzuzeigen, die zur Identifizierung anomalen Verhaltens führen. (B) Analysen der Merkmalswichtigkeit mithilfe des globalen SHAP-Ansatzes zur Identifizierung der einflussreichsten Merkmale. (C) Visualisierung der SHAP-Erklärbarkeit für eine fehlerhafte Maschine. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 10: Konfusionsmatrix für den MIMII-Datensatz. Die Konfusionsmatrix zeigt die Klassifikationsleistung des vorgeschlagenen CNN-Transformator-Algorithmus. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
| Parameter | Beschreibung |
| Datensatzname | MIMII (Malfunctioning Industrial Machine Investigation and Inspection) Datensatz |
| Typen von Maschinen | Armaturen, Pumpen, Gebläse, Schiebeschienen |
| Gesamtanzahl verwendeter Stichproben | Ungefähr 12.000 Audiodateien |
| Klassen | Normal, Anomal |
| Abtastrate | 16 kHz |
| Dauer pro Stichprobe | 10 Sekunden |
| Audioformat | WAV |
| Merkmalsdarstellung | Spektrogrammbilder auf STFT-Basis |
| Vorverarbeitungsschritte | Rauschfilterung, Normalisierung, Segmentierung, STFT-Transformation |
| Kennzeichnung der Klassen | Normale Aufnahmen sind als Klasse 0 gekennzeichnet; anormale Aufnahmen sind als Klasse 1 gekennzeichnet gemäß MIMII-Anmerkungen |
| Trainings-Test-Aufteilung | 80 % Training (≈9.600 Stichproben), 20 % Testen (≈2.400 Stichproben) |
| Strategie zur Datenaufteilung | Aufteilung auf Ebene der Audiodateien vor der Segmentierung, um Datenlecks zu vermeiden |
| Industrielle Schallbedingungen | Realistische industrielle akustische Umgebungen, enthalten in den MIMII-Aufnahmen |
| Anwendungsgebiet | Vorbeugende Wartung und Erkennung von Maschinenanomalien |
| Ziel der Industrie 5.0 | Erklärbare, menschenzentrierte und nachhaltige Fehlererkennung |
Tabelle 1: Datensatzbeschreibung. Beschreibung des MIMII-Datensatzes hinsichtlich Maschinentyp, Stichprobenverteilung, Signalcharakteristika, Merkmalsdarstellung und Anwendungsaspekten.
| Maschinentyp | Maschinen-IDs | Normale Aufnahmen | Anomale Aufnahmen | SNR-Bedingungen (dB) | Experimentelle Aufteilung |
| Lüfter | ID00–ID06 | 8.400 | 1.600 | −12 bis −9 | 80 % Training / 20 % Test |
| Getriebe | ID00–ID06 | 7.124 | 1.147 | −17 bis −15 | 80 % Training / 20 % Test |
| Pumpe | ID00–ID06 | 7.200 | 1.800 | −18 bis −9 | 80 % Training / 20 % Test |
| Schiebeschiene | ID00–ID05 | 7.000 | 1.800 | −14 bis −12 | 80 % Training / 20 % Test |
| Ventil | ID00–ID05 | 7.000 | 1.800 | −12 bis −9 | 80 % Training / 20 % Test |
Tabelle 2: Beschreibung des im vorliegenden Studium verwendeten MIMII-Datensatzes. Industrielle Maschinentypen, Maschinen-IDs, Anzahl der normalen und anomalen Aufnahmen, SNR-Bedingungen sowie die zur überwachten binären Klassifizierung verwendete Aufteilung in Trainings- und Testdatensatz.
| Modell | Genauigkeit (%) | Präzision (%) | Erkennungsrate (%) | F1-Wert (%) | Spezifität (%) |
| SVM | 88,6 | 86,9 | 85,4 | 86,1 | 90,2 |
| Random Forest (RF) | 91,2 | 89,8 | 88,5 | 89,1 | 92,6 |
| LSTM | 94,5 | 93,2 | 92,8 | 93 | 95,1 |
| CNN | 96,3 | 95,4 | 94,9 | 95,1 | 96,8 |
| CNN–Transformer (vorgeschlagen) | 98,5 | 98,06 | 99,02 | 98,54 | 97,96 |
Tabelle 3: Klassifizierungsergebnis für den MIMII-Datensatz. Vergleich der Leistung von maschinellen und tiefenlernbasierten Klassifikatoren anhand von Genauigkeit, Präzision, Sensitivität, F1-Score und Spezifität.
| Parameter | Wert |
| Datensatz | MIMII Dataset |
| Gesamtanzahl der Stichproben | ~12.000 |
| Aufteilung Trainings- und Testdaten | 80 % / 20 % |
| Eingabetyp | Spektrogrammbilder |
| CNN-Schichten | 3–5 Faltungsschichten |
| Transformer-Schichten | 2–4 Encoder-Schichten |
| Batch-Größe | 32 |
| Lernrate | 0,001 |
| Optimierer | Adam |
| Epochen | 100 |
| Hardware | Edge-Gerät + GPU (zum Training) |
| STFT-Fensterlänge | 1024 |
| Hubgröße | 512 |
| FFT-Größe | 1024 |
| Auflösung des Spektrogramms | 224 × 224 |
| Normalisierung | Z-Score |
| Datenaugmentierung | Zeitmaskierung, Frequenzmaskierung, Gaußsches Rauschen |
| CNN-Schichten | 4 |
| CNN-Filter | 32, 64, 128, 256 |
| Transformer-Schichten | 4 |
| Attention-Heads | 8 |
| Einbettungsdimension | 256 |
| Feed-Forward-Dimension | 1024 |
| Dropout | 0,3 |
| Zufallsstartwert | 42 |
| Geduld bei frühem Abbruch | 10 |
Tabelle 4: Umgebungssetup und Hyperparameter-Konfiguration. Parameter des Trainings und der Experimente, die zur Implementierung der vorgeschlagenen CNN-Transformer-Architektur verwendet wurden.
| Komponente | Spezifikation |
| Programmiersprache | Python |
| Framework | TensorFlow / PyTorch |
| Prozessor | Intel i7 / Ryzen 7 |
| GPU | NVIDIA GTX 1660 / RTX Series |
| Edge-Gerät | Raspberry Pi / NVIDIA Jetson |
| Visualisierung | Matplotlib, SHAP |
| Edge-Gerät | NVIDIA Jetson Xavier NX |
| CPU | 6-Kern-ARM v8.2 |
| Arbeitsspeicher (RAM) | 16 GB |
| Betriebssystem | Ubuntu 20.04 |
| Deep-Learning-Framework | TensorFlow/PyTorch |
| Datensatzgröße | 12.000 MIMII-Proben |
| Konnektivität | Ethernet/Wi-Fi |
| Modellparameter | 8,4 Millionen |
| Modellgröße | 32,7 MB |
| Framework | PyTorch + TensorRT |
| Batch-Größe | 1 |
| Quantisierung | FP16 |
| Pruning | Nein |
| Durchschnittliche Latenz | 17,8 ms |
| Median-Latenz | 17,2 ms |
| Latenz im 95. Perzentil | 19,6 ms |
| Throughput | 56 Proben/s |
| Speicherauslastung | 1,4 GB |
| Leistungsaufnahme | 11,3 W |
Tabelle 5: Simulationsumgebung. Hardware und Software wurden verwendet, um das vorgeschlagene Modell zu trainieren, bereitzustellen und zu evaluieren.
| Modell | ROC-AUC | PR-AUC | MCC |
| SVM | 0.913 | 0.765 | 0.73 |
| Zufallswald (RF) | 0.958 | 0.891 | 0.78 |
| LSTM | 0.97 | 0.912 | 0.86 |
| CNN | 0.98 | 0.95 | 0.9 |
| CNN–Transformer (vorgeschlagen) | 0.994 | 0.982 | 0.97 |
Tabelle 6: Vergleich der experimentellen Ergebnisse für ROC-AUC, PR-AUC und MCC. Vergleich zwischen verschiedenen Modellen basierend auf den Maßen ROC-AUC, PR-AUC und Matthews-Korrelationskoeffizient.
| Metrik | Mittelwert (%) + Std (%) | 95 % Konfidenz |
| Genauigkeit | 98,50 ± 0,19 | [98,1, 98,9] |
| Präzision | 98,06 ± 0,23 | [98,0, 98,2] |
| Recall | 99,22 ± 0,22 | [98,8, 99,4] |
| F1-Score | 98,54 ± 0,24 | [98,1, 98,9] |
| ROC-AUC | 0,984 | 0,004 |
| PR-AUC | 0,982 | 0,005 |
| MCC | 0,97 | 0,007 |
Tabelle 7: Ergebnis der statistischen Analyse bei mehrfacher Durchführung. Robustheitsstatistiken für die vorgeschlagene CNN-Transformer-Architektur über 10 Experimente hinweg, einschließlich Mittelwert, Standardabweichung, 95 %-Konfidenzintervall und Signifikanz für verschiedene Leistungskennzahlen.
| Metrik | Wert |
| Genauigkeit | 98,50 % |
| Präzision | 98,06 % |
| Erkennungsrate (Recall) | 99,02 % |
| Spezifität | 97,96 % |
| F1-Score | 98,54 % |
| MCC | 0,97 |
| ROC-AUC | 0,994 |
| PR-AUC | 0,982 |
Tabelle 8: Endgültiger Testdatensatz der Konfusionsmatrix. Die Konfusionsmatrix des vorgeschlagenen Architekturansatzes für den Testdatensatz, die die Klassifizierung normaler und abnormer Zustände der Maschinen zeigt, auf deren Grundlage die Evaluierungsmetriken berechnet wurden.
| Metrik | Cloud-Deployment | Edge-Deployment |
| Inferenz-Latenz (ms) | 85,6 | 18,7 |
| Durchsatz (Stichproben/s) | 22 | 53 |
| Netzwerkauslastung (MB/min) | 120 | 18 |
| Energie pro Vorhersage (J) | 0,52 | 0,17 |
| Echtzeitfähigkeit | Mäßig | Hoch |
Tabelle 9: Ergebnis der Inferenz-Latenz. Die Latenzleistung des Edge-Deployments des erklärbaren KI-fähigen CPHS-Frameworks, einschließlich Verarbeitungszeit, Durchsatz, Speicherverbrauch und weiterer Echtzeiteigenschaften.
| Metrik | Cloud-basiertes CPS | Vorgeschlagenes Edge-fähiges CPHS | Verbesserung |
| Datenübertragung pro Stunde | 100% | 35% | 65% Reduktion |
| Kommunikationslatenz | 75 ms | 20 ms | 73,3% Reduktion |
| Energieverbrauch | 100% | 68% | 32% Reduktion |
| Geschätzte Kohlenstoffemission | 100% | 70% | 30% Reduktion |
Tabelle 10: Nachhaltigkeitsbewertung des Edge-fähigen CPHS. Messungen des Kommunikationsaufwands, des Energieverbrauchs, der Ressourcennutzung und der Latenz werden zur Bewertung der Nachhaltigkeit des Edge-fähigen CPHS herangezogen.
| Modellvariante | Genauigkeit (%) | Präzision (%) | Erinnerung (%) | F1-Wert (%) |
| Nur CNN | 93,8 | 92,5 | 94,2 | 93,3 |
| Nur Transformer | 92,6 | 91,2 | 93,5 | 92,3 |
| CNN + Transformer | 98,5 | 98,06 | 99,02 | 98,54 |
Tabelle 11: Ergebnis der Ablationsstudie. Ergebnis der Ablationsstudie, das den Einfluss der CNN-, Transformer- und SHAP-Module auf die Leistung des vorgeschlagenen Modells verdeutlicht.