粒子滤波
假设系统的状态转移服从一阶马尔科夫模型
p(xk|xk−1)
ABCDEFGHIJKLMNOPQRSTUVWXYZ
ABC
ABCDEFGHIJKLMNOPQRSTUVWXYZABCDEFGHIJKLMNOPQRSTUVWXYZ
The following factorization
Rao-Blackwellized SIR(sampling importance resampling)
the observation model
p(zt|m,xt)
,即
p(zt|m(i)t−1,x(i)t)
the motion model
p(xt|xt−1,ut−1)
,即
p(x(i)t|x(i)t−1,ut−1)
Improved RBPF for Map Learning
Require:
St−1
,the sample set of the previous time step
zt
,the most recent laser scan
ut−1
,the most recent odometry measurement
Ensure:
St
,the new sample set
St={}
for all
s(i)t−1∈St−1
do
<x(i)t−1,w(i)t−1,m(i)t−1>=s(i)t−1
<script type="math/tex" id="MathJax-Element-928">
=s_{t-1}^{(i)}</script>
// scan matching
x′(i)t=x(i)t−1⊕ut−1
x^(i)t=argmaxxp(x|m(i)t−1,zt,x′(i)t)
其中,归一化常数
p(yk|yk−1)=∫p(yk|xk)p(xk|y1:k−1)dxk
计算里要用到积分,为了解决这个积分难的问题,可以用蒙特卡洛采样来代替计算后验概率
p^(i)n