exa121202.m

来自「功率谱估计的一些matlab例程」· M 代码 · 共 42 行

%-------------------------------------------------------------------------
% exa121202.m, for example 12.12.2;
% to test aryule,arburg,arcov,armcov;
%-------------------------------------------------------------------------
clear all;
load test x;

% 用自相关法求AR模型的系数和最小预测误差能量；
[a,E1]=aryule(x,10);
figure(1)
subplot(411)
stem(abs(a),'.')
ylabel(' for aryule');

% 用Burg法求AR模型的系数、反射系数和最小预测误差能量；
[a,E2,k]=arburg(x,10);
subplot(412)
stem(abs(a),'.')
ylabel(' for arburg');

% 用协方差法求AR模型的系数和最小预测误差能量；
[a,E3]=arcov(x,10);
subplot(413)
stem(abs(a),'.')
ylabel(' for cov');

% 用改进的协方差法求AR模型的系数和最小预测误差能量；
[a,E4]=armcov(x,10);
subplot(414)
stem(abs(a),'.')
ylabel(' for armcov');

figure(2)
stem(abs(k),'.');
ylabel(' for arburg');

E1
E2
E3
E4