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ń psychologicznych i zewnętrznych presji akademickich¹. Program nauczania, egzaminy, oczekiwania instytucjonalne oraz ocena społeczna nieustannie zaburzają stany emocjonalne studentów, wywołując wahania poziomu lęku, motywacji i zaangażowania². Reakcje emocjonalne te ewoluują w wielowymiarowej przestrzeni fazowej, którą sterują ocena poznawcza, strategie radzenia sobie oraz mechanizmy regulacji adaptacyjnej³. Klasyczne modele odpowiedzi na stres 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 i 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 conceptualizować jako otwarte układy dyssypatywne, które nieustannie wymieniają związaną ze stresem „energię” z otoczeniem akademickim poprzez procesy sprzężenia zwrotnego sterujące regulacją i adaptacją⁸.

Z perspektywy edukacyjnej i psychologicznej pojęcia wykorzystane w niniejszych ramach 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 postrzeganego stresu; aktywację emocjonalną można oceniać przy użyciu zwalidowanych narzędzi badania lęku, nastroju lub dobrostanu emocjonalnego; a zdolność regulacyjną można szacować za pomocą mierników rezyliencji, radzenia sobie lub samoregulacji. Mapowania te stanowią praktyczną podstawę dla przyszłej estymacji parametrów i kalibracji modelu. W tym kontekście pojęcia nieliniowe, takie jak stabilność, rezyliencja i przejścia krytyczne, mogą być rozumiane jako odrębne wzorce adaptacji psychologicznej. Stany stabilne reprezentują zrównoważone funkcjonowanie emocjonalne, rezyliencja odzwierciedla zdolność do regeneracji po zaburzeniach, a przejścia krytyczne odpowiadają gwałtownym przesunięciom w stronę wyników maladaptywnych, 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ść zaproponowanych ram dla odbiorców interdyscyplinarnych.

Aby przedstawić koncepcyjny przegląd tych oddziaływań, Rysunek 1 ilustruje strukturę dynamicznych interakcji stres–emocje w systemach szkolnictwa wyższego. Schemat przedstawia wejścia stresowe wynikające z wymagań akademickich, zmienne odpowiedzi emocjonalnej reprezentujące ewolucję stanu psychologicznego oraz adaptacyjne mechanizmy regulacyjne, które modulują system poprzez pętle sprzężenia zwrotnego. Komponenty te oddziałują 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. Choć ta koncepcyjna reprezentacja podkreśla dynamiczny charakter regulacji emocjonalnej, istniejące modele psychologiczne rzadko przekładają te relacje na jawne struktury matematyczne opisujące ewolucję czasową.

Rysunek 1: Schemat koncepcyjny dynamicznej interakcji stres–emocje w systemach szkolnictwa wyższego. Rysunek ten ilustruje koncepcyjną strukturę dynamiki stres–emocje w środowisku akademickim. Schemat przedstawia bodźce stresowe 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ę tego rysunku.

Ograniczenia obecnych podejść do modelowania zostały dodatkowo zilustrowane na Ryc. 2, która 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 skupionych jest 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ń sterujących opisujących ewolucję czasową 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 w stanie ilościowo określić odporność lub geometrię 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 sterowane nieliniowymi procesami sprzężenia zwrotnego¹². Wspólnie luki te podkreślają potrzebę stworzenia rygorystycznych ram modelowania dynamicznego, zdolnych do opisania zachowań oscylacyjnych, multistabilności i przejść fazowych w systemach stresu w szkolnictwie wyższym.

Rysunek 2: Schemat koncepcyjny ilustrujący luki badawcze w istniejącej literaturze dotyczącej modelowania stresu. Rysunek ten przedstawia koncepcyjny przegląd głównych ograniczeń w obecnych podejściach do modelowania stresu. W centrum diagramu znajduje się brak mechanistycznych ram dynamicznych opisujących 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 proponowaną architekturę nieliniowego modelowania dynamicznego opracowaną w niniejszym badaniu. Struktura ta integruje dynamikę akumulacji stresu, zmienne aktywacji emocjonalnej oraz adaptacyjne mechanizmy regulacji sprzężenia zwrotnego w ramach ujednoliconego nieliniowego sformułowania przestrzeni stanów. W ramach tej struktury 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ę interakcji między stresem a emocjami 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.

Rycina 3: Proponowana architektura nieliniowego modelowania dynamicznego dla regulacji stresu i emocji. Rycina ta przedstawia opracowane w niniejszym badaniu nieliniowe ramy dynamiczne służące do modelowania interakcji między stresem a emocjami w systemach szkolnictwa wyższego. Architektura ta integruje dynamikę akumulacji stresu, zmienne aktywacji emocjonalnej oraz adaptacyjne mechanizmy sprzężenia zwrotnego w ramach zunifikowanego sformułowania przestrzeni stanów. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

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 systemu dynamicznego zdolnego do wykazania wielostabilności, reżimów oscylacyjnych oraz przejść krytycznych. W przeciwieństwie do konwencjonalnych liniowych modeli indeksu stresu, proponowane ramy jawnie uwzględniają nieliniowe interakcje sprzężeń zwrotnych pomiędzy skumulowaną energią stresu, aktywacją emocjonalną a adaptacyjną zdolnością regulacji^17,18,19,20. Interakcje nieliniowe te 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 ramy te wprowadzają interpretację regulacji emocjonalnej opartą na krajobrazie energetycznym, w której funkcje potencjału typu Ljapunowa opisują strukturę stabilności stanów emocjonalnych^21,22,23. W tej reprezentacji rezyliencję można określić ilościowo za pomocą głębokości basenu przyciągania, natomiast progi niestabilności odpowiadają punktom bifurkacji, w których zanikają stabilne atraktory.

