声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2532|回复: 5

[HHT] 三维hht谱的说明

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

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

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

x
这几天将整个论坛的帖子过了一遍,受益匪浅,感谢各位前辈的无私奉献
看到很多人在问怎么做三维hht谱,其实论坛里已经有人给出了例子了,只是大家不会使用,这里我给出完整的例子,直接就可以运行
m=importdata('BYS.dat');
i=1;
N=1024;
f=500;
while i<1024
  
    data(i)=m(i);
    i=i+1;
end
imf=emd(data);
figure(1)
[m,n]=size(imf);

for i=1:m
   
    subplot(m,1,i);
    plot(imf(i,1:500));
   
end
figure(2)
dt=1/f;
h = nspab(imf',500,0,500,dt);
%h(1:5,:)=0;
subplot(211)
surf(h(1:200,11:1010))
shading interp
Xlabel('&Ecirc;±&frac14;&auml;(&micro;&atilde;&Ecirc;&yacute;)','FontSize',12);
Ylabel('&AElig;&micro;&Acirc;&Ecirc;(Hz)','FontSize',12);
Zlabel('&#8226;ù&Ouml;&micro;(Hz)','FontSize',12);
title('a','FontSize',12)
ms=mspc(h)
%ms(32)=1.2*ms(32)
subplot(212)
plot(1:250,ms(1:250))
Ylabel('&#8226;ù&Ouml;&micro;','FontSize',12);
Xlabel('&AElig;&micro;&Acirc;&Ecirc;(Hz)','FontSize',12)
title('c','FontSize',12)

运行该程序所需代码在http://www.veto.cn/down 看左边的hht


HHT QQ 群 11605343 资料共享

[ 本帖最后由 vfory 于 2008-9-17 21:58 编辑 ]
111.JPG

BYS.dat

14.75 KB, 下载次数: 79

回复
分享到:

使用道具 举报

发表于 2008-9-18 10:17 | 显示全部楼层
这是MATLAB 程序吧
有VC的吧
发表于 2016-1-17 09:33 | 显示全部楼层
谢谢分享
发表于 2016-1-19 20:09 | 显示全部楼层
谢谢分享
发表于 2016-2-22 15:33 | 显示全部楼层
谢谢分享!
发表于 2016-2-23 10:56 | 显示全部楼层
想问楼主mspc函数的程序代码在哪里能下到?楼主给的网站进不去呢
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-29 02:56 , Processed in 0.065665 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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