%-------------------------------------------------------------------------
% exa121201_pcov.m, for example 12.12.1 
% to test pcov.m;
%-------------------------------------------------------------------------
clear all;

% 調出數據；
load test x;

N=4096;tn=-0.5:1/N:0.5-1/N;
% 使用協方差算法得到功率譜估計；
xpsd=pcov(x,10,N);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
plot(tn,fftshift(xpsd));grid on;