間接法：
xn=[0,1,1,0,1,0];
nfft=128; %這個(gè)數(shù)可變,取決于您需要的頻率分辨率
cxn=xcorr(xn,'unbiased'); %計(jì)算序列的自相關(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;stem(k,plot_Pxx);
周期圖法：
Fs=2; %可設(shè)
xn=[0,1,1,0,1,0];
window=boxcar(length(xn)); %矩形窗
nfft=1024; %可設(shè)
[Pxx,f]=periodogram(xn,window,nfft,Fs); %直接法
figure;plot(f,10*log10(Pxx));
運(yùn)行會(huì)出圖的,可以看到兩種方法的圖類(lèi)似,關(guān)于功率譜的深入理解,還有一個(gè)非常好的網(wǎng)址：
比較了幾種功率譜估計(jì)的方法