aablackmantukeypsd.m

来自「Adaptive Filtering Primer with MATLAB」· M 代码 · 共 13 行

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);