bartlettse.m

来自「谱估计及阵列信号处理算法仿真库」· M 代码 · 共 24 行

function phi=bartlettse(y,M,L)
%
% The Bartlett method of spectra estimation.
%
% phi=bartlettse(y,M,L);
%
%      y -> the data vector
%      M -> the length of subsequences of y
%      L -> the number of psd samples
%
%    phi <- spectral estimates at L frequencies w=0, 2*pi/L, ..., 2*pi(L-1)/L

% Copyright 1996 by R. Moses


% check the lenth M
N=length(y);
if (M>N)
   error('M is greater than the data length.');
   return
end

phi=welchse(y,ones(M,1),M,L);   % bartlett is a special case of welch.