$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Wielokanałowe sygnały cykliczne i dystrybucja danych
W celu charakterystyki zachowania baterii litowo-jonowych podczas cykli przeprowadzono analizę reprezentatywnych wielokanałowych sygnałów cyklicznych oraz rozkładów wydajności. W badaniu wykorzystano pięć zbiorów danych baterii litowo-jonowych reprezentujących wiele układów chemicznych, w tym tlenek litu i kobaltu (LCO), tlenek niklu, kobaltu i aluminium (NCA), tlenek niklu, manganu i kobaltu (NMC) oraz LFP (Tabela 1). Podstawowe zmienne cykliczne, takie jak napięcie, prąd, pojemność, temperatura, rezystancja wewnętrzna oraz sprawność kulombowska lub energetyczna, zostały zintegrowane z cechami inżynieryjnymi, takimi jak charakterystyki pojemności przyrostowej i napięcia różnicowego, aby stworzyć wielokanałowy system charakterystyki. Profile pojemność-napięcie i rozkłady SOH dla zbioru danych CALCE (Rysunek 3A), średnie profile napięcia i rozkłady SOH dla zbioru danych SANYO (Rysunek 3B), trajektorie starzenia dla zbioru danych PANASONIC (Rysunek 3C), mapy ciepła ewolucji napięcia dla zbioru danych KOKAM (Rysunek 3D) oraz profile rozkładu napięcia dla zbioru danych GOTION (Rysunek 3E) wykazały odmienne zachowania degradacyjne w zależności od ocenianych składów chemicznych baterii. Na przykład zbiór danych GOTION IFP20100140A wykazał schodkowe charakterystyki napięcia typowe dla systemów LFP, podczas gdy zbiór danych PANASONIC NCR18650BD wykazał progresywną degradację pojemności połączoną ze przesunięciami profilu napięcia podczas długotrwałego cyklowania.
| ID zbioru danych | Typ / Model baterii | Chemia | Format | Liczba ogniw | Suma analizowanych cykli | Wykorzystane okna wczesnych cykli | Punkty próbkowania na cykl | Główne kanały cykliczne | Dodatkowe kanały stanu | Cechy pochodne | Cele predykcji |
| Zbiór danych #1 | Bateria CALCE LCO | LCO | Sakiewkowa (Pouch) | 38 | 4 216 | Pierwsze 50, pierwsze 100, pierwsze 20% życia | 128 | Napięcie, prąd, czas, pojemność | Temperatura, rezystancja wewnętrzna | dQ/dV, dV/dQ, czas ładowania, energia rozładowania | SOH, RUL, energia rozładowalna |
| Zbiór danych #2 | SANYO UR18650E | NCA | Walcowa 18650 | 24 | 2 487 | Pierwsze 50, pierwsze 100, pierwsze 20% życia | 128 | Napięcie, prąd, czas, pojemność | Temperatura, wydajność kulombowska | dQ/dV, nachylenie plateau napięcia, wydajność energetyczna | SOH, RUL, energia rozładowalna |
| Zbiór danych #3 | PANASONIC NCR18650BD | NCA | Walcowa 18650 | 42 | 5 134 | Pierwsze 50, pierwsze 100, pierwsze 20% życia | 128 | Napięcie, prąd, czas, pojemność | Temperatura, rezystancja wewnętrzna, wydajność kulombowska | dQ/dV, dV/dQ, czas ładowania, przepływ energii | SOH, RUL, energia rozładowalna |
| Zbiór danych #4 | KOKAM SLPB533459H4 | NMC | Sakiewkowa (Pouch) | 19 | 1 968 | Pierwsze 50, pierwsze 100, pierwsze 20% życia | 128 | Napięcie, prąd, czas, pojemność | Temperatura, rezystancja wewnętrzna | dQ/dV, czas rozładowania, nachylenie napięcia | SOH, RUL, energia rozładowalna |
| Zbiór danych #5 | GOTION IFP20100140A | LFP | Pryzmatyczna | 57 | 6 482 | Pierwsze 50, pierwsze 100, pierwsze 20% życia | 128 | Napięcie, prąd, czas, pojemność | Temperatura, rezystancja wewnętrzna, wydajność energetyczna | dQ/dV, czas trwania plateau, energia ładowania/rozładowania | SOH, RUL, energia rozładowalna |
Tabela 1: Charakterystyka zbiorów danych baterii litowo-jonowych, wielokanałowych zmiennych cyklicznych, cech inżynierii cech oraz celów predykcji wykorzystanych w reprodukowalnym przepływie pracy DL.Tabela podsumowuje zbiory danych baterii litowo-jonowych uwzględnione w niniejszym badaniu, w tym chemię baterii, format ogniwa, liczbę ogniw, całkowitą liczbę analizowanych cykli, okna obserwacji, gęstość próbkowania, wielokanałowe zmienne cykliczne, opracowane elektrochemiczne cechy oraz cele predykcji wykorzystane do opracowania i oceny modelu. Skróty: SOH = stan zdrowia (state of health); RUL = pozostały okres użytkowania (remaining useful life); LCO = tlenek litu i kobaltu; NCA = tlenek niklu, kobaltu i glinu; NMC = tlenek niklu, manganu i kobaltu; LFP = fosforan litu i żelaza.

