cepsd.m

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

% 直接法:未加窗，未分段平均
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)
var1=var(pxx1)
var2=var(pxx2)
var3=var(pxx3)
var4=var(pxx4)