MATLAB卡尔曼|卡尔曼滤波的公式【线性】

卡尔曼滤波

卡尔曼滤波(Kalman Filter) 一种用于估计系统状态的数学算法,不是类似于高通、低通滤波器那样的频域滤波。

卡尔曼滤波基于线性动态系统的假设,它将系统的状态表示为均值和协方差矩阵,通过递归地更新和预测这些值来实现对系统状态的估计。卡尔曼滤波有两个主要的步骤:预测和更新

卡尔曼滤波具有一些优点,例如对噪声和不确定性的鲁棒性较强,能够提供较为精确的估计结果,并且计算效率较高。然而,卡尔曼滤波的应用前提是系统满足线性动态系统的假设,对于非线性系统,需要通过扩展卡尔曼滤波(Extended Kalman Filter)或无迹卡尔曼滤波(Unscented Kalman Filter)等变种算法来进行处理。

滤波结构

  • 预测步骤中,卡尔曼滤波使用系统模型和上一时刻的状态估计来预测当前时刻的状态,并计算出预测状态的均值和协方差矩阵。

  • 更新步骤中,卡尔曼滤波使用当前时刻的测量数据和预测状态的均值和协方差矩阵,通过最小二乘法得到对当前状态的估计,并重新计算估计的状态的均值和协方差矩阵。

重要公式

下面从几个重要公式的角度来分析卡尔曼滤波:

  1. 观测的一步转移,根据当前时刻的X(k-1)、状态转移矩阵 F F F、输入向量 u k − 1 u_{k-1} uk−1计算下一时刻的X_{k}
    X ^ k − = F X k − 1 + G u k − 1 \hat X_{k}^- = FX_{k-1}+G u_{k-1} X^k−=FXk−1+Guk−1

  2. 协方差预测:
    P k − = F P k − 1 F T + Q k P_{k}^{-}=FP_{k-1}F^{T}+Q_{k} Pk−=FPk−1FT+Qk

  3. 计算增益:
    K k = P k − H T ( H P k − H T + R ) − 1 K_k=P_k^-H^T(HP_k^-H^T+R)^{-1} Kk=Pk−HT(HPk−HT+R)−1

  4. 计算估计值:
    X ^ t = X ^ t − + K t ( Z k − H X ^ k − ) \hat{X}{t}=\hat{X}{t}^{-}+K_{t}(Z_{k}-H\hat{X}_{k}^{-}) X^t=X^t−+Kt(Zk−HX^k−)

  5. 更新协方差:
    P k = ( I − K k H ) P k − P_k=(I-K_kH)P_k^- Pk=(I−KkH)Pk−

其他公式

表示观测预测的公式:
Z ^ k = H X k \hat{Z}{k} = HX{k} Z^k=HXk

相关推荐
2601_9669496522 分钟前
只拿到股票代码还不够:量化程序到底还需要哪些标的信息?——从数据建模开始理解股票元数据
开发语言·python·数据分析·pandas·量化交易·股票数据·quantdash
m0_7345717634 分钟前
深入理解5g <十八>传输块tbsize的计算
开发语言·5g
RAN11426760082 小时前
1.3第一个程序的介绍
开发语言·qt
blue_ice .3 小时前
Verilog DDS 数字信号发生器:DAC 输出模拟波形与 ADC 回采观测
开发语言·单片机·嵌入式硬件·fpga开发
远翔调光芯片^138287988723 小时前
ECP5702能芯科技PD取电芯片在市场上的优势有哪些?
开发语言·人工智能·单片机·嵌入式硬件·智能家居
泡海椒4 小时前
规则引擎对比:JQuick-Java vs Drools 轻量场景选型对比
java·开发语言
泡泡鱼(敲代码中)4 小时前
Python 容器类型学习笔记:列表、元组、字典、集合
开发语言·笔记·python·pycharm
2501_933670794 小时前
平台招商运营校招能力模型:数据分析、SQL与增长场景拆解
开发语言
菜鸟~noob2335 小时前
【Bonjour】华为Peerium架构深度解析:从冯·诺依曼到百万处理器“一台计算机”,附MATLAB性能仿真【含matlab代码,可直接运行】
开发语言·matlab·华为·peerium
白远山5 小时前
24小时自助健身房系统开发实战:从需求分析到完整指南
java·开发语言·数据库·数据挖掘·需求分析