三维hht谱的说明

vfory

这几天将整个论坛的帖子过了一遍，受益匪浅，感谢各位前辈的无私奉献
看到很多人在问怎么做三维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 编辑 ]

luoye919

这是MATLAB 程序吧
有VC的吧

1713573225

谢谢分享

hylqqq

谢谢分享

欣然的名单

谢谢分享！

欣然的名单

想问楼主mspc函数的程序代码在哪里能下到？楼主给的网站进不去呢