exa110701_pwelch.m

笔者在前人研究的基础上

%-------------------------------------------------------------------------
% exa1100701_pwelch.m, for example 11.7.1 and fig. 11.5.1(d)
% to estimate the PSD of x(n) by Welch average;
%--------------------------------------------------------------------------
clear;

% 调出数据；
load test x;
N=256;
Fn=-0.5:1/N:0.5-1/N;

plot(real(x));
figure;
plot(imag(x));

% 用 Welch 平均估计试验数据的功率谱；
xpsd=pwelch(x,hamming(33),16,N,'whole');        
figure;
plot(xpsd);