exa8_5.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 25 行

% exa080403_arspect.m ,for example 8.4.3 
% to test arorder.m,arrcest.m,and bispect.m 

clear;
% 装载数据
load ar1;

% 估计AR模型阶次
p=arorder(y);

% 计算AR模型参数
avec=arrcest(y,p)

% 计算基于AR模型的双谱，并画出其平面图
ma=1;
ar=avec;
subplot(211)
bisp=bispect(ma,ar,64);
title('基于AR模型估计的双谱平面图')

% 画出基于AR模型的双谱三维图
abs_bisp=abs(bisp);
subplot(212)
mesh(abs_bisp);
title('基于AR模型估计的双谱三维图')