%周期圖法是把隨機序列x(n)的N個觀測數據視為一能量有限的序列，
%直接計算x(n)的離散傅立葉變換，得X(k)，然后再取其幅值的平方，
%并除以N，作為序列x(n)真實功率譜的估計。
% author:閆朝星 2007.12.18
% version 0.1
%---------------------------------
Fs=600;%采樣頻率
n=0:1/Fs:1;%產生含有噪聲的序列
xn=cos(2*pi*40*n)+3*cos(2*pi*90*n)+0.1*randn(size(n));
window=boxcar(length(xn));%矩形窗
nfft=512;
[Pxx,f]=periodogram(xn,window,nfft,Fs);%直接法
plot(f,10*log10(Pxx));
window=boxcar(length(xn));%矩形窗
nfft=1024;
[Pxx,f]=periodogram(xn,window,nfft,Fs);%直接法
figure(1)
plot(f,10*log10(Pxx));