杭州锐达数字技术有限公司
查看: 23325|回复: 158

[HHT] 绘制HHT三维谱图例子

  [复制链接]
发表于 2011-4-13 21:25 | 显示全部楼层 |阅读模式

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

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

x
论坛中HHT时频图、边际谱图等绘图函数都已经介绍很多了,我看论坛上有不少人求助关于HHT三维图的绘制方法,实际很简单,要做的就是把利用计算出的瞬时频率将计算出的瞬时幅值映射到你所划分的频率轴上就可以了(映射后是个二维向量),把这个二维向量当成z,你所划分的频率轴当成y,采样点数当成x绘制三维图就可以了。我把程序(plot_hht_3d.m)及运行结果贴出来!PS:由于传不了m文件我用txt格式上传的,要安装到rilling的工具包,计算瞬时频率的时候用到了hhspectrum函数!
运行以下语句即可:
clc
clear
close all
fs = 1024;
n = 1024;
t = 0:1/fs:(n-1)/fs;
y = sin(2*pi*50*t) + 0.4*cos(2*pi*150*t);
figure,plot(y)
set(gca,'xlim',[0,1024])
imf = emd(y);
plot_hht_3d(imf,1024,1024,[-60,50])
函数参数分别为:
imf:emd分解后的结果
numfreq:划分的频率个数
fs:采样频率
ANGLE:视图角度

原始信号

原始信号

三维图

三维图

plot_hht_3d.txt

888 Bytes, 下载次数: 997

换成m文件即可

pack_emd.zip

98.17 KB, 下载次数: 55

评分

2

查看全部评分

回复
分享到:

使用道具 举报

发表于 2011-4-14 11:16 | 显示全部楼层
plot_HHT_3d.txt下来后打不开

点评

需要将文件名的后缀改为.m  发表于 2015-12-2 15:09
发表于 2011-4-14 21:51 | 显示全部楼层
发表于 2011-4-15 10:42 | 显示全部楼层
回复 3 # ben1984 的帖子

按照楼主说的下载后后缀名修改一下试试

评分

1

查看全部评分

发表于 2011-4-18 22:33 | 显示全部楼层
请问楼主,我想把采样点数改成时间,该修改哪里的程序呢?
 楼主| 发表于 2011-4-19 08:17 | 显示全部楼层
回复 7 # zt861217 的帖子

把那个计算出来的t改下就行了

点评

楼主讲下怎么改t,我也不懂这个问题,求解答  详情 回复 发表于 2016-6-4 14:52
 楼主| 发表于 2011-4-20 21:41 | 显示全部楼层
对于含有频率成分较多的信号,绘制出来的谱图最好平滑处理一下,以便于观察,下面两个图就是滚动轴承内圈故障信号HHT谱图平滑后的结果

3d

3d

时频图

时频图

评分

2

查看全部评分

发表于 2011-4-21 18:47 | 显示全部楼层
怎么下载不下来呢?下载后就是后缀为php的文件呢?
 楼主| 发表于 2011-4-21 19:01 | 显示全部楼层
回复 13 # cboboc 的帖子

用IE下载,别用迅雷等下载工具
 楼主| 发表于 2011-4-21 19:06 | 显示全部楼层
回复 13 # cboboc 的帖子

已经给你发过去了

点评

赞成: 5.0
赞成: 5
热心助人  发表于 2011-4-21 20:00
发表于 2011-4-21 19:07 | 显示全部楼层
回复 14 # Cena马 的帖子

恩,对,用IE 就可以下载下来,谢谢Cena马
发表于 2011-4-26 08:59 | 显示全部楼层
回复 18 # 凤鸣荆楚 的帖子

用IE浏览器直接点击下载来下。。。
再不行的话留言我发给你
发表于 2011-4-29 09:51 | 显示全部楼层
上面程序中的emd 什么 意思来呀,求教
发表于 2011-4-29 09:57 | 显示全部楼层
回复 24 # sliyuchun 的帖子

emd是经验模态分解。
是HHT中的Huang方法。
建议多看下书再看程序啊

评分

1

查看全部评分

发表于 2011-4-29 10:06 | 显示全部楼层
谢过,不过好像在上面的程序中,emd运行不了呀,是不是有什么工具箱?

点评

如何进行平滑处理  发表于 2013-4-9 11:47
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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