Nieliniowe modelowanie dynamiczne systemu regulacji stresu i zachowań emocjonalnych w szkolnictwie wyższym

Podsumowanie

Niniejszy protokół przedstawia inspirowane fizyką nieliniowe ramy przestrzeni stanów do modelowania dynamiki stresu i emocji w szkolnictwie wyższym, umożliwiając użytkownikom analizę progów wypalenia, oscylacyjnych cykli stresu, odporności oraz regulacji adaptacyjnej z wykorzystaniem równań różniczkowych, analizy stabilności i symulacji stochastycznych.

Streszczenie

Regulację stresu i zachowania emocjonalne w szkolnictwie wyższym powszechnie bada się przy użyciu opisowych modeli statystycznych, którym brakuje mechanistycznej struktury dynamicznej. Jednak stres akademicki wykazuje zachowania nieliniowe, w tym efekty progowe, wypalenie oraz odpowiedzi adaptacyjne, które nie są w sposób adekwatny ujmowane przez modele liniowe. Niniejsze badanie przyjmuje perspektywę nieliniowych systemów dynamicznych w celu modelowania regulacji stresu i zachowań emocjonalnych w szkolnictwie wyższym. Głównym celem badania jest opracowanie inspirowanych fizyką nieliniowych ram przestrzeni stanów, które rygorystycznie modelują sprzężenie stres–emocje oraz regulację adaptacyjną w szkolnictwie wyższym. Protokół ustanawia nieliniowy model przestrzeni stanów opisujący stres, aktywację emocjonalną i zdolność regulacyjną. Użytkownicy mogą wdrożyć ten protokół, aby symulować dynamikę systemu i oceniać stabilność w zmiennych warunkach akademickich. Ramy te obejmują nieliniowe analizy stabilności oraz symulacje stochastyczne w celu zbadania odporności, niestabilności i przejść związanych z wypaleniem. Wyniki symulacji wykazują nieliniową stabilność, efekty odporności oraz przejścia między reżimami, które nie są wychwytywane przez konwencjonalne modele liniowe. Wyniki wskazują, że regulacja stresu w szkolnictwie wyższym zachowuje się jak nieliniowy system dynamiczny. Poprzez integrację nieliniowej dynamiki, modelowania krajobrazu energetycznego i teorii regulacji adaptacyjnej, praca ta dostarcza zunifikowanych ram prognostycznych, które wyjaśniają procesy odporności, niestabilności i regeneracji w środowiskach akademickich. Model ten stanowi podstawę dla przyszłych badań nad dobrostanem studentów i ryzykiem wypalenia. Opracowane ramy są istotne dla badaczy i praktyków zainteresowanych dobrostanem studentów oraz regulacją stresu.

Wprowadzenie

Regulację stresu i zachowania emocjonalne w szkolnictwie wyższym można interpretować jako złożony proces dynamiczny, który ewoluuje w czasie pod wspólnym wpływem wewnętrznych sprzężeń zwrotnych o charakterze psychologicznym oraz zewnętrznych presji akademickich¹. Prace semestralne, egzaminy, oczekiwania instytucjonalne i ocena społeczna nieustannie zaburzają stany emocjonalne studentów, wywołując wahania poziomu lęku, motywacji i zaangażowania². Odpowiedzi emocjonalne te ewoluują w wielowymiarowej przestrzeni fazowej sterowanej przez ocenę poznawczą, strategie radzenia sobie oraz mechanizmy regulacji adaptacyjnej³. Klasyczne modele odpowiedzi stresowej w psychologii edukacyjnej zazwyczaj zakładają proporcjonalne zależności między bodźcami stresowymi a reakcjami emocjonalnymi, często opierając się na analizach regresji lub opisowych wskaźnikach statystycznych w celu scharakteryzowania tych relacji⁴. Jednak obserwacje empiryczne konsekwentnie wykazują cechy nieliniowe, takie jak efekty progowe, nasycenie, histereza i gwałtowne przejścia, w tym wypalenie oraz załamanie emocjonalne⁵. Małe stresory mogą powodować pomijalne przesunięcia emocjonalne w stabilnych warunkach; jednak w pobliżu progów krytycznych mogą one wywołać nieproporcjonalne reakcje emocjonalne⁶. Takie zachowanie przypomina oscylacje nieliniowe i zjawiska bifurkacji obserwowane w układach fizycznych utrzymywanych daleko od stanu równowagi⁷. W tej perspektywie studentów można konceptualizować jako otwarte układy dyssypatywne, które nieustannie wymieniają związaną ze stresem "energię" z ich środowiskiem akademickim poprzez procesy sprzężenia zwrotnego sterujące regulacją i adaptacją⁸.

Z perspektywy edukacyjnej i psychologicznej pojęcia wykorzystane w niniejszych ramach teoretycznych można interpretować intuicyjnie. Akumulacja stresu reprezentuje postępujące obciążenie generowane przez wymagania akademickie i presję środowiskową. Aktywacja emocjonalna odzwierciedla zmiany w stanach afektywnych, takich jak lęk, frustracja, motywacja i zaangażowanie, podczas gdy zdolność regulacyjna odpowiada mechanizmom radzenia sobie i samoregulacji, które pomagają studentom adaptować się do wyzwań. W przyszłej implementacji empirycznej akumulację stresu można reprezentować za pomocą skal stresu postrzeganego; aktywację emocjonalną można oceniać przy użyciu zwalidowanych narzędzi do pomiaru lęku, nastroju lub dobrostanu emocjonalnego; a zdolność regulacyjną można szacować za pomocą miar odporności psychicznej, radzenia sobie lub samoregulacji. Takie mapowania stanowią praktyczną podstawę dla przyszłej estymacji parametrów i kalibracji modelu. W tym kontekście pojęcia nieliniowe, takie jak stabilność, odporność i przejścia krytyczne, można rozumieć jako odrębne wzorce adaptacji psychologicznej. Stany stabilne reprezentują zrównoważone funkcjonowanie emocjonalne, odporność odzwierciedla zdolność do regeneracji po zaburzeniach, a przejścia krytyczne odpowiadają gwałtownym przesunięciom w stronę wyników maladaptacyjnych, takich jak stres chroniczny lub wypalenie. Interpretacje te budują wyraźniejszy pomost koncepcyjny między psychologią edukacyjną a teorią nieliniowych systemów dynamicznych, zwiększając tym samym dostępność proponowanych ram teoretycznych dla odbiorców interdyscyplinarnych.

Aby przedstawić koncepcyjny przegląd tych interakcji, Rysunek 1 ilustruje strukturę dynamicznych interakcji stres–emocje w systemach szkolnictwa wyższego. Schemat przedstawia sygnały wejściowe stresu wynikające z wymagań akademickich, zmienne reakcji emocjonalnej reprezentujące ewolucję stanu psychicznego oraz adaptacyjne mechanizmy regulacyjne, które modulują system poprzez pętle sprzężenia zwrotnego. Komponenty te oddziałują na siebie dynamicznie, tworząc trajektorie w przestrzeni fazowej stres–emocje, które mogą zbiegać się w stronę stabilnych stanów emocjonalnych lub rozbiegać w stronę reżimów nieadaptacyjnych. Chociaż ta reprezentacja koncepcyjna podkreśla dynamiczny charakter regulacji emocjonalnej, istniejące modele psychologiczne rzadko przekładają te zależności na jawne struktury matematyczne opisujące ewolucję czasową.

Rycina 1: Schemat koncepcyjny dynamicznej interakcji stres–emocje w systemach szkolnictwa wyższego. Rycina ta ilustruje koncepcyjną strukturę dynamiki stresu i emocji w środowiskach akademickich. Schemat przedstawia czynniki stresogenne wynikające z wymagań akademickich, zmienne odpowiedzi emocjonalnej reprezentujące ewolucję stanów psychologicznych oraz adaptacyjne mechanizmy regulacyjne, które modulują system poprzez pętle sprzężenia zwrotnego. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Ograniczenia obecnych podejść do modelowania zostały dodatkowo zilustrowane na Rysunku 2, który przedstawia koncepcyjne ramy identyfikujące główne luki badawcze w istniejącej literaturze dotyczącej modelowania stresu. W centrum schematu znajduje się brak mechanistycznych ram dynamicznych dla interakcji stres–emocje. Wokół tego ograniczenia skupia się kilka powiązanych ze sobą luk charakteryzujących obecne badania. Po pierwsze, wiele badań opiera się na wskaźnikach stresu opartych na korelacji, którym brakuje równań rządzących opisujących ewolucję w czasie i dynamikę sprzężeń zwrotnych⁹. Po drugie, wypalenie zawodowe i niestabilność emocjonalna są często opisywane jakościowo, a nie poprzez formalną analizę stabilności z wykorzystaniem macierzy Jacobiego lub kryteriów opartych na wartościach własnych¹⁰. Po trzecie, istniejącym modelom brakuje reprezentacji krajobrazu energetycznego, które byłyby zdolne do kwantyfikacji odporności lub geometrii basenu atrakcji w systemach regulacji emocjonalnej¹¹. Wreszcie, stres i reakcje emocjonalne są często traktowane jako zmienne niezależne, a nie jako dynamicznie sprzężone zmienne stanu rządzone nieliniowymi procesami sprzężenia zwrotnego¹². Wspólnie luki te podkreślają potrzebę opracowania rygorystycznych ram modelowania dynamicznego, zdolnych do opisu zachowań oscylacyjnych, multistabilności i przejść fazowych w systemach stresu w szkolnictwie wyższym.

