一次性学懂卡尔曼滤波

请看我的PPT,学了好久,也有幸得到了大佬指点。
在这里插入图片描述
在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

还有代码实现!

# -*- coding: utf-8 -*-
# Python源文件将以UTF-8格式保存
'''
    这个实践的已知量是导弹做水平运动,已知导弹每个时刻的速度,和速度测量仪器的标准差是v_std,
    还已知每个时刻用GPS测量出导弹位置position_noise以及GPS的方差。
'''
import numpy as np

t = np.linspace(1, 100, 100)  # 在1~100s内采样100次  等差数列
a = 0.5  # 加速度值
position = (a * t ** 2) / 2   # 这个位置是实际的位置。


position_noise = position + np.random.normal(0, 120, size=(t.shape[0]))  # 模拟生成GPS位置测量数据(带噪声)
# 生成一个正态分布的噪声,均值为0,标准差为120。
import matplotlib.pyplot as plt
# 写成一个函数   线  点  坐标。
plt.plot(t, position, label='truth position')
plt.plot(t, position_noise, label='only use measured position')
plt.plot()


# ---------------卡尔曼滤波----------------
# 初始的估计导弹的位置就直接用GPS测量的位置
predicts = [position[0]]   # 这个数组表示的是估计的导弹的位置。
position_predict = predicts[0]

predict_var = 0
pos_var = 120 ** 2  # 这是我们自己设定的位置测量仪器的方差,越大则测量值占比越低
v_std = 50  # 速度(加速度)测量仪器的标准差
for i in range(1, t.shape[0]):
    dv = (position[i] - position[i - 1]) + np.random.normal(0, 50)  # 模拟从IMU读取出的加速度
    position_predict = position_predict + dv  # 利用上个时刻的位置和速度预测当前位置
    predict_var += v_std ** 2  # 更新预测数据的方差   上一次最终的方差加上此时速度传感器的方差 此时是加法。
    # 下面是Kalman滤波
    position_predict = position_predict * pos_var / (predict_var + pos_var) + position_noise[i] * predict_var / (predict_var + pos_var)
    predict_var = (predict_var * pos_var) / (predict_var + pos_var)
    predicts.append(position_predict)

plt.plot(t, predicts, label='kalman filtered position')

plt.legend()  # 给图像加上图例。
plt.show()

谢谢阅读!应该对你会有帮助的!
欢迎关注我的公众号呀。“是编程还是玩”
在这里插入图片描述

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

指针的值是地址

觉得还不错,我会继续努力的。

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

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

打赏作者

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

抵扣说明:

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

余额充值