marry0539 发表于 2011-4-25 09:04

时间--小波能量谱对吗???

最近在做轴承的故障诊断,数据用的是美国西储大学的实验数据,轴承是外圈故障,特征频率是103HZ,我对其做时间小波能量谱,实验结果不太理想,请各位大虾赐教 程序哪里有问题啊??
load 200
fs=10000;
N=2048*4;
t=(0:1/fs:(N-1)/fs);
x=X200_DE_time ;
data=x(1:N,1);
=ddencmp('den','wv',data);%%消噪
data1=wdencmp('gbl',data,'sym4',2, thr,sorh,keepapp);
c=cwt(data1,1:48,'db5');
cc=abs(c).^2;
sum=0;
for a=1:48;
sum=sum+cc(a,:)/(a.^2);
end ;
figure(1);
subplot(211);
plot(t,sum);    %时间--小波能量谱
ff=fft(sum,N);
ayy=abs(ff);
ayy=ayy/(N/2);
ayy(1)=ayy(1)/N;
rr=((1:N)-1)*fs/N;
subplot(212);
plot(rr(1:N/2),ayy(1:N/2));   %%时间--小波能量谱的频谱
axis()
结果如下:上面图形为时间小波能量谱
               下面图形为时间--小波能量谱的频谱

marry0539 发表于 2011-4-25 12:02

是不是计算时间--小波能量的时候不正确啊 ?请各位大虾赐教!困扰我好久了!谢谢

marry0539 发表于 2011-5-2 21:54

比较困扰我的是:时间--小波能量谱是不是
sum=0;
for a=1:48;
sum=sum+cc(a,:)/(a.^2);
end ;
这样子求法啊??

uestcer03 发表于 2011-6-24 16:16

回复 3 # marry0539 的帖子

个人觉得你这个求法是不对的!

yangfan502 发表于 2011-7-18 11:04

我最近也在用时间小波能量谱做轴承的故障诊断。同样做不出来啊
{:{13}:}{:{41}:}

yangfan502 发表于 2011-7-18 11:34

x=X200_DE_time ;
这个是什么意思啊?谢谢

ziyou2050 发表于 2013-4-2 21:34

sum=sum+cc(a,:)/(a.^2);
看不懂,sum=sum+cc(a,:)可以吗
页: [1]
查看完整版本: 时间--小波能量谱对吗???