Rysunek 2: Rama koncepcyjna ilustrująca luki badawcze w istniejącej literaturze dotyczącej modelowania stresu. Rysunek ten przedstawia koncepcyjny przegląd głównych ograniczeń we współczesnych podejściach do modelowania stresu. W centrum schematu znajduje się brak mechanistycznej ramy dynamicznej opisującej interakcję między stresem a emocjami. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Aby rozwiązać te ograniczenia, na Rysunku 3 przedstawiono zaproponowaną w niniejszym badaniu architekturę nieliniowego modelowania dynamicznego. Ramy te integrują dynamikę akumulacji stresu, zmienne aktywacji emocjonalnej oraz adaptacyjne mechanizmy sprzężenia zwrotnego w ramach ujednoliconego nieliniowego sformułowania przestrzeni stanów. W tym ujęciu ewolucję systemu można opisać za pomocą nieliniowych równań różniczkowych, które oddają dynamikę regulacji wewnętrznej, efekty wymuszania zewnętrznego oraz perturbacje stochastyczne wynikające z nieoczekiwanych zdarzeń akademickich. Takie sformułowanie pozwala na analizę oddziaływań stres–emocje przy użyciu narzędzi z teorii nieliniowych systemów dynamicznych, w tym analizy stanu równowagi, analizy bifurkacji, teorii stabilności Lapunowa oraz metod stabilności stochastycznej^13,14,15,16. Poprzez wprowadzenie jawnych równań sterujących, model umożliwia rygorystyczną analizę tego, jak stany emocjonalne ewoluują w czasie pod wpływem ciągłego wymuszania środowiskowego i wewnętrznych procesów sprzężenia zwrotnego.

Rysunek 3: Proponowana architektura nieliniowego modelowania dynamicznego dla regulacji stresu i emocji. Niniejszy rysunek przedstawia opracowane w tej pracy nieliniowe ramy dynamiczne służące do modelowania interakcji między stresem a emocjami w systemach szkolnictwa wyższego. Architektura ta integruje dynamikę kumulacji stresu, zmienne aktywacji emocjonalnej oraz adaptacyjne mechanizmy sprzężenia zwrotnego w formie zunifikowanego sformułowania przestrzeni stanów. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Ogólnym celem metody przedstawionej w niniejszej pracy jest ustanowienie mechanistycznych ram modelowania nieliniowego, które umożliwiają analizę regulacji stresu w szkolnictwie wyższym jako układu dynamicznego zdolnego do wykazania multistabilności, reżimów oscylacyjnych oraz przejść krytycznych. W przeciwieństwie do konwencjonalnych liniowych modeli indeksu stresu, proponowane ramy jawnie uwzględniają nieliniowe oddziaływania sprzężenia zwrotnego między skumulowaną energią stresu, aktywacją emocjonalną a adaptacyjną zdolnością regulacyjną^17,18,19,20. Te nieliniowe oddziaływania umożliwiają identyfikację stanów równowagi, charakterystykę marginesów stabilności oraz wykrywanie punktów krytycznych związanych z przejściami w stan wypalenia. Ponadto w ramach tego podejścia wprowadza się interpretację regulacji emocjonalnej opartą na krajobrazie energetycznym, w której potencjały typu Ljapunowa opisują strukturę stabilności stanów emocjonalnych^21,22,23. W tej reprezentacji odporność psychiczną można określić za pomocą głębokości basenu przyciągania, podczas gdy progi niestabilności odpowiadają punktom bifurkacji, w których zanikają stabilne atraktory.

Aby zwiększyć realizm, model uwzględnia okresowe człony wymuszające reprezentujące cykle semestrów akademickich, a także perturbacje stochastyczne oddające nieoczekiwane zdarzenia akademickie lub fluktuacje środowiskowe^24,25,26. Rozszerzenia te pozwalają strukturze na symulowanie realistycznej dynamiki stresu akademickiego, w tym oscylacyjnych cykli stresu w okresach egzaminacyjnych oraz przejść między stanami emocjonalnymi napędzanych szumem. Takie możliwości modelowania zapewniają przewagę nad podejściami czysto statystycznymi, umożliwiając analizy prognostyczne progów niestabilności i mechanizmów odporności.

Nieliniowe podejścia dynamiczne były wcześniej stosowane w takich dziedzinach jak regulacja biologiczna, neuronauka i złożone systemy społeczne, w których procesy napędzane sprzężeniem zwrotnym generują zachowania emergentne, których nie można wyjaśnić za pomocą modeli liniowych^27,28,29. Zastosowanie podobnych zasad regulacji stresu w szkolnictwie wyższym łączy psychologię edukacyjną z teorią systemów dynamicznych, umożliwiając analizę procesów regulacji emocji w ramach rygorystycznych matematycznie ram. Poprzez połączenie nieliniowego modelowania dynamicznego, teorii stabilności oraz analizy krajobrazu energetycznego, opisana tutaj metoda zapewnia ujednoliconą platformę do badania dynamiki stresu, mechanizmów odporności i stabilności emocjonalnej w środowisku akademickim³⁰.

W ostatecznym rozrachunku ramy te umożliwiają badaczom analizę tego, w jaki sposób kumuluje się stres, jak ewoluują reakcje emocjonalne oraz jak adaptacyjne mechanizmy regulacyjne wpływają na długoterminową stabilność systemu. Podejście to stanowi zatem potężne narzędzie metodologiczne do badania odporności, przejść w stan wypalenia oraz procesów regeneracji w systemach szkolnictwa wyższego. Poza wkładem metodologicznym, proponowane ramy mają służyć jako interdyscyplinarny pomost między psychologią edukacyjną a nauką o nieliniowych systemach dynamicznych. Poprzez przełożenie procesów psychologicznych, takich jak kumulacja stresu, aktywacja emocjonalna, odporność i regulacja adaptacyjna, na matematycznie operacyjne zmienne stanu, ramy te zapewniają wspólny język analityczny do badania złożonych zjawisk edukacyjnych za pomocą rygorystycznych narzędzi dynamicznych. Integracja ta zwiększa klarowność koncepcyjną, ułatwia przyszłą kalibrację empiryczną z wykorzystaniem pomiarów behawioralnych i psychologicznych oraz wspiera rozwój modeli predykcyjnych i wyjaśniających regulację stresu w środowiskach akademickich³¹. Mówiąc prościej, proponowane ramy badają, w jaki sposób kumuluje się stres akademicki, jak reakcje emocjonalne zmieniają się w czasie oraz jak mechanizmy radzenia sobie wpływają na regenerację i odporność. Model ma na celu zidentyfikowanie warunków, w których studenci zachowują stabilność emocjonalną lub stają się podatni na wypalenie, łącząc tym samym teorię psychologiczną z ilościową analizą systemową.

Protokół

Niniejszy protokół nie obejmuje eksperymentów na ludziach ani zwierzętach. Badanie opiera się wyłącznie na modelowaniu matematycznym i symulacjach komputerowych dynamiki stresu i emocji w systemach szkolnictwa wyższego. Narzędzia obliczeniowe, środowiska programowe i zasoby numeryczne niezbędne do wdrożenia opisanego w protokole modelu nieliniowej dynamiki zostały wymienione w Tabeli Materiałów.

1. Inicjalizacja modelu i konfiguracja środowiska obliczeniowego

