gaojiepuindex.asv

内燃机转子仿真

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

by=by(:,i)-mean(by(:,i));

cy=iy4;
cy=aaii4(:,i)-mean(aaii4(:,i));

ar_order=5;
ma_order=3;
norder=3;
maxlag=18;
samp_seg=128;
%[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, 512);

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

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