ar.m

matlab环境下用AR谱估计算法实现AR模型谱估计

Fs = 1000 ;
N = 1024 ;
n = 0 : N-1 ;
t = n / Fs ;
xn = 2 * sin( 2 * pi * 100 * t ) + 6 * sin( 2 * pi * 105 * t ) ;
rn = randn ( 1 , N ) ;
x = xn + rn ;
[ P , f ] = pburg ( x , 100 , N , Fs) ;
Px = 10 * log10( P ) ;
figure
plot ( f , Px);
xlabel ( ' 频率 / Hz ' ) ;
ylabel ( ' 功率谱 / dB '  ) ;
title ( 'Burg法AR谱估计' ) ;
grid on ;