M=512;     %濾波器長度
L=M;      %每個數(shù)據(jù)塊的長度
N=2*M;    %FFT的長度
%% 需要添加u, y,d,e等數(shù)據(jù)，該程序就能給出類似于11.4節(jié)LMS算法的仿真圖
w=zeros(1,M);
wf=fft([w,zeros(1,M)]);       %頻域的抽頭系數(shù)
for k=0:blocknum-1;         %blocknum為數(shù)據(jù)塊的個數(shù)
    if(k==0)
        uf=fft([zeros(1,M),u(k*M+1:k*M+M)]);     %u為輸入信號，進行重疊保留法
    else
        uf=fft([u(k*M-M+1:k*M),u(k*M+1:k*M+M)]);
    end
    ytemp=ifft(uf.*wf);
    y(k*M+1:k*M+M)=ytemp(M+1:2*M);    %取后面的M個元素
    e(1,k*M+1:k*M+M)=y(k*M+1:k*M+M)-d(k*M+1:k*M+M);    %d為參考數(shù)據(jù)
    error(1,1:M)=e(1,k*M+1:k*M+M);
    ef=fft([zeros(1,M),error]);
    corrtemp=ifft(ef.*conj(uf));
    corr=corrtemp(1:M);%取前面的M個元素
    wf=wf+step*fft([corr,zeros(1,M)]);
end