%相關(guān)函數(shù)法：該方法先由序列估計出自相關(guān)函數(shù)R，
%然后對其進(jìn)行傅立葉變換，得到功率譜估計。
%當(dāng)延時和數(shù)據(jù)長度相比很小的時候可以有良好的估計精度。
% author:閆朝星 2007.12.18
% version 0.1
%---------------------------------
Fs=500;%采樣頻率
n=0:1/Fs:1;
%產(chǎn)生含有噪聲的序列
xn=cos(2*pi*40*n)+3*cos(2*pi*90*n)+randn(size(n));
nfft=512;
cxn=xcorr(xn,'unbiased');%計算序列的自相關(guān)函數(shù)
CXk=fft(cxn,nfft);
Pxx=abs(CXk);
index=0:round(nfft/2-1);
k=index*Fs/nfft;
plot_Pxx=10*log10(Pxx(index+1));
figure(1)
plot(k,plot_Pxx);