【Math】重要性采样 Importance sample推导【附带Python实现】

【Math】重要性采样 Importance sample推导【附带Python实现】

笔者在学习强化学习的过程中，经常遇到重要性采样的问题，这里简要记录一下推导过程方便以后查看。

1. Why need importance sample?

为什么需要重要性采样呢？假设我们目前有一个随机变量 $X\sim \mathcal{X}$, 并且该随机变量服从概率分布 $p_0(X)$，我们的目标是计算该随机变量的期望 ${\mathbb{E}}_{X\sim p_0}[X]$，我们同样知道随机变量 $X$ 的另一个易于获得的概率分布 $p_1(X)$，我们如何通过这个易于获得的概率分布来计算在概率分布 $p_0(X)$下的期望${\mathbb{E}}_{X\sim p_0}[X]$呢？这个技巧就被称作重要性采样（Importance Sample）。

2. Derivation of Discrete Distribution

假设随机变量$X$是服从离散的概率分布，并且我们在概率分布 $p_1(X)$下，有一致同分布的采样 $\{x_i{\}}_{i=1}^n$。在概率分布 $p_0(X)$下的期望${\mathbb{E}}_{X\sim p_0}[X]$就可以用下式来进行表示

$$\begin{gathered} {\mathbb{E}}_{X\sim p_0}[X]=\sum _{x\in \mathcal{X}}{p}_0(x)x=\sum _{x\in \mathcal{X}}{p}_1(x)\underset{f(x)}{\frac{p_0(x)}{p_1(x)}x}={\mathbb{E}}_{X\sim p_1}[f(X)] \\ {\mathbb{E}}_{X\sim p_0}[X]={\mathbb{E}}_{X\sim p_1}[f(X)]\approx \bar{f}=\frac{1}{n}\sum _{i=1}^{n}f(x_i)=\frac{1}{n}\sum _{i=1}^n\underset{\text{importance weight}}{\frac{p_0(x_i)}{p_1(x_i)}}{x}_i \end{gathered}$$

其中，$\frac{p_0(x_i)}{p_1(x_i)}$被称为重要性权重，那么通过这个重要性权重，我们就可以在概率分布$p_1(X)$下进行采样$\{x_i{\}}_{i=1}^n$，来计算期望${\mathbb{E}}_{X\sim p_0}[X]$了。

3. Derivation of Continuous Distribution

类似地，假设随机变量$X$是服从连续的概率分布，并且我们在概率密度函数 $p_1(x)$下，有一致同分布的采样 $\{x_i{\}}_{i=1}^n$。在概率函数 $p_0(x)$下的期望${\mathbb{E}}_{X\sim p_0}[X]$就可以用下式来进行表示

$${\mathbb{E}}_{X\sim p_0}[X]={\int }_{-\infty }^{\infty }x\times p_0(x)dx={\int }_{-\infty }^{\infty }{p}_1(x)\times \underset{f(x)}{\frac{p_0(x)}{p_1(x)}\times x}dx={\mathbb{E}}_{X\sim p_1}[f(X)]$$
然后我们使用大量的离散采样来估计连续的期望
$${\mathbb{E}}_{X\sim p_0}[X]={\mathbb{E}}_{X\sim p_1}[f(X)]\approx \bar{f}=\frac{1}{n}\sum _{i=1}^{n}f(x_i)=\frac{1}{n}\sum _{i=1}^n\underset{\text{importance weight}}{\frac{p_0(x_i)}{p_1(x_i)}}{x}_i$$

3. An Example

假设 $X\in \mathcal{X}=+1,-1$，概率分布$p_0(X)$满足

$$p_0(X=+1)=0.5,p_0(X=-1)=0.5$$

那么在概率分布$p_0$下的期望即为：

$${\mathbb{E}}_{X\sim p_0}[X]=(+1)\times 0.5+(-1)\times 0.5=0$$

假设另一个概率分布$p_1(X)$满足

$$p_0(X=+1)=0.8,p_0(X=-1)=0.2$$

那么在概率分布$p_1$下的期望即为：

$${\mathbb{E}}_{X\sim p_1}[X]=(+1)\times 0.8+(-1)\times 0.2=0.6$$

通过重要性采样，我们便可以通过在概率分布$p_1(X)$下的采样来计算在概率分布$p_0(X)$下的期望值，即

$${\mathbb{E}}_{X\sim p_0}[X]=\frac{1}{n}\sum _{i=1}^n\frac{p_0(x_i)}{p_1(x_i)}{x}_i$$

实现代码如下:

import numpy as np
import matplotlib.pyplot as plt
# reproducible
np.random.seed(0)

# 定义元素和对应的概率
elements = [1, -1]
probs1 = [0.5, 0.5]
probs2 = [0.8, 0.2]

# 重要性采样 importance sample
sample_times = 300
sample_list = []
i_sample_list = []
average_list = []
importance_list = []
for i in range(sample_times):
    sample = np.random.choice(elements, p=probs2)
    sample_list.append(sample)
    average_list.append(np.mean(sample_list))
    if sample == elements[0]:
        i_sample_list.append(probs1[0] / probs2[0] * sample)
    elif sample == elements[1]:
        i_sample_list.append(probs1[1] / probs2[1] * sample)
    importance_list.append(np.mean(i_sample_list))



plt.plot(range(len(sample_list)), sample_list, 'o', markerfacecolor='none', label='sample data')
plt.plot(range(len(average_list)), average_list, 'b--', label='average')
plt.plot(range(len(importance_list)), importance_list, 'g--', label='importance sampling')
plt.axhline(y=0.6, color='r', linestyle='--')
plt.axhline(y=0, color='r', linestyle='--')
plt.ylim(-1.5, 2.5) # 限制y轴显示范围
plt.xlim(0,sample_times) # 限制x轴显示范围
plt.legend(loc='upper right')
plt.show()

结果如下，可以看出之间用$p_1(X)$来进行期望的计算，随着样本数的增多，期望${\mathbb{E}}_{X\sim p_1}[X]$越来越逼近0.6，但是经过重要性采样，结果越来越逼近0，符合期望${\mathbb{E}}_{X\sim p_0}[X]=0$。