aablackmantukeypsd.m

自适应信号处理中LMS

function[s]=aablackmantukeyosd(x,lg,L)
%function[s]=aablackmantukeypsd(x,lg,L);
%the window used is the triangle (Bartlett) window;
%x=data;lg=lag number about 20-40% of length(x)=N;
%L=desired number of spectral points (bins);
[r]=aasamplebiasedautoc(x,lg);
n=-(lg-1):1:(lg-1);
w=1-(abs(n)/lg);
rc=[fliplr(r(1,2:lg)) r];
rcw=rc.*w;
s=fft(rcw,L);