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

hyhwqy

各位高手：
我用一下程序做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);
[m,n]=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);
[m,n]=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

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

qqvirile

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

hyhwqy

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

是我下载的hht_toolbox_20040808版本太低？谁有新版本的啊？可以传一个嘛？

qqvirile

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

haifeng1213

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

shants

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

yll

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

hyhwqy

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

我把我的mspc给你，你试试

abc198478

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

wains

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

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

lj2tt

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

misslfx

今天我也遇到了这个问题，个人认为是工具箱版本太低所致。
解决方法，上网搜索下载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

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

xiaokang

现在存在这样的问题，见下面网址
http://forum.vibunion.com/forum/ ... mp;extra=#pid457335