Non-Markovian dynamics: the governing equations of probability density functions for nonlinear dynamical systems with combined fractional Gaussian and Gaussian white noises

Bin Pei,Lifang Feng,Yunzhang Li,Yong Xu
2024-11-09
Abstract:This paper aims to investigate the non-Markovian dynamics. The governing equations are derived for probability density functions (PDFs) of non-Markovian stochastic responses to nonlinear dynamical systems (NDSs) excited by combined fractional Gaussian noises (FGNs) and Gaussian white noises (GWNs). The main difficulty here, is that, NDSs excited by FGNs cannot be augmented by filters excited by GWNs, leading to the inapplicability of Itô stochastic calculus theory. Thus, in the present work, based on fractional Wick Itô Skorohod integral and rough path theory, new non-Markovian probability density evolution methods (PDEMs) are established to derive theoretically the fractional-power time varying term-dependent governing equations (FPTV-GEs) for PDFs of non-Markovian stochastic responses to NDSs excited by FGNs and GWNs, which is a breakthrough to stochastic dynamics. Then, we extend an efficient algorithm, the local discontinuous Galerkin method, to numerically solve the FPTV-GEs. Remarkably, this proposed method attains a higher accuracy compared to the prevalent methods such as finite difference, path integral and Monte Carlo methods, and boasts a broader applicability than the path integral method, which fails to solve the FPTV-GEs for PDFs of non-Markovian stochastic responses. Finally, several numerical examples are illustrated to verify the proposed scheme.
Probability
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is how to derive the governing equations of the probability density functions (PDFs) for nonlinear dynamical systems (NDSs) jointly excited by fractional Gaussian noises (FGNs) and Gaussian white noises (GWNs). Specifically, due to the non - Markovian characteristics of these systems, the traditional Ito stochastic calculus theory is no longer applicable, so new theories and methods need to be developed to handle such problems. ### Main Challenges 1. **Non - Markovian Characteristics**: NDSs excited by FGNs do not possess the Markov property, causing the traditional Fokker - Planck - Kolmogorov (FPK) equation to be no longer applicable. 2. **Definition of Stochastic Integral**: A new form of stochastic integral needs to be defined to interpret stochastic differential equations (SDEs) driven by fractional Brownian motion (FBM), because FBM is not a semimartingale process and the traditional Ito integral is not applicable. ### Solutions The author establishes new non - Markovian probability density evolution methods (PDEMs) by combining the fractional Wick - Ito - Skorohod integral (FWIS integral) and rough path theory. This method allows solving SDEs in the path - wise sense without relying on Ito calculus. ### Key Results - **Derivation of Fractional - Power Time - Varying - Term - Dependent Governing Equations (FPTV - GEs)**: These equations describe the evolution of the probability density functions of NDSs jointly excited by FGNs and GWNs. - **Proposal of the Local Discontinuous Galerkin Method (LDG Method)**: It is used for numerically solving the above - mentioned FPTV - GEs. This method has higher accuracy and broader applicability, especially in handling the transient solutions of non - Markovian SDEs. ### Formula Example For example, for the mixed linear SDE: \[ X_t = X_0 + \int_0^t A_s X_s \, ds + \int_0^t B_s X_s \circ dW_s + \int_0^t C_s X_s \circ dB^H_s, \] where \( \circ \) represents the Stratonovich integral and \( B^H_t \) is fractional Brownian motion. Through the FWIS integral and rough path theory, the corresponding FPTV - GE can be derived: \[ \frac{\partial p(x,t)}{\partial t}=-\frac{\partial}{\partial x}\left((A_t x+\frac{1}{2} B_t^2 x+\hat{C}_t x) p(x,t)\right)+\frac{\partial^2}{\partial x^2}\left((\frac{1}{2} B_t^2 x^2+\hat{C}_t x^2) p(x,t)\right), \] where \( \hat{C}_t = C_t \int_0^t \phi(t,r) C_r \, dr \) and the initial condition is \( p(x,0)=\delta(x - x_0) \). In conclusion, this paper successfully solves the problem of probability density function evolution in non - Markovian dynamical systems by introducing new theories and numerical methods.