kalman滤波程序及其图像

发布时间:   来源:文档文库   
字号:
Forpersonaluseonlyinstudyandresearch;notforcommercialuse例子1:
clear
N=200;%随机序列长度Forpersonaluseonlyinstudyandresearch;notforcommercialuse
w(1=0;%初始化操作w=randn(1,N%产生随机序列x(1=0;
Forpersonaluseonlyinstudyandresearch;notforcommercialusea=1;fork=2:N;
x(k=a*x(k-1+w(k-1;%产生随机序列xnendsubplot(311t=1:N;
plot(t,x,'b'%绘制原始输入信号图形x(ntitle('原始输入信号图像'xlabel('x轴单位:t/s','color','b'ylabel('y轴单位:f/HZ','color','b'


V=randn(1,N;
q1=std(V;%计算随机列序V的标准差Rvv=q1.^2;%将标准差q1的平方赋予Rvvq2=std(w;%计算随机序列w的标准差Rww=q2.^2;c=0.2;
y=c*x+V;subplot(312t=1:N;
plot(t,y,'g'title('加噪后输入信号图像'xlabel('x轴单位:t/s','color','b'ylabel('y轴单位:f/HZ','color','b'
p(1=0;s(1=0;fort=2:N;
p1(t=a.^2*p(t-1+Rww;b(t=c*p1(t/(c.^2*p1(t+Rvv;
s(t=a*s(t-1+b(t*(y(t-a*c*s(t-1;p(t=p1(t-c*b(t*p1(t;end
%将标准差q2的平方赋予Rww%产生随机序列Y%绘制加噪后输入信号图像%初始化操作%初始化操作%输出预测值估计公式

subplot(313t=1:N;
plot(t,s,'r'%绘制经过卡尔曼滤波后图像title('经过卡尔曼滤波后图像'xlabel('x轴单位:t/s','color','b'ylabel('y轴单位:f/HZ','color','b'
出图结果为:

matlab下,kalman滤波程序:%卡尔曼滤波

本文来源:https://www.2haoxitong.net/k/doc/ed05519615791711cc7931b765ce0508763275c9.html

《kalman滤波程序及其图像.doc》
将本文的Word文档下载到电脑,方便收藏和打印
推荐度:
点击下载文档

文档为doc格式