hpb2.m

来自「谱估计及阵列信号处理算法仿真库」· M 代码 · 共 22 行

load('sa1y1.mat');
figure;plot(y1)


L = 1024;

phyper = fftshift(periodogramse(y1,ones(size(y1)),L));
L1 = L/8;
phyb = fftshift(bartlettse(y1,L1,L1));
phyw = fftshift(welchse(y1,hanning(L1),L1/2,L1));
w = bartlett(L);
phybt = fftshift(btse(y1,w(L/2:L),L));

omega0 = 2*pi/L;
omega1 = 2*pi/L1;

figure; 
semilogy(-pi:omega0:pi-omega0,phyper,'b',-pi:omega1:pi-omega1,phyb,'k',-pi:omega1:pi-omega1,phyw,'r',-pi:omega0:pi-omega0,phybt,'g');
legend('per','bart','welch','bt')