Sugerowana metodologia wykorzystuje systemowe podejście do cyberfizycznych systemów ludzkich (CPHS) w ramach Przemysłu 5.0, które łączy wyjaśnialne głębokie uczenie z ramami współpracy cyberfizycznej opartymi na przetwarzaniu brzegowym (edge computing) w celu wspierania zrównoważonych procesów produkcyjnych. W pierwszej kolejności problem zostaje zdefiniowany poprzez wskazanie pewnych ograniczeń obecnych rozwiązań CPHS, takich jak brak inteligencji opartej na AI, brak wyjaśnialności oraz brak współpracy skoncentrowanej na człowieku. Aby sprostać tym wyzwaniom, proponowane rozwiązanie wykorzystuje zbiór danych MIMII do monitorowania stanu maszyn przemysłowych przy użyciu przetwarzania sygnałów akustycznych.
Zaproponowany schemat rozpoczyna się od wprowadzenia wygenerowanych spektrogramów do modelu CNN, który wyodrębnia cechy przestrzenne poprzez rozpoznawanie istotnych wzorców częstotliwości powiązanych z różnymi stanami maszyn. Wyekstrahowane cechy przestrzenne zostaną przekształcone w sekwencje i wprowadzone do modułu kodującego transformera w celu uchwycenia zależności czasowych za pomocą mechanizmu self-attention. Przekształcona reprezentacja cech jest następnie przekazywana do w pełni połączonego modułu klasyfikatora w celu zaklasyfikowania stanów maszyn jako normalne lub nieprawidłowe. Ponadto, aby zapewnić przejrzystość i ułatwić podejmowanie decyzji zorientowanych na człowieka, zastosowano podejście do wyjaśnialności SHAP, aby wyjaśnić przewidywania poprzez wybór istotnych obszarów czasowo-częstotliwościowych w wyjściu modelu. Zaprojektowany model zostanie zaimplementowany w środowisku CPS z obsługą brzegu sieci (edge-enabled), co umożliwi wnioskowanie w czasie rzeczywistym przy niskim opóźnieniu obliczeniowym i niskim zużyciu energii. Dodatkowo zintegrowany zostanie schemat human-in-the-loop, aby umożliwić operatorom weryfikację wyników przy użyciu rezultatów wyjaśnialności. Wydajność jest walidowana za pomocą dokładności (accuracy), precyzji (precision), pełności (recall) i miary F1-score, a także innych metryk na poziomie systemu, takich jak opóźnienie i efektywność energetyczna.
Zaproponowany schemat stanowi mechanizm CPHS oparty na XAI dla zrównoważonej produkcji w kontekście Przemysłu 5.0, a jego przebieg przedstawiono jako potok procesowy – od pozyskiwania danych po inteligentne podejmowanie decyzji i ich realizację, co pokazano na Rysunku 1. W początkowej fazie procedury pozyskiwania danych sygnały dźwiękowe są zbierane z maszyn przemysłowych, takich jak silniki, pompy i łożyska, z wykorzystaniem zbiorów danych, w tym MIMII. Uzyskane w ten sposób sygnały danych są przekazywane do kolejnego etapu, zwanego wstępnym przetwarzaniem danych, gdzie podlegają filtrowaniu, normalizacji i segmentacji. Na warstwie reprezentacji cech posegmentowane sygnały danych są ostatecznie konwertowane na spektrogram przy użyciu krótkoczasowej transformaty Fouriera (STFT).
Następnie warstwa modelu AI wykorzystuje obrazy spektrogramów do ekstrakcji cech za pomocą modułu CNN, identyfikując wzorce przestrzenne i anomalie. Następnie wyekstrahowane cechy są przekształcane i wprowadzane do koderu transformera, który uczy się zależności czasowych i relacji długodystansowych poprzez mechanizmy self-attention. Wynikowe wektory cech są klasyfikowane jako normalne lub abnormalne przy użyciu warstwy klasyfikacji. Aby zapewnić przejrzystość procesu predykcji, warstwa wyjaśnialności wykorzystuje SHAP do wyróżnienia istotnych obszarów czasowo-częstotliwościowych na spektrogramie.
Dodatkowo istnieje proces z udziałem człowieka (human-in-the-loop), w którym operatorzy interpretują wyniki modelu i dostarczają na ich podstawie informacji zwrotnych. W warstwie wdrożeniowej model działa na infrastrukturze CPHS opartej na krawędzi (edge-based), co zapewnia szybką inferencję i niski pobór mocy, z opcjonalnym interfejsem łączącym z systemami chmurowymi w celu realizacji wymagań dotyczących przechowywania danych lub zaawansowanej analityki. W całości struktura ta zapewnia dokładne wykrywanie usterek, wyjaśnialność oraz efektywne działanie w czasie rzeczywistym, zgodnie z wymogami Przemysłu 5.0.
Definicja problemu i analiza wymagań
Cyberfizyczne systemy ludzkie (CPHS) w środowiskach Przemysłu 5.0 obiecują inteligencję, zrównoważony rozwój i koncentrację na człowieku. Niemniej jednak, współczesne ramy budowy CPHS są ograniczane przez brak zaawansowanych możliwości opartych na AI, wyjaśnialności oraz odpowiedniego wsparcia dla współpracy człowieka z maszyną. W rzeczywistości konwencjonalne CPHS opierają się na modelach regułowych lub nieinterpretowalnych, co utrudnia adaptację do zmieniających się warunków przemysłowych i ogranicza zaufanie operatorów do automatyzacji. Co więcej, konwencjonalne CPHS nie potrafią efektywnie analizować wysoce złożonych danych, takich jak przemysłowe sygnały akustyczne zawierające komponenty zarówno przestrzenne (tj. częstotliwościowe), jak i zależne od czasu, które są istotne dla konserwacji predykcyjnej. W związku z tym konieczne jest stworzenie ram, które umożliwią wykrywanie anomalii maszyn z wysoką precyzją, dostarczanie wyjaśnień oraz działanie na krawędzi sieci (edge computing), co w ten sposób umożliwi funkcjonowanie CPS.
Tutaj x(t) jest sygnałem wejściowym akustycznym, S(f,t) jest reprezentacją spektrogramową sygnału,
jest etykietą a θ jest zestawem parametrów architektury CNN-Transformer. Symbole użyte w kolejnych równaniach są wprowadzane bezpośrednio po ich pierwszym wystąpieniu, aby zachować poprawność matematyczną.
Z matematycznego punktu widzenia zadanie to można sformalizować jako problem uczenia nadzorowanego. Niech sygnał akustyczny przemysłowy będzie oznaczony jako x(t). STFT sygnału x(t) oblicza się zgodnie z wzorem:
(1)
gdzie ω(.) jest funkcją okna, a f jest częstotliwością. Wygenerowany spektrogram
służy jako dane wejściowe dla algorytmu głębokiego uczenia. Celem jest wytrenowanie funkcji odwzorowania
(2)
gdzie
to etykieta klasy (normalna lub nieprawidłowa), S to wejściowy spektrogram, y to przewidywana etykieta wyjściowa, a θ to zbiór parametrów możliwych do wyuczenia w architekturze CNN-Transformer. Problem optymalizacji sformułowano jako minimalizację straty klasyfikacji:
(3)
z zastrzeżeniem ograniczeń dotyczących minimalnej latencji
, efektywności energetycznej
oraz interpretowalności, gdzie wyjaśnienie ∅(SHAP) oddaje znaczenie cech w procesie predykcji,
(4)
gdzie
to wkład każdej poszczególnej cechy. Zgodnie z powyższym, wymagania projektowe obejmują następujące elementy: (i) niezawodną detekcję anomalii opartą na hybrydowym głębokim uczeniu, (ii) interpretowalne wyniki umożliwiające interwencję człowieka oraz (iii) implementację w ramach CPS wzmocnionego krawędziowo.
Pozyskiwanie i wstępne przetwarzanie zbioru danych
W związku z tym, aby ocenić proponowaną strukturę CPHS opartą na wyjaśnialnej sztucznej inteligencji (Explainable AI), wybrano zbiór danych MIMII13. Zbiór ten jest powszechnie stosowanym wzorcem do wykrywania anomalii w maszynach przemysłowych przy użyciu sygnałów akustycznych. Baza ta zawiera nagrania dźwiękowe różnych maszyn, takich jak prowadnice przesuwne, wentylatory, pompy i zawory, w normalnych oraz awaryjnych warunkach pracy. Dla każdej maszyny nagrania różnią się w zależności od ustawień i środowiska akustycznego. Dzięki temu zbiór ten nadaje się do trenowania i oceny modeli w rzeczywistych warunkach przemysłowych.
Podczas gdy zbiór danych MIMII jest powszechnie wykorzystywany w badaniach nad nienadzorowanym wykrywaniem anomalii, które podczas treningu opierają się wyłącznie na próbkach normalnych, zawiera on jawne adnotacje zarówno dla normalnych, jak i anomalnych stanów maszyn. Problem sformułowano jako nadzorowane zadanie klasyfikacji binarnej, aby przetestować zdolność zaproponowanej architektury CNN-Transformer opartej na wyjaśnialnej sztucznej inteligencji (Explainable AI) do rozróżniania zdrowych stanów maszyn od uszkodzonych. W szczególności zastosowano oficjalne etykiety MIMII, przypisując etykietę klasy 0 do próbek normalnych oraz etykietę klasy 1 do próbek anomalnych.
Tabela 1 przedstawia szczegółowy opis zbioru danych. W niniejszych badaniach rozważono około 12 000 nagrań sygnałów akustycznych, obejmujących zarówno prawidłową, jak i wadliwą pracę maszyn. Zgodnie z wymogami, wszystkie nagrania w tej studyjnej analizie zostały zarejestrowane z częstotliwością 16 kHz i trwają 10 s. Chociaż zbiór danych zawiera około 12 000 nagrań, wybrano zrównoważony pod względem klas podzbiór 2 000 nagrań testowych wyłącznie w celu wizualizacji macierzy pomyłek. Wszystkie raportowane metryki wydajności obliczono przy użyciu pełnego zestawu testowego (około 2 400 nagrań). Ponadto próbki treningowe i testowe zostały rozdzielone w stosunku 80:20 na potrzeby trenowania modelu. Baza ta jest dobrze dostosowana do zastosowań w Przemysle 5.0, ponieważ uwzględnia warunki hałasu w czasie rzeczywistym. Rozkład nagrań wykorzystanych w tej pracy przedstawiono w Tabeli 2. Z tego zbioru danych wybrano około 12 000 nagrań obejmujących wszystkie typy maszyn i warunki stosunku sygnału do szumu (SNR) dla eksperymentów nad nadzorowaną klasyfikacją binarną z podziałem na zbiór treningowy i testowy w stosunku 80:20; tabela ta przedstawia typy maszyn przemysłowych, odpowiadające im identyfikatory maszyn, całkowitą liczbę prawidłowych i anomalnych nagrań akustycznych, warunki SNR oraz eksperymentalny podział na zbiór treningowy/testowy zastosowany w nadzorowanej klasyfikacji binarnej.
Eksperyment przeprowadzono z wykorzystaniem zbioru danych MIMII (Malfunctioning Industrial Machine Investigation and Inspection) w wersji 1.0, który obejmuje cztery typy maszyn przemysłowych: wentylatory, pompy, zawory i szyny prowadzące. Dane zostały zebrane z różnych identyfikatorów maszyn (ID00–ID06, w zależności od typu maszyny) oraz w różnych warunkach. Baza danych MIMII zawiera dane pozyskane przy stosunkach sygnału do szumu (SNR) wynoszących −6 dB, 0 dB oraz +6 dB. Dane normalne odpowiadają sprawnemu stanowi maszyny przemysłowej, natomiast dane anomalne odpowiadają stanowi wadliwemu. Klasy zostały oznaczone na podstawie adnotacji w bazie danych, gdzie dane normalne oznaczono jako klasa 0, a dane anomalne jako klasa 1.
Rzetelna procedura oceny wymagała, aby podział na zbiory treningowy i testowy nastąpił na poziomie oryginalnych nagrań audio przed jakimikolwiek etapami dzielenia, augmentacji lub tworzenia spektrogramów. W rezultacie segmenty audio z tego samego nagrania nie mogły pojawić się jednocześnie w zbiorze treningowym i testowym. W eksperymencie oceniającym wykorzystano nagrania o różnych stosunkach sygnału do szumu (−6 dB, 0 dB i +6 dB), co umożliwiło ocenę proponowanego modelu w różnych warunkach przemysłowych.
Aby zapobiec wyciekowi danych, podział na zbiór treningowy i testowy przeprowadzono w oparciu o oryginalne pliki audio, przed jakimkolwiek segmentowaniem lub generowaniem spektrogramów. Segmenty z tego samego 10 s nagrania audio przypisano wyłącznie do zbioru treningowego lub testowego, tak aby żadne segmenty nie występowały w obu zbiorach jednocześnie. Następnie spektrogramy zostały wygenerowane i rozszerzone oddzielnie dla każdego zbioru danych. W ten sposób uniknięto problemu wycieku danych, co zapobiegło zbyt optymistycznej ocenie wydajności.
Przetwarzanie wstępne zbioru danych
Surowe sygnały akustyczne z zestawu danych MIMII są następnie przetwarzane w celu poprawy jakości sygnału i umożliwienia skutecznej ekstrakcji cech. Sygnał wejściowy jest oznaczony jako x(t), gdzie t oznacza czas. Pierwszym etapem jest filtrowanie szumów, podczas którego szumy tła są usuwane za pomocą funkcji filtrującej w następujący sposób:
(5)
gdzie ∗ reprezentuje proces splotu, a h(t) jest odpowiedzią impulsową filtra. Proces ten zwiększa klarowność sygnału akustycznego i minimalizuje skutki szumów, które często występują w środowiskach przemysłowych. Normalizacja sygnału jest stosowana w celu zapewnienia, aby amplitudy wszystkich sygnałów audio były porównywalne przed wprowadzeniem ich do modelu głębokiego uczenia. Sygnał znormalizowany, oznaczony jako xn(t), jest obliczany jako:
(6)
gdzie μ oraz σ to odpowiednio wartość średnia i odchylenie standardowe sygnału wejściowego. Dzięki tej procedurze wszystkie sygnały wejściowe mają zagwarantowaną wartość średnią równą zero oraz jednostkową wariancję. Na koniec, w celu przetwarzania, znormalizowany sygnał wejściowy zostaje podzielony na wiele małych fragmentów. Przyjmując, że długość całego sygnału wynosi T, zostaje on rozbity na N nakładających się ramek o długości L przy użyciu funkcji przesuwnego okna:
(7)
gdzie H to wielkość skoku (hop size). Spektrogram można zdefiniować za pomocą poniższego wzoru:
(8)
Spektrogramy uzyskane w ten sposób są wprowadzane do frameworka opartego na sieci CNN-transformer w celu wyodrębnienia cech przestrzennych i czasowych, co umożliwia skuteczną detekcję anomalii. Rysunek 2 przedstawia potok przetwarzania wstępnego.
Rysunek 2 przedstawia pionowy potok przetwarzania wstępnego, który można zastosować do danych akustycznych pozyskanych w przemyśle. W pierwszej kolejności surowe sygnały audio są wyodrębniane z zestawu danych MIMII i przekazywane do etapu usuwania szumów w celu redukcji zakłóceń. Następnie sygnał jest normalizowany, aby wyeliminować niespójności amplitudy i zapewnić jego jednolitość. Sygnał jest później dzielony na mniejsze, nakładające się sekcje. Po segmentacji do danych stosowane są transformacje STFT, aby przekształcić sygnały z dziedziny czasu w sygnały czasowo-częstotliwościowe. Następnie, w ramach generowania spektrogramu sygnału, obliczane jest widmo mocy. W ten sposób produktem końcowym procesu będą przetworzone wstępnie spektrogramy, które zostaną wykorzystane w proponowanym modelu CNN-transformer.
Generowanie spektrogramu
Po normalizacji i segmentacji danych sygnały akustyczne są przekształcane w reprezentację czasowo-częstotliwościową przy użyciu krótkoterminowej transformaty Fouriera (STFT). Powstający spektrogram stanowi graficzną reprezentację zmian w składzie częstotliwości w okresie analizy. W warunkach przemysłowych usterki maszyn zazwyczaj objawiają się jako wzorce częstotliwości, co sprawia, że spektrogramy są bardziej użyteczne w wykrywaniu anomalii maszynowych.
W rezultacie przetworzone dane są konwertowane na obrazy spektrogramów, aby odróżnić przypadki prawidłowego i nieprawidłowego działania maszyn. Dane audio zostaną przekształcone w wiele spektrogramów, co umożliwi modelowi uczenia się analizę lokalnych komponentów częstotliwości oraz zmian czasowych w procesie pracy maszyny. Obrazy te służą jako dane wejściowe dla modelu CNN-transformer, w którym sieć CNN obsługuje wzorce przestrzenne (częstotliwościowe), natomiast transformery odpowiadają za relacje czasowe.
Algorytm 1: Generowanie spektrogramu z sygnałów akustycznych
Wejście: Surowy sygnał audio x(t)
Wynik: Spektrogram S(t, f)
Krok 1: Załadowanie sygnału audio x(t)
Krok 2: Zastosuj filtrowanie szumów
Krok 3: Normalizacja sygnału xn(t)
Krok 4: Segmentacja sygnału na klatki:
dla k = 1 do N