1. Otwórz środowisko obliczeniowe i utwórz nowy katalog projektu dla badania symulacyjnego. Szczegółowe informacje na temat oprogramowania znajdują się w Pliku uzupełniającym 1.
2. Utwórz główny plik skryptu o nazwie StressEmotionModel.m w celu implementacji nieliniowego modelu stres–emocje–regulacja.
3. Zdefiniuj trzy zmienne stanu reprezentujące system dynamiczny.
   1. Zdefiniuj S(t) jako skumulowany stres akademicki. Zdefiniuj E(t) jako aktywację emocjonalną.
   2. Zdefiniuj R(t) jako adaptacyjną zdolność regulacyjną.
4. Utwórz plik inicjalizacji parametrów i przypisz wartości wszystkim podstawowym parametrom modelu.
5. Zdefiniuj podstawowy zestaw parametrów w następujący sposób.
   1. Ustaw parametr wymuszenia zewnętrznego F na 2.0. Ustaw parametr rozpraszania stresu β na 0.30.
   2. Ustaw parametr sprzężenia stres–emocje γ na 0.50. Ustaw parametr wzmocnienia regulacji κ na 0.40.
   3. Ustaw parametr wrażliwości emocjonalnej α na 0.80. Ustaw parametr relaksacji emocjonalnej δ na 0.50.
   4. Ustaw nieliniowy parametr nasycenia µ na 0.10. Ustaw parametr oddziaływania zapośredniczonego przez regulację η na 0.20. Ustaw parametr zaniku regulacji λ na 0.25.
6. Zdefiniuj wektor stanu początkowego.
   1. Ustaw S(0) = 1.0.
   2. Ustaw E(0) = 0.5.
   3. Ustaw R(0) = 0.8.
      UWAGA: Używaj identycznych wartości parametrów we wszystkich podstawowych symulacjach, chyba że przeprowadzana jest analiza wrażliwości.

2. implementacja nieliniowego modelu dynamicznego

1. Utwórz plik funkcji zawierający sprzężone równania różniczkowe opisujące dynamikę stresu, emocji i regulacji. Równania sterujące podano w Supplementary File 1 (Section S1.2), natomiast odpowiadającą im implementację w programie MATLAB przedstawiono w Supplementary File 1 (Section S1.8).
   UWAGA: Bazowe parametry modelu, warunki początkowe, konfigurację solvera numerycznego oraz sterujące nieliniowe równania różniczkowe zwyczajne podano w Supplementary File 1 (Sections S1.2–S1.6). Parametry te są stosowane konsekwentnie we wszystkich symulacjach numerycznych, chyba że zaznaczono inaczej.
2. Zaimplementuj równanie ewolucji stresu dla zmiennej stanu stresu skumulowanego S(t) zgodnie ze sterującym nieliniowym systemem dynamicznym opisanym w Supplementary File 1, Section S1.2. Odpowiednią implementację obliczeniową tego równania podano w Supplementary File 1, Section S1.8 (Function 1: Model Equations).
3. Zaimplementuj równanie aktywacji emocjonalnej opisane w Supplementary File 1, Section S1.2. Implementacja równania aktywacji emocjonalnej w programie MATLAB jest dostępna w Supplementary File 1, Section S1.8 (Function 1: Model Equations).
4. Zaimplementuj równanie adaptacyjnej regulacji sterujące zmienną stanu pojemności regulacyjnej R(t) zgodnie z nieliniowym modelem opisanym w Supplementary File 1, Section S1.2. Odpowiednią implementację w programie MATLAB podano w Supplementary File 1, Section S1.8 (Function 1: Model Equations).
5. Włącz opcjonalne perturbacje stochastyczne poprzez wprowadzenie niezależnych składników szumu gaussowskiego do każdej zmiennej stanu w celu symulacji losowych zakłóceń akademickich i zmienności środowiskowej.
6. Przeprowadź 30 niezależnych realizacji stochastycznych (n = 30), używając różnych sekwencji szumu losowego przy zachowaniu identycznych bazowych parametrów modelu i warunków początkowych.
7. Oblicz średnie trajektorie i odpowiadające im oszacowania zmienności; słupki błędów przedstawione w sekcji Wyniki reprezentują wariację statystyczną pomiędzy tymi niezależnymi realizacjami.
8. Zintegruj nieliniowy system dynamiczny za pomocą numerycznego solvera ODE45 z tolerancją względną 1 × 10⁻⁶, tolerancją bezwzględną 1 × 10⁻⁸, początkowym krokiem całkowania 0.001 i maksymalną wielkością kroku 0.05 w przedziale symulacji 0 ≤ t ≤ 500.
9. Szczegółową implementację obliczeniową, konfigurację solvera i schemat symulacji podano w Supplementary File 1, Sections S1.6–S1.11.
10. Zdefiniuj amplitudy stochastyczne w następujący sposób:
    1. Ustaw σS = 0.05.
    2. Ustaw σE = 0.05.
    3. Ustaw σR = 0.03.
11. Wymuś, aby wszystkie zmienne stanu pozostawały nieujemne podczas całkowania numerycznego.
12. Podczas postprocessingu zastąp ujemne wartości numeryczne zerem.

3. Symulacja numeryczna

1. Skonfiguruj solver numeryczny.
   1. Wybierz ODE45 dla symulacji deterministycznych.
   2. Wybierz równoważną procedurę całkowania stochastycznego w przypadku przeprowadzania symulacji stochastycznych.
   3. Wyznacz numerycznie stan równowagi, korzystając z funkcji fsolve w programie MATLAB z początkowym szacunkiem stanu [1, 1, 1]. Rozwiązanie iteracyjne uznaje się za zbieżne, gdy solver spełnia domyślne kryteria zbieżności nieliniowej przy określonych tolerancjach numerycznych.
      ​Uwaga: Pełna implementacja znajduje się w Pliku uzupełniającym 1, sekcji S1.8 (Solver równowagi).
2. Zdefiniuj ustawienia solvera numerycznego.
   1. Ustaw tolerancję względną na 10⁻⁶. Ustaw tolerancję bezwzględną na 10⁻⁸.
   2. Ustaw maksymalny rozmiar kroku na 0.05. Ustaw początkowy rozmiar kroku na 0.001.
3. Zdefiniuj przedział symulacji od t = 0 do t = 500.
4. Próbkuj trajektorie wyjściowe w odstępach Δt = 0.1.
5. Wykonaj całkowanie numeryczne.
6. Zapisz symulowane trajektorie S(t), E(t) i R(t).
7. Wyeksportuj wyniki symulacji w formacie wartości rozdzielanych przecinkami (CSV) do dalszej analizy.
   PUNKT PRZERWANIA: Na tym etapie wszystkie trajektorie symulacji i wyniki numeryczne zostały wygenerowane i wyeksportowane. Prace można przerwać i wznowić później, ponownie ładując zarchiwizowane pliki parametrów i wyniki symulacji przed przejściem do analizy równowagi, analizy stabilności, wizualizacji i analizy wrażliwości.

4. Analiza równowagi i stabilności

1. Oblicz rozwiązania stanu równowagi, przyjmując, że wszystkie pochodne stanu są równe zero. Rozwiąż wynikające z tego nieliniowe równania algebraiczne, korzystając z numerycznego algorytmu wyznaczania pierwiastków.
   UWAGA: Przyjmuje się, że rozwiązanie stanu ustalonego zostało osiągnięte, gdy różnica bezwzględna między dwoma kolejnymi wektorami stanu spełnia określone kryterium zbieżności. Szczegółowe kryterium zbieżności znajduje się w Pliku uzupełniającym 1, sekcji S1.7.
2. Oblicz macierz jakobianu dla każdego punktu równowagi. Wyznacz wartości własne macierzy jakobianu.
3. Klasyfikuj stany równowagi w zależności od znaku maksymalnej rzeczywistej wartości własnej. Oznacz punkty równowagi jako stabilne, gdy wszystkie wartości własne mają ujemne części rzeczywiste.
4. Wygeneruj mapy stabilności przedstawiające rozkład stabilnych i niestabilnych obszarów pracy.

5. Wizualizacja w przestrzeni fazowej i dziedzinie czasu

1. Wygeneruj wykresy zależności S(t), E(t) i R(t) w dziedzinie czasu. Utwórz dwuwymiarowe portrety fazowe dla płaszczyzny S–E.
2. Utwórz dwuwymiarowe portrety fazowe dla płaszczyzny E–R. Wygeneruj trójwymiarowe trajektorie w przestrzeni stanów S–E–R.
3. Nałóż pola wektorowe i wykresy linii prądu na portrety fazowe, tam gdzie jest to możliwe. Zapisz wszystkie rysunki w formacie obrazu o wysokiej rozdzielczości do celów dokumentacji i publikacji.

