利用小波变换检测轴承滚动体
轴承故障为滚动体局部剥落。轴承故障特征频率为外环通过频率及其倍频，理论值为219HZ%采样频率fs=10000;
%轴承滚动体故障信号
fid=fopen('bearingroll.dat','r');%故障
N=1024;
xdata=fread(fid,N,'int16');
fclose(fid);
xdata=(xdata-mean(xdata))/std(xdata,1);
%时域波形
figure(1);
plot(1:N,xdata);
xlabel('样本序号 n/个');
ylabel('电压 V/v');
%db10小波进行4层分解
%一维小波分解
= wavedec(xdata,4,'db10');
%重构第1-4层细节系数
d4 = wrcoef('d',c,l,'db10',4);
d3 = wrcoef('d',c,l,'db10',3);
d2 = wrcoef('d',c,l,'db10',2);
d1 = wrcoef('d',c,l,'db10',1);
%显示细节信号
figure(2)
subplot(4,1,1);
plot(d4,'LineWidth',2);
ylabel('d4');
subplot(4,1,2);
plot(d3,'LineWidth',2);
ylabel('d3');
subplot(4,1,3);
plot(d2,'LineWidth',2);
ylabel('d2');
subplot(4,1,4);
plot(d1,'LineWidth',2);
ylabel('d1');
xlabel('样本序号 n/个');
%第1层细节信号的包络谱
y=hilbert(d1);
ydata=abs(y);
y=y-mean(y);
nfft=1024;
p=abs(fft(ydata,nfft));
figure(3);
plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));
xlabel('频率 f/Hz');
ylabel('功率谱 P/W');
如果想将小波和双谱分析结合起来，该怎么做呢 'bearingroll.dat'
请楼主把这个文件传一传吧。
页:
[1]