随机过程——卡尔曼滤波学习笔记

随机过程——卡尔曼滤波学习笔记

news/2025/1/15 0:29:31/文章来源:https://blog.csdn.net/m0_46521579/article/details/135521570

一、均方预测和随机序列分解

考虑随机序列 $X(k),X(k-1),X(k-2),\cdots, X(2),X(1),X(0)$

使用 $X(k-1),X(k-2),\cdots,X(2),X(1),X(0)$ 预测 $X(k)$

定义 $X_{MS}(k) = E[X(k)|X(k-1),X(k-2),\cdots ,X(0)]$

称为 $X(k)$ 的均方可预测部分。

若 $X(k),X(k-1),\cdots ,X(1),X(0)$ 相互独立，则 $X(k)$ 是均方不可预测的。

$X_{MS}(k) = E[X(k)|X(k-1),X(k-2),\cdots ,X(0)] = E[X(k)]$

定义随机序列 $X(k)$ 的新息序列 $V(k)=X(k)-X_{MS}(k)$

V(k)基于样本观测的条件均值为0，即均方不可预测。

V(k)与 $X_{MS}(k)$ 是正交的，即 $E[V(k)X_{MS}(k)] = 0$ 。

二、卡尔曼滤波

输入观测量 $z(n)$ ，对 $x(n)$ 进行估计得到 $\hat{x}(n)$

1. 系统模型

状态方程

$x(n)=F(n,n-1)x(n-1)+v_1(n-1)$

观测方程

$z(n)=C(n)x(n)+v_2(n)$

其中，

$x(n)$ ：状态向量， $N \times 1$

$z(n)$ ：观测向量， $M \times 1$

$v_1(n)$ ：状态噪声， $N \times 1$ ，高斯白噪声

$v_2(n)$ ：观测噪声， $M \times 1$ ，高斯白噪声

$F(n, n-1)$ ：状态转移矩阵， $N \times N$

$C(n)$ ：观测矩阵， $M \times N$

相关性质：

（1）乘积率： $F(n+1, n-1) = F(n+1, n)F(n, n-1)$

（2）状态噪声自相关矩阵： $E\left[ v_1(n) v_1^H(n) \right ] = Q_1(n) \delta(n-k)$

（3）观测噪声自相关矩阵： $E\left[ v_2(n) v_2^H(n) \right ] = Q_2(n) \delta(n-k)$

（4）噪声独立性： $E\left[ v_1(n) v_2^H(n) \right ] = 0$

2. 新息过程

定义MMSE下的预测误差 $e(n)$ 为新息过程，记作 $a(n)$

$a(n) = z(n) - \hat{d}(n) = z(n) - W^Hz_{n-1}$

$\hat{d}(n)$ 记作 $\hat{z}\left( n | Z_{n-1} \right )$ 是用前n-1个观测值对z(n)进行MMSE估计。

其中，

观测向量： $z_{n-1} = \left[ z(1), z(2), z(3), \cdots , z(n-1) \right ]^T$

权向量： $W = \left[ w_1, w_2, \cdots, w_{n-1} \right ]^{T}$

性质

（1）当前新息与以前各观测量正交

$E\left[ a(n) z^{*}(n) \right ] = 0, k=0,1,2,\cdots,n-1$

（2）当前新息与以前各新息正交

$E\left[ a(n) a^{*}(n) \right ] = 0, k=0,1,2,\cdots,n-1$

（3）新息与观测量等价

3. 用观测量估计（MMSE估计）状态变量

$\hat{x}(n|Z_n) = W^{H} z_n = W^HL_{n}^{-1}a_n = b^{H}a_n$

其中，

$a_n = a(n) = \left[ a(1), a(2), \cdots, a(n) \right ]^T$

$b(n) = \left[ b(1), b(2), \cdots, b(n) \right ]^T$

$A(n) = E\left[ a_n a_n^{H} \right ]$

$p_a(n) = E\left[ a_n x^{*}(n) \right ]$

递推形式：

$\hat{x}(n|Z_n) = b^{H} a_n = \sum_{i=1}^{n} b^{*}(i)a(i) = \hat{x}(n|Z_{n-1}) + b^{*}(n)a(n)$

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.mzph.cn/news/615241.shtml

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

目标检测脚本之mmpose json转yolo txt格式

目标检测脚本之mmpose json转yolo txt格式

