hyhwqy 发表于 2009-3-5 21:36

关于HHT画三维图的时间点数的问题

各位高手：
我用一下程序做HHT，为什么画出的图形时间点数只有400点，怎么样得到与采样点数相同的图啊？非常急，希望各位高手能帮助解答，小女子感激不尽！！！
clear;
fs=2000;
N=1024;
t=0:1/fs:N/fs;
z=sin(20*pi*t)+sin(200*pi*t)+sin(400*pi*t);
imf=emd(z);
=size(imf);
for i=1:m
   
    subplot(m,1,i);
    plot(imf(i,:));
   
end
figure(3)
dt=1/fs;
clear;
fs=2000;
N=1024;
t=0:1/fs:N/fs;
z=sin(20*pi*t)+sin(200*pi*t)+sin(400*pi*t);
imf=emd(z);
=size(imf);
for i=1:m
   
    subplot(m,1,i);
    plot(imf(i,:));
   
end
figure(3)
dt=1/fs;

h = nspab(imf',500,0,500,0,N/fs);
subplot(211)
surf(h(1:200,50:end-20))
shading interp
xlabel('时间(点数)','FontSize',12);
Ylabel('频率(Hz)','FontSize',12);
zlabel('幅值','FontSize',12);
title('a','FontSize',12)
view([-75,25])
yt=subplot(223)
imagesc(h(1:200,:))
Ylabel('频率(Hz)','FontSize',12);
xlabel('时间(点数)','FontSize',12);
set(yt,'ydir','nor')
title('b','FontSize',12)
ms=mspc(h)
subplot(224)
plot((1:length(ms)),ms)
xlabel('频率(Hz)','FontSize',12);
ylabel('幅值','FontSize',12);
title('c','FontSize',12)
还有一个问题，有高手说用h = nspab(imf',500,0,500,dt);可以，但是我用这个命令时报错，必须用h = nspab(imf',500,0,500,t0,t1);格式才能识别，是不是我的hht_toolbox工具箱版本的问题？？实在搞不懂，希望各位高手能指点一下哈。




[ 本帖最后由 hyhwqy 于 2009-3-5 21:41 编辑 ]

hyhwqy 发表于 2009-3-6 09:43

怎么没有人回复啊？？？还望各位高手指点一二啊！:handshake

qqvirile 发表于 2009-3-6 17:29

回复 楼主 hyhwqy 的帖子

应该是工具箱的问题 我也照做了 是1000点的！但是我没有mspc.m可否给我传一个过来啊！！

hyhwqy 发表于 2009-3-6 19:50

回复 楼主 hyhwqy 的帖子
应该是工具箱的问题 我也照做了 是1000点的！但是我没有mspc.m可否给我传一个过来啊！！
是我下载的hht_toolbox_20040808版本太低？谁有新版本的啊？可以传一个嘛？

qqvirile 发表于 2009-3-6 20:08

回复 5楼 hyhwqy 的帖子

这个论坛上有的
把你qq给我吧 我正好也需要你的mspc.m文件

haifeng1213 发表于 2009-3-13 19:59

能问一下吗
surf(h(1:200,50:end-20))
为什么取1：200和50：end-20
新手，谢谢指教

shants 发表于 2009-3-25 18:24

我同样也遇到了这个问题，怎样才能使三维谱图的时间与采样点数一致呢，
哪位高手可以指点一下

yll 发表于 2009-3-26 09:28

回复 楼主 hyhwqy 的帖子

我运行那mspc时报错呢，你那ms=mspc(h)，只要一个h吗？我那还要输入f，ms=mspc(h，f)，f是什么呢？如何确定？还是我的mspc有问题，谢谢

hyhwqy 发表于 2009-4-16 14:22

本帖最后由 VibInfo 于 2016-10-21 15:25 编辑

原帖由 yll 于 2009-3-26 09:28 发表
我运行那mspc时报错呢，你那ms=mspc(h)，只要一个h吗？我那还要输入f，ms=mspc(h，f)，f是什么呢？如何确定？还是我的mspc有问题，谢谢
我把我的mspc给你，你试试

abc198478 发表于 2009-4-18 17:13

nspab是什么文件，我的工具箱里怎么没有啊？？你们的都是在哪里下载的？？

wains 发表于 2009-5-14 10:49

回复楼主

用h = nspab(imf',500,0,500,dt);效果见附图；
可能就是你用的hht_toolbox工具箱版本不是最新的。

[ 本帖最后由 wains 于 2009-5-14 10:51 编辑 ]

lj2tt 发表于 2009-7-30 12:54

我的也运行不对，难道也是工具箱版本太低了

misslfx 发表于 2009-8-17 16:09

今天我也遇到了这个问题，个人认为是工具箱版本太低所致。
解决方法，上网搜索下载hht_toobox工具箱，将nspab函数改成nspabz函数
% Z. Shen (JHU)                        July 2, 1995 Initial      -- nspab
% D. Xiang (JHU)                  March 27, 2002 Modified
% J. Marshak (NASA GSFC)      Feb. 11, 2004 Modified --nspabz
%                              (replaced zfam() by 'zfapanls()')
使用nspab时：

使用nspabz时：

当前的问题是第三幅边际谱仍然不对，而且200HZ处的频率也不对，不知道是什么原因。

qingqing5638 发表于 2009-10-19 20:31

回复 13楼 misslfx 的帖子

传给我此工具箱好吗我怎么搜索不到啊多谢qingqing5638@yahoo.com.cn

xiaokang 发表于 2010-1-17 20:37

现在存在这样的问题，见下面网址
http://forum.vibunion.com/forum/viewthread.php?tid=82724&pid=457335&page=1&extra=#pid457335
页: [1] 2
查看完整版本: 关于HHT画三维图的时间点数的问题