直接法.m

来自「现代谱估计的源程序」· M 代码 · 共 12 行

clc
clear
Fs=500;
N=128;  %采样点%
n=1:1/Fs:2;
F=0:Fs/N:1000;
Xn=(200.^0.5)*sin(2*pi*60*n)+(2.^0.5)*sin(2*pi*200*n)+randn(size(n));
nfft=1024;
Pw=(1/N)*(abs(fft(Xn)).^2);
plot(10*log10(Pw));
ylabel('PSD'); 
xlabel('Freq Hz');