Aby zwiększyć realizm, model uwzględnia okresowe wyrazy wymuszające reprezentujące cykle semestrów akademickich, a także perturbacje stochastyczne, które oddają nieoczekiwane zdarzenia akademickie lub fluktuacje środowiskowe^24,25,26. Rozszerzenia te pozwalają strukturze symulować realistyczną dynamikę stresu akademickiego, w tym oscylacyjne cykle stresu w okresach egzaminacyjnych oraz przejścia między stanami emocjonalnymi napędzane szumem. Takie możliwości modelowania zapewniają przewagę nad podejściami czysto statystycznymi, umożliwiając analizy predykcyjne progów niestabilności i mechanizmów odporności.

Nieliniowe podejścia dynamiczne były wcześniej stosowane w dziedzinach takich jak regulacja biologiczna, neuronauka i złożone systemy społeczne, w których procesy sterowane sprzężeniem zwrotnym wytwarzają 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 rygorystycznego matematycznie modelu. Dzięki połączeniu nieliniowego modelowania dynamicznego, teorii stabilności i analizy krajobrazu energetycznego, opisana tutaj metoda zapewnia zunifikowaną platformę do badania dynamiki stresu, mechanizmów odporności i stabilności emocjonalnej w środowiskach akademickich³⁰.

W ostateczności 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 kierunku 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 przy użyciu rygorystycznych narzędzi dynamiki. Integracja ta zwiększa jasność konceptualną, 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 środowisku akademickim³¹. 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 proponowanych nieliniowych ram dynamicznych oraz wykazania ich przewagi nad konwencjonalnymi sformułowaniami liniowymi przeprowadzono kompleksowe badanie symulacyjne (Supplementary File 2). O ile nie wskazano 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 miejscach, 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 zastosowano identyczne warunki bazowe dla obu modeli, podczas gdy pojedynczy parametr zmieniano w zdefiniowanym zakresie. Wynikowe stacjonarne stany emocjonalne, dynamika przejściowa oraz charakterystyka stabilności zostały zarejestrowane i zwizualizowane, aby uwydatnić różnice strukturalne w zachowaniu systemu. Prawidłowość wdrożenia protokołu potwierdzono, gdy system nieliniowy zbiegał do ograniczonych rozwiązań stanu stacjonarnego, wykazywał stabilne struktury atraktorów w przestrzeni fazowej i utrzymywał ujemne dominujące wartości własne macierzy Jacobiego 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 wykazują, czy odpowiedzi systemu wykazują proporcjonalne zachowanie skalowania, zgodnie z przewidywaniami modelu liniowego, czy też przejawiają zjawiska nieliniowe, takie jak nasycenie, buforowanie odporności i wielostabilność, zgodnie z przewidywaniami zaproponowanego sformułowania nieliniowego.

Wrażliwość równowagi emocjonalnej na siłę sprzężenia stres–emocje przedstawiono na Rysunku 4. Gdy parametr sprzężenia γ rośnie, model liniowy generuje niemal stałe reakcje emocjonalne, co wskazuje, że intensywność sprzężenia nie wpływa strukturalnie na wyniki równowagi. W przeciwieństwie do tego, sformułowanie nieliniowe wykazuje spadek stałej aktywacji emocjonalnej wraz ze wzrostem γ, co odzwierciedla wpływ nieliniowych członów oddziaływania, 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 (γ). Grupowe porównanie słupkowe 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 reprezentują model nieliniowy. Oś pozioma przedstawia siłę sprzężenia γ, a oś pionowa przedstawia stałą równowagę emocjonalną. Słupki reprezentują średnie wartości emocji w stanie ustalonym uzyskane z powtarzalnych przebiegów 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 natężenia sprzężenia. Słupki reprezentują średnie wartości równowagi emocjonalnej w stanie ustalonym uzyskane z 30 niezależnych przebiegów symulacji (n = 30), a słupki błędów wskazują ± odchylenie standardowe. Kliknij tutaj, aby wyświetlić powiększoną 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ą czułość parametrów. W przeciwieństwie do niego, model nieliniowy pozostaje względnie stabilny w tym samym zakresie parametrów ze względu na wewnętrzne tłumienie regulacyjne i nieliniowe mechanizmy sprzężenia zwrotnego.

