exa121201_pburg.m

来自「MATLAB程序关于数字信号处理的算法实现」· M 代码 · 共 19 行

%-------------------------------------------------------------------------
%exa121201_pburg.m, for example 12.12.1 
%to test pburg.m;
%-------------------------------------------------------------------------
clear all;

% 调出数据；
load test x;

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