为提高高维积分的计算速度,提出一种替换Monte Carlo积分方法.将积分区域以网格的形式离散化,再在网格上以相应的密度函数之值为权函数采用离散的Gibbs抽样算法抽样,对抽样得到的样本作均匀扰动后就可获得所需的新抽样序列,从而得到积分的近似估计值.模拟表明新算法计算速度较快.
To enhance the computing speed of high dimensional integral, an alternative Monte Carlo sampling algorithm is proposed in this paper. Firstly, the integral region is partitioned into net form. Secondly, grid points are sampled by using discrete Gibbs sampling method whose weighted functions are values of corresponding density function. Lastly, a new sampling sequence will be obtained by adding an uniform variable sequence to the original sequence correspondingly, and an estimation of the integral is given. The new sampling algorithm is as simple as the traditional numerical method. Simulating output showed that the new algorithm performs very well in computing speed.