Simulation of Highly Skewed Non-Gaussian Stochastic Processes
- 1 December 2001
- journal article
- Published by American Society of Civil Engineers (ASCE) in Journal of Engineering Mechanics
- Vol. 127 (12), 1284-1295
- https://doi.org/10.1061/(asce)0733-9399(2001)127:12(1284)
Abstract
In this paper, a simulation methodology is proposed to generate sample functions of a stationary, non-Gaussian stochastic process with prescribed spectral density function and prescribed marginal probability distribution. The proposed methodology is a modified version of the Yamazaki and Shinozuka iterative algorithm that has certain difficulties matching the prescribed marginal probability distribution. Although these difficulties are usually sufficiently small when simulating non-Gaussian stochastic processes with slightly skewed marginal probability distributions, they become more pronounced for highly skewed probability distributions (especially at the tails of such distributions). Two major modifications are introduced in the original Yamazaki and Shinozuka iterative algorithm to ensure a practically perfect match of the prescribed marginal probability distribution regardless of the skewness of the distribution considered. First, since the underlying “Gaussian” stochastic process from which the desired non-Gaussian process is obtained as a translation process becomes non-Gaussian after the first iteration, the empirical (non-Gaussian) marginal probability distribution of the underlying stochastic process is calculated at each iteration. This empirical non-Gaussian distribution is then used instead of the Gaussian to perform the nonlinear mapping of the underlying stochastic process to the desired non-Gaussian process. This modification ensures that at the end of the iterative scheme every generated non-Gaussian sample function will have the exact prescribed non-Gaussian marginal probability distribution. Second, before the start of the iterative scheme, a procedure named “spectral preconditioning” is carried out to check the compatibility between the prescribed spectral density function and prescribed marginal probability distribution. If these two quantities are found to be incompatible, then the spectral density function can be slightly modified to make it compatible with the prescribed marginal probability distribution. Finally, numerical examples (including a stochastic process with a highly skewed marginal probability distribution) are provided to demonstrate the capabilities of the proposed algorithm.Keywords
This publication has 15 references indexed in Scilit:
- Simulation of Stationary Non-Gaussian Translation ProcessesJournal of Engineering Mechanics, 1998
- Simulation of homogeneous nonGaussian stochastic vector fieldsProbabilistic Engineering Mechanics, 1998
- Simulation of Correlated Non-Gaussian Pressure FieldsMeccanica, 1998
- Effects of spatial variability on soil liquefaction: some design recommendationsGéotechnique, 1997
- Simulation of Stochastic Processes and Fields to Model Loading and Material UncertaintiesPublished by Springer Science and Business Media LLC ,1997
- Analysis and simulation tools for wind engineeringProbabilistic Engineering Mechanics, 1997
- Generation of non-Gaussian stationary stochastic processesPhysical Review E, 1996
- Simulation of Stochastic Processes and Fields for Monte Carlo Simulation Applications: Some Recent DevelopmentsPublished by ASME International ,1995
- Stochastic differential equations: an approach to the generation of continuous non-Gaussian processesIEEE Transactions on Signal Processing, 1995
- Crossings of Non‐Gaussian Translation ProcessesJournal of Engineering Mechanics, 1984