emd.m程序竟然分解不开跳频信号？？？

fay1014

>> t = linspace(0,20,500);
>> tt = linspace(20,40,500);
>> xt(1:500)=cos(2*pi*t*1);;
>> xt(500:1000)=cos(2*pi*tt*2);
>>imf = emd(xt);%新程序中的emd.m
竟然只分解出一个imf，而且imf 原信号信号一摸一样，这是怎么回事吗？请大家帮忙分析一下
分解得到以下图形：

eight

原帖由 fay1014 于 2007-4-24 09:35 发表
>> t = linspace(0,20,500);
>> tt = linspace(20,40,500);
>> xt(1:500)=cos(2*pi*t*1);;
>> xt(500:1000)=cos(2*pi*tt*2);
>>imf = emd(xt);%新程序中的emd.m
竟然只分解出一个imf，而且imf 原信号信号一 ...

请先理解好EMD的算法流程，特别是IMF的定义。然后，你会发现，分不开是正常的，也是正确的

qinle

同意老八的说法,没有什么算法是万能的.不能只拿来主义,最好先看看是怎么回事再用.

zhlong

nod, 楼主的跳频信号已经满足了IMF的定义

fay1014

哦，谢谢指教，这就是传说中的模式混频现象吧
那这样emd分解不适合奇异性点的检测？

eight

原帖由 fay1014 于 2007-4-24 10:32 发表
哦，谢谢指教，这就是传说中的模式混频现象吧
那这样emd分解不适合奇异性点的检测？

这跟模式混叠不同，自己好好理解一下

适合

fay1014

谢谢你，初学emd，见笑了
我的课题是用emd分解来实现奇异点的检测
可是我每次分解出来的imf分量都很难看出奇异点，用hht谱看效果也不好
你能给点建议吗？非常感谢

eight

原帖由 fay1014 于 2007-4-24 10:50 发表
谢谢你，初学emd，见笑了
我的课题是用emd分解来实现奇异点的检测
可是我每次分解出来的imf分量都很难看出奇异点，用hht谱看效果也不好
你能给点建议吗？非常感谢

不好意思，刚才想得太快，EMD在奇异点检测上可能不及小波有效，所以，跟其他方法结合试试吧。下面是我师兄在这方面发的一篇文章，看对你有否帮助：

李合龙, 基于EMD和小波分析的建筑结构损伤检测能力的探讨. 中山大学学报. 44(6): 20-23, 2005

fay1014

非常感谢，我在中国期刊网下了这边文章，但是很多图都没显示，你那有这篇文章吗，能麻烦你发一份给我吗？邮箱：fay1014gre@yahoo.com.cn
非常非常感谢

eight

原帖由 fay1014 于 2007-4-24 11:54 发表
非常感谢，我在中国期刊网下了这边文章，但是很多图都没显示，你那有这篇文章吗，能麻烦你发一份给我吗？邮箱：fay1014gre@yahoo.com.cn
非常非常感谢

刚下载看了一下，跟我手上的是一样的

fay1014

我下载的图1，3，5都没显示，你那个图片有显示马

eight

原帖由 fay1014 于 2007-4-24 12:33 发表
我下载的图1，3，5都没显示，你那个图片有显示马

有显示，我下载的是pdf格式，你可能要安装最新版本的Acrobat reader或者Acrobat Professional

fay1014

我重新下了pdf的，这次一张图都没有，能麻烦你把图1，3，5剪切一下发上来马？
非常非常感谢，麻烦你了

fay1014

谢谢你，下了个新版本的，可以看了，非常的感谢