间接法：
xn=[0,1,1,0,1,0];
nfft=128; %这个数可变,取决于您需要的频率分辨率
cxn=xcorr(xn,'unbiased'); %计算序列的自相关函数
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; %可设
xn=[0,1,1,0,1,0];
window=boxcar(length(xn)); %矩形窗
nfft=1024; %可设
[Pxx,f]=periodogram(xn,window,nfft,Fs); %直接法
figure;plot(f,10*log10(Pxx));
运行会出图的,可以看到两种方法的图类似,关于功率谱的深入理解,还有一个非常好的网址：
比较了几种功率谱估计的方法