The Maximal Almost Sure Lyapunov Exponent of Three-Dimensional Linear Stratonovich Stochastic Differential Equations

The sign of the maximal almost sure Lyapunov exponent determines the stability of stochastic systems, while its numerical computation for three-dimensional linear Stratonovich stochastic differential equations remains challenging due to the failure of classical two-dimensional strategies. The spherical angular motion of 3D systems produces a Fokker–Planck equation with intractable mixed partial derivatives, preventing conventional analytical solutions. This paper develops a unified computational framework for three-dimensional linear Stratonovich stochastic systems using analytical derivation for degenerate cases and physics-informed neural network (PINN) approximation for general non-degenerate scenarios. For degenerate systems, we reduce the coefficient matrix to a lower triangular form via orthogonal transformation and establish tight upper bounds based on the logarithmic growth property of the Wiener process, yielding closed-form expressions for the maximal almost sure Lyapunov exponent under all parameter sign configurations. For non-degenerate systems, we reformulate the Fokker–Planck equation in spherical coordinates and construct a customized PINN with trigonometric encoding to enforce periodic boundary conditions. The network is trained by joint loss functions of equation residuals, boundary constraints and normalization consistency, and the converged stationary density is substituted into the Furstenberg–Khasminskii formula to calculate the exponent via Gauss–Legendre quadrature. Monte Carlo simulations confirm the accuracy and robustness of the proposed method, which reliably identifies the sign of the maximal almost sure Lyapunov exponent even in near-critical regimes. Numerical experiments on a 3D stochastic Hopf bifurcation model show that noise negatively shifts the bifurcation point, with the offset linearly proportional to the squared noise intensity. This work extends Lyapunov stability analysis from two-dimensional to three-dimensional linear Stratonovich stochastic systems, offering an effective tool for stability evaluation of general three-dimensional stochastic dynamical models.

Authors

Institutions

Publication Details

Journal
Mathematics
Published
2026-06-19
DOI
https://doi.org/10.3390/math14122207
Primary Topic
Probabilistic and Robust Engineering Design
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

The Maximal Almost Sure Lyapunov Exponent of Three-Dimensional Linear Stratonovich Stochastic Differential Equations

Ziying He, Jianyue Su
Mathematics
Probabilistic and Robust Engineering Design
article

The Maximal Almost Sure Lyapunov Exponent of Three-Dimensional Linear Stratonovich Stochastic Differential Equations

Ziying He, Jianyue Su
article en

Abstract

The sign of the maximal almost sure Lyapunov exponent determines the stability of stochastic systems, while its numerical computation for three-dimensional linear Stratonovich stochastic differential equations remains challenging due to the failure of classical two-dimensional strategies. The spherical angular motion of 3D systems produces a Fokker–Planck equation with intractable mixed partial derivatives, preventing conventional analytical solutions. This paper develops a unified computational framework for three-dimensional linear Stratonovich stochastic systems using analytical derivation for degenerate cases and physics-informed neural network (PINN) approximation for general non-degenerate scenarios. For degenerate systems, we reduce the coefficient matrix to a lower triangular form via orthogonal transformation and establish tight upper bounds based on the logarithmic growth property of the Wiener process, yielding closed-form expressions for the maximal almost sure Lyapunov exponent under all parameter sign configurations. For non-degenerate systems, we reformulate the Fokker–Planck equation in spherical coordinates and construct a customized PINN with trigonometric encoding to enforce periodic boundary conditions. The network is trained by joint loss functions of equation residuals, boundary constraints and normalization consistency, and the converged stationary density is substituted into the Furstenberg–Khasminskii formula to calculate the exponent via Gauss–Legendre quadrature. Monte Carlo simulations confirm the accuracy and robustness of the proposed method, which reliably identifies the sign of the maximal almost sure Lyapunov exponent even in near-critical regimes. Numerical experiments on a 3D stochastic Hopf bifurcation model show that noise negatively shifts the bifurcation point, with the offset linearly proportional to the squared noise intensity. This work extends Lyapunov stability analysis from two-dimensional to three-dimensional linear Stratonovich stochastic systems, offering an effective tool for stability evaluation of general three-dimensional stochastic dynamical models.

MathematicsVol. 14(12)
Wuhan University of Technology (CN)
Openalex Percentile: Top 7%
Probabilistic and Robust Engineering Design
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.