Local Galerkin Method for the Approximate Solutions to General FPK Equations

In this paper, the method proposed recently by the author for the solution of probability density function (PDF) of nonlinear stochastic systems is presented in detail and extended for more general problems of stochastic differential equations (SDE), therefore the Fokker Planck Kolmogorov (FPK) equation is expressed in general form with no limitation on the degree of nonlinearity of the SDE, the type of δ correlated excitations, the existence of multiplicative excitations, and the dimension of SDE or FPK equation. Examples are given and numerical results are provided for comparing with known exact solution to show the effectiveness of the method.

A deep learning method based on prior knowledge with dual training for solving FPK equation

The evolution of the probability density function of a stochastic dynamical system over time can be described by a Fokker–Planck–Kolmogorov(FPK) equation, the solution of which determines the distribution of macroscopic variables in the stochastic dynamic system. Traditional methods for solving these equations often struggle with computational efficiency and scalability, particularly in high-dimensional contexts. To address these challenges, this paper proposes a novel deep learning method based on prior knowledge with dual training to solve the stationary FPK equations. Initially, the neural network is pre-trained through the prior knowledge obtained by Monte Carlo simulation(MCS). Subsequently, the second training phase incorporates the FPK differential operator into the loss function, while a supervisory term consisting of local maximum points is specifically included to mitigate the generation of zero solutions. This dual-training strategy not only expedites convergence but also enhances computational efficiency, making the method well-suited for high-dimensional systems. Numerical examples, including two different two-dimensional(2D), six-dimensional(6D), and eight-dimensional(8D) systems, are conducted to assess the efficacy of the proposed method. The results demonstrate robust performance in terms of both computational speed and accuracy for solving FPK equations in the first three systems. While the method is also applicable to high-dimensional systems, such as 8D, it should be noted that computational efficiency may be marginally compromised due to data volume constraints.

Dimension-reduction of FPK equation via equivalent drift coefficient

The Fokker–Planck–Kolmogorov（FPK） equation plays an essential role in nonlinear stochastic dynamics. However, neither analytical nor numerical solution is available as yet to FPK equations for high-dimensional systems. In the present paper, the dimension reduction of FPK equation for systems excited by additive white noise is studied. In the proposed method, probability density evolution method（PDEM）, in which a decoupled generalized density evolution equation is solved, is employed to reproduce the equivalent flux of probability for the marginalized FPK equation. A further step of constructing an equivalent coefficient finally completes the dimension-reduction of FPK equation. Examples are illustrated to verify the proposed method.

Semi-analytical steady-state response prediction for multi-dimensional quasi-Hamiltonian systems

The majority of nonlinear stochastic systems can be expressed as the quasi-Hamiltonian systems in science and engineering. Moreover, the corresponding Hamiltonian system offers two concepts of integrability and resonance that can fully describe the global relationship among the degrees-of-freedom(DOFs) of the system. In this work, an effective and promising approximate semi-analytical method is proposed for the steady-state response of multi-dimensional quasi-Hamiltonian systems. To be specific, the trial solution of the reduced Fokker–Plank–Kolmogorov(FPK) equation is obtained by using radial basis function(RBF) neural networks. Then, the residual generated by substituting the trial solution into the reduced FPK equation is considered, and a loss function is constructed by combining random sampling technique. The unknown weight coefficients are optimized by minimizing the loss function through the Lagrange multiplier method. Moreover, an efficient sampling strategy is employed to promote the implementation of algorithms. Finally, two numerical examples are studied in detail, and all the semi-analytical solutions are compared with Monte Carlo simulations(MCS) results. The results indicate that the complex nonlinear dynamic features of the system response can be captured through the proposed scheme accurately.

Stability analysis of the projectile based on random center manifold reduction

The center manifold method has been widely used in the field of stochastic dynamics as a dimensionality reduction method.This paper studied the angular motion stability of a projectile system under random disturbances.The random bifurcation of the projectile is studied using the idea of the Routh-Hurwitz stability criterion,the center manifold reduction,and the polar coordinates transformation.Then,an approximate analytical presentation for the stationary probability density function is found from the related Fokker–Planck equation.From the results,the random dynamical system of projectile generates three different dynamical behaviors with the changes of the bifurcation parameter and the noise strength,which can be a reference for projectile design.

A Novel PDF Shape Control Approach for Nonlinear Stochastic Systems

