窗函数法1.m

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

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;
Wn=0.54-0.46*cos((2*pi*n)./(N-1));
Rt=xcorr(Xn,'unbiased');
R=Rt(find((n>=1)&(n<=128)==1));
P=R.*Wn;
Pw=abs(fft(P));
plot(10*log10(Pw),'g');
ylabel('PSD'); 
xlabel('Freq Hz');