myautocorr.m

matlab 时频分析上面程序第八章的故障诊断的小程序

clear all
load bbb.m
%load myshuju02
%B=yy(4,:);

%load 齿轮台实验720_1_Waveform;
%B=a8;

%load mydata4xx1;
%B=xx1(3,:);
%B=0.2*xx1(1,:)+0.1*xx1(2,:)+xx1(3,:);

%load zh600_3_2_Waveform;
%B=a5(1:5000);

%load zh600_3_2_Waveform;
%%%%%%%%%%%%%%%%%%%%%%%%%%
%load mydatasvd e;
% load xie720-3Mpa-5_Waveform;
% x=a8(5001:11000);
x=bbb;
B=x-mean(x);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%$$load mydataxx1;
%$$B=0.5*xx1(1,:)+0.3*xx1(2,:)+xx1(3,:);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%load myxiedata41;
%B=xiedata41(4,:);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[ACF,Lags,Bounds]=autocorr(B,100,[],2);
figure(1);
subplot(211); 
plot(Lags,ACF);
grid;
%subplot(212);
%plot(e);
%grid;