Um den vorgeschlagenen nichtlinearen dynamischen Rahmen zu validieren und seine Vorteile gegenüber herkömmlichen linearen Formulierungen zu demonstrieren, wurde eine umfassende numerische Simulationsstudie durchgeführt (Supplementary File 2). Sofern nicht anders angegeben, werden alle Ergebnisse stochastischer Simulationen als Mittelwert ± Standardabweichung (SD) aus 30 unabhängigen Realisierungen dargestellt (n = 30). Statistische Vergleiche zwischen den Simulationsszenarien wurden mithilfe eines geeigneten statistischen Tests mit einem Signifikanzniveau von p < 0,05 durchgeführt. Exakte p-Werte werden dort angegeben, wo statistische Vergleiche durchgeführt wurden. Die Simulationen untersuchten, wie sich das nichtlineare Stress–Emotion–Regulationsmodell unter systematischer Variation der Parameter verhält, und verglichen dessen strukturelle Reaktionen mit denen des traditionellen linearen Stressmodells. Die numerische Integration der zugrundeliegenden Differentialgleichungen erfolgte über ausreichend lange Zeithorizonte, um eine Konvergenz hin zu stationären oder asymptotischen Zuständen sicherzustellen. Für jedes Experiment wurden identische Ausgangsbedingungen auf beide Modelle angewandt, während jeweils ein einzelner Parameter über einen vordefinierten Bereich variiert wurde. Die resultierenden stationären emotionalen Zustände, transiente Dynamiken und Stabilitätseigenschaften wurden aufgezeichnet und visualisiert, um strukturelle Unterschiede im Systemverhalten hervorzuheben. Eine erfolgreiche Durchführung des Protokolls wurde bestätigt, wenn das nichtlineare System zu beschränkten stationären Lösungen konvergierte, stabile Attraktorstrukturen im Phasenraum aufwies und unter den Ausgangsparameterbedingungen negative dominante Jacobische Eigenwerte beibehielt.
Die vergleichende Analyse konzentrierte sich auf mehrere zentrale dynamische Parameter, darunter die Kopplungsstärke zwischen Stress und Emotion (γ), die Stressdissipationsrate (β), die Amplitude der äußeren Anregung (F), die Regulierungsgewinnkonstante (κ) und der Nichtlinearitätskoeffizient der Emotion (µ). Diese Parameter beeinflussen direkt die Rückkopplungsintensität, Stabilitätsgrenzen und die Energieumverteilung innerhalb des modellierten akademischen Umfelds. Die Simulationen zeigen, ob die Systemreaktionen einem proportionalen Skalierungsverhalten folgen, wie vom linearen Modell vorhergesagt, oder nichtlineare Phänomene wie Sättigung, Pufferung der Resilienz und Multistabilität aufweisen, wie von der vorgeschlagenen nichtlinearen Formulierung vorhergesagt.
Die Empfindlichkeit des emotionalen Gleichgewichts gegenüber der Kopplungsstärke von Stress und Emotion ist in Abbildung 4 dargestellt. Wenn der Kopplungsparameter γ ansteigt, erzeugt das lineare Modell nahezu konstante emotionale Reaktionen, was darauf hinweist, dass die Kopplungsintensität die Gleichgewichtsergebnisse strukturell nicht beeinflusst. Im Gegensatz dazu zeigt die nichtlineare Formulierung eine Abnahme der stationären emotionalen Aktivierung mit steigendem γ, was den Einfluss nichtlinearer Wechselwirkungsterme widerspiegelt, die die Stress-Emotion-Rückkopplung dynamisch regulieren.

