自适应滤波器

用matlab编写的自适应滤波器，实现了自适应滤波器

clear all; close all;
%load s.wav and save
[s]=wavread（‘test.wav‘）; 
%add white gaussian noise to signal s and save
x=awgn（s30）;   

N=length（x）;
D=5;           %FIR滤波器长度为5
w=zeros（1D）;  %权矢量为W
te=zeros（1N）; %迭代过程中y输出
te1=zeros（1N）;%迭代过程中误差信号e的输出
te2=zeros（1N）;%迭代过程中均方误差信号的输出
   
%自适应滤波器迭代算法
for n=D+1:N
x1=x（n-1:-1:n-D）;
y=w*x1‘;
e=s（n）-y;       %误差e
mu=5e-3;         %收敛速度控制系数
w=w+mu*e*x1;     %调整权矢量
te（n-D）=y;
te1（n-D）=e;
te2（n-D）=e^2;
end

%plot s and x signal
figure（1）;
subplot（211）;
plot（s）; 
title（‘输入s信号‘）;
subplot（212）;
plot（x）; 
title（‘叠加噪声后的x信号‘）
figure（2）
plot（te）;
title（‘输出y信号‘）;

figure（3）
subplot（211）;
plot（te1）;
title（‘输出e信号‘）
subplot（212）
plot（te2）
title（‘输出均方误差信号e^2‘）

wavplay （s‘async‘）;
wavplay （x‘async‘）; 
wavplay（te‘sync‘）;

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

     文件        907  2010-11-17 18:19  自适应滤波器设计\zishiying.m

     文件     160768  2012-11-20 13:57  自适应滤波器设计\自适应滤器设计.doc

     目录          0  2012-11-20 13:57  自适应滤波器设计

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

               161675                    3