|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
大家好!我用HHT处理http://www.eecs.case.edu/laboratory/bearing/download.htm上的内圈故障数据,发现很多问题
具体程序如下:
close all
clear all
fs=12000;
dd=500;%新数据序列截取原数据中间一段,两端各500点防止边端效应
r=10;%重采样倍率
fid=fopen('inner_14.mat','r');
x=fread(fid);
status=fclose(fid);
s=decimate(x,r);
fs=fs/r;
N=length(s);
imf=emd(s);
emd_visu(s(1+dd:end-dd),1:length(s(1+dd:end-dd)),imf(:,1+dd:end-dd));
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A(:,dd:end-dd+1),f(:,dd:end-dd+1),tt(:,dd:end-dd+1),length(tt(:,dd:end-dd+1)));
f=(0:N-2*dd-1)/N*(fs/2);
tt=tt-dd-1;
figure
imagesc(tt,f,im);
set(gca,'YDir','normal')
colorbar
xlabel('时间(点数)','FontSize',12);
ylabel('频率(Hz)','FontSize',12);
IMF貌似很正常,为什么时频图作出来密布很多点,不知道错在哪里,请高手指教,谢谢!故障数据集见附件
大家一般是如何处理采样频率很高、数据量很大的数据的呢? |
|