gjpzxg.m

内燃机转子仿真

%gaojiepu(shuangpu)
clear all;
clc;
format long;
close all;
addpath('D:\Program Files\MATLAB71\toolbox\hosa');
which 'bispecd';
load 'bneiranji';
load 'bneiranji1';
load 'bneiranji2';
load 'bneiranji3';
load 'bneiranji4';
ip=7;

nfft=256;
wind=1;
nsamp=128;
overlap=30;
%n=length(xyz(7,140000:150000));
gxyz=xyz(ip,:);
gxyz1=xyz1(ip,:);
gxyz2=xyz2(ip,:);
gxyz3=xyz3(ip,:);
gxyz4=xyz4(ip,:);






[bic,waxis] = bicoher (gxyz,  nfft, wind, nsamp, overlap);


[bic1,waxis1] = bicoher (gxyz1,  nfft, wind, nsamp, overlap);


[bic2,waxis2] = bicoher (gxyz2,  nfft, wind, nsamp, overlap);


[bic3,waxis3] = bicoher (gxyz3,  nfft, wind, nsamp, overlap);


[bic4,waxis4] = bicoher (gxyz4,  nfft, wind, nsamp, overlap);


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

mesh(bic4);



end