ex_11_08_1.m

来自「清华大学胡广书教授< 数字信号处理理论算法实现>>第二版一书的」· M 代码 · 共 27 行

%  ex_11_08_1.m :  习题 11.8

clear all;
close all;
load test x;

xpsd=pwelch(x,triang(32),16,4096,1);        
mmax=max(xpsd);
xpsd=xpsd/mmax;
xpsd=10*log10(xpsd+0.000001);
psdx(1:2048)=xpsd(2049:4096);
psdx(2049:4096)=xpsd(1:2048);
t=-0.5:1/4096:0.5-1/4096;
figure('color','w')
subplot(211)
plot(t,psdx);grid on;

xpsd=pwelch(x,hamming(32),16,4096,1);        
mmax=max(xpsd);
xpsd=xpsd/mmax;
xpsd=10*log10(xpsd+0.000001);
psdx(1:2048)=xpsd(2049:4096);
psdx(2049:4096)=xpsd(1:2048);
subplot(212)
plot(t,psdx);grid on;