test.m

新的信号处理方法－希尔伯特黄变换的EMD分解MATLAB程序 把程序看懂 对希尔伯特黄变换的原理就理解了

clear
clc
t1=0:.001:2.5;
t2=2.5:.001:5;
t=0:.001:(5+.001);
a=sin(20*pi*t1+0.5*sin(5*pi*t1));
b=sin(10*pi*t2+.5*sin(5*pi*t2));
c=0.5*cos(3*pi*t);
x=[a(:,:) b(:,:)]+c;

stop = [0.1,0.1,0.1]; 
[imf,ort,nbits]=myemd(x,t,stop,2);


imfsize=size(imf);
columnNum=imfsize(1)/2+1;
for i=1:columnNum(1)
    subplot(columnNum,2,i);plot(imf(i,:));
end