load D:\105.mat;%内圈故障数据
x=X105_DE_time;%驱动计数段的内圈故障
fs=12000;%采样率
N=2048;%采样点数(100倍)
M=0;%采样数据段的起始位置
n=M:N-1;
t=n/fs;%信号时间序列
X=X105_DE_time(1:N);%装载 驱动计数端的内圈故障数据
%X=X107_FE_time(1:N);%装载 风扇计数端的内圈故障数据
%X=X107_BA_time(1:N);%装载 基础计数端的内圈故障数据
%X=X107_DE_time(1:N)-X107_BA_time(1:N);
y=X';%信号幅值序列
k_in=kurtosis(y);%峭度系数,正常轴承为3左右
figure(1)
plot(t,y);
ylabel('振幅 m/v');
xlabel('时间 t/s');
t=wpdec(y,3,'db10');
figure(2)
plot(t);
s130=wprcoef(t,[3,0]);
s131=wprcoef(t,[3,1]);
s132=wprcoef(t,[3,2]);
s133=wprcoef(t,[3,3]);
s134=wprcoef(t,[3,4]);
s135=wprcoef(t,[3,5]);
s136=wprcoef(t,[3,6]);
s137=wprcoef(t,[3,7]);<