exa8_7.m

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

% exa080405_armaspect.m ,for example 8.4.5 
% to test armarts.m,armaqs.m,and bispect.m 

clear;
% 装载数据
load arma1;

% 计算ARMA模型参数
ar_order=2;
ma_order=1;
norder=3;
maxlag=10;
samp_seg=128;
% 利用armarts函数计算ARMA模型参数
[bvec_rts avec_rts]=armarts(y,ar_order,ma_order,norder,maxlag,samp_seg)
% 利用armaqs函数计算ARMA模型参数
[bvec_qs avec_qs]=armaqs(y,ar_order,ma_order,norder,maxlag,samp_seg)

% 计算参数由armarts函数估计的ARMA模型的双谱，并画出该谱平面图
ma=bvec_rts;
ar=avec_rts;
subplot(221)
bisp_rts=bispect(ma,ar,128);
title('armarts函数估计的双谱平面图')

% 计算参数由armaqs函数估计的ARMA模型的双谱，并画出该谱平面图
ma=bvec_qs;
ar=avec_qs;
subplot(222)
bisp_qs=bispect(ma,ar,128);
title('armaqs函数估计的双谱平面图')

% 画出参数由armarts函数估计的ARMA模型的双谱三维图
abs_bisp=abs(bisp_rts);
subplot(223)
mesh(abs_bisp);
title('armarts函数估计的双谱三维图')

% 画出参数由armaqs函数估计的ARMA模型的双谱三维图
abs_bisp=abs(bisp_qs);
subplot(224)
mesh(abs_bisp);
title('armaqs函数估计的双谱三维图')