目标检测脚本之mmpose json转yolo txt格式一、需求分析在使用yolopose及yolov8-pose 网络进行人体姿态检测任务时，有时需要标注一些特定场景的中的人型目标数据，用来扩充训练集，提升自己训练模型的效果。因为单纯的人工标注耗时费力&…

阅读更多...

响应式Web开发项目教程（HTML5+CSS3+Bootstrap）第2版例3-4 CSS 立方体

响应式Web开发项目教程（HTML5+CSS3+Bootstrap）第2版例3-4 CSS 立方体

代码 <!doctype html> <html> <head> <meta charset"utf-8"> <title>CSS 立方体</title> <link href"CSS/style.css" rel"stylesheet" type"text/css"> <style> .box {width: 200px…

阅读更多...

Spring上下文之support模块MessageSourceAccessor

Spring上下文之support模块MessageSourceAccessor

博主介绍：✌全网粉丝5W+，全栈开发工程师，从事多年软件开发，在大厂呆过。持有软件中级、六级等证书。可提供微服务项目搭建与毕业项目实战，博主也曾写过优秀论文，查重率极低，在这方面有丰富的经验✌ 博主作品：《Java项目案例》主要基于SpringBoot+MyBatis/MyBatis-plus+…

阅读更多...

「JavaSE」类和对象1

「JavaSE」类和对象1

🎇个人主页：Ice_Sugar_7 🎇所属专栏：快来卷Java啦 🎇欢迎点赞收藏加关注哦！ 类和对象 🍉类的定义🍌类的实例化 🍉this引用🍉对象的构造及初始化🍌…

阅读更多...

【JAVA】Java 中什么叫单例设计模式？请用 Java 写出线程安全的单例模式

【JAVA】Java 中什么叫单例设计模式？请用 Java 写出线程安全的单例模式

🍎个人博客：个人主页 🏆个人专栏：JAVA ⛳️ 功不唐捐，玉汝于成目录前言正文懒汉式（Lazy Initialization）： 双重检查锁定（Double-Checked Locking）…

阅读更多...

MySQL篇—自带物理克隆数据工具Clone插件介绍（第一篇，总共三篇）

MySQL篇—自带物理克隆数据工具Clone插件介绍（第一篇，总共三篇）

各位小伙伴，今天我为大家介绍一下MySQL Clone Plugin这个插件，简单来说，就是MySQL 8.0.17版本之后的一个物理克隆数据工具，它能够帮助我们快速、高效地克隆或复制数据库，极大地简化了数据库迁移、备份和恢复的过程&…

阅读更多...

5分钟了解股票交易！上海股票开户交易佣金最低是多少？怎么开户费用最低？

5分钟了解股票交易！上海股票开户交易佣金最低是多少？怎么开户费用最低？

股票交易是指通过证券市场买卖股票的活动。以下是股票交易的基本步骤： 开立证券账户：首先需要选择一家证券公司，向其提交相关材料开立证券账户，并完成账户开立手续。研究和选择股票：在决定购买股票之前，建…

阅读更多...

GPCR蛋白一般残基编号(Generic residue numbering)

GPCR蛋白一般残基编号(Generic residue numbering)

文章目录前言定义特殊情况参考连接前言在相应的文章中看到对于对于描述GPCR中的序列位置，往往在除了在当前蛋白的氨基酸序列序号意外，会在右上角标注一个类似于6 x 49的编号。经查这个编号有一个统一名称：Generic residue numbering。本文…

阅读更多...

【AT 指令开发】软件框架与接口

【AT 指令开发】软件框架与接口

目录 1 软件逻辑2.代码2.1 at_command.h2.2 at_command.c本文主要用于记录纯MCU无OS下，AT 指令开发软件框架 1 软件逻辑 2.代码 2.1 at_command.h #ifndef AT_COMMAND_H #define AT_COMMAND_Hvoid AT_CMD_Process(uint8_t *uartBuffer, uint8_t dataLen);/*描述AT指令返回值…

阅读更多...

Vue+element-china-area-data实现省市区三级联动

Vue+element-china-area-data实现省市区三级联动

安装依赖 npm install element-china-area-data -S cnpm install element-china-area-data -S 引用 import { provinceAndCityData, regionData, provinceAndCityDataPlus, regionDataPlus, CodeToText, TextToCode } from element-china-area-data； provinceAnd…

阅读更多...

写点东西《Docker入门（上）》

写点东西《Docker入门（上）》

写点东西《Docker入门（上）》环境变量 Docker 镜像 Docker CMD 与 ENTRYPOINT 有什么区别 Docker 中的网络： Docker 存储： Docker 是一个工具，允许开发人员将他们的应用程序及其所有依赖项打包到一个容器中。然后&…

阅读更多...

【Redis】Redis 进阶

【Redis】Redis 进阶

文章目录 1. BigKey1.1 MoreKey1.2 BigKey 2. 缓存双写一致性更新策略2.1 读缓存数据2.2 数据库和缓存一致性的更新策略2.3 canal 实现双写一致性 3. 进阶应用3.1 统计应用3.2 hyperloglog3.3 GEO3.4 bitmap 4. 布隆过滤器5. Redis 经典问题5.1 缓存预热5.2 缓存穿透5.3 缓存击…

阅读更多...

Android Persistent自启机制

Android Persistent自启机制

1.persistent属性的使用在开发系统级的App时，很有可能就会用persistent属性。当在AndroidManifest.xml中将persistent属性设置为true时，那么该App就会具有如下两个特性： 在系统刚起来的时候，该App也会被启动起来该App被强制杀…

阅读更多...

华为完全自研之后，中国的手机会变得更便宜好用吗？

华为完全自研之后，中国的手机会变得更便宜好用吗？

华为完全自研之后，中国的手机是否会变得更便宜好用，这是一个复杂的问题，涉及到多个因素。首先，华为完全自研意味着公司需要自主研发和生产手机的所有组件，包括处理器、摄像头、屏幕等。这将有助于降低成本&#xff0c…

阅读更多...

antd时间选择器，设置显示中文

antd时间选择器，设置显示中文

需求在实现react，里面引入antd时间选择器，默认显示为英文思路入口处使用ConfigProvider全局化配置，设置 locale 属性为中文来实现。官方文档介绍全局化配置 ConfigProvider - Ant Design 代码 import React from react; import { Prov…

阅读更多...

（Arcgis）matlab编程批量处理hdf5格式转换为tif格式

（Arcgis）matlab编程批量处理hdf5格式转换为tif格式

国家青藏高原科学数据中心全球逐日0.05时空连续地表温度数据集（2002-2022） 此代码仅用于该数据集处理版本：arcgis10.2 matlab2020 参考：MATLAB hdf(h5)文件转成tif图片格式（批量处理） 一、遇到问题 h5…

阅读更多...

MidTool的AIGC与NFT的结合-艺术创作和版权保护的革新

MidTool的AIGC与NFT的结合-艺术创作和版权保护的革新

在数字艺术和区块链技术的交汇点上，NFT（非同质化代币）正以其独特的方式重塑艺术品的收藏与交易。将MidTool（https://www.aimidtool.com/）的AIGC（人工智能生成内容）创作的图片转为NFT&#xff0c…

阅读更多...

GEE数据集——2000 年至 2022 年与传感器无关的 MODIS 和 VIIRS LAI/FPAR CDR

GEE数据集——2000 年至 2022 年与传感器无关的 MODIS 和 VIIRS LAI/FPAR CDR

2000 年至 2022 年与传感器无关的 MODIS 和 VIIRS LAI/FPAR CDR 该地理空间数据集包含关键的生物物理参数，即叶面积指数（LAI）和光合有效辐射分量（FPAR），是描述陆地生态系统特征不可或缺的参数。该数据集解…

阅读更多...

uniapp自定义封装只有时分秒的组件，时分秒范围选择

uniapp自定义封装只有时分秒的组件，时分秒范围选择

说实话，uniapp和uview的关于只有时分秒的组件实在是不行。全是日历，但是实际根本就不需要日历这玩意。百度了下，终于看到了一个只有时分秒的组件。原地址：原地址，如若侵犯请联系我删除 <template><view clas…

阅读更多...

使用Eexcl调换txt文本中的两列数据

使用Eexcl调换txt文本中的两列数据

问题描述本方法使用对txt存储的数据格式有特别要求。需要数据每行具有相同个数数据，且具有统一的间隔符号隔开。（常见的间隔符号有tab键、空格、逗号、分号等） 对于一个有空格间隔每行只有三列数据的txt文件，对调第二列和第三列…

阅读更多...

推荐文章

最新文章