Rysunek 5: Porównawcza analiza pola w funkcji szybkości rozpraszania stresu (β). Porównanie oparte na polu powierzchni dla stanu ustalonej 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 ustaloną równowagę emocjonalną. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Zależność między aktywacją emocjonalną a zewnętrznym wymuszaniem akademickim przedstawiono na Rysunku 6. Formuła liniowa 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 suboptymalnych warunkach parametrów, takich jak nadmierne sprzężenie stres–emocje lub niewystarczające wzmocnienie regulacyjne, system wykazywał niestabilne trajektorie, powię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 psychicznym oraz zwiększoną podatnością na przejścia typu wypalenie.

Rycina 6: Porównanie odpowiedzi równowagi emocjonalnej pod wpływem wymuszenia zewnętrznego (F). Stan stabilny odpowiedzi emocjonalnej w funkcji 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 stabilną równowagę emocjonalną. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Wytrzymałość 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.

Rysunek 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) w warunkach perturbacji stochastycznych. Oś pozioma przedstawia kroki czasowe symulacji. Obie zmienne pozostają ograniczone pomimo ciągłych zaburzeń wywołanych szumem. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Analiza krajobrazu energetycznego dostarcza dodatkowych informacji na temat stabilności systemu. Rysunek 8 przedstawia reprezentację konturową funkcji energii opartej na metodzie Ljapunowa w płaszczyźnie fazowej regulacji emocji, w której trajektoria systemu zbiega się do stabilnego basenu przyciągania. Trójwymiarowa reprezentacja tej struktury energetycznej została pokazana na Rysunku 9, ujawniając wiele studni potencjału, co sugeruje możliwość istnienia 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 przyciągania. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 9: Trójwymiarowy nieliniowy krajobraz energetyczny ilustrujący potencjał alternatywnych stanów regulacji emocjonalnej. Trójwymiarowa reprezentacja krajobrazu energetycznego Lapunowa w przestrzeni faz regulacji emocji (E–R). Liczne studnie potencjału sugerują możliwość występowania 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ę tej ryciny.

Właściwości stabilności układu nieliniowego zostały dodatkowo przeanalizowane przy użyciu metod opartych na wartościach własnych. Rysunek 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 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 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 tę wartość wskazują progi niestabilności i potencjalne przejścia między reżimami.

Rycina 10: Mapa ciepła stabilności nieliniowego systemu stres–emocje w oparciu o maksymalną rzeczywistą wartość własną. 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 oznaczają zmniejszoną stabilność i większe prawdopodobieństwo niestabilności. Kliknij tutaj, aby zobaczyć powiększoną wersję tej ryciny.

Dynamikę przestrzeni fazowej przedstawiono na Rysunku 11, który obrazuje pole wektorowe i 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. Ta zbieżność potwierdza pomyślną implementację protokołu i dowodzi, że nieliniowe ramy teoretyczne spójnie odwzorowują 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 warunków początkowych zbiegają się w stronę wspólnego atraktora w płaszczyźnie stres–emocje.

Rysunek 11: Reprezentacja pola wektorowego w płaszczyźnie fazowej i linii prądu dla dynamiki stres–emocje. Reprezentacja pola wektorowego i linii prądu dla nieliniowego systemu stres–emocje w 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ć większą 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 zobaczyć powiększoną wersję tej ryciny.

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

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) oraz 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 stanów emocjonalnych 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 zobaczyć powiększoną wersję tego rysunku.

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

Rysunek 15: Dwuosiowa analiza porównawcza przy zmiennym wzmocnieniu regulacji (κ). Porównanie stanu ustalonego równowagi emocjonalnej przy różnym 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 Rysunek 16 przedstawia porównanie czułości w formie wykresu rozrzutu w zależności od zmian współczynnika nieliniowości emocjonalnej (µ). W ramach modelu liniowego równowaga emocjonalna pozostaje niezmieniona ze względu na brak składników nieliniowych. W przeciwieństwie do tego, model nieliniowy wykazuje spadek równowagi emocjonalnej wraz ze wzrostem µ, co demonstruje stabilizujący wpływ nasycenia sześciennego na dynamikę emocjonalną.

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

Poza ich 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ą buforować skutki nieoczekiwanych stresorów akademickich, wspierając tym samym odporność i stabilność emocjonalną. Podobnie, istnienie stabilnych atraktorów można interpretować jako psychologicznie 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 maladaptacyjne reakcje stresowe. Analizy wrażliwości wskazują ponadto, ż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. 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 obejmują kilka właściwości strukturalnych nieobecnych 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 zmniejszone 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 powtarzalnoś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 przepływ obliczeniowy niezbędny do powtórzenia wszystkich symulacji, analiz stabilności, trajektorii w przestrzeni fazowej, analiz wrażliwości, map ciepła wartości własnych i krajobrazów energii Lapunowa 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 pełny przebieg symulacji; mapowanie matematycznych zmiennych stanu i parametrów modelu na mierzalne konstrukty psychologiczne; zalecane narzędzia oceny psychologicznej; wytyczne dotyczące szacowania parametrów; proponowane ramy dla przyszłej walidacji empirycznej i kalibracji z wykorzystaniem danych studenckich; oraz potencjalne strategie wdrażania w celu walidacji podłużnej i zastosowań edukacyjnych.Kliknij 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.