声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1756|回复: 4

[HHT] 频谱显示时出现的问题

[复制链接]
发表于 2008-4-2 15:27 | 显示全部楼层 |阅读模式

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

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

x
我在运行HHT程序时,没有使用toimage和disp_hhs这两个函数而是直接用plot输出瞬时频率,得到的图如 第一张图所示,而使用toimage和disp_hhs这两个函数得到的图如第二张所示,请问在第一张图中端点处出现一些幅度很大的线,该怎么理解呢?主程序如下

fs=200;
N=1001;
n=0:1:N-1;dt=1/fs;t=n*dt;
x1=sin(2*pi*10*t);
x2=sin(2*pi*6*t);
x3=cos(2*pi*3*t);
x=x1+x2+x3;

%经验模态分解
%imf= emd(x,t,[0.05,0.5,0.05],2);
imf=emd(x)
imf=emd_visu(x,t,imf)%EMD分解
[A,f,tt]=hhspectrum(imf);
ff=f*fs;
figure
plot(tt,ff,'k')
xlabel('时间点数')
ylabel('实际频率')
title('希尔伯特-黄谱')

%[A,f,tt]=hhspectrum(imf);
%[im,tt]=toimage(A,f);
%disp_hhs(im);
%colormap(flipud(gray))

[localimg=514,406]5[/localimg]

[localimg=506,406]4[/localimg]
HHT2.jpg
HHT.jpg
回复
分享到:

使用道具 举报

 楼主| 发表于 2008-4-2 15:30 | 显示全部楼层
不好意思上面的图我发反了,第一张图是使用toimage跟disp_hhs函数得到的图,第二张是没有使用这两个函数直接用plot 输出的
发表于 2008-4-2 22:32 | 显示全部楼层

回复 2楼 的帖子

建议楼主修改disp_hhs这个函数就可以 得到实际频率的图形。
如果直接用[A,f,tt]=hhspectrum(imf);这个函数,请看清楚输出变量的说明 tt是被截断的时间,不能直接拿来用
发表于 2011-12-2 16:19 | 显示全部楼层
顶起,期望大牛出现
发表于 2013-3-20 17:37 | 显示全部楼层
能把dis_hhs的原程序发给我吗?我下了个错的ludaheu@163.com
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-25 09:37 , Processed in 0.066877 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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