我的目标是使用蒙特卡洛积分获得五维积分的数值。我使用平均值方法得到了很好的结果,但我想尝试使用重要性抽样来获得更好的结果。
我的主要资料给了我一个如何使用一维积分来执行此方法的示例,但它甚至没有提到如何将它复制到多变量积分中。我的主要问题是从多变量概率分布函数生成一组随机数。
如果我理解正确,为单个变量执行此操作的过程是建立方程:
然后求解 (是生成的均匀分布的随机数)。
如何为任意 5 维 pdf复制此过程?
我找不到很好的资料,我发现的大多数资料都是专门处理 n 维高斯分布的。
我的目标是使用蒙特卡洛积分获得五维积分的数值。我使用平均值方法得到了很好的结果,但我想尝试使用重要性抽样来获得更好的结果。
我的主要资料给了我一个如何使用一维积分来执行此方法的示例,但它甚至没有提到如何将它复制到多变量积分中。我的主要问题是从多变量概率分布函数生成一组随机数。
如果我理解正确,为单个变量执行此操作的过程是建立方程:
然后求解 (是生成的均匀分布的随机数)。
如何为任意 5 维 pdf复制此过程?
我找不到很好的资料,我发现的大多数资料都是专门处理 n 维高斯分布的。
让我们回顾一下你想做的事情:你有一些集合并且想要近似一些函数的积分:
“均值法”听起来像蒙特卡罗类型的近似,形式为
其中是独立同分布的。上具有均匀分布的随机变量是 V 的。
这种普通蒙特卡罗方法的问题在于,根据函数,大多数样本不会对积分做出贡献,因为。在重要性采样中,想法是从适应形状的分布中采样,而不是从中均匀采样。
为此,设上概率分布的密度。我们近似
其中是独立同分布的。从中抽取的随机变量。
正如您所提到的,我们可以通过使用反转采样从 P 中精确地从。这只有在我们可以反转分布的累积分布函数时才有效。
在多个维度中,有几种方法可以从给定的分布中提取:
为了选择一种方法,我们需要更多关于实际分布的信息。是否有一些可以利用的结构(例如,高斯混合)?