一、概述
对于卡尔曼滤波,可以通过高斯分布的性质直接解得其概率分布。但是对于粒子滤波,由于其状态转移概率和发射概率是任意的,所以没办法得到其概率分布,只能通过采样的方法来进行估计,通常我们在实际应用中更加关心的并非概率分布而是概率分布关于某函数的期望,因此本文主要介绍如何通过采样的方法来求某个函数依概率分布的期望。
二、序列重要性采样
所以有:
三、重要性重采样
重要性重采样(Sampling Importance Resampling,SIR)是在上述SIS基础上进行改进的算法,与SIS相比其添加了两个部分,一个是重采样,另一个是特定的概率分布。
- 重采样
SIS算法会出现权值退化的现象,在一定时间后可能会出现大部分权重都逼近于的情况。简单解释一下,这主要还是由于维度灾难的问题,在高维空间中需要大量的样本。使用重采样可以缓解这一现象。
下面以为例来阐述重采样的过程,在下表中展示了3个采样粒子的权重(归一化后的权重可以理解为pdf),并且计算了pdf的cdf:
cdf
在SIS的基础上加上重采样,就是基本的粒子滤波算法(Basic Particle Filter)。
- 选择合适的提议分布
选择恰当的提议分布也是一种解决权重衰减的方法。这里一个常用的选择就是:
则此时的权重的计算公式就变为了:
- 算法
在SIS的基础上加上重采样和上述特定的提议分布就得到了SIR算法,其迭代的流程如下: