Characterized by an outer integral connected to an inner integral through a nonlinear function, nested integration is a challenging problem in various fields, such as engineering and mathematical finance. The available numerical methods for nested integration based on Monte Carlo (MC) methods can be prohibitively expensive owing to the error propagating from the inner to the outer integral. Attempts to enhance the efficiency of these approximations using the quasi-MC (QMC) or randomized QMC (rQMC) method have focused on either the inner or outer integral approximation. This work introduces a novel nested rQMC method that simultaneously addresses the approximation of the inner and outer integrals. The method leverages the unique nested integral structure to offer a more efficient approximation mechanism. Incorporating Owen's scrambling techniques, we address integrands exhibiting infinite variation in the Hardy--Krause sense, enabling theoretically sound error estimates. As the primary contribution, we derive asymptotic error bounds for the bias and variance of our estimator, along with the regularity conditions under which these bounds can be attained. Moreover, we derive a truncation scheme for applications in the context of expected information gain estimation and indicate how to use importance sampling to remedy the measure concentration arising in the inner integral. We verify the estimator quality through numerical experiments by comparing the computational efficiency of the nested rQMC method against standard nested MC estimation for two case studies: one in thermomechanics and the other in pharmacokinetics. These examples highlight the computational savings and enhanced applicability of the proposed approach.
翻译:嵌套积分问题在工程学和数理金融等多个领域中具有重要应用,其特点是通过非线性函数将外层积分与内层积分相连接。现有的基于蒙特卡洛(MC)方法的嵌套积分数值方法,由于误差从内层积分传播至外层积分,往往计算成本极高。此前尝试使用拟蒙特卡洛(QMC)或随机拟蒙特卡洛(rQMC)方法提升近似效率的研究,主要集中于对内层或外层积分的单独近似。本文提出了一种新颖的嵌套rQMC方法,能够同时处理内层和外层积分的近似问题。该方法利用嵌套积分特有的结构,提供了一种更高效的近似机制。通过结合Owen的置乱技术,我们处理了在Hardy--Krause意义下具有无穷变差的被积函数,从而获得了理论严密的误差估计。作为主要贡献,我们推导了估计器偏差和方差的渐近误差界,以及达到这些误差界所需的正则性条件。此外,我们针对期望信息增益估计的应用提出了一种截断方案,并说明了如何使用重要性采样来缓解内层积分中出现的测度集中问题。通过两个案例研究(热力学和药代动力学),我们将嵌套rQMC方法与标准嵌套MC估计的计算效率进行比较,从而验证了估计器的质量。这些算例凸显了所提方法在计算效率上的优势及其更广泛的应用潜力。