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