webpsd1.m

基于matlab的功率谱分析以及短时傅立叶变换

Fs=1000;
n=0:1/Fs:1;
xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));
window=boxcar(length(xn));
nfft=1024;
[Pxx,f]=periodogram(xn,window,nfft,Fs);
plot(f,10*log10(Pxx));