clear?all;
M=64;
L=M;
N=2*M;
g=100;?%?統計仿真次數為g
Z=2560;?%?輸入信號抽樣點數N
step=0.04;
lambda=0.9;
block=11;
t=1:Z;
?a=1;
?s=a*sin（0.05*pi*t）;?%?輸入單頻信號s
?figure（1）;
??subplot（311）
??plot（treal（s））;?%?信號s時域波形
??title（‘信號s時域波形‘）;
??xlabel（‘n‘）;
??ylabel（‘s‘）;
??axis（[0Z-a-1a+1]）;
??u=awgn（s5）;?%?加入均值為零的高斯白噪聲，信噪比為3dB
????%?設置初值
??y=zeros（1Z）;?%?輸出信號y
??e=zeros（1Z）;
???estimated_power=?ones（12*M）/100;
???%?y（1:k）=u（1:k）;?%?將輸入信號xn的前k個值作為輸出y的前k個值
???w=zeros（1M）;
???wf=fft（[wzeros（1M）]）;???????%頻域的抽頭系數
???blocknum=fix（Z/M）;
???for?k=0:blocknum-1;?????????%blocknum為數據塊的個數
????if（k==0）
????????uf=fft（[zeros（1M）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=real（ifft（uf.*wf））;??%矩陣點乘
????y（k*M+1:k*M+M）=ytemp（M+1:2*M）;????%取后面的M個元素
????e（1k*M+1:k*M+M）=s（k*M+1:k*M+M）-y（k*M+1:k*M+M