cepsd.asv

来自「功率谱估计matlab防真程序 从介绍功率谱的估计原理入手,分析了经典谱估计和现」· ASV 代码 · 共 27 行

%直接法:未加窗，未分段平均
x=wgn(100000,1,0);
[pxx1,f1]=pwelch(x,rectwin(100000),0,1024,1000);
subplot(2,2,1);
plot(f1,10*log10(pxx1));
 xlabel('Frequency'), ylabel('dBW / Hz');
grid;
%Bartlett法:加矩形窗，分段平均但不重叠
[pxx2,f2]=pwelch(x,rectwin(256),0,1024,1000);
subplot(2,2,2);
plot(f2,10*log10(pxx2));
grid;
%Welch法:加hann窗，分段平均且重叠1/2
[pxx3,f3]=pwelch(x,hann(256),128,1024,1000);
subplot(2,2,3);
plot(f3,10*log10(pxx3));
grid;
% Welch法:nuttall（加权余弦）窗，分段平均且重叠1/2
[pxx4,f4]=pwelch(x,nuttallwin(256),128,1024,1000);
subplot(2,2,4);
plot(f4,10*log10(pxx4));
grid;
% 估计的均值和方差
mean1=mean(pxx1)
mean2=mean(pxx2)
mean3=mean(pxx3)
mean4=mean(pxx4)