声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4626|回复: 15

[综合] 斜拉桥实测信号去噪处理_平均法、小波变换、SVD去噪

  [复制链接]
发表于 2010-12-5 20:27 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 secondye 于 2010-12-5 20:32 编辑

        最近在对一斜拉桥主梁的振动信号进行去噪处理,尝试过五点三次滑动平均法、小波变换去噪、SVD去噪,但好像效果都不怎么好,去噪后的FFT图和去噪前FFT图的几乎一样。

        五点三次滑动平均法,就是用王济老师的《Matlab在振动信号处理中的应用》这本书提供的代码,也尝试过matlab的“smooth”函数,好像都是一样的吧;我先对信号去趋势项detrend、再smooth、然后FFT。但平不平均,好像对于噪声没什么作用,光是从FFT图上,不能识别真实的峰值。
                原始信号                                                                           
signal-smooth.jpg FFT FFT-smooth.jpg

        至于小波变换去噪,本人也是刚接触,看过有软阀值、硬阀值、缺省值、抑制细节等去噪方法,但不知道该用哪种方法好,看网上一些代码提供的都是仿真的信号,去噪结果似乎都很好。但是我用之对斜拉桥主梁实测信号进行处理,看有人做过,先对信号去均值、然后用Battle-Lemarie正交小波进行5层分解,再对低频部分进行单支重构。
        我照着做了,但没找到Battle-Lemarie正交小波,只能用db4来代替,但好像去噪前后的FFT图都不能很好地识别真实的模态峰值;而且还出现了去噪前没有出现的虚假模态,有限元计算的频率值为0.337HZ,0.431HZ,0.59HZ,而小波去噪后还出现了0.16HZ,0.26HZ,而且峰值很明显。
                             原始信号                                                                     
signal-wavelet.jpg FFT FFT-wavelet.jpg

        SVD去噪,也是最近看东方振动与噪声研究所的一篇论文才知道的,昨天在网上搜来一段程序,该程序去噪有选取阈值为特征值的平均值,也有选择阈值为特征值的中值,但去噪效果也不怎么好。SVD去噪昨天刚看,不怎么熟悉,请大家指教下
                              原始信号                                                                           
signal-SVD.jpg FFT FFT SVD.jpg

         denoise program.rar (33.6 KB, 下载次数: 178)   把数据和三种去噪的代码都贴上了,采样频率为20HZ,欢迎各位高手过来指教。小弟感激不尽
FFT-wavelet.jpg
回复
分享到:

使用道具 举报

发表于 2010-12-7 09:04 | 显示全部楼层
不知道你研究过EMD吗?好像有用EMD进行去噪的 没有深入研究过 你可以试试
 楼主| 发表于 2010-12-7 14:45 | 显示全部楼层
回复 2 # hitdely 的帖子

哦,有看到过。
好像也是像小波变换那样,先分解再重构的。

对于各种去噪的方法,有什么体会呢?说来跟大家分享一下
发表于 2010-12-7 19:21 | 显示全部楼层
回复 3 # secondye 的帖子

我没做这方面 我做的是hht的
发表于 2010-12-8 21:58 | 显示全部楼层
EMD的效果会更好,之前有一篇《对一实测信号的处理》上有谈到楼主的问题。还有楼主可以试试用vondark滤波试试。
发表于 2011-5-6 15:01 | 显示全部楼层
发表于 2011-5-9 10:48 | 显示全部楼层
本帖最后由 zhumeichen 于 2011-5-9 10:50 编辑

EMD效果是比较好,重大的汤宝平老师研究的很深,最近我也在做这方面的,可刚入手不久。认识很肤浅
发表于 2012-4-13 22:14 | 显示全部楼层
回复 1 # secondye 的帖子

个人认为利用SVD去噪的原理是认为后几阶奇异值对应到的就是噪声张成空间的特征值,在恢复信号的时候去掉那几个奇异值就行。
发表于 2012-4-14 01:22 | 显示全部楼层
smooth不是滤时域的,是光滑频域的,调用方法是x=smooth(x,3);频域显示纵坐标用对数坐标
发表于 2012-7-15 10:37 | 显示全部楼层
关注,求高手指点
发表于 2012-7-17 21:59 | 显示全部楼层
本帖最后由 吃书的老虎 于 2012-7-17 23:15 编辑

这是利用EMD分解到的IMF,再将包含0.337HZ,0.431HZ,0.59HZ的IMF合并时域图,及其FFT。不知是否符合你的要求。不过0.59频率成分几乎看不到。即使原始图的FFT也几乎看不到。
lv.jpg
ori.jpg
发表于 2012-9-29 19:50 | 显示全部楼层
最近也在学习去噪
发表于 2014-4-17 21:04 | 显示全部楼层
学习了,谢谢楼主
发表于 2014-5-7 15:28 | 显示全部楼层
谢谢你 正好需要  一起交流下
发表于 2014-5-11 11:56 | 显示全部楼层
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-5-26 23:30 , Processed in 0.080164 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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