如何估计一组数据的真实值

月下冯

于 2022-09-22 20:06:49 发布

阅读量430

点赞数

分类专栏： matlab 文章标签：算法机器学习人工智能

本文链接：https://blog.csdn.net/ilovewold/article/details/126998542

版权

matlab 专栏收录该内容

12 篇文章 8 订阅

订阅专栏

$\huge \begin{matrix} \\ k_k=\frac{e_{estk-1}}{e_{estk-1}+e_{meak}} \\ \hat{x}=\hat{x_{k-1}}+k_k(Z_k-\hat{x}_{k-1}) \\e_{estk}=(1-k_k)e_{estk-1} \end{matrix}$

clear all;
close;
clc;

%给定一组数据，测量误差为3；
%Z=10-3*randn(100,1);
Z=normrnd(10,3,100,1);
%测量误差
E_mea=3;
%给定数据初始估计值和估计误差
X_hat0=8;
E_est0=4;
%数据估计误差和估计值
E_est=E_est0;
X_hat=X_hat0;
Z(1)=X_hat0;
for i=1:length(Z)
    t(i)=i;
%计算卡尔曼增益
Kk=E_est(i)/(E_est(i)+E_mea);
%计算估计值
X_hat(i+1)=X_hat(i)+Kk*(Z(i)-X_hat(i));
%更新估计误差
E_est(i+1)=(1-Kk)*E_est(i);
end
plot(Z);hold on
plot(X_hat);hold on
scatter(t,Z,'o');