马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
对齿轮箱故障数据分析,采样频率在8192,齿轮转速在700RPM,齿轮啮合频率在500HZ,得到得频率和频谱图如下:
不知道为什么在频率为100HZ左右的振幅特别大,在四中状态下都存在,不知道这是啥原因该如何去除。打算对信号小波去噪处理,但不知道阀值如何设?
附上代码,希望大家给点意见。
clear all;
clc;
load testdata.mat;
ttrain=t;
% train1=normal;
% train2=unbalance;
% train3=loosness;
% train4=misalignment;
fs=8192;N=8192;
n=0:N-1;t1=n/fs;
y1=fft(train1,N);
y2=fft(train2,N);
y3=fft(train3,N);
y4=fft(train4,N);
mag1=abs(y1);
mag2=abs(y2);
mag3=abs(y3);
mag4=abs(y4);
f=n*fs/N;
figure(1);
subplot(4,3,1);
plot(ttrain,train1);
axis([0 1 -0.2 0.2]);
ylabel('正常信号');
subplot(4,3,2);
plot(f(1:N/2),mag1(1:N/2));
py1=y1.*conj(y1)/N;
subplot(4,3,3);
plot(f(1:N/2),py1(1:N/2));
subplot(4,3,4);
plot(ttrain,train2);
axis([0 1 -0.2 0.2]);
ylabel('轴不平衡');
subplot(4,3,5);
plot(f(1:N/2),mag2(1:N/2));
py2=y2.*conj(y2)/N;
subplot(4,3,6);
plot(f(1:N/2),py2(1:N/2));
subplot(4,3,7);
plot(ttrain,train3);
axis([0 1 -0.2 0.2]);
ylabel('螺钉松动');
subplot(4,3,8);
plot(f(1:N/2),mag3(1:N/2));
py3=y3.*conj(y3)/N;
subplot(4,3,9);
plot(f(1:N/2),py3(1:N/2));
subplot(4,3,10);
axis([0 1 -0.2 0.2]);
plot(ttrain,train4);
ylabel('轴不对中');
xlabel('时间');
subplot(4,3,11);
plot(f(1:N/2),mag4(1:N/2));
xlabel('频率');
py4=y4.*conj(y4)/N;
subplot(4,3,12);
plot(f(1:N/2),py4(1:N/2));
xlabel('功率'); |