6. Analiza wrażliwości

1. Stwórz skrypt do analizy wrażliwości. Zmieniaj jeden parametr na raz, utrzymując pozostałe parametry na ich wartościach bazowych. Przeprowadź analizę wrażliwości dla parametru sprzężenia stres–emocje.
2. Przeprowadź analizę wrażliwości dla parametru rozpraszania stresu. Przeprowadź analizę wrażliwości dla parametru wzmocnienia regulacji. Przeprowadź analizę wrażliwości dla parametru wymuszenia zewnętrznego.
3. Przeprowadź analizę wrażliwości dla nieliniowego parametru nasycenia emocjonalnego. Uruchom symulację dla każdej wartości parametru.
4. Oblicz stan stacjonarny równowagi emocjonalnej po każdej symulacji. Zapisz wszystkie wartości równowagi w bazie danych wyników.
5. Przeprowadź jednoparametrowe analizy wrażliwości, niezależnie zmieniając główne parametry modelu w predefined zakresach wymienionych w Supplementary File 1, Section S1.10. Konkretnie: zmieniaj γ od 0.1 do 2.0 w krokach co 0.1, β od 0.1 do 1.0 w krokach co 0.05, κ od 0.1 do 2.0 w krokach co 0.1, F od 0.5 do 5.0 w krokach co 0.25 oraz µ od 0.01 do 1.0 w krokach co 0.05.
   UWAGA: Dla każdej wartości parametru przeprowadź całkowanie numeryczne, oblicz rozwiązanie w stanie stacjonarnym, zapisz wynikowe zmienne stanu oraz wygeneruj odpowiadające im krzywe wrażliwości i wykresy porównawcze.
6. Wygeneruj krzywe wrażliwości. Wygeneruj porównawcze wykresy punktowe.
7. Wygeneruj porównawcze wykresy słupkowe. Wygeneruj mapy ciepła odpowiedzi parametrów.

7. Generowanie i weryfikacja wyników

1. Zweryfikuj zbieżność wszystkich symulacji. Potwierdź, że trajektorie stanów pozostają ograniczone przez cały przedział symulacji.
2. Zweryfikuj, czy rozwiązanie numeryczne osiągnęło stan ustalony przed przejściem do kolejnych analiz.
   UWAGA: Zbieżność przyjmuje się, gdy następujące po sobie zmiany wszystkich zmiennych stanu spełniają zdefiniowane kryterium zbieżności opisane w Pliku uzupełniającym 1, Sekcji S1.7 (Kryteria zbieżności). Do analiz równowagi, stabilności i czułości wykorzystuje się wyłącznie rozwiązania zbieżne.
3. Potwierdź, że rozwiązania równowagowe spełniają kryterium zbieżności. Eksportuj wszystkie wyniki symulacji, mapy stabilności, portrety fazowe, wykresy czułości oraz zbiory danych numerycznych.
   UWAGA: O ile nie określono inaczej, wszystkie symulacje wykorzystują bazowy zestaw parametrów przedstawiony w Pliku uzupełniającym 1, a wszystkie analizy numeryczne przebiegają zgodnie z przepływem obliczeniowym podsumowanym w Pliku uzupełniającym 1, Sekcji S1.9. Pełny kod źródłowy MATLAB, pliki parametrów, konfiguracje solverów oraz procedury wizualizacji niezbędne do powtórzenia opisanych symulacji znajdują się w Pliku uzupełniającym 1.
4. Zarchiwizuj wszystkie skrypty, pliki parametrów i pliki wyjściowe w celu zapewnienia odtwarzalności i przyszłych badań walidacyjnych.
   UWAGA: Dostarcz cały kod źródłowy MATLAB, pliki parametrów, procedury wykresowania i przepływy symulacji jako materiały uzupełniające, aby umożliwić niezależne powtórzenie przedstawionych wyników.

Wyniki

W celu walidacji proponowanego nieliniowego modelu dynamicznego oraz wykazania jego przewag nad konwencjonalnymi sformułowaniami liniowymi przeprowadzono kompleksowe badanie symulacji numerycznych (Supplementary File 2). O ile nie zaznaczono inaczej, wszystkie wyniki symulacji stochastycznych przedstawiono jako średnią ± odchylenie standardowe (SD) uzyskane z 30 niezależnych realizacji (n = 30). Porównania statystyczne między scenariuszami symulacji przeprowadzono przy użyciu odpowiedniego testu statystycznego z poziomem istotności p < 0.05. W przypadkach, w których przeprowadzono porównania statystyczne, podano dokładne wartości p. Symulacje oceniały zachowanie nieliniowego modelu stres–emocje–regulacja przy systematycznej zmianie parametrów oraz porównywały jego odpowiedzi strukturalne z odpowiedziami tradycyjnego liniowego modelu stresu. Całkowanie numeryczne rządzących równań różniczkowych przeprowadzono w wystarczająco długich horyzontach czasowych, aby zapewnić zbieżność do stanów stacjonarnych lub reżimów asymptotycznych. W każdym eksperymencie do obu modeli zastosowano identyczne warunki bazowe, zmieniając jeden parametr w predefinedowanym zakresie. Wynikowe stacjonarne stany emocjonalne, dynamika przejściowa oraz charakterystyka stabilności zostały zarejestrowane i zwizualizowane, aby uwypuklić różnice strukturalne w zachowaniu systemu. Prawidłową implementację protokołu potwierdzono w momencie, gdy system nieliniowy zbiegał do ograniczonych rozwiązań stanu stacjonarnego, wykazywał stabilne struktury atraktorów w przestrzeni fazowej oraz utrzymywał ujemne dominujące wartości własne macierzy Jakobiego w bazowych warunkach parametrów.

Analiza porównawcza skupiła się na kilku kluczowych parametrach dynamicznych, w tym na sile sprzężenia stres–emocje (γ), szybkości rozpraszania stresu (β), amplitudzie wymuszenia zewnętrznego (F), wzmocnieniu regulacji (κ) oraz współczynniku nieliniowości emocjonalnej (µ). Parametry te bezpośrednio wpływają na intensywność sprzężenia zwrotnego, marginesy stabilności oraz redystrybucję energii w modelowanym środowisku akademickim. Symulacje ujawniają, czy odpowiedzi systemu wykazują zachowanie proporcjonalnego skalowania, zgodnie z przewidywaniami modelu liniowego, czy też przejawiają zjawiska nieliniowe, takie jak nasycenie, buforowanie odporności oraz wielostabilność, zgodnie z przewidywaniami proponowanej formulacji nieliniowej.

Wrażliwość równowagi emocjonalnej na siłę sprzężenia stres-emocje przedstawiono na Rysunku 4. W przypadku zwiększenia parametru sprzężenia γ model liniowy generuje niemal stałe reakcje emocjonalne, co wskazuje, że intensywność sprzężenia nie wpływa strukturalnie na wyniki równowagi. Z kolei sformułowanie nieliniowe wykazuje spadek stałej aktywacji emocjonalnej wraz ze wzrostem γ, co odzwierciedla wpływ nieliniowych składników interakcji, które dynamicznie regulują sprzężenie zwrotne stres-emocje.

Rysunek 4: Porównawcza analiza czułości w zależności od siły sprzężenia stres–emocje (γ). Porównanie zgrupowane słupkowo odpowiedzi emocjonalnej w stanie ustalonym wraz ze wzrostem siły sprzężenia stres–emocje γ. Niebieskie słupki reprezentują model liniowy, a pomarańczowe słupki model nieliniowy. Oś pozioma przedstawia siłę sprzężenia γ, a oś pionowa przedstawia ustaloną równowagę emocjonalną. Słupki reprezentują średnie wartości emocji w stanie ustalonym uzyskane z powtórzonych serii symulacji, natomiast słupki błędów wskazują ± odchylenie standardowe (SD) od średniej. Wyniki wykazują, że model nieliniowy wykazuje zachowanie nasycenia i ograniczone odpowiedzi emocjonalne wraz ze wzrostem siły sprzężenia, podczas gdy model liniowy pozostaje stosunkowo niewrażliwy na zmiany intensywności sprzężenia. Słupki reprezentują średnie wartości ustalonej równowagi emocjonalnej uzyskane z 30 niezależnych serii symulacji (n = 30), a słupki błędów wskazują ± odchylenie standardowe. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Wpływ szybkości rozpraszania stresu β na równowagę emocjonalną przedstawiono na Rysunku 5. Model liniowy przewiduje gwałtowny spadek aktywacji emocjonalnej wraz ze wzrostem rozpraszania, co wykazuje silną wrażliwość parametrów. W przeciwieństwie do niego model nieliniowy pozostaje stosunkowo stabilny w tym samym zakresie parametrów dzięki wewnętrznemu tłumieniu regulacyjnemu oraz nieliniowym mechanizmom sprzężenia zwrotnego.

