%?變步長頻域?qū)ο?br/>function?[e?W?mu]?=?blmsfft_vss（xdLNmu0mu_maxmu_minlampdabetaPianYiW）
%?e?對消輸出
%?x?input?reference?signal?
%?d?回波
%?L?filter?length
%?N?FFT?length
%?mu0?初始步長
%?mu_max最大步長
%?mu_min最小步長
%lambda=?0.99;?系數(shù)遺忘因子
%beta=?0.01;?步長調(diào)整系數(shù)
%?PianYi??直達(dá)波提前?點(diǎn)數(shù)
%?W?filter?coeff
ntr?=?length（x）;????????????%??temporary?number?of?iterations?
ntrB?=?floor（（ntr-N）/（N-L））;??%??temporary?number?of?block?iterations
e?=?zeros（size（x））;??????????????????????%??initialize?error?signal?vector
ee?=?zeros（1N）;
power=sum（abs（?x（1+L:N）.^2?））;???%?能量歸一化用，power可設(shè)為常數(shù)
mu（1）=mu0;
%

%??Main?loop?
for?k=0:ntrB?
????xx?=?x（1+k*（N-L）+PianYi:N+k*（N-L）+PianYi）;
????X=fft（xx）;???????%?直達(dá)波信號組成N點(diǎn)向量，每次步進(jìn)N-L點(diǎn)N-L即為塊長度
????y=ifft（X.*W）;???????????????????????????%?濾波預(yù)測值，取y的后N-L點(diǎn)
????ee（1+L:N）=d（1+L+k*（N-L）:N+k*（N-L））-y（1+L:N）;??%?回波-?預(yù)測值=對消中間向量ee，回波也每次步進(jìn)N-L點(diǎn)
????e（1+L+k*（N-L）:N+k*（N-L））=ee（1+L:N）;????%?ee的后N-L（塊長度）點(diǎn)作為對消輸出
????xe=ifft（conj（X）.*fft（ee））;???????????????%?xe的作用是更新頻域?yàn)V波器系數(shù)
????xe（1+L:N）=0;?????????????????????????????%?xe前L點(diǎn)有用，后N-L點(diǎn)置為0
????W=W+fft（xe）*mu（k+1）/power;???????????????%?歸一化LMS即NLMS更新公式，power做為能量歸一化用，power可設(shè)為常數(shù)以減小運(yùn)算量
????mu（k+2）=lampda*mu（k+1）+beta*sum（abs（?ee（1+L:N）））/sum（abs（?d（1+L+k*（N-L）:N+k*（N-L））））;
????if?mu（k+2）>=mu_max
???????mu（k+2）=mu_max;
????elseif?mu（k+2）<=mu_min
???????mu（k+2）?=?mu_min;
????end;
end
%?coeff?=?W（1:L）.‘;???????%??save?final?coefficient?vector

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-03-21?09:23??micro\
?????文件????????1346??2011-01-13?19:14??micro\CAF_CIC_FIR.m
?????文件????????4533??2013-05-16?08:49??micro\CAF_mDgen.m
?????文件????????5140??2012-07-24?17:03??micro\DBF_cancellation_CAF.m
?????文件????????1605??2011-04-08?21:44??micro\blmsfft_vss.m
?????文件????????1230??2013-05-15?21:22??micro\mD_CAF_CIC_FIR.m
?????文件????????4373??2013-05-07?09:52??micro\tfrspwv.m
?????文件????????3069??2013-05-07?09:52??micro\tfrstft.m
?????文件????????5356??2005-10-13?23:15??micro\tftb_window.m