Abbildung 4: Vergleichende Sensitivitätsanalyse in Abhängigkeit der Stärke der Stress-Emotions-Kopplung (γ). Vergleichende Balkendarstellung der stationären emotionalen Reaktion mit zunehmender Stärke γ der Stress-Emotions-Kopplung. Die blauen Balken repräsentieren das lineare Modell, die orangefarbenen Balken das nichtlineare Modell. Die horizontale Achse zeigt die Kopplungsstärke γ, die vertikale Achse das stabile emotionale Gleichgewicht. Die Balken zeigen die mittleren Werte der stationären Emotionszustände aus wiederholten Simulationen, während die Fehlerbalken die ± Standardabweichung (SD) um den Mittelwert angeben. Die Ergebnisse zeigen, dass das nichtlineare Modell eine Sättigungsdynamik und begrenzte emotionale Reaktionen bei zunehmender Kopplungsstärke aufweist, während das lineare Modell relativ unempfindlich gegenüber Änderungen der Kopplungsintensität bleibt. Die Balken repräsentieren die mittleren Werte des stationären emotionalen Gleichgewichts aus 30 unabhängigen Simulationen (n = 30), und die Fehlerbalken zeigen die ± Standardabweichung an. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Der Einfluss der Stressdissipationsrate β auf das emotionale Gleichgewicht ist in Abbildung 5 dargestellt. Das lineare Modell sagt einen steilen Rückgang der emotionalen Aktivierung mit zunehmender Dissipation voraus, was auf eine starke Parametersensitivität hinweist. Im Gegensatz dazu bleibt das nichtlineare Modell im gleichen Parameterbereich relativ stabil, aufgrund intrinsischer Regelungsdämpfung und nichtlinearer Rückkopplungsmechanismen.

Abbildung 5: Vergleichende Flächenanalyse in Abhängigkeit der Stressdissipationsrate (β). Flächenbasierter Vergleich des stationären emotionalen Gleichgewichts in Abhängigkeit von der Stressdissipationsrate β. Der blau schraffierte Bereich stellt das lineare Modell dar, der orange schraffierte Bereich das nichtlineare Modell. Die horizontale Achse zeigt β, die vertikale Achse das stationäre emotionale Gleichgewicht. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Der Zusammenhang zwischen emotionaler Aktivierung und äußerem akademischem Zwang ist in Abbildung 6 dargestellt. Die lineare Formulierung zeigt ein proportionales Ansteigen der emotionalen Aktivierung mit zunehmender Amplitude der äußeren Einwirkung. Das nichtlineare Modell hingegen weist eine Sättigungsantwort auf, bei der die emotionale Aktivierung anfangs ansteigt, sich jedoch allmählich aufgrund nichtlinearer Dämpfung und adaptiver regulatorischer Effekte stabilisiert. Unter suboptimalen Parameterbedingungen, wie einer übermäßigen Kopplung von Stress und Emotion oder unzureichender regulatorischer Verstärkung, zeigte das System instabile Trajektorien, vergrößerte Oszillationen oder einen Verlust der Gleichgewichtsstabilität, was auf eine verringerte Systemstabilität hinweist und möglicherweise Zustände repräsentiert, die mit erhöhter psychischer Belastung und einer gesteigerten Anfälligkeit für Burnout-ähnliche Übergänge assoziiert sind.

Abbildung 6: Vergleichende Reaktion des emotionalen Gleichgewichts unter äußerer Anregung (F). Die stationäre emotionale Reaktion in Abhängigkeit von der Amplitude F der äußeren Anregung. Die rote gestrichelte Kurve stellt das lineare Modell dar, während die blaue durchgezogene Kurve das nichtlineare Modell darstellt. Die horizontale Achse zeigt die Amplitude F der Anregung, und die vertikale Achse zeigt das stabile emotionale Gleichgewicht. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Die Robustheit des nichtlinearen Systems unter stochastischen Störungen wird in Abbildung 7 veranschaulicht, die die zeitliche Entwicklung der emotionalen Aktivierung und der Stressenergie unter zufälligen Störungen zeigt. Beide Variablen schwanken trotz kontinuierlicher Rauscheinleitung innerhalb begrenzter Bereiche, was darauf hinweist, dass nichtlineare Rückkopplungsmechanismen die Stabilität unter Umweltvariabilität aufrechterhalten.

