1、数学模型如下:

2、simulink建模

3、滤波结果如下图

4、系统仿真程序与滤波程序

%Plant
function [u]=kalman_sim(u1,u2,u3,u4)
%u1:输入加速度
dv=u1;
%u2:系统噪声W
W=u2;
%u3:测量噪声V
V=u3;
%u4:时间t
t=u4;persistent v s Ts
if t==0  Ts=0.1;s=0.0;v=0.0;z=s+V;
end
if t>0s=s+v*Ts+0.5*Ts^2*(dv+W);v=v+Ts*dv+Ts*W;z=s+V;
end
u(1)=s;
u(2)=v;
u(3)=t;
u(4)=z;
%Plant
function [u]=kalman_sim(u1,u2,u3)
%u1:输入加速度
dv=u1;
%u2:测量位移
z=u2;
%u3:时间t
t=u3;persistent x Ts Q R P C Jf B
if t==0Ts=0.1;%系统方程:C=[1.0,0.0];  %Covariances of w:B=[0.5*Ts^2;Ts];Q=B*[0.1]*B';%Q=diag([Ts*0.1*Ts',0,0]);%Covariances of v:R=[0.1];%初始估计值:x=[0.0;0.0];%初始估计误差协方差:P=diag([0.0,0.0]);%P=diag([Ts*0.1*Ts',9000000000,0.000000]);    end%后验,Measurement update:
%根据估计误差协方差和测量噪声协方差计算卡尔曼增益:
Kk=P*C'/(C*P*C'+R);%计算最优估计值:
%x=A*x+Mn*(yv-C*A*x);
x=x+Kk*(z-C*x);%计算最优估计值和真实值之间的误差协方差矩阵,为下次递推做准备:
P=(eye(2)-Kk*C)*P;%ye=C*x;           %Filtered value
u=x;    %Filtered signalerrcov=C*P*C';      %Covariance of estimation error%先验,Time update:
%根据系统状态方程计算下一状态预测值:
%x=A*x+B*u1;
Jf=eye(2);
Jf(1,1)=1;
Jf(1,2)=Ts;s=x(1);
v=x(2);
s=s+v*Ts+0.5*Ts^2*dv;
v=v+Ts*dv;
x(1)=s;
x(2)=v;%预测误差协方差:
P=Jf*P*Jf'+Q;

卡尔曼滤波simulink例子,位移和速度2变量估计相关推荐

  1. 使用加速度传感器与GPS位移求速度

    根据卡尔曼滤波来进行计算,刚开始时,不清楚为什么卡尔曼滤波能够预测物体的状态,看了几天后发现,采样率是很重要的一个点. 假设车辆做匀加速直线运动,要预测车辆的速度vt,在采样率确定的情况下,很明显时间 ...

  2. 细说振动位移、速度、加速度

    轴承的状态监测和故障诊断中,经常会使用振动信号进行分析.振动分析常用的量包括振动的位移.振动的速度和振动的加速度. 在以前的文章中我们介绍过振动位移.速度.加速度信号在进行振动分析中的选择及其原则.事 ...

  3. 卡尔曼滤波-Simulink仿真

    卡尔曼滤波-Simulink仿真 ​ 描述:卡尔曼滤波是一种优化估计算法(an optimal estimation algorithm) ​ 卡尔曼滤波器也称传感器融合算法 卡尔曼滤波器最早应用在阿 ...

  4. 0025-中点位移的速度

    题目 中点位移的速度 难度级别:A: 运行时间限制:1000ms: 运行空间限制:256000KB: 代码长度限制:2000000B 试题描述 很多物理运动非常有趣,有一种运动叫做匀加速直线运动:其速 ...

  5. 51单片机基本刷屏测试实验_基于单片机的发动机振动速度、位移和加速度测量方法...

    Single-chip microcomputer-based measuring of engine vibration speed.displacement and acceleration Ab ...

  6. 矢量求导之位移与速度及加速度

    矢量求导之位移与速度及加速度

  7. 首师大附中OJ系统 0025 中点位移的速度

    中点位移的速度 这里可以快速跳转: 试题要求 试题分析 试题代码 试题要求 难度级别:A: 运行时间限制:1000ms: 运行空间限制:256000KB: 代码长度限制:2000000B 试题描述 很 ...

  8. 位移到速度的变化 matlab,单摆运动过程的matlabsimulink建模与仿真.doc

    单摆运动过程的matlabsimulink建模与仿真.doc 1单摆运动过程的MATLA2周期与摆球质量无关3单摆振动的周期与摆长有关单摆周期的平方与摆长成正比以上结论是在理想条件下得到的结论,现对这 ...

  9. 如何加快Simulink模型的仿真速度

    有几种原因可能导致仿真运行速度很慢: 1. 模型中有一个 MATLAB 函数块 :当存在 MATLAB Fcn 模块时,在每个采样时间都会调用 MATLAB 解释器.这会大大降低仿真速度.所以,应尽可 ...

最新文章

  1. jquery和css3打造超梦幻的三维动画背景
  2. 视觉Transformer最新综述
  3. 每日一皮:完美解释 UDP 为什么会丢包...
  4. VC++ 多文档模板(添加新文档模板)编程实例
  5. 大厂JVM GC面试题
  6. linux命令之有关网络的操作命令
  7. python一行包含if else需要注意的点
  8. Android 滑动定位+吸附悬停效果实现
  9. Bootstrap优秀模板-INSPINIA.2.9.2
  10. 【STM32】中断和事件
  11. ios android 图文混排,[分享]iOS开发-50行代码实现图文混排
  12. (二)Tushare Pro教程:上市公司财务数据接口
  13. buuctf pwn wp(第四波)格式化字符串漏洞系列
  14. 中国科学技术大学计算机科学夏令营,中国科学技术大学计算机科学与技术学院(专业学位)软件工程保研夏令营...
  15. 跟着铁头干混淆4.1 ollvm控制流平坦化基本概念
  16. 找不到Microsoft Access Driver(*.mdb)ODBC驱动程序的安装例程。请重新安装驱动
  17. 运筹学实验8 最短路的求解
  18. MPC和PID实现自动驾驶车道保持原理
  19. 计算机桌面打不开原因,电脑桌面图标打不开怎么办?
  20. 云原生系列1 pod基础

热门文章

  1. 分布式数据库中间件 MyCat 安装及使用
  2. 厦门大学c语言模拟考试题,厦门大学《C语言》模拟试卷(10级).doc
  3. leetcode 买卖股票的最佳时机 II
  4. C语言/找出任意两整数之间的素数以及他们的和
  5. 【记录】jenkins 安装及环境配置(二)
  6. 一个列表包揽所有你需要的Windows应用
  7. C 语言实例 - 斐波那契数列
  8. 用于Elasticsearch数据可视化和分析的强大工具
  9. java 调试_我最喜欢的Java调试技术
  10. 建造模式实际应用_建造多少,购买多少:为聊天和消息传递应用程序提供支持...