绘制HHT三维谱图例子

Cena马 · 发表于 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：视图角度

ben1984 · 发表于 2011-4-14 11:16

plot_HHT_3d.txt下来后打不开

chenlu1986 · 发表于 2011-4-14 21:51

和http://forum.vibunion.com/thread-28698-1-1.html有异曲同工之妙

chenlu1986 · 发表于 2011-4-15 10:42

回复 3 # ben1984 的帖子

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

zt861217 · 发表于 2011-4-18 22:33

请问楼主，我想把采样点数改成时间，该修改哪里的程序呢？

Cena马 · 发表于 2011-4-19 08:17

回复 7 # zt861217 的帖子

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

Cena马 · 发表于 2011-4-20 21:41

对于含有频率成分较多的信号，绘制出来的谱图最好平滑处理一下，以便于观察，下面两个图就是滚动轴承内圈故障信号HHT谱图平滑后的结果

cboboc · 发表于 2011-4-21 18:47

怎么下载不下来呢？下载后就是后缀为php的文件呢？

Cena马 · 发表于 2011-4-21 19:01

回复 13 # cboboc 的帖子

用IE下载，别用迅雷等下载工具

Cena马 · 发表于 2011-4-21 19:06

回复 13 # cboboc 的帖子

已经给你发过去了

cboboc · 发表于 2011-4-21 19:07

回复 14 # Cena马的帖子

恩，对，用IE 就可以下载下来，谢谢Cena马

chenlu1986 · 发表于 2011-4-26 08:59

回复 18 # 凤鸣荆楚的帖子

用IE浏览器直接点击下载来下。。。
再不行的话留言我发给你

sliyuchun · 发表于 2011-4-29 09:51

上面程序中的emd 什么意思来呀，求教

chenlu1986 · 发表于 2011-4-29 09:57

回复 24 # sliyuchun 的帖子

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

sliyuchun · 发表于 2011-4-29 10:06

谢过，不过好像在上面的程序中，emd运行不了呀，是不是有什么工具箱？

账号		自动登录	找回密码
密码			我要加入

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

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

评分

点评

评分

点评

评分

点评

评分

点评