关于HHT画三维图的问题

svca_0

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

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

谢谢大家啊

zhangnan3509

Undefined function or variable 'dt'.
dt 没做设置吧

svca_0

谢谢你的回复，很弱的问问dt的具体含义是什么，应该怎么设置啊  
谢谢

jinnian

dt是采样时间间隔，即采样率的倒数。

svca_0

你好，我刚才试了试，设置了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);

zhangnan3509

dt=1:2048? 这和采样频率的倒数有关系吗

svca_0

不好意思，手误了，应该是dt=1/2048，运行的时候还是一样报错。

zhangnan3509

nspab函数中的500 也要变一下吧 不能还是500，我现在电脑上没有MATLAB，没办法验证

svca_0

好的，谢谢回复啊。我先试试

prideheart

应该是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
我试了一下设置了一个终止时间就可以了

svca_0

你好，谢谢你的回复，不过能解释一下nyy，minw，maxw这三个参数的具体含义么，我刚才试了一下，还是报错

ndefined function or variable 'nyy'

谢谢啊

翟瑞光

nyy是频率分辨率，minw，maxw是最小频率和最大频率

svca_0

谢谢楼上的回复啊，受教了

yll

我想请教一下，nspab中的500如何根据现实信号设置，谢谢

yll

h = nspab(imf',500,0,500,dt);
surf(h(1:200,50:end-20))
imagesc(h(1:200,:))
第一句那个500如何设置，还有第二句和第三句是什么意思？万分感谢