声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2295|回复: 6

[滤波] 谐波小波滤波

[复制链接]
发表于 2009-6-4 23:22 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我在使用谐波小波对信号进行滤波时出现了一个奇怪的问题,
现将方法域代码贴出,请高手指点,不胜感激。

使用的图片和公式在 2 楼

方法一:李舜酩教授《振动信号的现代分析技术与应用》一书中方法
图1
谐波小波从频域返回到时域取共轭再到频域滤波

方法二:何正嘉《机械设备非平稳信号的故障诊断原理及应用》书中方法
图2
直接在频域利用谐波小波进行滤波


实现代码如下: 其中delg为欲滤波信号

M=length(delg);  %信号长度
ggg=fft(delg);
for i=1:M;
   fai(i)=1/0.006;
   if (i/M>0.006)
      fai(i)=0;
   end
end



%利用方法1
kkk=ifft(fai);
kkk=kkk';
xx=conj(kkk);
xxx=fft(xx);
yy=ggg.*xxx;
aa=ifft(yy);
figure1
plot(aa)

%使用方法2
yy1=ggg.*fai';
figure2
plot(ifft(yy1))



最后返回时域值为复数,不知道为什么?
图形如下:
figure1 :方法一的图形
图3
figure2  方法二得到的图形
图4

如上所示得到的数据为复数,敬请高手指点,谢谢。


[ 本帖最后由 yanftian 于 2009-6-4 23:31 编辑 ]
回复
分享到:

使用道具 举报

 楼主| 发表于 2009-6-4 23:27 | 显示全部楼层
图片刚才没发上去,不好意思

图1

图1

图2

图2

图3

图3

图4

图4
发表于 2009-8-26 15:13 | 显示全部楼层
两个图是一样的吗,怎么判断出故障呢,请明示啊,谢谢了
发表于 2009-8-31 20:50 | 显示全部楼层

回复 楼主 yanftian 的帖子

aa是复数 作图取摸
abs(aa)
发表于 2010-1-30 10:07 | 显示全部楼层
没看明白:@L
发表于 2010-2-1 11:35 | 显示全部楼层
滤波的时域结果一般是复数序列的,只要取实数部分的一半即可。
你的代码就是相当于把振动信号用理想型的低通滤波器滤波了,效果不好。
发表于 2013-10-28 21:45 | 显示全部楼层
楼上说得对,不过和普通低通滤波器不同的是,谐波小波在频域的虚部也有取值,因此不能直接用ones函数构造,正逆傅立叶变换的目的正是获得谐波小波的复包络
而当有了复包络以后,就是简单的两个复数的频域相乘了

令,楼主这个问题,考虑到正负频率特性,李老师的是更严格的说法
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 00:48 , Processed in 0.072856 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表