代码示例:РРclear;РРFs=1000;РРn=0:1/Fs:1;РРxn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));РРnfft=1024;РРwindow=boxcar(100); %?矩形窗Р[matlab实现经典功率谱估计]matlab功率谱估计Р[matlab实现经典功率谱估计]matlab功率谱估计Р8 / 4Р[matlab实现经典功率谱估计]matlab功率谱估计РРРРРРРwindow1=hamming(100); % 海明窗РРwindow2=blackman(100); %blackman?窗РРnoverlap=20; %?数据无重叠РРrange= “half ” ; % 频率间隔为 [0 Fs/2]?,只计算一半的频РРР率РР[Pxx,f]=pwelch(xn,window,noverlap,nfft,Fs,range);РР[Pxx1,f]=pwelch(xn,window1,noverlap,nfft,Fs,range);РР[Pxx2,f]=pwelch(xn,window2,noverlap,nfft,Fs,range);РРplot_Pxx=10*log10(Pxx);РРplot_Pxx1=10*log10(Pxx1);РРplot_Pxx2=10*log10(Pxx2);РРfigure(1)РРplot(f,plot_Pxx);РРpause;РРfigure(2)РРplot(f,plot_Pxx1);РРpause;РРfigure(3)РРplot(f,plot_Pxx2);Р[matlab实现经典功率谱估计]matlab功率谱估计Р[matlab实现经典功率谱估计]matlab功率谱估计Р9 / 4Р[matlab实现经典功率谱估计]matlab功率谱估计