Rysunek 3. Reprezentatywne wielokanałowe profile cykliczne i rozkłady stanu zdrowia baterii w zbiorach danych dotyczących akumulatorów litowo-jonowych. (A) Profile pojemność-napięcie i rozkłady SOH dla zbioru danych CALCE dla tlenku litu i kobaltu. (B) Średnie profile napięcia i rozkłady SOH dla zbioru danych SANYO dla tlenku niklu, kobaltu i glinu. (C) Trajektorie starzenia i skumulowane rozkłady SOH dla zbioru danych PANASONIC dla tlenku niklu, kobaltu i glinu. (D) Mapy ciepła ewolucji napięcia i rozkłady próbek SOH dla zbioru danych KOKAM dla tlenku niklu, manganu i kobaltu. (E) Profile rozkładu napięcia i zależności SOH od rezystancji wewnętrznej dla zbioru danych GOTION dla fosforanu litu i żelaza. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Ogólna wydajność prognostyczna dla magazynowania energii
Zaproponowany schemat postępowania wykazał niskie błędy predykcji oraz spójną wydajność prognostyczną we wszystkich zadaniach ewaluacyjnych (Tabela 2). W przypadku szacowania SOH model osiągnął wartość MAE wynoszącą 0,021 dla zbioru walidacyjnego oraz 0,024 dla zbioru testowego, przy odpowiadających im wartościach R2 wynoszących odpowiednio 0,962 i 0,955. Wyniki te wskazują na stabilną wydajność predykcyjną w ciągłym szacowaniu stanu zdrowia baterii w ocenianych warunkach. Schemat wykazał również niezawodną zdolność predykcyjną w szacowaniu RUL oraz energii rozładowywalnej, osiągając dla zbioru testowego wartości R2 odpowiednio 0,927 i 0,941. Aby dodatkowo ocenić stabilność predykcji, wszystkie parametry wydajności obliczono dla pięciu niezależnych inicjalizacji modelu przy użyciu identycznych partycji zbiorów danych. Uzyskane odchylenia standardowe oraz wartości współczynnika zmienności pozostały niskie w powtórzonych ewaluacjach, co wskazuje na ograniczoną wrażliwość na stochastyczną inicjalizację i stabilne zachowanie predykcyjne w powtarzalnych warunkach treningu. Ponadto 95% przedziały ufności wykazały wąską zmienność w powtórzonych uruchomieniach, co potwierdza odporność i reprodukowalność zaproponowanego wielokanałowego schematu postępowania.
| Zadanie predykcji | Podział zbioru danych | MAE | RMSE | MAPE (%) | R² |
| Szacowanie SOH | Walidacja | 0.021 | 0.029 | 2.47 | 0.962 |
| Szacowanie SOH | Test | 0.024 | 0.033 | 2.83 | 0.955 |
| Predykcja RUL | Walidacja | 18.6 | 25.9 | 8.74 | 0.938 |
| Predykcja RUL | Test | 21.4 | 29.7 | 9.63 | 0.927 |
| Predykcja energii rozładowalnej | Walidacja | 0.087 | 0.121 | 3.18 | 0.948 |
| Predykcja energii rozładowalnej | Test | 0.094 | 0.129 | 3.56 | 0.941 |
Tabela 2: Wydajność predykcyjna powtarzalnego przepływu pracy DL dla szacowania parametrów baterii litowo-jonowych w zestawach walidacyjnych i testowych.Tabela podsumowuje wydajność predykcyjną proponowanego przepływu pracy dla szacowania SOH, predykcji RUL oraz predykcji energii rozładowywalnej z wykorzystaniem zestawów walidacyjnych i testowych. Metryki wydajności obejmują średni błąd absolutny (MAE), pierwiastek błędu średniokwadratowego (RMSE), średni procentowy błąd absolutny (MAPE) oraz współczynnik determinacji (R2). Skróty: SOH = stan techniczny (state of health); RUL = pozostały okres użytkowania (remaining useful life); MAE = średni błąd absolutny; RMSE = pierwiastek błędu średniokwadratowego; MAPE = średni procentowy błąd absolutny.
Porównanie wydajności różnych architektur DL
Ocena różnych architektur sieci i konfiguracji wejściowych wykazała znaczne zróżnicowanie w wydajności predykcyjnej (Tabela 3). Rozkłady RMSE przed i po wyborze architektury w oparciu o powtarzalność przedstawiono na Rysunku 4A. Architektura głębokiej sieci neuronowej wykazała porównawczo wyższe wartości MAE i RMSE w stosunku do architektur opartych na sekwencjach oraz architektur hybrydowych. Porównania różnych wielokanałowych konfiguracji wejściowych podsumowano na Rysunku 4B. Modele łączące lokalną ekstrakcję cech z czasowym modelowaniem międzycyklicznym wykazały poprawioną dokładność predykcyjną w ocenianych warunkach. Modele wykazujące mniejszą zmienność predykcji grupowały się zazwyczaj przy wyższych średnich wynikach wydajności i mniejszych odchyleniach standardowych (Rysunek 4C). Zależności między średnią wydajnością a odchyleniem standardowym dla ocenianych konfiguracji wejściowych zilustrowano na Rysunku 4D. Spośród testowanych konfiguracji, architektura CNN–LSTM osiągnęła najniższą wartość MAE (0,024) oraz najwyższą wartość R2 (0,955) w ramach ocenianych zbiorów danych i ustawień eksperymentalnych. Ponadto integracja cech wygenerowanych z pełnymi wejściami wielokanałowymi obniżyła MAE do 0,022 przy wartości R2 wynoszącej 0,961, podczas gdy wejścia oparte wyłącznie na pojemności wykazały niższą wydajność predykcyjną.
| Kategoria | Model / Konfiguracja wejściowa | MAE | RMSE | MAPE (%) | R² |
| Architektura głębokiego uczenia | DNN | 0.031 | 0.042 | 3.58 | 0.928 |
| CNN | 0.027 | 0.037 | 3.14 | 0.941 |
| LSTM | 0.026 | 0.036 | 3.05 | 0.944 |
| CNN–LSTM | 0.024 | 0.033 | 2.81 | 0.955 |
| Transformer | 0.025 | 0.034 | 2.89 | 0.951 |
| Konfiguracja wejściowa | Tylko pojemność | 0.034 | 0.046 | 3.92 | 0.916 |
| Napięcie + prąd | 0.029 | 0.039 | 3.29 | 0.936 |
| Pełne wejście wielokanałowe | 0.024 | 0.033 | 2.83 | 0.955 |
| Wielokanałowe + cechy inżynieryjne | 0.022 | 0.031 | 2.61 | 0.961 |
Tabela 3: Porównanie wydajności predykcyjnej dla różnych architektur DL i konfiguracji wejściowych.Tabela porównuje wydajność predykcyjną różnych architektur DL oraz konfiguracji wejściowych zastosowanych w powtarzalnym przepływie pracy wielokanałowej. Metryki wydajności obejmują średni błąd absolutny (MAE), pierwiastkowy średni błąd kwadratowy (RMSE), średni bezwzględny błąd procentowy (MAPE) oraz współczynnik determinacji (R2). Skróty: DNN = głęboka sieć neuronowa; CNN = splotowa sieć neuronowa; LSTM = długa krótkotrwała pamięć; MAE = średni błąd absolutny; RMSE = pierwiastkowy średni błąd kwadratowy; MAPE = średni bezwzględny błąd procentowy.