Rysunek 5: Porównawcza analiza obszarów w zależności od szybkości rozpraszania stresu (β). Porównanie oparte na obszarach dla stanu stacjonarnej równowagi emocjonalnej przy zmiennej szybkości rozpraszania stresu β. Niebieski obszar zacieniowany reprezentuje model liniowy, a pomarańczowy obszar zacieniowany reprezentuje model nieliniowy. Oś pozioma przedstawia β, a oś pionowa przedstawia stacjonarną równowagę emocjonalną. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Zależność między aktywacją emocjonalną a zewnętrznym wymuszaniem akademickim przedstawiono na Rysunku 6. Sformułowanie liniowe wykazuje proporcjonalny wzrost aktywacji emocjonalnej wraz ze zwiększaniem amplitudy wymuszania. Model nieliniowy wykazuje natomiast odpowiedź nasycającą, w której aktywacja emocjonalna początkowo rośnie, a następnie stopniowo stabilizuje się ze względu na nieliniowe tłumienie i adaptacyjne efekty regulacyjne. W podoptymalnych warunkach parametrów, takich jak nadmierne sprzężenie stres–emocje lub niewystarczające wzmocnienie regulacyjne, system wykazywał niestabilne trajektorie, zwiększone oscylacje lub utratę stabilności równowagi, co wskazuje na obniżoną stabilność systemu i może reprezentować warunki związane z podwyższonym napięciem psychologicznym oraz zwiększoną podatnością na przejścia typu wypaleniowego.

Rysunek 6: Porównanie odpowiedzi równowagi emocjonalnej pod wpływem wymuszenia zewnętrznego (F). Stacjonarna odpowiedź emocjonalna jako funkcja amplitudy wymuszenia zewnętrznego F. Czerwona przerywana krzywa reprezentuje model liniowy, natomiast niebieska ciągła krzywa reprezentuje model nieliniowy. Oś pozioma przedstawia amplitudę wymuszenia F, a oś pionowa stacjonarną równowagę emocjonalną. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Odporność układu nieliniowego na perturbacje stochastyczne przedstawiono na Rysunku 7, który obrazuje ewolucję czasową aktywacji emocjonalnej i energii stresu w warunkach losowych zakłóceń. Obie zmienne fluktuują w ograniczonych zakresach pomimo ciągłego wprowadzania szumu, co wskazuje, że nieliniowe mechanizmy sprzężenia zwrotnego utrzymują stabilność w obliczu zmienności środowiskowej.

Rycina 7: Stochastyczna dynamika stresu i emocji pod wpływem perturbacji indukowanych szumem. Ewolucja czasowa aktywacji emocjonalnej E(t) (krzywa niebieska) oraz energii stresu S(t) (krzywa pomarańczowa) pod wpływem perturbacji stochastycznych. Oś pozioma reprezentuje kroki czasowe symulacji. Obie zmienne pozostają ograniczone pomimo ciągłych zaburzeń wywołanych szumem. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Analiza krajobrazu energetycznego dostarcza dodatkowych informacji na temat stabilności systemu. Rycina 8 przedstawia reprezentację konturową funkcji energii opartej na metodzie Ljapunowa w płaszczyźnie fazowej regulacji emocji, w której trajektoria systemu zbiega się w stronę stabilnego basenu przyciągania. Trójwymiarowa reprezentacja tej struktury energetycznej została pokazana na Rycini 9, ujawniając wiele potencjalnych studni, co sugeruje możliwość wystąpienia alternatywnych stabilnych stanów regulacji emocjonalnej.

Rysunek 8: Mapa konturowa krajobrazu energetycznego z trajektorią dynamiczną w płaszczyźnie fazowej regulacji emocji. Przedstawienie konturowe krajobrazu energetycznego opartego na funkcji Lapunowa w płaszczyźnie fazowej aktywacji emocjonalnej (E) i zdolności regulacyjnej (R). Czarna trajektoria ilustruje ewolucję systemu w kierunku stabilnego basenu atrakcji. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 9: Trójwymiarowy nieliniowy krajobraz energetyczny ilustrujący potencjał dla alternatywnych stanów regulacji emocjonalnej. Trójwymiarowa reprezentacja krajobrazu energetycznego Ljapunowa w przestrzeni faz regulacji emocji (E–R). Wiele studni potencjału sugeruje możliwość wystąpienia alternatywnych stabilnych stanów regulacji emocjonalnej w różnych warunkach systemowych. Krajobraz energetyczny zapewnia jakościową wizualizację struktury stabilności systemu; jednak bezpośrednie potwierdzenie wielostabilności wymaga dodatkowych dowodów dynamicznych, takich jak przełączanie trajektorii lub analiza bifurkacji. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Właściwości stabilności systemu nieliniowego zostały poddane dalszej analizie przy użyciu metod opartych na wartościach własnych. Rycina 10 przedstawia mapę ciepła maksymalnej rzeczywistej wartości własnej macierzy Jakobiego dla różnych wartości siły sprzężenia stres–emocje (γ) oraz wzmocnienia regulacji (κ). Zwiększenie siły regulacji prowadzi do uzyskania bardziej ujemnych wartości własnych, co wskazuje na silniejszą stabilność asymptotyczną, podczas gdy nadmierne sprzężenie może obniżyć stabilność, jeśli nie zostanie zrównoważone przez wystarczającą regulację. Obszary charakteryzujące się ujemnymi maksymalnymi rzeczywistymi wartościami własnymi odpowiadają stabilnym warunkom pracy, natomiast obszary zbliżające się do zera lub przekraczające go wskazują progi niestabilności i potencjalne przejścia między reżimami.

Rycina 10: Mapa ciepła stabilności nieliniowego układu stres–emocje na podstawie maksymalnej rzeczywistej wartości własnej. Mapa ciepła przedstawiająca maksymalną część rzeczywistą wartości własnych macierzy Jakobiego w zależności od siły sprzężenia stres–emocje γ oraz wzmocnienia regulacji κ. Bardziej ujemne wartości wskazują na silniejszą stabilność asymptotyczną, natomiast wartości zbliżające się do zera wskazują na zmniejszoną stabilność i większe prawdopodobieństwo niestabilności. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Dynamika przestrzeni fazowej została przedstawiona na Rysunku 11, który obrazuje pole wektorowe oraz linie prądu systemu stres–emocje w płaszczyźnie fazowej S–E. Trajektorie zbiegają się w stronę stabilnego obszaru równowagi, co wykazuje zachowanie atraktora. Taka zbieżność potwierdza prawidłową implementację protokołu i dowodzi, że nieliniowe ramy teoretyczne spójnie odtwarzają stabilną dynamikę regulacji stres–emocje w określonych warunkach symulacji. Globalna struktura stabilności systemu została dodatkowo zilustrowana na Rysunku 12, gdzie trajektorie rozpoczynające się z wielu różnych warunków początkowych zbiegają się w stronę wspólnego atraktora w płaszczyźnie stres–emocje.

Rysunek 11: Reprezentacja pola wektorowego na płaszczyźnie fazowej oraz linii prądu dynamiki stres–emocje. Reprezentacja pola wektorowego i linii prądu nieliniowego układu stres–emocje na płaszczyźnie fazowej energii stresu (S) i aktywacji emocjonalnej (E). Linie prądu zbiegają się w stronę stabilnego obszaru równowagi, co wskazuje na zachowanie atraktora. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 12: Gęsty portret fazowy dla wielu warunków początkowych w płaszczyźnie stres–emocje. Portret fazowy wygenerowany z wielu warunków początkowych w płaszczyźnie fazowej S–E. Trajektorie zbiegają się do wspólnego atraktora, co wykazuje silną stabilność w różnych stanach początkowych. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Trójwymiarowa struktura atraktora układu nieliniowego została przedstawiona na Rysunku 13, gdzie trajektorie ewoluują w połączonej przestrzeni stanów stres–emocje–regulacja i dążą do stabilnego atraktora. Dynamika w dziedzinie czasu sprzężonych zmiennych została zilustrowana na Rysunku 14, gdzie aktywacja emocjonalna i energia stresu wykazują przejściowe dostrojenia przed zbieżnością do wartości stanu równowagi.

