gaojiepuindex1.m

内燃机转子仿真

%gaojiepu(shuangpu)
clear all;
clc;
format long;
close all;
load 'bneiranji1';
load 'cneiranji11d2e7';
%load 'arma1';
i=7;


xyz1=xyz1';
cxyz11d2e7=cxyz11d2e7';

by=xyz1(:,i)-mean(xyz1(:,i));
by=by/max(by)+randn(nn,1);
stdby=std(by);
%by=y;
cy=cxyz11d2e7(:,i)-mean(cxyz11d2e7(:,i));
cy=cy/max(cy)+randn(nn,1);
stdcy=std(cy);
%cy=y;
norder=3;

bma_order= maorder(by,0,16,0.05, '1');
bpmax=bma_order+1;
bqmax=bma_order+1;
bar_order= arorder(by,norder, bpmax,bqmax, '1');

cma_order= maorder(cy,0,16,0.05, '1');
cpmax=cma_order+1;
cqmax=cma_order+1;
car_order= arorder(cy,norder, cpmax,cqmax, '1');


bmaxlag=bar_order+bma_order+20;
cmaxlag=car_order+cma_order+20;
samp_seg=512;
%[bavec,bbvec] = armarts(by,ar_order,ma_order, norder,maxlag,samp_seg);
[bavec,bbvec] = armaqs(by,bar_order,bma_order, norder,bmaxlag,samp_seg);

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

bpq=max(bar_order,bma_order);
%bnfft=2^nextpow2(bpq);
bnfft=1024;
bbisp = bispect (ar, ma, bnfft);

ar=cavec;
ma=cbvec;
cpq=max(car_order,cma_order);
%cnfft=2^nextpow2(cpq);
cnfft=1024;
cbisp = bispect (ar, ma, cnfft);

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