Rycina 4. Porównanie architektur modeli i konfiguracji danych wejściowych w oparciu o powtarzalność. (A) Rozkłady błędu średniokwadratowego przed i po wyborze architektury w oparciu o powtarzalność. (B) Rozkłady błędu średniokwadratowego dla różnych konfiguracji danych wejściowych przed i po wyborze stabilnego modelu. (C) Porównanie średniej wydajności i odchylenia standardowego dla kandydujących architektur modeli, gdzie kolor oznacza błąd średniokwadratowy, a puste okręgi oznaczają wybrane stabilne modele. (D) Porównanie średniej wydajności i odchylenia standardowego dla konfiguracji danych wejściowych, gdzie kolor oznacza błąd średniokwadratowy, a puste okręgi oznaczają wybrane optymalne konfiguracje. Aby wyświetlić powiększoną wersję tej ryciny, kliknij tutaj.
Rozkład błędów i ocena spójności
Rozkłady błędów bezwzględnych przeanalizowano w celu oceny stabilności predykcji w odniesieniu do konwencjonalnych modeli bazowych ML oraz architektur DL. Tradycyjne modele bazowe ML wykazywały szersze rozkłady błędów i wydłużone ogony błędów wysokich, co wskazuje na zmniejszoną stabilność w zadaniach predykcji złożonych wzorców degradacji. Spośród ocenianych modeli DL, architektura CNN–LSTM wykazała porównawczo węższe rozkłady błędów i krótsze ogony błędów wysokich, co sugeruje lepszą integrację ekstrakcji cech lokalnych z modelowaniem sekwencji czasowych (Rysunek 5). Włączenie zmiennych związanych z temperaturą, rezystancją wewnętrzną i sprawnością wiązało się z sukcesywnym zmniejszaniem dyspersji błędu predykcji w analizowanych zbiorach danych.

