转载:https://blog.csdn.net/wangpeng138375/article/details/74645637
问题:
g(x)已知,为x的函数,p(x)为随机变量x的概率密度函数,p(x)解析式已知,计算g(x)在p(x)分布下的均值(数字特征)。
方法:
Ep(x)(g(x))=∫xg(x)p(x)dx
Ep(x)(g(x))=∫xg(x)p(x)dx
但是现在的问题是此定积分很难计算,所以需要使用采样的方法。参照p(x)采样N个样本,按照如下方法计算:
Ep(x)(g(x))≈E^p(x)(g(x))=1N∑i=1Ng(xi)
Ep(x)(g(x))≈E^p(x)(g(x))=1N∑i=1Ng(xi)
但是,问题又来了,很难从p(x)中采样!
方法:
选择一个容易采样的分布q(x),把计算均值的公式做如下转换:
Ep(x)(g(x))=∫xg(x)p(x)dx=∫xg(x)p(x)q(x)q(x)dx=Eq(x)(g(x)p(x)q(x))
Ep(x)(g(x))=∫xg(x)p(x)dx=∫xg(x)p(x)q(x)q(x)dx=Eq(x)(g(x)p(x)q(x))
现在参照q(x)采样N个样本(x1,x2,……,xN),计算如下:
Ep(x)(g(x))≈E^p(x)(g(x))=E^q(x)(g(x)p(x)q(x))=1N∑i=1N[g(xi)p(xi)q(xi)]