koniec pętli
Krok 5: Zastosuj funkcję okna w(t) do każdego segmentu
Krok 6: Obliczenie STFT:
(9)
Krok 7: Oblicz wartość amplitudy:
(10)
Krok 8: Zapisz spektrogram S(t, f)
Krok 9: Zwrócenie zbioru danych spektrogramu
Proces generowania spektrogramów rozpoczyna się od pozyskania oryginalnego sygnału audio, po czym następują etapy przetwarzania wstępnego, takie jak redukcja szumów i normalizacja sygnału, aby zapewnić jego integralność i jakość. Następnie znormalizowany sygnał audio jest dzielony na wiele nachodzących na siebie okien, aby uchwycić charakterystykę czasową zarówno okresów długoterminowych, jak i krótkoterminowych. Aby zapobiec wyciekowi widma, każdy z fragmentów jest mnożony przez funkcję okna przed przekazaniem do krótkoczasowej transformaty Fouriera (STFT). Następnie przetworzony sygnał jest reprezentowany jako funkcja o wartościach zespolonych w dziedzinie częstotliwości, a następnie podnosi się do kwadratu wartości ich amplitud. W ten sposób energia częstotliwości w danym czasie jest reprezentowana przez spektrogramy. Na koniec są one zapisywane jako dane w formie obrazów i przekazywane do sieci neuronowej w celu przeprowadzenia treningu.
Ekstrakcja cech przestrzennych przy użyciu CNN
Po procesie tworzenia spektrogramów uzyskane obrazy czasowo-częstotliwościowe są przyjmowane jako dane wejściowe i przekazywane do konwolucyjnej sieci neuronowej (CNN). Przyjmując spektrogram wejściowy jako
, gdzie H oznacza domenę częstotliwości, a W domenę czasu. Celem zastosowania sieci CNN w tym przypadku jest nauczenie się dyskryminatywnych cech przestrzennych ze spektrogramu wejściowego. Cechy te pomogą odróżnić normalne od nieprawidłowego zachowania maszyn w środowisku przemysłowym.
Podstawową operacją wykonywaną przez splotową sieć neuronową (CNN) jest splot. Splot z określonym filtrem
definiuje się jako:
(11)
Opisana powyżej operacja umożliwia sieci neuronowej identyfikację lokalnych cech przestrzennych, takich jak krawędzie, częstotliwości i zawartość harmoniczna, w spektrogramie. Odbywa się to poprzez wielokrotne operacje filtrowania, które generują różne reprezentacje cech, tworząc kilka map cech. Po splotach następuje nieliniowa funkcja aktywacji, jak pokazano poniżej:
(12)
Ta nieliniowość umożliwia sieci wychwytywanie i reprezentowanie złożonych zależności między zmiennymi. Kolejnym etapem w konwolucyjnej sieci neuronowej jest zastosowanie operacji poolingu. Powszechnie stosowana technika max poolingu jest wyrażona matematycznie w następujący sposób:
(13)
gdzie Ω oznacza okno poolingu. Pooling pomaga w osiągnięciu niezmienniczości translacyjnej i zachowaniu jedynie istotnych informacji przy jednoczesnym odrzuceniu redundancji. Wykorzystując splotową sieć neuronową, można hierarchicznie wyodrębnić cechy wyższego poziomu – od zmiennych niskopoziomowych, takich jak częstotliwości, aż po złożone zachowania maszyn. W konsekwencji mapy cech mogą być wykorzystane do przedstawienia końcowego wyniku sieci:
(14)
gdzie reprezentuje wyuczalne parametry sieci CNN. Mapy cech zawierają bogate informacje przestrzenne i są przekształcane w sekwencje, aby służyć jako wejścia dla koderach transformera w kolejnej fazie. Opisany powyżej proces ekstrakcji cech przestrzennych jest kluczowy dla umożliwienia modelowi wykrywania drobnych anomalii w przemysłowych sygnałach akustycznych.
Rycina 3 przedstawia sposób horyzontalnej ekstrakcji cech przestrzennych z wejściowego spektrogramu przy użyciu modelu CNN. Proces rozpoczyna się od spektrogramu wejściowego, który koduje właściwości czasowo-częstotliwościowe sygnału audio w formie obrazowej. Spektrogram wejściowy jest następnie przekazywany do warstwy splotowej, w której różne rodzaje filtrów ekstrahują cechy przestrzenne, takie jak pasma częstotliwości, krawędzie i tekstury, generując w ten sposób mapy cech. Następnie wyekstrahowane mapy cech przechodzą przez nieliniową funkcję aktywacji — jednostkę liniową z prostowaniem (ReLU) — co pozwala sieci neuronowej wykrywać zależności w danych wejściowych. Mapy cech są następnie przekazywane do warstwy poolingu, w której rozmiar mapy zostaje zredukowany przy zachowaniu jej najważniejszych charakterystyk. Ten proces warstw splotowych, aktywacji i poolingu powtarza się wielokrotnie, tworząc hierarchię reprezentacji cech, która zaczyna się od prostych elementów i ewoluuje w stronę bardziej złożonych cech związanych z pracą maszyn.
Modelowanie czasowe z wykorzystaniem koder transformer
Po wyekstrahowaniu cech przestrzennych za pomocą CNN, kolejnym krokiem jest wyprowadzenie map cech wysokiego poziomu, które przechwytują wzorce oparte na częstotliwościach wyuczone ze spektrogramu. Powstałe mapy cech są następnie spłaszczane i przekształcane w sekwencje w celu uchwycenia informacji czasowych w sygnale akustycznym. Przy założeniu, że mapy cech są reprezentowane przez
, gdzie H', W', i C to odpowiednio wysokość, szerokość i liczba kanałów, wyprowadzana jest następująca sekwencja:
(15)
gdzie T oznacza długość sekwencji, a d jest wymiarem cechy. Następnie Transformer wykorzystuje powyższą sekwencję do modelowania długoterminowych zależności przy użyciu mechanizmu self-attention. Model transformer, w przeciwieństwie do konwencjonalnych modeli rekurencyjnych, przetwarza każdy element sekwencji równolegle i wykorzystuje wagi uwagi do określenia względnej istotności każdego fragmentu. Mechanizm uwagi działa poprzez wykorzystanie wektorów Zapytania Q (Query), Klucza K (Key) i Wartości V (Value) generowanych dla każdego wejścia xt zgodnie z:
(16)
podczas gdy WQ, WK, WV są wagami trenowalnymi. Podobieństwa między elementami są wyznaczane za pomocą:
(17)
dk to istotny wymiar wektorów. Ta procedura ułatwia modelowi naukę zależności z odległych kroków czasowych w sygnale i pozwala skupić się wyłącznie na istotnych segmentach czasowych. Aby zwiększyć zdolność sieci do uczenia się, zastosowano mechanizm uwagi wielogłowicowej (multihead attention), który polega na równoległym uczeniu się z wejścia przez kilka głowic uwagi:
(18)
Każda głowica odpowiada za odmienne cechy sekwencji wejściowej w dziedzinie czasu. Następnie zastosowano sieci neuronowe z przesuwniem w przód oraz połączenia rezydualne:
(19)
Wynikowa zakodowana sekwencja rejestruje zarówno cechy przestrzenne, jak i czasowe danych wejściowych. Sekwencja ta jest przesyłana do warstwy klasyfikatora w celu wykrycia anomalii.
Rysunek 4 przedstawia ogólny proces modelowania czasowego z wykorzystaniem koderen transformer na cechach wyekstrahowanych z CNN. Najpierw wykorzystuje się mapy cech wyekstrahowane z CNN, które zawierają informacje przestrzenne ze spektrogramu. Mapy cech są spłaszczane do postaci sekwencji, w której każdy element odpowiada określonej znacznikowi czasu. Następnie sekwencja przechodzi przez transformację liniową, która odwzorowuje dane na wektory o określonym wymiarze. Ponieważ modele transformer nie posiadają wrodzonej zdolności do przechwytywania zależności czasowych, do wektorów zanurzeń dodawane jest kodowanie pozycyjne w celu zachowania struktury czasowej sekwencji. Na koniec sekwencja przechodzi przez koder transformer, składający się ze stosu wielu warstw, w których zastosowano mechanizm wielogłowicowej samouwag (multi-head self-attention), aby umożliwić modelowi naukę powiązań między znacznikami czasu w sekwencji. Sieć neuronowa typu feed-forward dodatkowo udoskonala te zanurzenia, podczas gdy połączenia rezydualne zapewniają efektywną propagację gradientów. Proces ten jest powtarzany w kolejnych warstwach stosu kodera transformer, które przechwytują progresywnie bardziej złożone interakcje czasowe. Ostatecznie zanurzenia sekwencji są łączone w reprezentację o stałej długości poprzez pooling (np. average pooling lub token klasyfikacji (CLS)), co daje globalną reprezentację cech czasowych.
Warstwa klasyfikacji
Po kodowaniu czasowym za pomocą koder-transformera (transformer encoder), wyjściem będzie sekwencja zakodowanych wektorów, z których każdy zawiera cechy związane zarówno z przestrzennym, jak i czasowym aspektem sygnału wejściowego. Pierwszym krokiem w kierunku klasyfikacji na podstawie zakodowanych danych wyjściowych jest zagregowanie wszystkich tych cech w jeden wektor. Osiąga się toy poprzez zastosowanie uśredniania (mean pooling) sekwencji wektorów cech lub poprzez wykorzystanie tzw. „tokenu klasyfikacyjnego” (classification token). Przyjmijmy, że wyjściem Transformera jest
, gdzie każdy
. Zagregowana cecha h jest obliczana jako:
(20)
Następnie zagregowany wektor cech h reprezentuje dane wejściowe w formie skondensowanej i służy jako wejście klasyfikatora. Część klasyfikacyjna składa się z jednej lub więcej warstw gęstych, które przekształcają wyuczoną reprezentację w przestrzeń wyjściową. Warstwa gęsta realizuje przekształcenie zdefiniowane w następujący sposób:
(21)
gdzie o to wyjście (logity), b to wektor obciążenia, a W to macierz wag. Te ostatnie stanowią nieznormowane wyniki klas. Aby wprowadzić nieliniowość do modelu, często stosuje się funkcje aktywacji, takie jak ReLU. Na koniec do logitów stosowana jest funkcja aktywacji softmax w celu uzyskania prawdopodobieństwa dla każdej klasy:
(22)
gdzie C oznacza liczbę klas (C = 2; normalna i nieprawidłowa). Predykcja klasy będzie oparta na maksymalnej wartości prawdopodobieństwa. Podczas trenowania optymalne wartości parametrów modelu są uzyskiwane poprzez minimalizację funkcji straty entropii krzyżowej między etykietami rzeczywistymi a przewidywanymi, co pozwala na dokładną klasyfikację stanów maszyn. Podsumowując, warstwa klasyfikacji pełni kluczową funkcję w przekształcaniu wyekstrahowanych cech czasowych, uzyskanych przy użyciu architektury Transformer, w praktyczne decyzje. Dzięki odpowiedniemu rozróżnianiu między zdrowymi a uszkodzonymi stanami maszyn, w paradygmacie Industry 5.0 możliwe staje się wykrywanie anomalii i inteligentne podejmowanie decyzji.
Rycina 5 przedstawia zasadę działania warstwy klasyfikacji, która następuje po enkoderze transformera. Proces rozpoczyna się od wyjściowej zakodowanej sekwencji, w której każdy wektor zawiera informacje czasowe wyuczone z danych wejściowych. Wektory te są łączone w jeden wektor, znany jako globalny wektor cech. Następnie globalny wektor cech jest przekazywany przez kilka gęstych warstw sieci neuronowej. W każdej gęstej warstwie wykonywane są operacje mnożenia macierzy w celu stworzenia kombinacji liniowych danych wejściowych. Po tym etapie można zastosować nieliniową funkcję aktywacji, taką jak ReLU, aby zwiększyć elastyczność w uczeniu granic decyzyjnych. Na końcu ostatniej gęstej warstwy znajdują się wyniki logitów dla każdej klasy. Wyniki logitów te mogą zostać wykorzystane w funkcji Softmax do obliczenia wyników prawdopodobieństwa dla każdej klasy. Wyniki prawdopodobieństwa te służą do określenia, czy dane wejściowe należą do Klasy 0 (Normalna) czy Klasy 1 (Nieprawidłowa).
Wyjaśnialność z wykorzystaniem SHAP
Po zaklasyfikowaniu sygnałów wejściowych przez model, generowane są predykcje określające, czy maszyna pracuje prawidłowo, czy w sposób nieprawidłowy. Niemniej jednak w środowisku Przemysłu 5.0 samo tworzenie predykcji jest niewystarczające; system musi również zapewnić przejrzystość i wyjaśnialność, aby wspierać proces podejmowania decyzji przez człowieka. Poniższy schemat wykorzystuje SHAP – model oparty na teorii gier, który kwantyfikuje znaczenie każdej cechy w określaniu predykcji. W ten sposób SHAP przypisuje wagi cechom, co pozwala ocenić względną istotność różnych sekcji spektrogramu wejściowego dla końcowej predykcji.
Matematyczna formulacja SHAP opiera się na dekompozycji modelu f(x):
(23)
gdzie ∅0 jest wartością bazową (średnim wynikiem modelu), ∅i jest efektem i-tej cechy, a M to liczba cech wejściowych. Na potrzeby niniejszego badania cechami są przedziały czasowo-częstotliwościowe ze spektrogramu. Wartości SHAP są dodatnie, jeśli przyczyniają się do prawdopodobieństwa wystąpienia stanu nieprawidłowego, a ujemne w przeciwnym razie. Ta formulacja gwarantuje unikalne, lokalnie spójne wyjaśnienie dla każdej predykcji. Metodę SHAP można zastosować do wytrenowanego modelu CNN-transformer, aby uzyskać wyjaśnienie oparte na istotności cech. Reprezentacja wizualna wskaże najważniejsze przedziały czasowo-częstotliwościowe w spektrogramie, które wpływają na proces podejmowania decyzji. Na przykład nieprawidłowa praca urządzenia może skutkować wystąpieniem określonych pasm częstotliwości o wysokiej energii, które framework SHAP zidentyfikuje jako istotne czynniki wpływające na proces decyzyjny.
Dodatkowo, SHAP może sprzyjać współpracy między ludźmi a maszynami w środowisku cyber-fizycznym. Wynika to z faktu, że przewidywania mogą być wyjaśnione za pomocą intuicyjnych wizualizacji, co pozwala operatorom i systemom decyzyjnym weryfikować prognozy, analizować problem i podejmować odpowiednie kroki. Nie tylko zwiększa to zaufanie, ale także sprawia, że procesy produkcyjne oparte na AI są niezawodne i rozliczalne. Dzięki integracji algorytmów głębokiego uczenia, spełniających aktualne standardy wydajności, z SHAP, który zapewnia wyjaśnialność, proponowana architektura spełnia wymagania Przemysłu 5.0.
Algorytm 2: Wyjaśnialność predykcji modelu w oparciu o SHAP
Dane wejściowe: Wytrenowany model f(x), wejściowy spektrogram S
Wynik: wartości SHAP
oraz mapa objaśnień
Krok 1: Wprowadzenie spektrogramu S do wytrenowanego modelu
Krok 2: Obliczenie prognozy: y = f(S)
Krok 3: Inicjalizacja wyjaśniacza SHAP zorientowanego na głębokie uczenie
Krok 4: Oblicz wartości SHAP:

