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);