EMD得到的波形如何进行HHT变换
楼主才学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)
=hhspectrum(imf);
=toimage(A,f);
disp_hhs(im);
请问得到的结果如何再次进行HHT变换?
楼主的意思好像不太明白。。。 {:{39}:},谢谢回复!求高人! 看看好像不太明白。 hfutql 发表于 2014-3-14 09:51
,谢谢回复!求高人!
关键程序就是这几行:
imf=emd(z); %EMD分解
emd_visu(z,t,imf) %EMD分解结果展示
=hhspectrum(imf); %求出IMF(本征模态函数)的瞬时频率与幅值
=toimage(A,f); %将一维函数谱转化为图像
disp_hhs(im); %HHT谱图
页:
[1]