armoxing.m

内燃机转子仿真

clear all;
clc;
format long;
close all;
load 'arma1';
ar_order=2;
ma_order=1;
norder=3;
maxlag=10;
samp_seg=128;
[avec_rts,bvec_rts] = armarts(y,ar_order,ma_order, norder,maxlag,samp_seg);
[avec_qs, bvec_qs] = armaqs(y,ar_order,ma_order, norder,maxlag,samp_seg);
ar=avec_rts;
ma=bvec_rts;
subplot(2,2,1);
bisp_rts=bispect(ar,ma,128);
title('armarts函数估计的双谱平面图');
ar=avec_qs;
ma=bvec_qs;
subplot(2,2,2);
bisp_qs=bispect(ar,ma,128);
title('armarts函数估计的双谱平面图');
abs_bisp=abs(bisp_rts);
subplot(2,2,3);
mesh(abs_bisp);
title('armarts')
abs_bisp=abs(bisp_qs);
subplot(2,2,4);
mesh(abs_bisp);
[avec,bvec]=armaqs(y,2,1,3,10,128);