Determination of the stationary state densities of the stochastic nonlinear dynamical systems

The stationary state probability densities appear not only in the study of dynamical systems with random vector fields, but also in the deterministic dynamical systems exhibiting chaotic behavior when the uncertainties in the initial conditions are represented with the probability densities. But since it is very hard problem to determine these densities, in this paper the new efficient method to obtain an approximate solution of Fokker-Planck-Kolmogorov equation which arises in the determination of the stationary state probability densities has been given by representing the densities with compactly supported functions. With specific choice of the compactly supported functions as piecewise multivariable polynomials which are supported on the ellipsoidal regions, the parameters to be calculated for determining the densities can be considerably decreased compared to Multi-Gaussian Closure scheme, in which the stationary densities are assumed to be the weighted average of the Gaussian densities. The main motivation to choose the compactly supported functions is that, in the chaotic dynamics the states are trapped in a specific compact subspace of the state space. The stationary state densities of two basic examples commonly considered in the literature have been estimated using the Parzen's estimator, and the densities obtained using the newly proposed method have been compared with these estimated densities and the densities obtained with the Multi-Gaussian Closure scheme. The results indicate that the presented compactly supported piecewise polynomial scheme can be successful compared to Multi-Gaussian scheme, when the system is highly nonlinear. (c) 2006 Elsevier Ltd. All rights reserved