论文部分内容阅读
在数学、物理、工程等领域经常需要计算多重积分,而传统数值积分方法对这些积分并不适用.一个主要的原因就是传统数值积分方法在高维度问题会遇上可怕的”维数诅咒”(Curse of Dimensionality).根据统计理论,这些积分可以转化为计算服从某一分布的随机变量的均值,可以使用蒙特卡洛方法(Monte Carlo Methods)或者拟蒙特卡洛方法(Quasi-Monte Carlo Methods)来估计其均值的置信区间,从而得到这些积分的可靠近似结果.这两种方法的主要思想是当独立同分布的样本数量趋向于无穷大时,样本序列的均值收敛到随机变量均值.但是这两种方法都存在缺陷,得到的结果有可能是不准确的.首先蒙特卡洛方法采用随机抽样,基于中心极限定理来估计置信区间,但是这个理论是渐近的,并没有给出如何确定样本大小的方法,并且在该理论中含有未知的信息,比如母体方差,所以不能从理论上保证随机变量的均值以给定置信水平落入估计的置信区间.另一方面,根据复杂性理论(Complexity theory),拟蒙特卡洛方法采用确定性序列估计多重积分,虽然可以得到精确的误差分析,但是往往依赖于问题先验信息的大小,比如有界变分,其计算非常困难.本文主要研究多重积分的蒙特卡洛计算问题.为了克服现有方法的不足,我们提出一种新的自适应蒙特卡洛方法来估计积分.该方法适合具有有限峰度的随机变量均值估计,在一定概率控制误差的条件下,构建出随机变量均值准确的置信区间(非渐近),同时给出相应的蒙特卡洛随机取样的样本容量计算算法.这一方法具体分为两个阶段,在第一个阶段中,利用Cantelli不等式得到的一个样本容量,这个样本容量自适应性地依赖被积函数的方差估计.第二个阶段利用前一个阶段中计算出来的方差估计和样本容量,结合Berry-Esseen不等式建立出随机变量均值的置信区间.本文方法主要有以下优点:首先,这一方法能够保证在一定置信水平下得到随机变量的准确置信区间,提高了精度.其次,将此方法向拟蒙特卡洛理论的推广是非常有意义的.在被积函数峰度有限的假设下,估计其方差的上限,以及通过确定性序列分析标准误差的技术,能有效提高收敛速度,同时降低先验信息的需求,从而极大地提高计算效率,降低了计算复杂性.第三,方法假设所受限制很少,对于四阶矩存在的随机变量其峰度一定存在,更进一步,该假设下的函数空间对于数乘运算是封闭的.本文包含四章.第一章主要介绍研究背景和研究内容以及意义.第二章为两阶段自适应蒙特卡洛方法和理论阐述.第三章给出算例和仿真研究结果.最后一章对全文进行总结,并提出进一步研究.