Krok 5: Dla każdej cechy i w S:
Obliczenie wkładu 
Krok 6: Wygeneruj mapę wyjaśnień:
Wyróżnij obszary o wysokim 
Krok 7: Wizualizacja istotności cech (mapa ciepła)
Krok 8: Zwrócenie wartości SHAP i mapy wyjaśnień
Algorytm wyjaśnialności oparty na metodzie SHAP rozpoczyna pracę od pobrania przetworzonego obrazu spektrogramu jako danych wejściowych i wprowadzenia go do przeszkolonej architektury CNN-transformer w celu dokonania predykcji (np. normalny lub nieprawidłowy stan maszyny). Po uzyskaniu predykcji za pomocą przeszkolonej architektury, metoda SHAP jest wykorzystywana do interpretacji modelu. W szczególności obliczane są wartości SHAP (ϕ) dla każdej cechy wejściowej, aby ustalić, które komponenty czasowo-częstotliwościowe sygnału wejściowego mają największy wpływ na predykcję. Aby obliczyć te wartości, należy porównać wyniki modelu dla danych wejściowych z określonymi cechami oraz bez nich.
Po obliczeniu wartości SHAP dla każdej cechy, jej wpływ na prognozę można ocenić na podstawie znaku i wielkości wartości SHAP. Na przykład cechy o wyższych wartościach SHAP z większym prawdopodobieństwem wpływają na wynik, a znak dodatni lub ujemny wskazuje, że dana cecha przesuwa prognozę odpowiednio w stronę klasy prawidłowej lub nieprawidłowej. Na podstawie uzyskanych wyników można stworzyć mapę wyjaśnień w formie mapy ciepła.
Na koniec algorytm generuje zarówno wartości SHAP, jak i wizualizacje, co pozwala zrozumieć, w jaki sposób ludzie mogą analizować decyzje modelu. Technika ta zwiększa przejrzystość procesu i umożliwia zaangażowanie człowieka w podejmowanie decyzji, zapewniając tym samym zarówno dokładność przewidywań algorytmu, jak i jego niezawodność w przemyśle opartym na AI. W celu zwiększenia wiarygodności analizy wyjaśnialności wybrano reprezentatywne próbki testowe z dokładnie zdiagnozowanych anomalnych przypadków dla różnych typów maszyn. Podczas gdy globalna analiza istotności cech SHAP została przeprowadzona poprzez agregację wartości SHAP w całym zestawie danych testowych w celu zidentyfikowania spójnie wpływowych lokalizacji czasowo-częstotliwościowych, dla tych reprezentatywnych przypadków wygenerowano lokalne wyjaśnienia SHAP. Bez konieczności wprowadzania adnotacji eksperckich, ta zintegrowana analiza lokalna i globalna oferuje wiarygodną interpretację przewidywań modelu.
Integracja CPS z obsługą krawędziową
Proponowane rozwiązanie jest wykorzystywane w ramach CPS z obsługą krawędziową, którego celem jest osiągnięcie wydajności, zrównoważenia oraz operacji przemysłowych w czasie rzeczywistym. Po przeszkoleniu modelu, transformator CNN zostanie wdrożony za pomocą podejścia obliczeń krawędziowych. Wspomniana warstwa obliczeń krawędziowych znajduje się w pobliżu urządzeń przemysłowego Internetu Rzeczy (IIoT), w tym czujników i maszyn. Zamiast przesyłać wszystkie zebrane dane do chmury, urządzenie krawędziowe wykonuje wnioskowanie na przychodzącym sygnale akustycznym bez konieczności jego wcześniejszego transferu. Dzięki temu system może niemal natychmiast wykryć nieprawidłowości w działaniu maszyn i podjąć niezbędne kroki, zanim sytuacja ulegnie dalszemu pogorszeniu. Inną zaletą tego podejścia z obsługą krawędziową jest poprawiona latencja i wydajność, ponieważ dane są analizowane w pobliżu miejsca ich generowania. Innymi słowy, w przeciwieństwie do podejścia chmurowego, gdzie czas wnioskowania Tinf
(24)
Opóźnienie, w ujęciu ogólnym, zostanie znacznie zredukowane dzięki zmniejszeniu Tcommunication. Inną korzyścią jest zmniejszone zużycie przepustowości oraz energii wydatkowanej podczas ciągłej transmisji danych. System ten jest zatem skalowalny i zrównoważony w przypadku dużych zakładów przemysłowych z licznymi podłączonymi urządzeniami. Ponadto CPS umożliwia interakcję z warstwą fizyczną (maszynami), AI (warstwą obliczeniową) oraz warstwą ludzką.
Współpraca z człowiekiem w pętli (human-in-the-loop)
Kluczowym elementem sugerowanego systemu CPS w ramach Przemysłu 5.0, który obejmuje aktywny udział operatorów ludzkich w procesie podejmowania decyzji obok modeli AI, jest współpraca typu human-in-the-loop. Mianowicie, zgodnie z zaprojektowaną strukturą, system nie staje się autonomiczną „czarną skrzynką”, lecz umożliwia operatorom interakcję z nim oraz obserwację wyjaśnialnych wyników predykcji w formie objaśnień SHAP. Przedmiotowe predykcje odnoszą się do wyników (normalny/nienormalny), którym towarzyszą wizualizacje fragmentów spektrogramów wykorzystanych do ich sporządzenia.
Taka interakcja pomaga operatorom lepiej zrozumieć, w jaki sposób dokładnie tworzona jest prognoza oraz które czynniki mają na nią wpływ. Udział człowieka polega na weryfikacji wyników systemu i umożliwieniu operatorom ich zmiany w razie potrzeby, w oparciu o ich profesjonalny osąd. Współpraca z operatorami sprawia, że system jest bardziej godny zaufania, wydajny, bezpieczny i niezawodny, co pomaga zapobiegać błędom.