间接法.m

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

clc
clear
Fs=500;
N=128;  %采样点%
n=1:1/Fs:2;
F=0:Fs/N:1000;  %分辨率是Fs/N%
Xn=(200.^0.5)*sin(2*pi*60*n)+(2.^0.5)*sin(2*pi*200*n)+randn(size(n));
nfft=512;
R=xcorr(Xn,'unbiased');
Rk=fft(R);
Pw=abs(Rk);
plot(10*log10(Pw),'r');
ylabel('PSD'); 
xlabel('Freq Hz');