EMD得到的波形如何进行HHT变换

hfutql

楼主才学EMD,按照论坛上的牛人装了EMD,得到了IMF1,IMF2,....可是接下来得到的波形如何进行HHT变化啊？附原程序：
  clc
clear all
close all
% 原始数据
fs = 1000;
ts = 1/fs;
t=0:ts:2;
z=sin(2*pi*6*t)+0.5.*(1+sin(2*pi*40*t)).*sin(2*pi*200*t)+(1+sin(2*pi*40*t)).*sin

(2*pi*400*t)+n(t);
figure
plot(t, z)
title('原始信号')
% EMD
imf=emd(z);
emd_visu(z,t,imf)
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
请问得到的结果如何再次进行HHT变换？         

仍小呆

楼主的意思好像不太明白。。。

hfutql

，谢谢回复！求高人！

h798471649

看看好像不太明白。

鱼音rao梁

hfutql 发表于 2014-3-14 09:51
，谢谢回复！求高人！

关键程序就是这几行：
imf=emd(z);                           %EMD分解
emd_visu(z,t,imf)                   %EMD分解结果展示
[A,f,tt]=hhspectrum(imf);       %求出IMF（本征模态函数）的瞬时频率与幅值
[im,tt]=toimage(A,f);             %将一维函数谱转化为图像
disp_hhs(im);                        %HHT谱图