论文部分内容阅读
正态分布作为自然界和科学领域最常见的分布,广泛应用于工程技术的各个领域。当前,一维正态分布积分数值计算的理论和算法已趋于成熟。近年来,由于理论和实际的需要,针对多维正态分布积分的计算成为一个热门的研究方向。由于多维正态分布的维数较大,而传统的计算一维或者二维的数值算法大多都是随着维数的增大,计算量呈几何倍的变化,这种算法并不适用于高维的情况。蒙特卡罗和拟蒙特卡罗算法,作为解决高维数值积分的常用算法,走进了人们的视野。 本文首先简单地介绍了计算多维正态分布积分的研究现状,阐述了计算多维正态分布积分面临的两个问题:变换形式和计算方法。变换形式方面,介绍了三种不同的变换形式,通过对比并确定了本文所采用的Genz的方法。计算方法方面,详细介绍了蒙特卡罗算法及拟蒙特卡罗算法。通过仿真实验得出,拟蒙特卡罗算法无论在精度上还是时间上较蒙特卡罗算法都有比较大的改进。最后针对现有的拟蒙特卡罗算法,提出了两方面的改进:一方面,当相关系数较大,协方差矩阵接近奇异时,求解出来的值会有较大的误差。通过对积分区域进行有效的排序,可以在一定程度上减少精度的损失;另一方面,在选点的问题上,除了选取已经确定的点列后,还选取了它们在这个空间的对偶点,保证了结果的可靠性。