您现在的位置:首页 >> 项目讲座 >> 通信/信号处理 >> 内容

卡尔曼滤波

时间:2022/9/7 16:09:10 点击:

  核心提示:卡尔曼滤波...

function data_out2 = func_kalman(data);

 

T             = 0.7;

LL            = length(data);

data_out      = zeros(2,LL);  %产生2*LL的全零矩阵

Y0            = [0;1];

data_out(:,1) = Y0;           %Y的第一列等于Y0

A             = [1 T

                 0 1];   

      

B      = [1/2*(T)^2 T]';

H      = [1 0];

 

P0     = [0 0

          0 1];

      

P      = [P0 zeros(2,2*(LL-1))];

Q      = (0.1)^2; 

R      = (0.1)^2; 

X      = zeros(1,LL);

 

% kalman

for n               = 1:LL

    i               = (n-1)*2+1;

    K               = P(:,i:i+1)*H'*inv(H*P(:,i:i+1)*H'+R);%滤波增益

    data_out(:,n)   = data_out(:,n)+K*(data(:,n)-H*data_out(:,n));  %估计

    data_out(:,n+1) = A*data_out(:,n);                     %预测

    P(:,i:i+1)      = (eye(2,2)-K*H)*P(:,i:i+1);  %误差

    P(:,i+2:i+3)    = A*P(:,i:i+1)*A'+B*Q*B';   %kalman滤波

end

 

data_out2 = data_out(1,2:end);

联系:highspeedlogic

QQ :1224848052

微信:HuangL1121

邮箱:1224848052@qq.com

网站:http://www.mat7lab.com/

网站:http://www.hslogic.com/

微信扫一扫:

作者:卡尔曼滤波 来源:卡尔曼滤波
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做|FPGA项目合作(www.hslogic.com) © 2023 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科课题、Coursework、Eassy、Assignment