Rysunek 13: Trójwymiarowy nieliniowy atraktor w przestrzeni stanów stres–emocje–regulacja. Trójwymiarowa trajektoria układu nieliniowego w przestrzeni stanów zdefiniowanej przez stres (S), aktywację emocjonalną (E) i zdolność regulacyjną (R). Trajektoria zbiega się do stabilnego atraktora reprezentującego długoterminową równowagę układu. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 14: Ewolucja w dziedzinie czasu sprzężonych stanów stresu i emocji pod wpływem regulacji nieliniowej. Ewolucja czasowa aktywacji emocjonalnej E(t) (oś lewa) oraz energii stresu S(t) (oś prawa). Obie zmienne wykazują przejściowe dostrojenie, po którym następuje zbieżność do stanu równowagi ustalonej. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Wpływ siły regulacji adaptacyjnej na równowagę emocjonalną został przedstawiony na Rysunku 15. Wraz ze wzrostem wzmocnienia regulacyjnego κ, model liniowy przewiduje znaczne zmniejszenie pobudzenia emocjonalnego, podczas gdy model nieliniowy utrzymuje niemal stałe wartości równowagi dzięki mechanizmom saturacji adaptacyjnej.

Rysunek 15: Dwuosiowa analiza porównawcza przy zmiennym wzmocnieniu regulacji (κ). Porównanie stanu ustalonego równowagi emocjonalnej przy zmiennym wzmocnieniu regulacji κ. Krzywa niebieska reprezentuje model liniowy, a czerwona krzywa przerywana reprezentuje model nieliniowy. Oś pozioma przedstawia wzmocnienie regulacji κ. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Na koniec Rycina 16 przedstawia porównanie czułości w formie wykresu rozrzutu w zależności od zmienności współczynnika nieliniowości emocjonalnej (µ). W modelu liniowym równowaga emocjonalna pozostaje niezmieniona ze względu na brak członów nieliniowych. W przeciwieństwie do niego, model nieliniowy wykazuje spadek równowagi emocjonalnej wraz ze wzrostem µ, co demonstruje stabilizujący wpływ saturacji sześciennej na dynamikę emocjonalną.

Rysunek 16: Porównanie czułości oparte na rozproszeniu w zależności od parametru nieliniowości emocjonalnej (µ). Porównanie rozproszenia stanu ustalonego równowagi emocjonalnej przy zmiennym współczynniku nieliniowości emocjonalnej µ. Czerwone znaczniki reprezentują model liniowy, a niebieskie znaczniki reprezentują model nieliniowy. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Poza znaczeniem obliczeniowym, zaobserwowane zachowania dynamiczne mają istotne interpretacje w kontekstach edukacyjnych i psychologicznych. Ograniczone odpowiedzi obserwowane przy perturbacjach stochastycznych sugerują, że adaptacyjne mechanizmy regulacyjne mogą amortyzować skutki nieoczekiwanych stresorów akademickich, wspierając tym samym odporność i stabilność emocjonalną. Podobnie, istnienie stabilnych atraktorów można interpretować jako psychicznie zrównoważone stany, w których studenci skutecznie regulują naciski akademickie, podczas gdy obszary niestabilności i progi bifurkacji mogą odpowiadać warunkom, w których zasoby radzenia sobie stają się niewystarczające, co zwiększa podatność na wypalenie, wyczerpanie emocjonalne lub nieadaptacyjne reakcje stresowe. Analizy wrażliwości dodatkowo wskazują, że wzmocnienie zdolności regulacyjnych może rozszerzyć obszary stabilności i zmniejszyć podatność na destabilizujące przejścia, co podkreśla potencjalne implikacje dla interwencji mających na celu poprawę umiejętności radzenia sobie, regulacji emocjonalnej i dobrostanu studentów. Te wyniki oparte na symulacjach sugerują, że proponowane ramy mogą służyć jako użyteczna podstawa koncepcyjna dla przyszłych badań nad trajektoriami stresu, mechanizmami odporności i potencjalnymi strategiami interwencyjnymi. Niemniej jednak, przed wdrożeniem praktycznym wymagana jest walidacja z wykorzystaniem empirycznych danych studenckich.

Wspólnie wyniki te wykazują, że proponowane ramy nieliniowe oddają kilka właściwości strukturalnych, których brakuje w tradycyjnych liniowych modelach stresu, w tym zachowanie nasycenia, buforowanie stabilności podczas perturbacji, wielostabilne krajobrazy energetyczne oraz odporność dzięki regulacji adaptacyjnej. W testowanych zakresach parametrów model nieliniowy konsekwentnie utrzymywał ograniczoną aktywację emocjonalną i stabilne zachowanie atraktora, podczas gdy model liniowy wykazywał znacznie większą wrażliwość na zmienność parametrów i zredukowane marginesy stabilności. Właściwości te zapewniają bardziej realistyczną reprezentację dynamiki stres–emocje w złożonych środowiskach akademickich.

Plik uzupełniający 1: Kod źródłowy MATLAB, równania rządzące, implementacja numeryczna i dokumentacja odtwarzalności.Ten plik uzupełniający zawiera rządzące nieliniowe równania różniczkowe, bazowe parametry modelu, warunki początkowe, konfigurację solvera numerycznego, kryteria zbieżności, kod źródłowy MATLAB, solver stanu równowagi, procedury analizy wrażliwości, rutyny wizualizacji, specyfikacje oprogramowania, zakresy parametrów oraz schemat obliczeniowy niezbędny do odtworzenia wszystkich symulacji, analiz stabilności, trajektorii w przestrzeni fazowej, analiz wrażliwości, map ciepła wartości własnych oraz krajobrazów energii Ljapunowa przedstawionych w manuskrypcie.Prosimy kliknąć tutaj, aby pobrać ten plik.

Plik uzupełniający 2: Dokumentacja przepływu pracy symulacji, mapowanie zmiennych psychologicznych i ramy walidacji empirycznej.Ten dodatek zawiera kompletny przepływ pracy symulacji; mapowanie matematycznych zmiennych stanu i parametrów modelu na mierzalne konstrukty psychologiczne; zalecane instrumenty oceny psychologicznej; wytyczne dotyczące szacowania parametrów; proponowane ramy dla przyszłej walidacji empirycznej i kalibracji z wykorzystaniem danych uczniów; oraz potencjalne strategie wdrożenia dla walidacji podłużnej i zastosowań edukacyjnych.Prosimy kliknąć tutaj, aby pobrać ten plik.

Dyskusja

Przedstawiony tutaj protokół obliczeniowy demonstruje, w jaki sposób można zaimplementować nieliniowe ramy dynamiczne w celu symulacji interakcji między stresem a emocjami poprzez całkowanie numeryczne, analizę stabilności oraz analizę wrażliwości. Krytycznym krokiem w implementacji tej metody jest prawidłowe całkowanie numeryczne rządzących równań różniczkowych w wystarczająco długich odstępach czasu, aby zapewnić zbieżność do stanu stacjonarnego lub dynamiki atraktora. Analizy wrażliwości wykazują, że sformułowanie nieliniowe reaguje na zmiany parametrów inaczej niż tradycyjny model liniowy. W szczególności siła sprzężenia (γ) bezpośrednio wpływa na równowagę emocjonalną w układzie nieliniowym, podczas gdy szybkość rozpraszania stresu (β) wywołuje przesadzone reakcje głównie w sformułowaniu liniowym. Podobnie, zwiększenie wymuszenia zewnętrznego prowadzi do nieograniczonego wzmocnienia w modelu liniowym, ale do ograniczonego nasycenia w ramach nieliniowych. Obserwacje te podkreślają znaczenie starannej konfiguracji parametrów podczas stosowania protokołu w celu zachowania realistycznej dynamiki układu³¹.

Protokół ten można dostosować do symulacji różnych warunków środowiskowych i scenariuszy behawioralnych. Na przykład do struktury symulacji można włączyć perturbacje stochastyczne, reprezentujące nieoczekiwane zdarzenia akademickie lub stresory społeczne, aby zbadać odporność systemu na losowe zaburzenia. W takich warunkach układ nieliniowy wykazuje ograniczone fluktuacje, a nie niekontrolowaną dywergencję, co wskazuje na odporność na perturbacje wywołane szumem. Dodatkowo analiza krajobrazu energetycznego dostarcza geometrycznego wglądu w stabilność systemu. Wynikowe struktury energetyczne ilustrują, w jaki sposób trajektorie systemu zbiegają się w stronę stabilnych basenów atrakcji oraz jak wiele potencjalnych studni może reprezentować odrębne stany regulacji emocjonalnej. Te narzędzia analityczne stanowią cenne wsparcie w diagnozowaniu niestabilnych symulacji i identyfikowaniu zakresów parametrów, które wspierają stabilną dynamikę regulacyjną. Do powszechnych przyczyn niestabilności symulacji należą nadmiernie duże parametry sprzężenia, niewystarczające wzmocnienie regulacyjne lub niewłaściwe ustawienia całkowania numerycznego. Problemy te można często rozwiązać poprzez rekalibrację parametrów i weryfikację stabilności solvera numerycznego.

