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