exa121201_pcov.m

来自「《数字信号处理-理论、算法与实现》（第二版）源代码」· M 代码 · 共 18 行

%-------------------------------------------------------------------------
% exa121201_pcov.m, for example 12.12.1 
% to test pcov.m;
%-------------------------------------------------------------------------
clear all;

% 调出数据；
load test x;

N=4096;tn=-0.5:1/N:0.5-1/N;
% 使用协方差算法得到功率谱估计；
xpsd=pcov(x,10,N);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
plot(tn,fftshift(xpsd));grid on;