|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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('ʱ¼ä(µãÊý)','FontSize',12);
Ylabel('ƵÂÊ(Hz)','FontSize',12);
Zlabel('•ùÖµ(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('•ùÖµ','FontSize',12);
Xlabel('ƵÂÊ(Hz)','FontSize',12)
title('c','FontSize',12)
运行该程序所需代码在http://www.veto.cn/down 看左边的hht
HHT QQ 群 11605343 资料共享
[ 本帖最后由 vfory 于 2008-9-17 21:58 编辑 ] |
|