gaojiepuindex.m

内燃机转子仿真

%gaojiepu(shuangpu)
clear all;
clc;
format long;
close all;
%load 'bneiranji';
%load 'bneiranji4';
load 'hfneiranji';
%i=7;
by=real(iy1);

by=by-mean(by);

cy=real(iy4);
cy=cy-mean(cy);

ar_order=2;
ma_order=1;
norder=3;
maxlag=18;
samp_seg=1024;
%[bavec,bbvec] = armarts(by,ar_order,ma_order, norder,maxlag,samp_seg);
[bavec,bbvec] = armaqs(by,ar_order,ma_order, norder,maxlag,samp_seg);

%[cavec,cbvec] = armarts(cy,ar_order,ma_order, norder,maxlag,samp_seg);
[cavec,cbvec] = armaqs(cy,ar_order,ma_order, norder,maxlag,samp_seg);
ar=bavec;
ma=bbvec;

bbisp = bispect (ar, ma, 2048);

ar=cavec;
ma=cbvec;
cbisp = bispect (ar, ma, 2048);

%for j=1:nfft
    % bBspec(j,:)= bBspec(j,:)/m3Bspec;
   % cBspec(j,:)= cBspec(j,:)/m3Bspec;
%end

abs_bbisp=abs(bbisp);
subplot(2,1,1);
%mesh(bwaxis*(1/h),bwaxis*(1/h),bBspec);

mesh(abs_bbisp);
abs_cbisp=abs(cbisp);
subplot(2,1,2);
mesh(abs_cbisp);
%mesh(cwaxis*(1/h),cwaxis*(1/h),cBspec);
%subplot(2,1,2);