clear?all
clc
N=100;
M=2;?%濾波器階數
del=0.001;
a1=1.558;
a2=-0.81;
v=randn（1N）;
u（1）=v（1）;
for?n=2:N-1
????u（2）=a1*u（1）+v（2）;
????u（n+1）=a1*u（n）+a2*u（n-1）+v（n+1）;?
????d（2）=0.2*u（2）+0.7*u（1）;
????d（n+1）=0.2*u（n+1）+0.7*u（n）;
end
%初始化
am=zeros（NM）;
bm=zeros（NM）;
wm=zeros（NM）;
gm=zeros（NM）;
k=0;
r（1:）=ones（1M+1）;
r（2:）=ones（1M+1）;
epsf=del;
epsb=del;
km=zeros（1M）;??
for?n=3:N??
????uf=[u（n-1）?u（n-2）];
????ub=[u（n）?u（n-1）];
????ef（nn-1）=u（n）-uf*am（n-1:）‘;
????ef（nn）=r（n-1M）*ef（nn-1）;%ef是標量
????epsf0=epsf;
????epsf=epsf+ef（nn-1）;
????am（n:）=am（n-1:）+ef（nn-1）*gm（n-1:）;
????r（nM+1）=epsf0*r（n-1M）/epsf;
????km1=[0?gm（n-1:）]+ef（nn）/epsf*[1?-am（n:）];
???for?m=1:M
???????km（m）=km1（m）;
???end
???k=km1（end）;
????eb（2n-1）=ub*bm（n-1:）‘;
????eb（n+1n-1）=u（n+1-M）-ub*bm（n-1:）‘;
????r（nM）=1/（1-k*eb（nn-1））*r（nM+1）;
????eb（nn）=r（nM）*eb（nn-1）;?
????epsb=epsb+eb（nn）*eb（nn-1）;
????gm（n:）=（km+k*bm（n-1:））*r（nM）/r（nM+1）;
????bm（n:）=bm（n-1:）+gm（n:）*eb（nn-1）;
????e（nn-1）=d（n）-ub*wm（n-1:）‘;
????wm（n:）=wm（n-1:）+gm（n:）*e（nn-1）;??
end
plot（wm）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????17447??2008-06-04?20:11??FTF\1.jpg

?????文件???????2163??2008-05-28?20:43??FTF\FTF.fig

?????文件???????1135??2008-05-28?20:30??FTF\FTF.m

?????文件???????1089??2013-11-11?16:32??FTF\FTF2.asv

?????文件???????1089??2013-11-11?16:35??FTF\FTF2.m

?????目錄??????????0??2013-11-11?16:35??FTF

-----------?---------??----------?-----??----

????????????????23141????????????????????7