故障诊断实验——转子实验台振动信号的时域、频域分析
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
故障诊断第一次实验报告
——转子实验台振动信号的时域、频域分析1.由时域信号得出震动的平均值、方差、峭度信息
程序代码如下:
y=Data(1:6000,1);
y=y';
[a,b]=xcorr(y,'unbiased');
figure(1);
plot(b,a);
grid;
xlabel('位移信号自相关')
figure(2);
x=0:1/2560:1;
plot(x,mean(y));
gtext('平均值');
hold on;
figure(3);
plot(x,var(y));
gtext('方差');
hold on;
figure(4);
plot(x,kurtosis(y));
gtext('峭度');
位移信号
当转速为600r/min时
当转速为1000r/min时
当转速为1500r/min时
速度信号速度为600
速度1000
速度1500
1.由频域信号得出信号的傅里叶变换、功率谱密度信息
程序代码如下:
t=0:1/2560:1;
f=Data;
f=f';
y=fft(f,82944);
m=abs(y);
f1=(0:length(y)/2-1)'*2560/length(y); %计算变换后不同点对应的幅值figure(1);
plot(f1,m(1:length(y)/2));
ylabel('幅值的模');
xlabel('时间(s)');
title('原始信号傅里叶变换');
grid;
p=y.*conj(y)/82944; %计算功率谱密度
ff=(0:length(y)/2-1)'*2560/length(y); %计算变换后不同点对应的频率值figure(2);
plot(ff,p(1:length(y)/2));
ylabel('幅值');
xlabel('频率(Hz)');
title('功率谱密度');
grid;
位移信号速度600
速度1000
速度1500
速度信号速度600
速度1000
速度1500