Study of weighted Monte Carlo algorithms with branching |
| |
Authors: | I. N. Medvedev G. A. Mikhailov |
| |
Affiliation: | (1) Institute of Computational Mathematics and Mathematical Geophysics, Siberian Branch, Russian Academy of Sciences, pr. Akademika Lavrent’eva 6, Novosibirsk, 630090, Russia |
| |
Abstract: | Various weighted algorithms for numerical statistical simulation are formulated and studied. The trajectory of an algorithm branches when the current weighting factor exceeds unity. As a result, the weight of an individual branch does not exceed unity and the variance of the estimate for the computed functional is finite. The unbiasedness and finiteness of the variance of estimates are analyzed using the recurrence “partial“ averaging method formulated in this study. The estimation of the particle reproduction factor and solutions to the Helmholtz equation are considered as applications. The comparative complexity of the algorithms is examined using a test problem. The variances of weighted algorithms with branching as applied to integral equations with power nonlinearity are analyzed. |
| |
Keywords: | Monte Carlo method variance of weighted estimates branching trajectory complexity reduction numerical solution to the Helmholtz equation integral equations with power singularity |
本文献已被 SpringerLink 等数据库收录! |