exa6_13.m

新版Matlab 7辅助信号处理技术与应用

% exa060302_Burg.m, for example 6.3.2 
% to estimate AR model parameters of x(n) by Burg function;

clear;
randn('seed',0)
a=[1 0.12 0.33 0.65 0.94];
xn=impz(1,a,20)+randn(20,1)/20;

% 计算AR模型的参数
Ba=arburg(xn,4)