|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
s=X130_DE_time(2048:6143);
s1=s*2^7;
nfft=4096;
fs=12000;
s2=round(s1);
%s2=s1;
[cA,cD]=dwt(s2,'db13','mode','sym');
cD1=idwt([],cD,'db13');
cA1=idwt(cA,[],'db13');
y=abs(cD1);
y1=fft(y,4096);
y2=abs(y1);
figure(1);
plot((0:nfft/2-1)/nfft*fs,y2(1:nfft/2));
xlabel('频率 f/Hz');
ylabel('功率谱 P/W');
经过希尔伯特变换
我用此程序对凯斯西储大学轴承外圈数据进行处理,就是小波分解重构后取模,然后FFT,得到的故障频率也非常清楚(右图)。
我有一个地方不明白。一般的做法是小波分解重构之后还要进行希尔伯特变换取模后FFT,也能很清楚的看到故障频率(左图),第一种方法怎么解释y=abs(cD1)?求大神指导
|
|