zperiodogram.m

来自「用MATLAB实现数字信号处理的程序」· M 代码 · 共 19 行

%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);