Właściwości stabilności układu nieliniowego można również zbadać za pomocą analizy opartej na wartościach własnych oraz technik wizualizacji przestrzeni fazowej. Analizy te wykazują, w jaki sposób oddziaływanie między sprzężeniem stres–emocje a wzmocnieniem regulacyjnym determinuje ogólną stabilność układu. Zwiększenie siły regulacji zazwyczaj zwiększa stabilność asymptotyczną, podczas gdy nadmiernie silne sprzężenie może zdestabilizować układ, jeśli regulacja jest niewystarczająca. Reprezentacje w przestrzeni fazowej dodatkowo ujawniają, że trajektorie wychodzące z różnych warunków początkowych mają tendencję do zbiegania się w stronę wspólnego atraktora, co wskazuje na wysoką stabilność w ramach modelu nieliniowego. Skoordynowana ewolucja stresu, aktywacji emocjonalnej i zdolności regulacyjnych obrazuje, w jaki sposób adaptacyjne mechanizmy sprzężenia zwrotnego prowadzą układ do stanów równowagi³².

Analiza w dziedzinie czasu dostarcza dodatkowych informacji na temat przejściowych zachowań systemu oraz formowania się stanu równowagi. Zarówno aktywacja emocjonalna, jak i skumulowany stres zazwyczaj wykazują początkowe fluktuacje przejściowe, po których następuje stopniowa zbieżność do poziomów równowagi stanu ustalonego. Analiza wrażliwości parametrów regulacji adaptacyjnej wykazuje, że zwiększona siła regulacji silnie tłumi aktywację emocjonalną w modelach liniowych, ale powoduje bardziej umiarkowane zmiany w sformułowaniu nieliniowym ze względu na mechanizmy nasycenia i buforowania sprzężeń zwrotnych. Podobnie, uwzględnienie nieliniowych członów nasycenia emocjonalnego redukuje nadmierną amplifikację i stabilizuje odpowiedzi emocjonalne. Wspólnie wyniki te dowodzą, że ramy nieliniowe oddają dynamikę regulacji i mechanizmy stabilizacji, które nie są reprezentowane w konwencjonalnych liniowych modelach stresu. W porównaniu z konwencjonalnymi modelami stresu opartymi na liniowości lub korelacji, proponowane ramy umożliwiają analizę granic stabilności, dynamiki atraktorów oraz przejść nieliniowych, których nie udaje się łatwo uchwycić za pomocą opisowych podejść statystycznych.

Przyszłe badania powinny przeprowadzić analizę porównawczą proponowanego modelu z uznanymi podejściami, w tym z transakcyjnymi modelami stresu, sieciowymi modelami psychologicznymi, podłużnymi modelami latentnego wzrostu oraz ramami predykcyjnymi opartymi na uczeniu maszynowym. Takie porównania pozwoliłyby na bardziej rygorystyczną ocenę zalet wyjaśniających i predykcyjnych oferowanych przez nieliniowe modelowanie dynamiczne. Poza wkładem metodologicznym, proponowany model ma kilka szerszych implikacji dla psychologii edukacyjnej i praktyki szkolnictwa wyższego. Z perspektywy doradztwa, identyfikacja obszarów stabilności, mechanizmów odporności oraz progów niestabilności może pomóc w konceptualizacji sposobu, w jaki studenci reagują na przedłużającą się presję akademicką oraz dlaczego niektóre osoby skutecznie regenerują się po stresie, podczas gdy inne doświadczają narastających trudności emocjonalnych. Model sugeruje, że odporność psychiczna nie jest jedynie statyczną cechą osobniczą, lecz może wynikać z dynamicznych interakcji między ekspozycją na stres, reakcjami emocjonalnymi a adaptacyjnymi procesami regulacyjnymi³³.

Opracowane ramy mogą również stanowić fundament teoretyczny dla przyszłych systemów wsparcia studentów oraz inicjatyw monitorowania dobrostanu. Na przykład identyfikacja reżimów parametrów związanych ze zmniejszoną stabilnością mogłaby wesprzeć opracowanie wskaźników wczesnego ostrzegania przed podwyższonym ryzykiem wypalenia, podczas gdy analiza mechanizmów sprzężenia zwrotnego regulacji mogłaby pomóc w opracowaniu interwencji mających na celu wzmocnienie umiejętności radzenia sobie oraz emocjonalnej samoregulacji. Na poziomie instytucjonalnym ramy te oferują perspektywę zorientowaną na systemy, która może pomóc w ocenie, w jaki sposób obciążenie akademickie, harmonogramy oceniania i zasoby wsparcia wspólnie wpływają na dobrostan studentów. Proponowane podejście ma na celu uzupełnienie istniejących teorii psychologicznych dotyczących stresu i regulacji, a nie ich zastąpienie. Tradycyjne modele, w tym transakcyjne teorie stresu, ramy oparte na radzeniu sobie, teorie regulacji emocjonalnej oraz modele odporności psychicznej, dostarczają cennych wyjaśnień procesów adaptacji poznawczej, behawioralnej i emocjonalnej. Nieliniowe ramy dynamiczne rozszerzają te perspektywy, wprowadzając ilościową reprezentację tego, jak procesy te oddziałują na siebie w czasie, potencjalnie generując stabilność, zachowania oscylacyjne, odporność lub gwałtowne przejścia w kierunku wyników nieadaptacyjnych. Taka integracja pozwala osadzić opracowane ramy w szerszej literaturze dotyczącej stresu u studentów i adaptacji psychologicznej, podkreślając jednocześnie ich potencjalną wartość jako pomostu między psychologią edukacyjną a nauką o systemach dynamicznych.

Należy również wziąć pod uwagę szerszą wartość edukacyjną ram opartych na symulacjach. Choć niniejsze badanie koncentruje się na nieliniowej dynamice stres–emocje, zastosowane podejście modelowe może wspierać refleksyjne uczenie się i adaptacyjne podejmowanie decyzji, umożliwiając studentom, pedagogom i specjalistom wspierającym badanie konsekwencji różnych ścieżek regulacji stresu w symulowanych środowiskach akademickich. Takie symulacje mogą sprawić, że inaczej abstrakcyjne procesy psychologiczne staną się bardziej widoczne i zrozumiałe, co w konsekwencji sprzyja większej samoświadomości i refleksji nad strategiami radzenia sobie, reakcjami emocjonalnymi oraz rozwojem odporności psychicznej.

Perspektywa ta jest zgodna z aktualnymi badaniami wykazującymi korzyści edukacyjne płynące z doświadczeń opartych na symulacjach. Na przykład Shi et al. (2025)³⁴ wykazali, że symulowane doświadczenia z transportu publicznego, zaprojektowane w celu odwzorowania wyzwań, z jakimi mierzą się osoby z lekkimi zaburzeniami widzenia, sprzyjały rozwijaniu empatii, myślenia refleksyjnego oraz poprawiały proces podejmowania decyzji projektowych. Podobne zasady mogą mieć zastosowanie w regulacji stresu w szkolnictwie wyższym, gdzie symulacyjne reprezentacje kumulacji stresu, adaptacji emocjonalnej i dynamiki odporności psychicznej mogłyby pomóc studentom lepiej zrozumieć ich reakcje behawioralne, wspierając jednocześnie wykładowców i doradców w projektowaniu skuteczniejszych interwencji. W konsekwencji proponowane ramy mogą służyć nie tylko jako narzędzie analityczne do badania dynamiki stresu, ale także jako podstawa dla przyszłych symulacji edukacyjnych, które wspierają empatię, autorefleksję i adaptacyjne podejmowanie decyzji w środowisku akademickim.

