gaojiepufenxi.m

内燃机转子仿真

%gaojiepu(shuangpu)
clear all;
clc;
format long;
close all;
addpath('D:\Program Files\MATLAB71\toolbox\hosa');
which 'bispecd';
load 'b2neiranji';
load 'b2neiranji1';
load 'b2neiranji2';
load 'b2neiranji3';
load 'b2neiranji4';
load 'b2neiranji5';
ip=10;

nfft=256;
wind=5;
nsamp=16;
overlap=30;
dl=20000;
dh=40000;
%n=length(xyz(7,140000:150000));
gxyz=xyz(ip,dl:dh);
gxyz1=xyz1(ip,dl:dh);
gxyz2=xyz2(ip,dl:dh);
gxyz3=xyz3(ip,dl:dh);
gxyz4=xyz4(ip,dl:dh);
gxyz5=xyz5(ip,dl:dh);


gxyz_mean=gxyz-mean(gxyz);
gxyz_std=std(gxyz_mean);
gy=gxyz_mean/gxyz_std;



gxyz1_mean=gxyz1-mean(gxyz1);
gxyz1_std=std(gxyz1_mean);
gy1=gxyz1_mean/gxyz1_std;



gxyz2_mean=gxyz2-mean(gxyz2);
gxyz2_std=std(gxyz2_mean);
gy2=gxyz2_mean/gxyz2_std;



gxyz3_mean=gxyz3-mean(gxyz3);
gxyz3_std=std(gxyz3_mean);
gy3=gxyz3_mean/gxyz3_std;



gxyz4_mean=gxyz4-mean(gxyz4);
gxyz4_std=std(gxyz4_mean);
gy4=gxyz4_mean/gxyz4_std;


gxyz5_mean=gxyz5-mean(gxyz5);
gxyz5_std=std(gxyz5_mean);
gy5=gxyz5_mean/gxyz5_std;



[bBspec,bwaxis] = bispecd (gy,  nfft, wind, nsamp, overlap);
bBspec_abs=abs(bBspec);
[bic,waxis] = bicoher (gy,  nfft, wind, nsamp, overlap);

[bBspec1,bwaxis1] = bispecd (gy1,  nfft, wind, nsamp, overlap);
bBspec1_abs=abs(bBspec1);
[bic1,waxis1] = bicoher (gy1,  nfft, wind, nsamp, overlap);

[bBspec2,bwaxis2] = bispecd (gy2,  nfft, wind, nsamp, overlap);
bBspec2_abs=abs(bBspec2);
[bic2,waxis2] = bicoher (gy2,  nfft, wind, nsamp, overlap);

[bBspec3,bwaxis3] = bispecd (gy3,  nfft, wind, nsamp, overlap);
bBspec3_abs=abs(bBspec3);
[bic3,waxis3] = bicoher (gy3,  nfft, wind, nsamp, overlap);

[bBspec4,bwaxis4] = bispecd (gy4,  nfft, wind, nsamp, overlap);
bBspec4_abs=abs(bBspec4);
[bic4,waxis4] = bicoher (gy4,  nfft, wind, nsamp, overlap);

[bBspec5,bwaxis5] = bispecd (gy5,  nfft, wind, nsamp, overlap);
bBspec5_abs=abs(bBspec5);
[bic5,waxis5] = bicoher (gy5,  nfft, wind, nsamp, overlap);


%for j=1:nfft
    % bBspec(j,:)= bBspec(j,:)/m3Bspec;
   % cBspec(j,:)= cBspec(j,:)/m3Bspec;
%end
tu=0;
if tu==1
figure(1)
subplot(3,2,1)
mesh(bic);

title('无裂纹第9轴段相对扭角归一化双谱');

subplot(3,2,2)
mesh(bic1);

title('e=0.1第9轴段相对扭角归一化双谱');
subplot(3,2,3)
mesh(bic2);

title('e=0.2第9轴段相对扭角归一化双谱');
subplot(3,2,4)
mesh(bic3);

title('e=0.3第9轴段相对扭角归一化双谱');
subplot(3,2,5)

mesh(bic4);

title('e=0.4第9轴段相对扭角归一化双谱');
subplot(3,2,6)

mesh(bic5);

title('e=0.45第9轴段相对扭角归一化双谱');

figure(2);
subplot(3,2,1);
mesh(bBspec_abs);

title('无裂纹第9轴段相对扭角双谱');
%subplot(2,1,2);
%surf(bf1,bf2,bpBspec);
%zlim([0,1]*10e-10);

subplot(3,2,2);
mesh(bBspec1_abs);

title('e=0.1第9轴段相对扭角双谱');
%subplot(2,1,2);
%surf(cf1,cf2,cpBspec);
%zlim([0,1]*10e-10);
subplot(3,2,3);
mesh(bBspec2_abs);

title('e=0.2第9轴段相对扭角双谱');

subplot(3,2,4);
mesh(bBspec3_abs);

title('e=0.3第9轴段相对扭角双谱');

subplot(3,2,5);
mesh(bBspec4_abs);

title('e=0.4第9轴段相对扭角双谱');

subplot(3,2,6);
mesh(bBspec5_abs);

title('e=0.45第9轴段相对扭角双谱');
end