频域自适应滤波源程序

function e = FDLMS（ x d M beta mu ）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Frequency domain LMS algorithm with self-orthogonalization
% Assume the block length equals the filter length
% x:    input column vector
% d:    desired response.
% M:    filter length
% beta: forgetting factor
% mu:   step size
% e:    output error
% program by Xilin Li
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Z = 2*M;                   % padding Z zeros
x = [zeros（Z1）; x];        
d = [zeros（Z1）; d];
Len = length（d）;
e=zeros（Len1）;

W = zeros（2*M1）;  
X = zeros（2*M1）;   
E = zeros（2*M1）;