yghit08
发表于 2012-6-16 17:07
提供一个实测信号的的处理结果试验件:980*50*9mm的悬臂梁(钢)
锤击实验得到的实验数据
处理结果见下图,低阶频率分离的还可以(高阶我看了,一般,几乎处于没有被分离出来的状态)。大家可以利用我提供的程序处理,我也用的是那个。欢迎大家交流自己写的程序,我求移动平均法的处理程序和有理样条的处理程序
yghit08
发表于 2012-6-16 19:13
同一组数据,利用EMD做出来的结果,貌似要比LMD做的要好。大家看看啊
雨夜1011
发表于 2012-6-16 21:48
我用的是labvIEW来实现的,但是实现的效果真的不咋地~~~也不知道哪里有问题
解析表达式为
雨夜1011
发表于 2012-6-16 21:53
分解出来后,给我的感觉就最后的一个幅值是对的~~我估计第二个的信号特征出现了泄露
yghit08
发表于 2012-6-17 12:06
回复 19 # 雨夜1011 的帖子
你实现的结果要比我做的要好。你用的是移动平均法做的平滑处理?有程序或者伪代码否?
昨天我给LMD的作者发了邮件,针对楼上的实测数据做的(EMD和LMD),他肯定了我的程序处理结果没什么问题(间接肯定了我的处理程序没什么问题),但是从结果上来看的话好像是EMD做出来的要好一些
yghit08
发表于 2012-6-17 19:42
本帖最后由 yghit08 于 2012-6-17 20:00 编辑
利用分别利用上下包络做局部均值和包络,也就是在LMD原始文献中不区分局部极大值和局部极小值(统称为局部极值),利用局部极值做均值和包络分析。在EMD中区分局部极大/小值,分别利用极大/小值做包络,然后利用上下包络做均值和包络分析。在这一楼传的内容就是利用上下包络取代原LMD算法中相应的内容(不同于前面传的文件,这里不需要CombinMaxMin和LocalMeanEn这两个函数),处理信号的结果和前面传的那个LMD是一样的。本楼的图片是处理楼上给的那个信号的分离图(调频部分,因为我更关心频率,利用的是本楼的这个程序处理的)
雨夜1011
发表于 2012-6-18 09:33
回复 21 # yghit08 的帖子
你第一层的终止条件和原文上的是一样的吗?最后求出来的包络估计函数接近1。还有就是你最后去减的时候是怎么处理的,最后的终止条件是什么?我是直接去减的,但是应用于实际信号的时候会出现畸变,就是有的值已经超出了他应有的范围
yghit08
发表于 2012-6-18 12:43
回复 22 # 雨夜1011 的帖子
目前都没有做终止条件(指的是完全分离),针对你说的第一层的终止条件貌似不是和原文一样(完全小于1或者大于-1可能导致死循环,因此选用的小于1+delata或者大于-1-delta,delta=1e-6)。我也是直接去掉了,实际信号的用EMD处理了,见楼上的图。我认为没有必要做最后的终止条件,因为分离是从高频到低频(当然做也没问题,判断一下极点数就行)。实际信号的LMD见楼上,效果不如EMD,这个结果我发给了LMD的作者,他肯定了我做的程序没什么问题,估计他也应用了一下,结果和我给他的一样。你用的移动平均做的???是的话,分享一下啊!
雨夜1011
发表于 2012-6-18 16:16
我做的程序太乱了,估计会扰乱你的~~如果你想知道怎么做我可以告诉你~~加我QQ!!!117220306
cbjia
发表于 2012-6-18 16:59
LMD三次样条改造很容易实现,但是真正用到实际信号处理时,原始信号的信噪比对分离结果影响很大,建议在仿真信号中加入噪声,再观察分离结果
yghit08
发表于 2012-6-18 21:03
回复 24 # 雨夜1011 的帖子
加了。这样,方便联系一些
yghit08
发表于 2012-6-18 21:04
回复 25 # cbjia 的帖子
可以参考楼上利用实测信号做的结果,EMD和LMD做的比较的那一楼。你也是这样做的??
雨夜1011
发表于 2012-6-19 10:02
回复 25 # cbjia 的帖子
你也是做这方面的??你做出来的是什么样的~~~给个参考呗
yghit08
发表于 2012-6-19 11:30
回复 25 # cbjia 的帖子
图片是用三次样条替换移动平均法做的结果(楼上提及的悬臂梁的实测数据)的前三阶固有频率,结果和EMD相比几乎是一样的(用的程序是楼上提供的lmd_test)
李清志
发表于 2012-6-19 14:55
回复 1 # yghit08 的帖子
我试了人家给的改进的LMD程序,应用于仿真还可以,但是用于实际就不行了。改进的地方:三次样条插值,镜像延拓。
页:
1
[2]
3
4
5
6
7
8
9
10
11