% 離散隨機信號處理 Project I
% 功率譜估計
% 馮偉  2005310367  fengw@wireless.mdc.tsinghua.edu.cn
% 功能: 周期圖法進行譜估計   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;