19-4.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函数估计的双谱三维图') [avec,bvec]=armaqs(y,2,1,3,10,128);