particle filter中important sampling能否用gibbs算法,以及importance distribution怎样选取?
假设我们的未知后验是normal 分布的, 则对于 important sampling, 你要propose一个尽可能的与后验概率相似的分布(比如t-distributino)作为important function然后对这个分布进行sampling, 然后计算weight 来distort 我们的sample。 propose的distribution与posterior distribution越接近,sample的效率越好。当然,我们不知道posterior distribution是什么,但是我们可以用优化的方法大概知道其mode 与方差,这样我们就能propose一个至少mean 和 variance 差不太多的import function。important sample的一大缺点就是我们必须要知道posterior distribution的normalized constant(即其marginal likelyhood), 所以,最好能用rejection sampling. 当然,important sampling , rejection sampling 与MCMC比起类,还是MCMC比较好用:一个重要的原因就是高维灾难,在高纬度问题中,大多数地方是分布稀疏的,在选important function 的mean 或者 选rejection sampling 中要用到的 envelope function 的时候, 往往确定不了posterior distribution的"bump"在什么地方. 而是用MCMC类的方法,比如H-M算法,会自动找到bump的位置进行sampling (当然初始位置没选好也会有很长时间的burn-in,不过可以用multiple chain来解决).
推荐阅读
- 有哪些基于UKF(Unscented Kalman Filter)的改进算法其中哪些适用于对电池SOC的估计
- Unbiased Finite Impulse Response UFIR Filter 到底是啥
- 想问问关于责任链的几个问题?
- 请问这段删除list中空字符串的python代码中的filter是咋起作用的
- 怎样给ABAQUS中的SPH particle施加一个surface tension呢
- 鲁棒卡尔曼滤波 Robust Kalman Filter 到底是啥
- 有哪些类似于BloomFilter或者GeoHash这样原理容易理解但设计很巧妙的算法可以推荐
- Kalman Filter在非匀加速情况下还可以应用吗
- 卡尔曼滤波 Kalman Filter 有哪些推导方式
- python 中用filter求素数的时候,这里的筛选函数是啥意思啊
