卡尔曼滤波(kalman)

卡尔曼(kalman)滤波原理

        kalman滤波器可以看做状态变量在由观测生成的线性空间上的射影。

        如下状态空间模型描述的动态系统:

                                \mathit{X\left ( k+1 \right )= \Phi X\left ( k \right )+\Gamma W\left ( k \right )}            (1)

                        ​​​​​​​        ​​​​​​​        \mathit{Y(k)=HX(k)+V(k)}                  (2)

式中,k为离散时间,系统在时刻k的状态为X(k)\in R^n\boldsymbol{\mathit{}}Y(k)\in R^m\boldsymbol{\mathit{}}为对应状态的观测信号;\mathit{W(k)\in R^r}为输入的白噪声;V(k)\in R^m\mathit{}为观测噪声。

        称式(1)为状态方程,称式(2)为观测方程。称\mathit{\Phi}为状态转移矩阵,\Gamma \mathit{}为噪声驱动矩阵,H为观测矩阵。

        假设W(k)和V(k)是均值为0,方差各为Q、R的不相关白噪声,初始状态X(0)不相关与W(k)和V(k) 。

由式(1)、(2)在假设条件下,递推出kalman滤波器如下:

状态一步预测:                        \hat{X}(k+1|k)=\Phi \hat{X}(k|k)

状态更新:                        \hat{X}(k+1|k+1)=\hat{X}(k+1|k)+K(k+1)\varepsilon (k+1)

                                                      \varepsilon (k+1)=Y(k+1)-H\hat{X}(k+1|k)

滤波增益矩阵:                           K(k+1)=P(k+1|k)H^T[HP(k+1|k)H^T+R]^{-1}

一步预测协方差阵:        ​​​​​​​         P(k+1|k)=\Phi P(k|k)\Phi^{T}+\Gamma Q \Gamma ^T

协方差更新:        ​​​​​​​        ​​​​​​​     P(k+1|k+1)=[I_n-K(k+1)H]P(k+1|k)

                                                \hat{X}(0|0)=\mu _0, P(0|0)=P_0

  • 1
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Skywalker...

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值