贝叶斯滤波器

贝叶斯滤波器原理与推导
本文详细介绍了贝叶斯滤波的基本概念,包括运动和观测模型,以及非线性运动和观测函数。通过数学推导,展示了如何利用贝叶斯公式和一阶马尔可夫性来更新和预测系统状态的概率分布,从而实现对未知状态的估计。该过程涉及到后验置信度、先验置信度和利用观测模型进行更新等关键步骤。

贝叶斯滤波器

定义以下运动和观测模型:

运动方程:xk=f(xk−1,vk,wk),k=1,⋯ ,K观测方程:yk=g(xk,nk),k=0,⋯ ,K \begin{aligned}&\text{运动方程:}\boldsymbol{x}_{k}=f\left(\boldsymbol{x}_{k-1}, \boldsymbol{v}_{k},\boldsymbol{w}_{k}\right), \quad k=1, \cdots, K \\&\text{观测方程:}\boldsymbol{y}_{k}=g\left(\boldsymbol{x}_{k}, \boldsymbol{n}_{k}\right), \quad k=0, \cdots, K\end{aligned} 运动方程:xk=f(xk1,vk,wk),k=1,,K观测方程:yk=g(xk,nk),k=0,,K

变量和符号定义如下:

  • xk∈RN\boldsymbol{x}_{k} \in \mathbb{R}^{N}xkRN 表示系统状态
  • x0∈RN\boldsymbol{x}_{0} \in \mathbb{R}^{N}x0RN 表示初始状态
  • vk∈RN\boldsymbol{v}_{k} \in \mathbb{R}^{N}vkRN 表示输入
  • wk∈RN\boldsymbol{w}_{k} \in \mathbb{R}^{N}wkRN 表示过程噪声
  • yk∈RM\boldsymbol{y}_{k} \in \mathbb{R}^{M}ykRM 表示测量
  • nk∈RM\boldsymbol{n}_{k} \in \mathbb{R}^{M}nkRM 表示测量噪声
  • kkk 为时间下标,最大值为 KKK
  • 函数 f(⋅)f(\cdot)f() 为非线性的运动模型
  • 函数 g(⋅)g(\cdot)g() 为非线性的观测模型

贝叶斯滤波仅使用过去以及当前的测量,通过构造一个完整的概率密度函数PDF,来计算当前状态 xk\boldsymbol{x}_{k}xk 的置信度:

P(xk∣xˇ0,v1:k,y0:k)(1) P\left(\boldsymbol{x}_{k} \mid \check{\boldsymbol{x}}_{0}, \boldsymbol{v}_{1: k}, \boldsymbol{y}_{0: k}\right) \tag{1} P(xkxˇ0,v1:k,y0:k)(1)

其中, xˇ0\check{\boldsymbol{x}}_{0}xˇ0 表示 0 时刻系统状态的先验。

贝叶斯滤波器可以用下面公式表示:

p(xk∣xˇ0,v1:k,y0:k)⏟后验置信度 =ηp(yk∣xk)⏟利用 g(⋅) 进行更新 ∫p(xk∣xk−1,vk)⏟利用 f(⋅) 进行预测 p(xk−1∣xˇ0,v1:k−1,y0:k−1)⏟先验置信度 dxk−1(2) \begin{aligned}& \underbrace{p\left(\boldsymbol{x}_{k} \mid \check{\boldsymbol{x}}_{0}, \boldsymbol{v}_{1: k}, \boldsymbol{y}_{0: k}\right)}_{\text {后验置信度 }} \\=& \eta \underbrace{p\left(\boldsymbol{y}_{k} \mid \boldsymbol{x}_{k}\right)}_{\text {利用 } \boldsymbol{g}(\cdot) \text { 进行更新 }} \int \underbrace{p\left(\boldsymbol{x}_{k} \mid \boldsymbol{x}_{k-1}, \boldsymbol{v}_{k}\right)}_{\text {利用 } \boldsymbol{f}(\cdot) \text { 进行预测 }} \underbrace{p\left(\boldsymbol{x}_{k-1} \mid \check{\boldsymbol{x}}_{0}, \boldsymbol{v}_{1: k-1}, \boldsymbol{y}_{0: k-1}\right)}_{\text {先验置信度 }} \mathrm{d} \boldsymbol{x}_{k-1}\end{aligned}\tag{2} =后验置信度  p(xkxˇ0,v1:k,y0:k)η利用 g() 进行更新 

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值