齿轮故障诊断的程序和例程

用高阶统计量中的三阶累积量的切片谱实现对齿轮故障诊断的程序和例程.

clear all
close all
fs=1024
K=64;  %段数
nfft=512;
f=fs*（0:nfft/2-1）/nfft; %归一化频率
 
load E:\datadecompose\LMD\gear\gearlw\CQLB11c1.dat
 
LB111=[CQLB11c1CQLB11c1CQLB11c1CQLB11c1CQLB11c1CQLB11c1CQLB11c1CQLB11c1];
% LB11=[CQLB11c1]
 x=LB111;
%  xx=LB11
%  ft（xx1024）
% %  
%  [pFXbl]=myblp（xx1024）
% % x=FXbl;
N=length（x）
  
M=fix（N/K） % 每段的数据长度  每段数据长128 
for i=0:K-1
y=x（i*M+1:（i+1）*M）;
y=y-mean（y）;
xk（i+1:）=y;
end
%计算累积量
 
for t=-（nfft/2-1）:nfft/2     
c1（t+nfft/2）=third_slice_cumulant（xkt）;  %三阶累积量对角切片
end
 
% FFT变换
% subplot（211） 
C1=c1.*hann（length（c1））‘;  
Y=fft（C1nfft）;   %对三阶累积量对角切片求FFT
P1=Y.*conj（Y）/nfft;   %%%  P1
% f=fs*（0:nfft/2-1）/nfft; %归一化频率
% plot（fP1（1:nfft/2））;    %画1.5维谱图
% xlabel（‘频率‘）
% ylabel（‘幅值‘）
 
 
load E:\datadecompose\LMD\gear\gearlw\CQLB11c2.dat
LB112=[CQLB11c2CQLB11c2CQLB11c2CQLB11c2CQLB11c2CQLB11c2CQLB11c2CQLB11c2];
% LB12=[CQLB11c2];
 
x=LB112;
% N=length（x）
% K=64      %段数
% M=fix（N/K） % 每段的数据长度  每段数据长128 
 
for i=0:K-1
y=x（i*M+1:（i+1）*M）;
y=y-mean（y）;
xk（i+1:）=y;
end
         
%计算累积量
 
for t=-（nfft/2-1）:nfft/2     
c2（t+nfft/2）=third_slice_cumulant（xkt）;  %三阶累积量对角切片
end
 
% FFT变换
% subplot（211） 
C2=c2.*hann（length（c2））‘;  
Y=fft（C2nfft）;  %对三阶累积量对角切片求FFT
P2=Y.*conj（Y）/nfft;   %%%  P1
% f=fs*（0:nfft/2-1）/nfft; %归一化频率
% plot（fP2（1:nfft/2））;    %画1.5维谱图
% xlabel（‘频率‘）
% ylabel（‘幅值‘）
 
 
load E:\datadecompose\LMD\gear\gearlw\CQLB11c3.dat
LB113=[CQLB11c3CQLB11c3CQLB11c3CQLB11c3CQLB11c3CQLB11c3CQLB11c3CQLB11c3];
% LB13=[CQLB11c3];
x=LB113;
% N=length（x）
% K=64      %段数
% M=fix（N/K） % 每段的数据长度  每段数据长128 
 
for i=0:K-1
y=x（i*M+1:（i+1）*M）;
y=y-mean（y）;
xk（i+1:）=y;
end
         
%计算累积量
for t=-（nfft/2-1）:nfft/2     
c3（t+nfft/2）=third_slice_cumulant（xkt）;  %三阶累积量对角切片
end
 
% FFT变换
% subplot（211） 
C3=c3.*hann（length（c3））‘;  
Y=fft（C3nfft）;   %对三阶累积量对角切片求FFT
P3=Y.*conj（Y）/nfft;   %%%  P1
% f=fs*（0:nfft/2-1）/nfft; %归一化频率
% plot（fP3（1:nfft/2））;    %画1.5维谱图
% xlabel（‘频率‘）
% ylabel（‘幅值‘）
 
 
 
load E:\datadecompose\LMD\gear\gearlw\CQLB11c4.dat
LB114=[CQLB11c4CQLB11c4CQLB11c4CQLB11c4CQLB11c4CQLB11c4CQLB11c4CQLB11c4];
% LB14=[CQLB11c4];
 
x=LB114;
% N=length（x）
% K=64      %段数
% M=fix（N/K） % 每段的数据长度  每段数据长128 
 
for i=0:K-1
y=x（i*M+1:（i+1）*M）;
y=y-mean（y）;
xk（i+1:）=y;
end
         
%计算累积量
 
for t=-（nfft/2-1）:nfft/2     
c4（t+nfft/2）=third_slice_cumulant（xkt）;  %三阶累积量对角切片
end
 
% FFT变换
% subplot（211） 
C4=c4.*hann（length（c4））‘;  
Y=fft（C4nfft）;   %对三阶累积量对角切片求FFT
P4=Y.*conj（Y）/nfft;   %%%  P1
 
X=[c1;c2;c3;c4]
SLB11=svd（X）;


% clear all
% close all
% fs=1024
% K=64;  %段数
% nfft=512;
% f=fs*（0:nfft/2-1）/nfft; %归一化频率
 
load E:\datadecompose\LMD\gear\gearlw\CQLB12c1.dat
 
LB121=[CQLB12c1CQLB12c1CQLB12c1CQLB12c1CQLB12c1CQLB12c1CQLB12c1CQLB12c1];
% LB12=[CQLB12c1]
 x=LB121;
%  xx=LB12
%  ft（xx1024）
% %  
%  [pFXbl]=myblp（xx1024）
% % x

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      42213  2013-04-11 11:10  hospeciPRG\QDqiepian.m

     文件        683  2013-04-11 11:39  hospeciPRG\SQD.mat

     文件      22174  2013-04-11 12:46  hospeciPRG\ZCqiepian.m

     文件        211  2013-03-29 21:45  hospeciPRG\third_slice_cumulant.m

     文件        619  2013-04-11 11:39  hospeciPRG\SLB.mat

     文件        415  2013-04-11 14:35  hospeciPRG\SZC.mat

     文件      37202  2013-04-11 11:39  hospeciPRG\LBqiepian.m

     目录          0  2013-04-11 13:24  hospeciPRG

----------- ---------  ---------- -----  ----

               103517                    8