Rysunek 5. Rozkłady błędu bezwzględnego dla konwencjonalnych baz uczenia maszynowego (ML) i wielokanałowych modeli głębokiego uczenia (DL). Wykresy skrzypcowe porównujące rozkład bezwzględnych błędów predykcji dla konwencjonalnych baz ML, architektur DL oraz wielokanałowych konfiguracji wejściowych. Czarne linie oznaczają średni błąd bezwzględny, a kolorowe linie oznaczają średnią ± odchylenie standardowe. Aby wyświetlić powiększoną wersję tego rysunku, kliknij tutaj.
Powtarzalność i odporność w różnych scenariuszach walidacji
Przepływ pracy (workflow) wykazał spójną wydajność predykcyjną w powtórzonych cyklach szkolenia oraz przy zastosowaniu wielu strategii walidacji (Tabela 4). Powtórne uruchomienia przeprowadzone z wykorzystaniem stałych partycji danych i stałej inicjalizacji ziarna losowego (random-seed) pozwoliły uzyskać stabilne wskaźniki wydajności, przy czym wartości R2 konsekwentnie pozostawały powyżej 0,95, a odchylenia standardowe MAE były niskie. Dodatkowa ocena z wykorzystaniem pięciokrotnej walidacji krzyżowej oraz zewnętrznych zbiorów walidacyjnych wykazała ograniczoną zmienność wydajności w powtórzonych warunkach eksperymentalnych. Faza walidacji zewnętrznej opierała się wyłącznie na zbiorze danych GOTION, który został całkowicie wykluczony z procedur szkolenia modelu i optymalizacji hiperparametrów. Ponieważ zbiór ten reprezentuje odmienną chemię LFP oraz konfigurację ogniw pryzmatycznych, zewnętrzna ocena zapewniła rygorystyczną weryfikację zdolności do generalizacji między zbiorami danych w heterogenicznych warunkach elektrochemicznych i strukturalnych. Ponadto analizy współczynnika zmienności oraz oceny powtórnych uruchomień wykazały niską dyspersję predykcyjną w niezależnych eksperymentach, co potwierdza solidność i powtarzalność proponowanego wielokanałowego przepływu pracy.
| Ustawienia walidacji | Metryka wydajności | Średnia | SD | Minimum | Maksimum | Współczynnik zmienności (%) |
| Powtórzone uruchomienia ze stałym podziałem danych (n = 5) | MAE | 0.023 | 0.0021 | 0.021 | 0.026 | 9.13 |
| RMSE | 0.032 | 0.0028 | 0.029 | 0.036 | 8.75 |
| MAPE (%) | 2.76 | 0.24 | 2.48 | 3.11 | 8.7 |
| R² | 0.957 | 0.006 | 0.949 | 0.964 | 0.63 |
| Powtórzone uruchomienia z różnymi ziarnami losowości (n = 5) | MAE | 0.024 | 0.0024 | 0.021 | 0.027 | 10 |
| RMSE | 0.033 | 0.0031 | 0.029 | 0.037 | 9.39 |
| MAPE (%) | 2.84 | 0.27 | 2.52 | 3.19 | 9.51 |
| R² | 0.954 | 0.007 | 0.946 | 0.962 | 0.73 |
| Pięciokrotna walidacja krzyżowa | MAE | 0.025 | 0.003 | 0.022 | 0.029 | 12 |
| RMSE | 0.034 | 0.0034 | 0.03 | 0.039 | 10 |
| MAPE (%) | 2.95 | 0.31 | 2.57 | 3.41 | 10.51 |
| R² | 0.951 | 0.009 | 0.94 | 0.961 | 0.95 |
| Walidacja zewnętrzna | MAE | 0.027 | 0.0032 | 0.023 | 0.031 | 11.85 |
| RMSE | 0.037 | 0.0038 | 0.032 | 0.042 | 10.27 |
| MAPE (%) | 3.18 | 0.35 | 2.76 | 3.67 | 11.01 |
| R² | 0.944 | 0.011 | 0.931 | 0.956 | 1.17 |
Tabela 4: Analiza powtarzalności i odporności proponowanego wielokanałowego przepływu pracy DL w ramach wielokrotnych uruchomień treningowych i strategii walidacji.Tabela podsumowuje wydajność w zakresie powtarzalności i odporności proponowanego przepływu pracy podczas wielokrotnych uruchomień z ustalonymi partycjami danych, wielokrotnych uruchomień z różnymi ziarnami losowymi, pięciokrotnej walidacji krzyżowej oraz walidacji zewnętrznej. Metryki wydajności obejmują średni błąd bezwzględny (MAE), pierwiastkowy błąd średniokwadratowy (RMSE), średni bezwzględny błąd procentowy (MAPE), współczynnik determinacji (R2), odchylenie standardowe (SD) oraz współczynnik zmienności. Skróty: MAE = średni błąd bezwzględny; RMSE = pierwiastkowy błąd średniokwadratowy; MAPE = średni bezwzględny błąd procentowy; SD = odchylenie standardowe.
Wrażliwość hiperparametrów i ablacja cech
Analiza wrażliwości hiperparametrów wykazała, że wydajność predykcji różniła się w zależności od rozmiaru okna obserwacyjnego, funkcji aktywacji, rozmiaru wymiarów ukrytych oraz głębokości modelu sekwencyjnego. Wpływ rozmiaru okna obserwacyjnego na rozkłady błędów predykcji przedstawiono na Rysunku 6A. W badanych warunkach okna obserwacyjne od 100 do 200 cykli wiązały się z niższymi błędami predykcji, podczas gdy dłuższe okna obserwacyjne wykazywały malejące przyrosty dokładności predykcyjnej. Zmienność predykcji zależna od funkcji aktywacji została podsumowana na Rysunku 6B. Funkcje aktywacji ReLU (Rectified Linear Unit) oraz GELU (Gaussian Error Linear Unit) wykazały niższe błędy predykcji w porównaniu z alternatywnymi funkcjami aktywacji w badanych warunkach. Wrażliwość wydajności predykcji na rozmiar wymiarów ukrytych oraz głębokość modelu sekwencyjnego przedstawiono na Rysunku 6C. Wymiary ukryte równe 128 oraz dwie do trzech warstw modelowania sekwencyjnego wiązały się z poprawioną wydajnością predykcyjną, podczas gdy płytkie architektury i zredukowane wymiary ukryte wiązały się ze zwiększonym błędem predykcji. Analiza ablacji cech dodatkowo wykazała, że usunięcie cech wygenerowanych lub zmiennych związanych z wydajnością zwiększyło wartość MAE w stosunku do pełnej konfiguracji wielokanałowej (Tabela 5), co sugeruje, że połączenie wejść wielokanałowych i wygenerowanych cech elektrochemicznych przyczyniło się do poprawy wydajności predykcyjnej.

