随机数由 LCG(seed=42 固定)和 Box-Muller 方法生成,因此相同设置下结果可重现。
暂停时,拖动滑块即可即时更新结果。
上半部:蓝=N(0,1) 目标分布/绿=建议分布 q(x)/红点线=尾阈值 t/红填充=尾区域 x>t 下半部:蓝=直接MC估计值/红=IS估计值/绿破线=真值
目标分布 $p$ 中的期望值可用另一个建议分布 $q$ 的样本加权估计:
$$I = \mathbb{E}_p[h(X)] = \int h(x)\,p(x)\,dx = \mathbb{E}_q\!\left[h(X)\,\frac{p(X)}{q(X)}\right]$$从 $q$ 中抽取 $N$ 个样本 $X_i$,用权重 $w(X_i)=p(X_i)/q(X_i)$ 补正的重要抽样估计量:
$$\hat I_{\text{IS}} = \frac{1}{N}\sum_{i=1}^{N} h(X_i)\,w(X_i),\quad X_i \sim q$$尾概率 $P(X\gt t)$ 中 $h(x)=\mathbf 1(x\gt t)$、$p=\mathcal N(0,1)$、$q=\mathcal N(\mu,\sigma^2)$。权重为
$$w(x)=\frac{\sigma\exp(-x^2/2)}{\exp(-(x-\mu)^2/(2\sigma^2))}$$当 $\mu\approx t$ 时,几乎所有样本都落在尾部,方差降低比 $\operatorname{Var}(\hat I_{\text{MC}})/\operatorname{Var}(\hat I_{\text{IS}})$ 可达数倍至 1000 倍。