声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2656|回复: 11

[HHT] EMD分解问题!!!!!

[复制链接]
发表于 2012-9-24 15:57 | 显示全部楼层 |阅读模式

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

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

x
各位老师:
       你们好!
       我是一个初学EMD的菜鸟,应用xray老师的程序,对一调制信号进行分析,但是还是有些问题不太明白,希望各位老师给予指导!
%%%%%程序
clear;clc;close all;
fs=100;T=10;N=T*fs;
t=0:1/fs:T;                   % 采样率100HZ
z=(1+sin(2*pi*t)).*sin(2*pi*10*t);%调幅信号
imf = emd(z); %%%%%%%%%%EMD
[A, fa, tt] = hhspectrum(imf);
if size(imf,1) > 1
    [A,fa,tt] = hhspectrum(imf(1:end-1, :));
else
    [A,fa,tt] = hhspectrum(imf);
end
[E, tt1] = toimage(A,fa,tt,length(tt));   
for k = 1:size(E,1)
    bjp(k) = sum(E(k,:))*1/fs*1/T;  
end
f = (0:N-2)/N*(fs/2);
figure(1)
plot(t,z);
figure(2)
imagesc(tt1,[0,0.5*fs],E);
set(gca,'YDir','normal')
% 使用灰度图显示
% colormap(flipud(gray))
figure(3)
plot(f, bjp);
分析结果如图:
我有几个问题不明白:
1。EMD分解后,各阶的IMF是存在什么位置的?工作空间里的矩阵A,和fa都是9行的矩阵,不知道哪个才是?
2。边际谱中只包含10HZ的频率成分,而我的主要目的是获得调制信号中的包络信号的频率(1HZ),应该怎么获得?
3。应用EMD分解,我是否可以直接从分解中的IMF中提取1HZ的频率信号,再进行FFT变化,从而,获得1HZ的频率?
以上是我的问题,希望各位老师给予指导!!!




2.jpg
3.jpg
1.jpg

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

 楼主| 发表于 2012-9-25 11:39 | 显示全部楼层
怎么没人帮我回答啊!各位大侠,指导一下啊!!
发表于 2012-9-25 21:32 | 显示全部楼层
imf <10x1001>:表明有10个imf,只不过第2个至第10个imf的幅度都很小。第1个imf与原信号z几乎相等。这是因为,imf本来就是一个AM-FM函数。

求信号包络的最简单的办法是先做Hilbert变换,得到解释信号,再对解释信号(复信号)求模值,这个模值即是你想要的包络。
 楼主| 发表于 2012-9-25 22:16 | 显示全部楼层
谢谢你的回答。
我想知道的是:是否可以通过EMD分解获得包络信号,比如某一个IMF就是包络信号。
另外,通过边际谱的计算,是否可以获得包络信号的频率,即1HZ的频率。
 楼主| 发表于 2012-9-27 22:47 | 显示全部楼层
有木有人帮我回答一下啊!!!!
发表于 2012-9-29 17:00 | 显示全部楼层
如果是纯调幅信号的话,我想你通过EMD做不出来它包络(被认为是某一阶IMF),但是你做出来后通过hhspectrum这个函数处理某一阶IMFs的话,能做出来包络,这样做出来的结果正如二楼说的是直接通过HIlbert变换做出来的(hhspectrum的功能去除求解瞬时频率的部分就是二楼说的那样的做法)。看来楼主对EMD没了解啊!
发表于 2012-11-6 14:46 | 显示全部楼层

是用解析信号的模的正负当做上下包络吗?还是另有算法?
发表于 2012-11-7 12:58 | 显示全部楼层
用Hilbert法求信号包络,前提条件是该信号已经是单分量信号。

若是多分量信号,则可用EMD中的方法来求上下包络。
发表于 2012-11-8 17:31 | 显示全部楼层
dsp2008 发表于 2012-11-7 12:58
用Hilbert法求信号包络,前提条件是该信号已经是单分量信号。

若是多分量信号,则可用EMD中的方法来求上 ...

单分量信号意思是IMF吗?哦,这样的话最大最小值差不多就是关于横坐标轴对称的是吧
发表于 2012-11-8 19:04 | 显示全部楼层
回答正确!加十分。
发表于 2012-11-9 15:32 | 显示全部楼层
这个问题问得非常好,我也在思考。
发表于 2014-12-22 20:36 | 显示全部楼层
怎么又看不到???能不能不设置门槛?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-24 14:16 , Processed in 0.089437 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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