%Periodogram函數：功率譜估計
%用FFT直接法計算功率譜
Fs=2000;
NFFT=1024;
n=0:1/Fs:1;
x=sin(2*pi*100*n)+4*sin(2*pi*500*n)+randn(size(n));
figure(1)
window=boxcar(length(x));
periodogram(x,window,NFFT,Fs);

%對噪聲序列進行計算
X=fft(x,NFFT);
Pxx=abs(X).^2/length(n);        %求解PSD
t=0:round(NFFT/2-1);
k=t*Fs/NFFT;
P=10*log10(Pxx(t+1));           %縱坐標的單位為dB
figure(2)
plot(k,P);