Rysunek 6. Wrażliwość wydajności predykcji na zmienność hiperparametrów. (A) Wpływ rozmiaru okna wejściowego na rozkład błędów predykcji. (B) Wpływ wyboru funkcji aktywacji na rozkład błędów predykcji. (C) Mapa ciepła przedstawiająca wrażliwość predykcji na rozmiar wymiaru ukrytego oraz liczbę warstw modelu sekwencyjnego. Czarne linie oznaczają średni błąd bezwzględny, a kolorowe linie oznaczają średnią ± odchylenie standardowe, gdzie dotyczy. Aby wyświetlić większą wersję tego rysunku, kliknij tutaj.
| Kategoria | Ustawienie | MAE | RMSE | MAPE (%) | R² |
| Rozmiar okna wejściowego | 20 cykli | 0.031 | 0.043 | 3.74 | 0.927 |
| 50 cykli | 0.026 | 0.036 | 3.08 | 0.944 |
| 100 cykli | 0.023 | 0.032 | 2.71 | 0.957 |
| 150 cykli | 0.022 | 0.031 | 2.63 | 0.96 |
| 200 cykli | 0.022 | 0.031 | 2.59 | 0.961 |
| 250 cykli | 0.023 | 0.032 | 2.67 | 0.958 |
| 300 cykli | 0.024 | 0.033 | 2.79 | 0.955 |
| Funkcja aktywacji | ReLU | 0.023 | 0.032 | 2.74 | 0.957 |
| GELU | 0.022 | 0.031 | 2.61 | 0.961 |
| Tanh | 0.026 | 0.036 | 3.05 | 0.946 |
| Sigmoid | 0.034 | 0.047 | 3.96 | 0.918 |
| Wymiar ukryty / liczba kanałów na warstwę | 32 | 0.028 | 0.039 | 3.22 | 0.939 |
| 64 | 0.024 | 0.033 | 2.81 | 0.954 |
| 128 | 0.022 | 0.031 | 2.58 | 0.961 |
| 256 | 0.023 | 0.032 | 2.69 | 0.958 |
| Liczba warstw modelowania sekwencji | 1 | 0.027 | 0.038 | 3.14 | 0.942 |
| 2 | 0.023 | 0.032 | 2.73 | 0.956 |
| 3 | 0.022 | 0.031 | 2.6 | 0.961 |
| 4 | 0.024 | 0.034 | 2.86 | 0.952 |
| Ablacja cech | Pełny wielokanałowy zestaw + cechy generowane | 0.022 | 0.031 | 2.61 | 0.961 |
| Bez cech generowanych | 0.024 | 0.033 | 2.84 | 0.955 |
| Bez temperatury | 0.025 | 0.034 | 2.93 | 0.951 |
| Bez rezystancji wewnętrznej | 0.026 | 0.035 | 3.01 | 0.949 |
| Bez zmiennych związanych z wydajnością | 0.025 | 0.034 | 2.95 | 0.95 |
| Tylko napięcie + prąd | 0.029 | 0.039 | 3.29 | 0.936 |
| Tylko pojemność | 0.034 | 0.046 | 3.92 | 0.916 |
Tabela 5: Analiza wrażliwości wydajności modelu przy różnych ustawieniach hiperparametrów i konfiguracjach ablacji cech. Tabela podsumowuje wpływ rozmiaru okna wejściowego, funkcji aktywacji, rozmiaru wymiaru ukrytego, głębokości modelu sekwencyjnego oraz strategii ablacji cech na wydajność predykcji w ramach proponowanego wielokanałowego przepływu pracy DL. Metryki wydajności obejmują średni błąd bezwzględny (MAE), pierwiastkowy błąd średniokwadratowy (RMSE), średni procentowy błąd bezwzględny (MAPE) oraz współczynnik determinacji (R2). Skróty: MAE = średni błąd bezwzględny; RMSE = pierwiastkowy błąd średniokwadratowy; MAPE = średni procentowy błąd bezwzględny; ReLU = jednostka liniowa z prostowaniem; GELU = jednostka liniowa z błędem Gaussa.