振动信号的频谱分析

zjheme

对齿轮箱故障数据分析，采样频率在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('功率');

michaelpeng

:@Q :@L 难啊 不会

zjheme

高手给点意见

华电机械

运行了你的程序，怎么出现错误啊
Undefined function or variable 't'.

Error in fft1 (line 4)
ttrain=t;

白兔糖

华电机械 发表于 2013-5-8 09:43

                               
登录/注册后可看大图

运行了你的程序，怎么出现错误啊
Undefined function or variable 't'.

正常信号 频谱分析 该怎么分析啊 有什么特点？

雅文

围观一下

PXMAX

白兔糖 发表于 2013-5-28 15:07

                               
登录/注册后可看大图

正常信号 频谱分析 该怎么分析啊 有什么特点？

不同状态下的波形不同，自然频谱分析的结果也不同，根据你得到的波形，通过频率的分布，观测特征频率，对特征频率进行比对，可以推测出问题的所在，例如：转速为1200rpm的旋转体，它的失衡特征频率为20HZ，那么观察20HZ这个区间的频谱，就可以知道它是否失衡。

hucanboy1

个人认为只要分析小频率的频谱。楼主如果不介意，可以分享一下你的testdata.mat.呀

cdut2013

{:{03}:}{:{03}:}