声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3004|回复: 15

[HHT] 关于HHT画三维图的问题

[复制链接]
发表于 2008-5-5 10:46 | 显示全部楼层 |阅读模式

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

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

x
用的是小议emd的使用里面的例子,结果总是报这个错,实在是不明白怎么回事

??? Undefined function or variable 'dt'.
Error in ==> nspab at 45
omg=abs(diff(unwrap(angle(data))))/(2*pi*dt);

谢谢大家啊

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2008-5-5 11:19 | 显示全部楼层

回复 楼主 的帖子

Undefined function or variable 'dt'.
dt 没做设置吧
 楼主| 发表于 2008-5-5 15:22 | 显示全部楼层
谢谢你的回复,很弱的问问dt的具体含义是什么,应该怎么设置啊  
谢谢
发表于 2008-5-5 16:35 | 显示全部楼层
dt是采样时间间隔,即采样率的倒数。
 楼主| 发表于 2008-5-6 09:25 | 显示全部楼层
你好,我刚才试了试,设置了dt可还是不行啊,快急死了,我用的例子就是前几天论坛上发的那个,结果还是和原来一样报错,盼望各位大牛帮忙看看啊
>> imf=emd(s1);
[m,n]=size(imf);
for i=1:m
   
    subplot(m,1,i);
    plot(imf(i,:));
   
end
figure(3)
dt=1:2048;
h = nspab(imf',500,0,500,dt);
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)
??? Undefined function or variable 'dt'.
Error in ==> nspab at 45
omg=abs(diff(unwrap(angle(data))))/(2*pi*dt);

发表于 2008-5-6 09:29 | 显示全部楼层
dt=1:2048? 这和采样频率的倒数有关系吗
 楼主| 发表于 2008-5-6 09:50 | 显示全部楼层
不好意思,手误了,应该是dt=1/2048,运行的时候还是一样报错。
发表于 2008-5-6 10:56 | 显示全部楼层

回复 7楼 的帖子

nspab函数中的500 也要变一下吧 不能还是500,我现在电脑上没有MATLAB,没办法验证
 楼主| 发表于 2008-5-6 10:57 | 显示全部楼层
好的,谢谢回复啊。我先试试
发表于 2008-5-11 16:12 | 显示全部楼层
应该是nspab函数的问题吧?
h = nspab(imf',500,0,500,dt);这行好象调用的不对吧?
function [h,xs,w] = nspab(data,nyy,minw,maxw,t0,t1)
%        t0        - the start time
%        t1        - the end time
我试了一下设置了一个终止时间就可以了
 楼主| 发表于 2008-5-11 16:35 | 显示全部楼层
你好,谢谢你的回复,不过能解释一下nyy,minw,maxw这三个参数的具体含义么,我刚才试了一下,还是报错

ndefined function or variable 'nyy'

谢谢啊
发表于 2008-5-17 09:42 | 显示全部楼层
nyy是频率分辨率,minw,maxw是最小频率和最大频率
 楼主| 发表于 2008-5-18 12:39 | 显示全部楼层
谢谢楼上的回复啊,受教了
发表于 2009-3-25 09:46 | 显示全部楼层

回复 8楼 zhangnan3509 的帖子

我想请教一下,nspab中的500如何根据现实信号设置,谢谢
发表于 2009-3-25 09:49 | 显示全部楼层

回复 8楼 zhangnan3509 的帖子

h = nspab(imf',500,0,500,dt);
surf(h(1:200,50:end-20))
imagesc(h(1:200,:))
第一句那个500如何设置,还有第二句和第三句是什么意思?万分感谢
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-17 00:40 , Processed in 0.077008 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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