声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1211|回复: 5

[HHT] 请教信号预处理问题

[复制链接]
发表于 2008-12-19 21:10 | 显示全部楼层 |阅读模式

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

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

x
经验模式分解(EMD)方法在处理imf分量间能量差很远的情况时,不能很好的辨识各imf分量,有文献说,可以用基于傅里叶变换的emd方法来处理这类问题。
不知有没有现成的matlab程序?多谢
比如
t = 0:1/3200:0.24
x = 10*sin(2*pi*50*t) + sin(2*pi*150*t);
对于这个信号,做emd处理,无法得到频率为50和150的两个imf分量啊。。
请大家指教~
回复
分享到:

使用道具 举报

发表于 2008-12-19 22:35 | 显示全部楼层
我没有做过这个方面的尝试
但是EMD本身就是存在这样的问题
模态混叠
端点效应

期待有人来解答楼主的问题
发表于 2008-12-22 15:52 | 显示全部楼层
楼主的采样频率是多少呀?
发表于 2008-12-30 11:32 | 显示全部楼层

看本论坛的精贴

clear;
fs=1000;    %fs为采样频率;
N=1000;     %采样点数
t=1/fs:1/fs:1;
y1=10*sin(2*pi*50*t);
y2=sin(2*pi*150*t);
y=[y1;y2;zeros(size(y1))]; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);
[E,tt1]=toimage(A,fa,tt,length(tt));
%E=flipud(E);
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;  
end
f=(0:N-3)/N*fs/2;
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
QQ截图未命名.jpg

评分

1

查看全部评分

发表于 2008-12-30 11:34 | 显示全部楼层

回复 地板 qqvirile 的帖子

我也是昨天刚看的 今天正好拿来用一下 这个程序没有emd分解 而是直接 y1=10*sin(2*pi*50*t);y2=sin(2*pi*150*t);作为单分量函数 大家一起学习吧!我也是初学!
发表于 2008-12-30 19:42 | 显示全部楼层
这个程序没有emd分解 而是直接 y1=10*sin(2*pi*50*t);y2=sin(2*pi*150*t);作为单分量函数
正弦信号本身可以作为一个IMF
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 23:19 , Processed in 0.065599 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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