In this work,a novel shape control approach of the probability density function(PDF)for nonlinear stochastic systems is presented.First,we provide the formula for the PDF shape controller without devising the control law of the controller.Then,based on the exact analytical solution of the Fokker-PlanckKolmogorov(FPK)equation,the product function of the polynomial and the exponential polynomial is regarded as the stationary PDF of the state response.To validate the performance of the proposed control approach,we compared it with the exponential polynomial method and the multi-Gaussian closure method by implementing comparative simulation experiments.The results show that the novel PDF shape control approach is effective and feasible.Using an equal number of parameters,our method can achieve a similar or better control effect as the exponential polynomial method.By comparison with the multiGaussian closure method,our method has clear advantages in PDF shape control performance.For all cases,the integral of squared error and the errors of first four moments of our proposed method were very small,indicating superior performance and promising good overall control effects of our method.The approach presented in this study provides an alternative for PDF shape control in nonlinear stochastic systems.

ON THE MAXIMAL LYAPUNOV EXPONENT FOR A REAL NOISE PARAMETRICALLY EXCITED CO-DIMENSION TWO BIFURCATION SYSTEM (Ⅱ)

For a co_dimension two bifurcation system on a three_dimensional central manifold, which is parametrically excited by a real noise, a rather general model is obtained by assuming that the real noise is an output of a linear filter system_a zeromean stationary Gaussian diffusion process which satisfies detailed balance condition. By means of the asymptotic analysis approach given by L. Arnold and the expression of the eigenvalue spectrum of Fokker_Planck operator, the asymptotic expansions of invariant measure and maximal Lyapunov exponent for the relevant system are obtained.

ON THE MAXIMAL LYAPUNOV EXPONENT FOR A REAL NOISE PARAMETRICALLY EXCITED CO_DIMENSION TWO BIFURCATION SYSTEM (Ⅰ)

For a real noise parametrically excited co_dimension two bifurcation system on a three_dimensional central manifold, a model of enhanced generality is developed in the present paper by assuming the real noise to be an output of a linear filter system, namely,a zero_mean stationary Gaussian diffusion process that satisfies the detailed balance condition. On such basis, asymptotic expansions of invariant measure and maximal Lyapunov exponent for the relevant system are established by use of Arnold asymptotic analysis approach in parallel with the eigenvalue spectrum of Fokker_Planck operator.

A highly-efficient method for stationary response of multi-degree-of-freedom nonlinear stochastic systems

Analytical and numerical studies of multi-degree-of-freedom(MDOF) nonlinear stochastic or deterministic dynamic systems have long been a technical challenge.This paper presents a highly-efficient method for determining the stationary probability density functions(PDFs) of MDOF nonlinear systems subjected to both additive and multiplicative Gaussian white noises. The proposed method takes advantages of the sufficient conditions of the reduced Fokker-Planck-Kolmogorov(FPK) equation when constructing the trial solution. The assumed solution consists of the analytically constructed trial solutions satisfying the sufficient conditions and an exponential polynomial of the state variables, and delivers a high accuracy of the solution because the analytically constructed trial solutions capture the main characteristics of the nonlinear system. We also make use of the concept from the data-science and propose a symbolic integration over a hypercube to replace the numerical integrations in a higher-dimensional space, which has been regarded as the insurmountable difficulty in the classical method of weighted residuals or stochastic averaging for high-dimensional dynamic systems. Three illustrative examples of MDOF nonlinear systems are analyzed in detail. The accuracy of the numerical results is validated by comparison with the Monte Carlo simulation(MCS) or the available exact solution. Furthermore, we also show the substantial gain in the computational efficiency of the proposed method compared with the MCS.

Technique of Probability Density Function Shape Control for Nonlinear Stochastic Systems

The shape control of probability density function(PDF) of the system state is an important topic in stochastic systems. In this paper, we propose a control technique for PDF shape of the state variable in nonlinear stochastic systems. Firstly, we derive and prove the form of the controller by investigating the Fokker-PlanckKolmogorov(FPK) equation arising from the stochastic system. Secondly, an approach for getting approximate solution of the FPK equation is provided. A special function including some parameters is taken as the approximate stationary solution of the FPK equation. We use nonlinear least square method to solve the parameters in the function, and capture the approximate solution of the FPK equation. Substituting the approximate solution into the form of the controller, we can acquire the PDF shape controller. Lastly, some example simulations are conducted to verify the algorithm.