exa110701_pwelch.m

来自「数字信号处理书中的程序,程序序分别由FORTRAN语言和C语言编写.」· M 代码 · 共 20 行

%-------------------------------------------------------------------------
% 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=4096;
Fn=-0.5:1/N:0.5-1/N;

% 用 Welch 平均估计试验数据的功率谱；
% xpsd=pwelch(x,N,1,hamming(33),16,'whole');     for MATLAB 5.3
xpsd=pwelch(x,hamming(33),16,N,'whole');        %for MATLAB 6.1
mmax=max(xpsd);
xpsd=xpsd/mmax;
xpsd=10*log10(xpsd+0.000001);
plot(Fn,fftshift(xpsd));grid on;