Abbildung 7: Stochastische Stress-Emotions-Dynamik unter durch Rauschen verursachten Störungen. Zeitliche Entwicklung der emotionalen Aktivierung E(t) (blaue Kurve) und der Stressenergie S(t) (orangefarbene Kurve) unter stochastischen Störungen. Die horizontale Achse stellt die Zeitschritte der Simulation dar. Beide Variablen bleiben trotz kontinuierlicher, durch Rauschen verursachter Störungen begrenzt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Die Analyse des Energielandschafts bietet zusätzliche Einblicke in die Systemstabilität. Abbildung 8 zeigt die Konturdarstellung der auf Lyapunov basierenden Energie-Funktion in der Emotions-Regulations-Phasenebene, wobei sich die Systemtrajektorie einem stabilen Attraktionsbecken nähert. Die dreidimensionale Darstellung dieser Energiestruktur ist in Abbildung 9 dargestellt und zeigt mehrere Potentialtöpfe, die die Möglichkeit alternativer stabiler emotional-regulatorischer Zustände nahelegen.

Abbildung 8: Höhenliniendarstellung des Energielandes mit dynamischer Trajektorie in der Emotions–Regulations-Phasenebene. Höhenliniendarstellung des auf Lyapunov basierenden Energielandes in der Phasenebene aus emotionaler Aktivierung (E) und regulatorischer Kapazität (R). Die schwarze Trajektorie veranschaulicht die Systementwicklung in Richtung eines stabilen Attraktors. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 9: Dreidimensionale nichtlineare Energielandschaft, die das Vorhandensein alternativer emotional-regulatorischer Zustände veranschaulicht. Dreidimensionale Darstellung der Lyapunov-Energielandschaft im Emotions-Regulations- (E–R-)Phasenraum. Mehrere potenzielle Energiemulden deuten auf die Möglichkeit alternativer stabiler emotional-regulatorischer Zustände unter unterschiedlichen Systembedingungen hin. Die Energielandschaft bietet eine qualitative Visualisierung der Stabilitätsstruktur des Systems; eine direkte Bestätigung der Multistabilität erfordert jedoch zusätzliche dynamische Nachweise, wie beispielsweise Trajektorienwechsel oder eine Bifurkationsanalyse. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Die Stabilitätseigenschaften des nichtlinearen Systems wurden weiterhin mit eigenwertbasierten Methoden analysiert. Abbildung 10 zeigt eine Wärmekarte des maximalen reellen Eigenwerts der Jacobi-Matrix über verschiedene Werte der Kopplungsstärke zwischen Stress und Emotion (γ) und der Regulierungsverstärkung (κ). Eine zunehmende Regulierungsstärke erzeugt stärker negative Eigenwerte, was auf eine stärkere asymptotische Stabilität hinweist, während eine übermäßige Kopplung die Stabilität verringern kann, wenn sie nicht durch ausreichende Regulierung ausgeglichen wird. Bereiche, die durch negative maximale reelle Eigenwerte gekennzeichnet sind, entsprechen stabilen Betriebsbedingungen, während Bereiche, die sich null annähern oder diesen überschreiten, Instabilitätsschwellen und mögliche Regimeübergänge anzeigen.

Abbildung 10: Stabilitäts-Heatmap des nichtlinearen Stress-Emotions-Systems basierend auf dem maximalen Realteil des Eigenwerts. Heatmap, die den maximalen Realteil der Eigenwerte der Jacobi-Matrix über die Kopplungsstärke γ zwischen Stress und Emotion sowie die Regulierungsverstärkung κ darstellt. Je negativer die Werte sind, desto stärker ist die asymptotische Stabilität; Werte, die sich Null nähern, weisen auf eine geringere Stabilität und eine höhere Wahrscheinlichkeit von Instabilität hin. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Die Phasenraum-Dynamik ist in Abbildung 11 dargestellt, die das Vektorfeld und die Stromlinien des Stress-Emotions-Systems in der S–E-Phasenebene zeigt. Die Trajektorien konvergieren in Richtung eines stabilen Gleichgewichtsbereichs, was ein Attraktorverhalten belegt. Dieses Konvergenzverhalten bestätigt die erfolgreiche Implementierung des Protokolls und zeigt, dass das nichtlineare Modell unter den angegebenen Simulationsbedingungen konsistent stabile Dynamiken der Stress-Emotions-Regulation reproduziert. Die globale Stabilitätsstruktur des Systems wird weiterhin in Abbildung 12 veranschaulicht, wobei Trajektorien, die von mehreren Anfangsbedingungen ausgehen, sich einem gemeinsamen Attraktor in der Stress-Emotions-Ebene nähern.