pomimo tych zalet, proponowany protokół posiada kilka ograniczeń. Model przedstawia regulację emocji za pomocą uproszczonego zestawu zmiennych dynamicznych opisujących stres, aktywację emocjonalną i zdolność regulacyjną. Rzeczywista dynamika emocjonalna zależy od szerszego zakresu czynników, w tym procesów oceny poznawczej, interakcji społecznych i indywidualnych strategii radzenia sobie. W związku z tym obecne ramy należy interpretować przede wszystkim jako koncepcyjne i obliczeniowe narzędzie modelowania, a nie jako model predykcyjny dla indywidualnych zachowań psychologicznych. Niemniej jednak, podejście to stanowi wartościową platformę metodologiczną do badania dynamiki odporności psychicznej, warunków stabilności i nieliniowych mechanizmów sprzężenia zwrotnego w systemach stresu akademickiego. Przyszłe badania mogą rozszerzyć te ramy poprzez włączenie empirycznych zbiorów danych, spersonalizowaną estymację parametrów oraz metody kalibracji oparte na danych, aby poprawić zdolność predykcyjną i przydatność w rzeczywistych warunkach. Przyszłe badanie walidacyjne mogłoby zostać przeprowadzone poprzez gromadzenie podłużnych danych studenckich przy użyciu zwalidowanych instrumentów psychologicznych, takich jak Perceived Stress Scale (PSS), Maslach Burnout Inventory–Student Survey (MBI-SS), Connor–Davidson Resilience Scale (CD-RISC) oraz miary dobrostanu emocjonalnego. Obserwacje szeregów czasowych uzyskane w ciągu semestru akademickiego mogłyby posłużyć do estymacji parametrów modelu i porównania przewidywanych trajektorii z obserwowanymi wynikami w zakresie stresu i emocji. Taka walidacja umożliwiłaby ilościową ocenę dokładności modelu i jego wydajności predykcyjnej w rzeczywistych warunkach edukacyjnych³⁵.

Dodatkowym ograniczeniem jest brak walidacji empirycznej z wykorzystaniem rzeczywistych danych studenckich. Niniejsze badanie koncentruje się na opracowaniu teoretycznych i obliczeniowych ram modelu, a zatem opiera się przede wszystkim na analizie matematycznej i symulacjach numerycznych. Chociaż badania symulacyjne są cenne dla eksploracji zachowania systemu i identyfikacji potencjalnych mechanizmów, nie mogą one same w sobie potwierdzić, że modelowana dynamika dokładnie reprezentuje procesy regulacji stresu w rzeczywistych środowiskach edukacyjnych. Ponadto szacowanie parametrów modelu stanowi istotne wyzwanie, ponieważ zmienne takie jak akumulacja stresu, aktywacja emocjonalna i zdolność regulacyjna nie są bezpośrednio obserwowalne i mogą różnić się w zależności od jednostki, instytucji oraz kontekstu kulturowego. Dodatkowo obecny model zakłada uogólnioną populację studentów i nie uwzględnia w sposób jawny różnic demograficznych, kulturowych, instytucjonalnych lub dyscyplinarnych, które mogą wpływać na dynamikę regulacji stresu. Przyszłe badania powinny ocenić wydajność modelu w zróżnicowanych kontekstach edukacyjnych.

Ważnym kierunkiem przyszłych badań jest integracja zaproponowanego nieliniowego modelu stres–emocje–regulacja z rozwijającymi się imersyjnymi technologiami edukacyjnymi i zdrowia psychicznego. Niedawne prace nad wirtualnymi krajobrazami terapeutycznymi wykazały, że imersyjne środowiska cyfrowe mogą redukować lęk, poprawiać dobrostan emocjonalny i wspierać regenerację psychiczną studentów poprzez starannie zaprojektowane doświadczenia regeneracyjne. Wyniki te sugerują, że środowiska oparte na symulacjach mogą stanowić praktyczne platformy do przekształcania teoretycznych modeli regulacji stresu w aplikacje zorientowane na interwencję. W tym kontekście obszary stabilności, progi odporności, adaptacyjne mechanizmy regeneracji oraz stochastyczne wzorce odpowiedzi zidentyfikowane w ramach niniejszego modelu mogłyby potencjalnie wspomóc projektowanie responsywnych psychologicznie środowisk uczenia się, zdolnych do monitorowania, symulowania i wspierania dobrostanu studentów³⁶.

Przyszłe badania powinny zatem skoncentrować się na empirycznej kalibracji i walidacji proponowanych ram teoretycznych. Obiecującym podejściem byłoby gromadzenie danych podłużnych przy użyciu zwalidowanych instrumentów psychologicznych mierzących postrzegany stres, dobrostan emocjonalny, skuteczność radzenia sobie, rezyliencję, objawy wypalenia oraz zaangażowanie akademickie. Takie zbiory danych mogłyby zostać wykorzystane do oszacowania parametrów modelu, oceny wydajności predykcyjnej, zidentyfikowania różnic indywidualnych w zachowaniach dynamicznych oraz oceny stopnia, w jakim proponowane mechanizmy stabilności i rezyliencji odpowiadają obserwowanym wynikom studentów. Ponadto przyszłe badania mogą porównać ten model z uznanymi podejściami psychologicznymi, statystycznymi i opartymi na uczeniu maszynowym, aby określić jego względną wartość wyjaśniającą i predykcyjną. Rozwiązania te stanowiłyby kluczowy krok w stronę przekształcenia obecnych ram teoretycznych w empirycznie zwalidowane narzędzie służące do rozumienia i wspierania dobrostanu studentów w szkolnictwie wyższym³⁷.

Materiały

Lista materiałów użytych w tym artykule

Nazwa	Firma	Numer katalogowy	Komentarze
Rozwiązywanie równań różniczkowych (ODE45)	MathWorks	MATLAB R2024a	Solver numeryczny użyty do całkowania sprzężonych nieliniowych równań różniczkowych. Tolerancja względna = 1×10-6; Tolerancja bezwzględna = 1×10-8.
Oprogramowanie MATLAB	MathWorks, Natick, MA, USA	RRID: SCR_001622	Podstawowe oprogramowanie do obliczeń numerycznych wykorzystane do implementacji modelu, analizy stabilności, analizy wrażliwości, wizualizacji przestrzeni fazowej, generowania krajobrazu energetycznego Ljapunowa oraz przygotowania rycin.
Pakiet dokumentacji metodologicznej	Plik uzupełniający S2	Nie dotyczy	Zawiera dokumentację ram teoretycznych, interpretację parametrów, strategię walidacji, ograniczenia oraz opisy przepływu obliczeń.
Biblioteka do przetwarzania tablic numerycznych (NumPy)	Deweloperzy NumPy	NumPy 2.0	Biblioteka wykorzystywana do operacji macierzowych, obliczeń numerycznych, przemiatań parametrów i przetwarzania danych w implementacji w języku Python.
Środowisko symulacji numerycznych	Środowisko MATLAB R2024a	N/A	Środowisko obliczeniowe wykorzystane do wykonania nieliniowej analizy naprężeń–emocje–ramy regulacyjne, analizy zakresowe parametrów, analiza stanu równowagi oraz symulacje stochastyczne.
Optimization Toolbox	MathWorks	MATLAB Optimization Toolbox	Stosowane do obliczeń stanu równowagi za pomocą nieliniowych metod poszukiwania pierwiastków (fsolve).
Środowisko programistyczne Python	Python Software Foundation	Python 3.12	Alternatywna platforma obliczeniowa wykorzystana do weryfikacji wyników numerycznych oraz odtworzenia wybranych symulacji i wizualizacji.
Biblioteka Obliczeń Naukowych (SciPy)	Społeczność SciPy	SciPy 1.14	Biblioteka wykorzystywana do obliczeń numerycznych, rozwiązywania równań nieliniowych oraz analizy równań różniczkowych w implementacji w języku Python.
Pakiet kodu źródłowego	Plik uzupełniający S1	Nie dotyczy	Zawiera skrypty MATLAB, konfiguracje solvera, pliki parametrów, rutyny analizy wrażliwości, kody analizy stabilności oraz przepływy pracy wizualizacji niezbędne do pełnej reprodukowalności.
Oprogramowanie arkusz kalkulacyjny	Microsoft Corporation	Microsoft Excel 365	Służy do przechowywania, organizacji, eksportu i inspekcji wyników symulacji oraz uzupełniających zbiorów danych.
Narzędzie do wizualizacji (Matplotlib)	Zespół programistyczny Matplotlib	Matplotlib 3.9	Biblioteka do wykresów naukowych wykorzystywana do generowania wykresów czułości, diagramów stabilności, trajektorii w przestrzeni fazowej, map ciepła oraz rycin o jakości publikacyjnej.