rfb.m

来自「现代信号谱估计算法的matlab代码 调用函数代码」· M 代码 · 共 24 行

function phi=rfb(y,K,L)
%
% The statically stable REFIL spectral estimator.
%
% phi=rfb(y,K,L);
%
%    y   <- the data vector (length N)
%    K   <- the ratio the baseband filter bandwidth to 1/N
%    L   <- the number of estimated spectral samples
%    phi -> the estimated spectrum
%

% Copyright 1996 by R. Moses

y=y(:);
N=length(y);       % data length

h=slepian(N,K,K);  % get the first K Slepian filters

phi=abs(fft(flipud(h).*(y*ones(1,K)),L)).^2;
if (K > 1)
   phi=mean(phi')';
end