阿尔法-贝塔滤波器

线性滤波器中的alpha——β滤波器设计，同卡尔曼滤波器有同样的效果

clc
clear all
close all

x0=150;
v0=0.3;
w=randn（11000）/1000;
v=randn（11000）;
t=0:999;
fai=[1 1
     0 1];
X=zeros（21000）;
X（:1）=[x0
        v0];
Xk=zeros（21000）;
Xkk=zeros（21000）;
Y=zeros（11000）;
h=[1 0];
G=[0.5
    1];
Xkk（:1）=X（:1）;
k=[0.5
    0.3];
for i=1:999
    X（:i+1）=fai*X（:i）+w（1i）*G;
    Y（1i）=h*X（:i）+v（1i）;
    Xk（:i+1）=fai*Xkk（:i）;
    Xkk（:i+1）=Xk（:i+1）+k*（Y（1i）-h*Xk（:i+1））;
end

 figure;
 plot（tX（1:））;
 grid on
 hold on
 plot（tXkk（1:）‘r--‘）;
 title（‘跟踪对比图‘）;
 xlabel（‘时间（s）‘）;
 ylabel（‘距离（Km）‘）;

 E=zeros（11000）;
 E=X-Xkk;
 figure
 plot（tE（1:））;
 grid on
 title（‘跟踪误差‘）;
 xlabel（‘时间（s）‘）;
 ylabel（‘误差（Km）‘）;
 
 figure
 plot（tE（2:））;
 grid on
 title（‘速度误差‘）;
 xlabel（‘时间（s）‘）;
 ylabel（‘误差（Km/s）‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        830  2013-09-16 00:02  α-β滤波\ab.m

     文件    1568518  2013-09-16 00:04  α-β滤波\跟踪对比图.bmp

     文件    1568518  2013-09-16 00:04  α-β滤波\跟踪误差.bmp

     文件    1568518  2013-09-16 00:04  α-β滤波\速度误差.bmp

     目录          0  2013-09-16 00:05  α-β滤波

----------- ---------  ---------- -----  ----

              4706384                    5