Kalman Filter-3放弃
1.卡尔曼公式详解 - 重点!!!
预测当前的值也叫先验估计\(\hat x_t^- = F\hat x_{t-1}+B\mu_{t-1}\)
基于之前的最优估计来推出当阶段的先验估计
协方差公式
\(P_t^-=FP_{t-1}F^T+Q\)
\(P_t^- 是\hat x_t^时刻的协方差矩阵\)
\(K_t ,卡尔曼增益, 很重要!!!\)
\(\hat x_t 最优估计-\color{red}{这是我们最后要得到的东西}\)
\(\hat x_t = \hat x_t^- + K_t (z_t - H \hat x_t^-)\)
\(z_t是观测值,\hat x_t是最优估计 = \hat x_t^-(估计值) + K_t (z_t - H \hat x_t^-)(观测的修正)\)
协方差矩阵的更新
\(P_t=(I-K_tH)P_t^-\)
2.举例说明
预测
有一个二维矢量,一个维度是p-position,小车的位置,还有一个是v,小车的速度\(\begin{bmatrix} p\\ v\\ \end{bmatrix}\)
假设小车是匀加速直线运动,则根据高中的物理模型,有如下预测模型
\(这里的表达式就是预测中的\color{red}{第一个公式},先验估计公式\)
然后计算先验估计对应的协方差矩阵
这里应用了几个协方差公式
\(cov(x,x)=Var(x)\)
\(cov(Ax,Ax)=Acov(x,x)A^T\)
\(cov(Ax+k,Ax+k)=Acov(x,x)A^T\)
\(所以有\)
这里可以看出和原来的2式就差了后面的Q,是因为上面的1式中的先验估计少了噪声\(w_t\)
测量
给小车加个gps观察
现在有测量模型
\(\Delta p是观测器的噪声,误差\)
状态更新
这三个公式的推导没有讲,就讲了应用,其实就是一个循环迭代的过程,类似梯度下降
这里H举了一个一维的例子,主要为了告诉我们,K是关于Q和R的,可以通过Q,R这两个超参数来调整K,告知模型是更加信任预测值还是更加信任观测值
最后就是一个循环迭代的过程
2.调整超参数
\(图片里面的化简指的是F,H取一维的化简\)