Abbildung 11: Vektorfeld- und Stromliniendarstellung der Stress-Emotions-Dynamik in der Phasenebene. Vektorfeld- und Stromliniendarstellung des nichtlinearen Stress-Emotions-Systems in der Phasenebene aus Stressenergie (S) und emotionaler Aktivierung (E). Die Stromlinien konvergieren in Richtung eines stabilen Gleichgewichtsbereichs, was auf ein Attraktorverhalten hinweist. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 12: Dichtes Phasenporträt unter mehreren Anfangsbedingungen in der Stress-Emotions-Ebene. Phasenporträt, berechnet aus mehreren Anfangsbedingungen in der S–E-Phasenebene. Die Trajektorien konvergieren zu einem gemeinsamen Attraktor, was eine robuste Stabilität über verschiedene Anfangszustände hinweg zeigt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Die dreidimensionale Attraktorstruktur des nichtlinearen Systems ist in Abbildung 13 dargestellt, wobei Trajektorien im kombinierten Zustandsraum aus Stress–Emotion–Regulation verlaufen und sich einem stabilen Attraktor nähern. Die Dynamik der gekoppelten Variablen im Zeitbereich ist in Abbildung 14 veranschaulicht, wobei emotionale Aktivierung und Stressenergie vor dem Konvergieren zu stationären Gleichgewichtswerten vorübergehende Anpassungen aufweisen.

Abbildung 13: Dreidimensionales nichtlineares Attraktor im Zustandsraum von Stress–Emotion–Regulation. Dreidimensionale Trajektorie des nichtlinearen Systems im Zustandsraum, definiert durch Stress (S), emotionale Aktivierung (E) und regulatorische Kapazität (R). Die Trajektorie konvergiert zu einem stabilen Attraktor, der das langfristige Gleichgewicht des Systems darstellt. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 14: Zeitliche Entwicklung gekoppelter Stress- und emotionaler Zustände unter nichtlinearer Regulation. Zeitliche Entwicklung der emotionalen Aktivierung E(t) (linke Achse) und der Stressenergie S(t) (rechte Achse). Beide Variablen zeigen eine vorübergehende Anpassung, gefolgt von einer Konvergenz in Richtung eines stationären Gleichgewichtszustands. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Der Einfluss der Stärke der adaptiven Regelung auf das emotionale Gleichgewicht wird in Abbildung 15 untersucht. Während das lineare Modell mit steigender Regelverstärkung κ deutliche Verringerungen der emotionalen Aktivierung vorhersagt, hält das nichtlineare Modell aufgrund adaptiver Sättigungsmechanismen nahezu konstante Gleichgewichtswerte aufrecht.

Abbildung 15: Zweiachsiges Vergleichsanalyse bei Variation der Regelverstärkung (κ). Vergleich des stationären emotionalen Gleichgewichts bei variierender Regelverstärkung κ. Die blaue Kurve stellt das lineare Modell dar, und die rote gestrichelte Kurve stellt das nichtlineare Modell dar. Die horizontale Achse zeigt die Regelverstärkung κ. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Schließlich zeigt Abbildung 16 einen streudiagrammbasierten Vergleich der Empfindlichkeit, während sich der Koeffizient der emotionalen Nichtlinearität (µ) ändert. Im linearen Modell bleibt das emotionale Gleichgewicht unverändert, da nichtlineare Terme fehlen. Im Gegensatz dazu weist das nichtlineare Modell ein abnehmendes emotionales Gleichgewicht mit steigendem µ auf, was den stabilisierenden Einfluss der kubischen Sättigung auf die emotionale Dynamik verdeutlicht.

