xiangguan.m

自适应滤波算法

function out=xiangguan(D,r,lenth,JuLi,PinYi)

%……直达波加窗……
window=hanning(lenth);
WWin=D.*window';
% WWin=D; %不加窗

x1=[r zeros(1,JuLi)];      
for p=1:JuLi		
    p
   CorrDr=x1(1+p-1:lenth+p-1).*conj(WWin(1:lenth));              %加窗 
   FFTDr=(abs(fftshift(fft(CorrDr))));  
%    TufftDr(p,1:PinYi)=FFTDr(lenth-PinYi+1:lenth);
%    TufftDr(p,PinYi+1:PinYi+PinYi)=FFTDr(1:PinYi);
TufftDr(p,1:2*PinYi)=FFTDr(lenth/2-PinYi+1:lenth/2+PinYi);
end

out=TufftDr;