多变量 pdf 中多维积分和随机数的重要性采样

计算科学 数字 一体化 蒙特卡洛 可能性
2021-12-09 09:56:32

我的目标是使用蒙特卡洛积分获得五维积分的数值。我使用平均值方法得到了很好的结果,但我想尝试使用重要性抽样来获得更好的结果。

我的主要资料给了我一个如何使用一维积分来执行此方法的示例,但它甚至没有提到如何将它复制到多变量积分中。我的主要问题是从多变量概率分布函数生成一组随机数。

如果我理解正确,为单个变量执行此操作的过程是建立方程:

0xdxp(x)=z

然后求解 (是生成的均匀分布的随机数)。xz

如何为任意 5 维 pdf复制此过程?p(x1,x2,x3,x4,x5)

我找不到很好的资料,我发现的大多数资料都是专门处理 n 维高斯分布的。

1个回答

让我们回顾一下你想做的事情:你有一些集合并且想要近似一些函数的积分: “均值法”听起来像蒙特卡罗类型的近似,形式为 其中是独立同分布的。上具有均匀分布的随机变量是 V 的VR5f:VR

Vf(x)dx
Vf(x)dxvol(V)Ni=1Nf(Xi),
XiVvol(V)V

这种普通蒙特卡罗方法的问题在于,根据函数,大多数样本不会对积分做出贡献,因为在重要性采样中,想法是从适应形状的分布中采样,而不是从中均匀采样。ff(Xi)0fV

为此,设概率分布的密度。我们近似 其中是独立同分布的。中抽取的随机变量。p:V[0,)PV

Vf(x)dx=Vf(x)p(x)p(x)dx=Vf(x)p(x)dP(x)1Ni=1Nf(Yi)p(Yi),
YiP

正如您所提到的,我们可以通过使用反转采样从 P 中精确地从这只有在我们可以反转分布的累积分布函数时才有效。YiP

在多个维度中,有几种方法可以从给定的分布中提取:

  • 拒绝抽样,
  • 吉布斯采样(如果我们可以有效地从一维条件中提取),
  • 马尔可夫链蒙特卡罗(一般只提供渐近正确的样本,但可以直接应用于积分问题而无需重要性采样绕道),
  • 多一点。

为了选择一种方法,我们需要更多关于实际分布的信息。是否有一些可以利用的结构(例如,高斯混合)?