Abbildung 16: Streudiagramm zum Vergleich der Empfindlichkeit unter dem Parameter der emotionalen Nichtlinearität (µ). Vergleich des dynamischen emotionalen Gleichgewichts im stationären Zustand, während sich der Koeffizient der emotionalen Nichtlinearität µ ändert. Rote Markierungen stellen das lineare Modell dar, blaue Markierungen das nichtlineare Modell. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Über ihre rechnerische Bedeutung hinaus lassen die beobachteten dynamischen Verhaltensweisen sinnvolle Interpretationen in pädagogischen und psychologischen Kontexten zu. Die unter stochastischen Störungen beobachteten begrenzten Reaktionen deuten darauf hin, dass adaptive Regulierungsmechanismen die Auswirkungen unerwarteter akademischer Belastungen abfedern können, wodurch Resilienz und emotionale Stabilität unterstützt werden. Ebenso können die Existenz stabiler Attraktoren als psychologisch ausgeglichene Zustände interpretiert werden, in denen Studierende akademischen Druck erfolgreich regulieren, während Instabilitätsbereiche und Bifurkationsschwellen Bedingungen entsprechen könnten, unter denen Bewältigungsressourcen unzureichend werden und die Anfälligkeit für Burnout, emotionale Erschöpfung oder maladaptive Stressreaktionen zunimmt. Die Sensitivitätsanalysen zeigen zudem, dass eine Stärkung der Regulierungskapazität die Stabilitätsbereiche vergrößern und die Anfälligkeit für disruptive Übergänge verringern kann, was mögliche Implikationen für Interventionen zur Verbesserung von Bewältigungsstrategien, emotionaler Regulation und dem Wohlbefinden von Studierenden aufzeigt. Diese simulationsbasierten Ergebnisse legen nahe, dass der vorgeschlagene Rahmen als nützliche konzeptionelle Grundlage für zukünftige Untersuchungen von Stressverläufen, Resilienzmechanismen und potenziellen Interventionsstrategien dienen könnte. Eine Validierung anhand empirischer Daten von Studierenden ist jedoch erforderlich, bevor eine praktische Anwendung etabliert werden kann.
Zusammen zeigen diese Ergebnisse, dass der vorgeschlagene nichtlineare Rahmen mehrere strukturelle Eigenschaften erfasst, die in traditionellen linearen Spannungsmodellen fehlen, darunter Sättigungsverhalten, Stabilitätspufferung unter Störungen, multistabile Energielandschaften und Resilienz durch adaptive Regulation. Über die getesteten Parameterbereiche hinweg hielt das nichtlineare Modell durchgängig begrenzte emotionale Aktivierung und stabiles Attraktorverhalten aufrecht, während das lineare Modell erheblich größere Empfindlichkeit gegenüber Parameteränderungen und reduzierte Stabilitätsgrenzen aufwies. Diese Eigenschaften ermöglichen eine realistischere Darstellung der Dynamik von Stress und Emotion in komplexen akademischen Umgebungen.
Zusatzdatei 1: MATLAB-Quellcode, beschreibende Gleichungen, numerische Implementierung und Dokumentation zur Reproduzierbarkeit. Diese Zusatzdatei enthält die beschreibenden nichtlinearen Differentialgleichungen, die Basisparametersätze des Modells, Anfangsbedingungen, Konfiguration des numerischen Lösers, Konvergenzkriterien, MATLAB-Quellcode, Gleichgewichtslöser, Verfahren zur Sensitivitätsanalyse, Visualisierungsroutinen, Softwarespezifikationen, Parameterbereiche und den erforderlichen Rechenworkflow, um alle Simulationen, Stabilitätsanalysen, Trajektorien im Phasenraum, Sensitivitätsanalysen, Heatmaps der Eigenwerte sowie Lyapunov-Energielandschaften, die in der Arbeit präsentiert werden, reproduzieren zu können.Bitte klicken Sie hier, um diese Datei herunterzuladen.
Zusatzdatei 2: Dokumentation des Simulationsworkflows, Abbildung psychologischer Variablen und Rahmenwerk zur empirischen Validierung. Diese Zusatzdatei enthält den vollständigen Simulationsworkflow; die Zuordnung mathematischer Zustandsvariablen und Modellparameter zu messbaren psychologischen Konstrukten; empfohlene psychologische Assessments; Anleitungen zur Parameterschätzung; einen vorgeschlagenen Rahmen für die zukünftige empirische Validierung und Kalibrierung mithilfe von Schülerdaten; sowie mögliche Umsetzungsstrategien für die longitudinale Validierung und pädagogische Anwendungen.Bitte klicken Sie hier, um diese Datei herunterzuladen.