%-------------------------------------------------------------------------
%exa121201_pmusic.m, for example 12.12.1 
%to test pmusic.m;
%-------------------------------------------------------------------------
clear all;

%調(diào)出數(shù)據(jù)
load test x;

N=4096;
fn=-0.5:1/N:0.5-1/N;

% 使用自相關矩陣分解的 MUSIC 算法得到功率譜估計；
xpsd=pmusic(x',10,N);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
for i=1:N
    xxpsd(i)=xpsd(N+1-i);
end
plot(fn,fftshift(xxpsd));grid on;