Comparison of PDEM and MCS: Accuracy and efficiency

Abstract

The accuracy and efficiency of two methods for stochastic analysis, the probability density evolution method (PDEM) and the Monte Carlo simulation (MCS) method, are compared in terms of how well they reflect the physical properties of stochastic systems. The basic principle and the numerical implementation details of PDEM and MCS are revisited. The analytical solutions of generalized probability density evolution equation (GDEE) for three typical stochastic systems are given and are to be used as the basis for comparing the two methods. It is verified that, with the rational partition of the probability space, the PDEM provides a continuous and complete reflection of physical properties over the whole probability space. Meanwhile, with the help of the numerical solution of GDEE, PDEM is efficient and accurate to describe the process of the probability density evolution of stochastic systems. In contrast, the random samples in the MCS may not reflect the physical properties of a stochastic system adequately, and the local cluster of sample points may cause redundant calculation, which leads to lower computational efficiency. Through three typical numerical examples, the paper compares the accuracy and efficiency of PDEM and MCS specifically. It is shown that, as the numerical approaches for the stochastic response of a system, the PDEM could get much higher numerical accuracy than MCS with the same number of samples. To achieve the same level of calculation accuracy, MCS needs a much higher number of samples than PDEM.