psd_zqt.m

周期图法进行功率谱估计的函数

% 离散随机信号处理 Project I
% 功率谱估计
% 冯伟  2005310367  [e-mail]
% 功能: 周期图法进行谱估计   2007年5月18日

function S = PSD_zqt(x,M)     
N = length(x);
pi = 3.1415926;
S0 = zeros(1,M+1);
for mm = 1:M+1
    kk = 0:N-1;
    x_dft = x*exp(-j*pi/M*(mm-1)*kk).';    
    S0(mm) = 1/N*abs(x_dft)^2;      %周期图的计算    
end
S = S0;