求助HHT作图问题

cliffking

我还想请教你一个问题，第二副图形 f2c是代表什么意思呢？

[ 本帖最后由 zhlong 于 2007-6-8 21:13 编辑 ]

zhangnan3509

在小议EMD程序里我解释过了，你好好观察一下就明白了

cliffking

请教对于http://forum.vibunion.com/forum/thread-44797-3-2.html 40楼的例子，我用什么办法能变成时间---频率---幅值得三维图呢 谢谢指导！

zhangnan3509

在小议EMD程序里面有

cliffking

恩 我看到EMD小议里面你写的这个程序。是用这个程序吗？
m=importdata();
i=1;
N=1024;
f=500;
while i<1024
  
    data(i)=m(i);
    i=i+1;
end
imf=emd(data);
figure(1)
[m,n]=size(imf);

for i=1:m
   
    subplot(m,1,i);
    plot(imf(i,1:500));
   
end
figure(2)
dt=1/f;
h = nspab(imf',500,0,500,dt);
%h(1:5,:)=0;
subplot(211)
surf(h(1:200,11:1010))
shading interp
Xlabel('&Ecirc;±&frac14;&auml;(&micro;&atilde;&Ecirc;&yacute;)','FontSize',12);
Ylabel('&AElig;&micro;&Acirc;&Ecirc;(Hz)','FontSize',12);
Zlabel('&#8226;ù&Ouml;&micro;(Hz)','FontSize',12);
title('a','FontSize',12)
ms=mspc(h)
%ms(32)=1.2*ms(32)
subplot(212)
plot(1:250,ms(1:250))
Ylabel('&#8226;ù&Ouml;&micro;','FontSize',12);
Xlabel('&AElig;&micro;&Acirc;&Ecirc;(Hz)','FontSize',12)
title('c','FontSize',12)

cliffking

是不是要把这个函数也和主函数放在一个文件夹内呢。

zhangnan3509

嗯！但是你要有数据，这不是仿真的程序

cliffking

请教 能告诉我那个例子分解出来的IFM，我想把它画出时间--频率---幅值得三维图应该怎么弄呢，小弟是新手，有很多都没有弄明白， 恳请指教 谢谢

cliffking

就是那个例子中的两个正弦相加的信号不能做出 时间--频率---幅值得三维图的三维图谱吗？ 请教

zhangnan3509

程序都给你了，再画不出来我就没办法了:'(

zhangnan3509

一开始你拿到的程序就是8楼zhlong的，和这个小议里面的结合一下就可以了，把前面的导入数据部分，变成你要处理的仿真信号选好点就可以了，主体框架已经有了，细节上你变一下就行了

cliffking

能不能把你的http://forum.vibunion.com/forum/thread-44797-3-2.html  26楼的那几幅图形和程序发给我
thankyou！！ cliffking2006@yahoo.com.cn

zhangnan3509

小议里面的，因为这个仿真信号对我没什么实际价值，原程序我没保存。

[ 本帖最后由 zhlong 于 2007-6-10 14:48 编辑 ]

cliffking

能不能帮我用这个例子作一个时间---频率---幅值得三维图谱，我真是弄不出来了:'( ，我弄了很长时间，:@Q
t=1:500;
t=t*1/400;
x1=sin(2*pi*20*t);
x2=0.4*sin(2*pi*40*t+140);
z=x1+x2
plot(t,z)
y= linspace(1,2,500);%添加趋势项
z=z+y
imf=emd(z);
emd_visu(z,1:length(z),imf)
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))

zhangnan3509

f=400;
y=linspace(1,2,500);

t=1:500;
t=t*1/400;
x1=sin(2*pi*20*t);
x2=0.4*sin(2*pi*40*t+140);
z=x1+x2+y
%在后面添加一个趋势向可以把三个分量完全分解出来
imf=emd(z);

emd_visu(z,1:length(z),imf)
figure (4)
dt=1/f;
h = nspab(imf',500,0,500,dt);
h(1:2,:)=0;

subplot(211)
surf(h(1:200,1: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)
简单写了一下，你先看看。（这几天没怎么在线，不好意思）